Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  344 1350 (N35E50 L= 65)   326     42        0 ( 317/1335- 373/1365)   0.0 E53-E46 N33-N35
      2 1180 1255 (N18W13 L=128)   343    162        0 (1158/1190-1221/1327) -83.7 W12-W14 N13-N22
      3 1261 1288 (N19W19 L=134)   355    169        1 (1192/1228-1307/1326) -54.5 W15-W21 N14-N21
      4 1686 1024 (S05W44 L=158)   522    119        1 (1639/ 995-1732/1058) -35.5 W40-W47 S07-S01
      5  930  602 (S18E12 L=102)   594    203        3 ( 885/ 528- 986/ 681)  45.0 E16-E07 S22-S13
      6 1186  520 (S28W03 L=118)   750    328        1 (1088/ 412-1284/ 630)  38.7 E05-W10 S34-S21
      7  612  489 (S22E37 L= 77)   649    202        0 ( 538/ 427- 675/ 546)  23.6 E46-E30 S26-S18
      8  733  403 (S30E32 L= 83)  1029    270        2 ( 629/ 340- 829/ 478)  20.0 E43-E21 S33-S25
      9 1196  396 (S38W02 L=116)   486    157        1 (1128/ 362-1259/ 439)  21.8 E03-W07 S39-S35
  

    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