Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1685 1331 (N17W51 L=199)   482    146        0 (1627/1316-1751/1356) -20.0 W46-W58 N15-N19
      2  914 1039 (N07E07 L=141)   381    157        0 ( 890/ 966- 943/1099) -90.0 E07-E04 N02-N10
      3 1730 1063 (N00W49 L=197)   377     55        2 (1705/1040-1753/1093) -90.0 W46-W51 S01-N01
      4 1200  946 (S01W10 L=158)   251     86        1 (1158/ 917-1235/ 966) -21.8 W07-W12 S03-N00
      5  901  818 (S07E09 L=139)   492    209        0 ( 836/ 760- 986/ 876) -41.2 E12-E04 S10-S02
      6 1591  701 (S21W37 L=185)  1574    278        4 (1485/ 619-1675/ 801)  39.8 W28-W43 S25-S14
      7  932  716 (S13E08 L=140)   290    122        1 ( 873/ 695- 967/ 751)   9.5 E11-E05 S14-S11
      8 1391  490 (S33W22 L=170)   499     69        0 (1355/ 462-1424/ 523)  31.0 W19-W24 S35-S31
      9 1210  293 (S48W08 L=156)   330     35        0 (1190/ 276-1232/ 314) -45.0 W06-W10 S50-S45
  

    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