Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  621 1725 (N64E23 L= 76)   377     57        0 ( 594/1706- 654/1737)   0.0 E28-E17 N62-N66
      2  782 1411 (N35E08 L= 91)   421    144        0 ( 717/1369- 844/1436) -36.9 E12-E04 N31-N37
      3 1831 1137 (S06W59 L=158)   480     88        0 (1815/1100-1849/1177)  68.2 W57-W61 S08-S03
      4 1853 1041 (S12W60 L=158)  1676    307        4 (1779/ 912-1930/1142)  17.7 W49-W75 S17-S07
      5  755 1095 (N16E15 L= 84)   245    111        0 ( 718/1066- 801/1133) -56.3 E16-E12 N13-N19
      6  194 1043 (N21E63 L= 35)   476     61        1 ( 165/1015- 213/1073) -23.2 E67-E60 N19-N22
      7 1315  890 (S07W14 L=113)   306     73        1 (1298/ 862-1345/ 928)  76.0 W13-W16 S08-S04
      8 1736  857 (S19W44 L=143)   453    124        0 (1682/ 839-1772/ 897)  -6.3 W39-W48 S20-S17
      9 1427  725 (S19W19 L=117)   541    223        0 (1378/ 648-1448/ 823)  55.0 W14-W21 S23-S13
     10 1755  554 (S42W51 L=150)   720    100        2 (1713/ 525-1791/ 593)  -4.8 W45-W57 S44-S39
     11  580  466 (S19E43 L= 56)   518    167        0 ( 523/ 442- 635/ 502) -20.6 E46-E37 S21-S16
     12  693  259 (S37E48 L= 51)  1804    214        1 ( 657/ 161- 743/ 351)  25.3 E72-E36 S47-S29
  

    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