Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1489 1587 (N25W50 L=227)  1221    168        0 (1440/1540-1558/1619)   0.0 W44-W57 N22-N27
      2 1681 1491 (N14W60 L=237)  1488    272        0 (1585/1396-1773/1556) -61.9 W54-W64 N05-N20
      3 1607 1408 (N12W48 L=226)   316     84        0 (1557/1397-1646/1419) -23.2 W44-W51 N10-N13
      4 1189 1403 (N24W21 L=198)   277     95        0 (1140/1386-1232/1416) -36.9 W17-W23 N22-N26
      5  226 1032 (N24E58 L=119)  1308    279        1 ( 190/ 942- 290/1166) -60.3 E62-E51 N18-N32
      6 1824  993 (S18W53 L=231)   451     51        0 (1799/ 965-1856/1021)   9.5 W50-W57 S18-S16
      7 1523  941 (S11W27 L=205)   429    113        0 (1494/ 886-1543/ 983) -90.0 W26-W28 S14-S07
      8  197  805 (N11E66 L=111)  3288    408        3 ( 140/ 701- 266/ 896)  76.0 E76-E59 N04-N16
      9 1179  660 (S18E00 L=177)   215     58        0 (1165/ 627-1196/ 691)  63.4 E01-W01 S19-S15
     10 1812  633 (S42W54 L=231)   854    122        2 (1773/ 613-1847/ 659) -20.6 W46-W62 S44-S37
     11 1699  545 (S44W39 L=216)   517     39        1 (1677/ 524-1725/ 563)  11.3 W35-W42 S45-S41
  

    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