Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1690 1612 (N28W64 L=314)  2161    169        3 (1629/1554-1739/1673) -90.0 W59-W74 N24-N33
      2 1333 1585 (N31W29 L=278)   610    116        1 (1280/1529-1375/1620) -60.3 W26-W31 N26-N33
      3 1674 1460 (N19W52 L=301)   458     92        0 (1637/1430-1711/1489) -45.0 W49-W54 N16-N20
      4  964 1457 (N26W01 L=250)   284     82        0 ( 927/1437-1012/1474) -31.0 E01-W04 N24-N27
      5 1632  908 (S14W38 L=287)   305     84        0 (1603/ 881-1650/ 950)  53.1 W36-W39 S15-S11
      6  723  864 (S07E20 L=229)   215     63        0 ( 708/ 830- 740/ 891)  56.3 E20-E18 S08-S05
      7 1786  619 (S34W59 L=308)   559     56        1 (1766/ 594-1821/ 647) -15.3 W54-W65 S35-S32
      8 1540  586 (S34W32 L=282)   373     66        1 (1511/ 559-1557/ 622)  90.0 W29-W34 S35-S31
      9  966  471 (S34E11 L=238)   467    125        1 ( 942/ 421- 997/ 528) -90.0 E12-E09 S38-S30
     10  963  248 (S52E20 L=229)   311     24        1 ( 945/ 226- 984/ 269)  76.0 E21-E17 S54-S50
     11  806  171 (S57E45 L=204)   804     88        0 ( 771/ 141- 833/ 207)   9.5 E56-E38 S60-S53
  

    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