Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1627 1545 (N21W59 L=187)   579    114        0 (1579/1527-1690/1565) -29.1 W54-W63 N17-N22
      2 1254 1490 (N28W28 L=156)   397     90        3 (1218/1460-1285/1528)  15.9 W24-W31 N26-N30
      3 1634 1468 (N16W54 L=182)   394     71        0 (1614/1438-1657/1499) -90.0 W52-W55 N14-N18
      4  199 1045 (N23E62 L= 66)  2089    226        3 ( 151/ 941- 247/1117) -50.7 E70-E56 N16-N27
      5 1854  943 (S20W58 L=186)   530     72        1 (1834/ 913-1877/ 974)  56.3 W54-W60 S21-S18
      6 1633  863 (S17W35 L=163)   679    165        1 (1586/ 796-1687/ 922)  29.1 W30-W39 S20-S14
      7 1711  643 (S35W41 L=169)   353     41        1 (1686/ 624-1733/ 660)   0.0 W37-W43 S36-S33
      8  735  435 (S22E34 L= 95)  3504    797        4 ( 537/ 298- 939/ 637)  18.0 E53-E13 S31-S13
      9  953  465 (S25E18 L=110)   495    173        0 ( 892/ 416-1005/ 497) -66.8 E23-E16 S29-S22
     10 1082  439 (S30E11 L=117)   307     80        1 (1040/ 420-1123/ 466)   8.1 E14-E07 S30-S28
  

    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