Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  858 1808 (N44E39 L= 53)  3003    318        4 ( 804/1690- 952/1941) -23.1 E72-E24 N36-N56
      2  416 1606 (N19E60 L= 32)   985    113        3 ( 387/1563- 440/1658) -36.9 E64-E55 N16-N22
      3  601 1537 (N18E40 L= 52)   610    120        2 ( 579/1482- 639/1602) -58.0 E43-E36 N14-N22
      4  581 1441 (N12E37 L= 54)   272     79        1 ( 558/1418- 617/1467) -21.8 E39-E34 N11-N13
      5  285 1316 (S01E57 L= 35)  1114    165        2 ( 253/1247- 319/1395)  90.0 E59-E54 S05-N04
      6  167 1172 (S11E67 L= 25)  1248    115        7 ( 129/1137- 202/1230) -29.1 E74-E62 S13-S07
      7 1016 1054 (S01E01 L= 91)   314     77        1 ( 996/1020-1040/1095)  76.0 E02-W00 S03-N01
      8 1420  934 (N01W25 L=117)   955    211        3 (1370/ 833-1468/1014) -37.9 W20-W30 S03-N05
      9 1495  583 (S17W41 L=133)   401     71        1 (1474/ 556-1524/ 618) -18.4 W38-W44 S18-S15
  

    OPENR, u, binfile, /COMP, /GET_LUN
    head = FLTARR(6) ;[xsize,ysize,radius,hh,mm,ss]
    READU, u, head
    data = INTARR(FIX(head[0]),FIX(head[1]))
    READU, u, data
    FREE_LUN, u
    

• each filament has the value of Fil_ID
• all values are in pixel

  c = WHERE(data EQ Fil_ID, area)
  f = '(2(A,I6))'
  PRINT, FORM=f, "Area of Fil_ID", Fil_ID," =", area 

  data = ROTATE(data,7)
  Solar_y = 2048-Solar_y