Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1181 1673 (N36E10 L= 99)   621    133        1 (1149/1603-1225/1724) -84.3 E11-E07 N30-N40
      2 1218 1221 (N09W05 L=115)   378    126        2 (1174/1187-1289/1252)  56.3 W03-W08 N06-N12
      3  140 1142 (S19E70 L= 40)  1216    114        1 (  78/1119- 200/1165)   4.8 E86-E61 S20-S17
      4 1847 1101 (N22W53 L=163)   537    114        3 (1806/1049-1884/1144) -16.7 W48-W58 N20-N23
      5  769 1038 (S13E15 L= 95)   273     95        0 ( 733/1005- 795/1087) -31.0 E17-E12 S14-S10
      6 1824  981 (N14W53 L=162)   439     48        1 (1801/ 962-1851/1007)  11.3 W50-W56 N13-N15
      7 1840  918 (N12W56 L=166)   421     47        1 (1818/ 889-1859/ 944) -45.0 W53-W58 N10-N13
      8 1872  847 (N09W63 L=172)   551     74        2 (1845/ 813-1903/ 873)  76.0 W59-W67 N07-N11
      9 1630  476 (S18W59 L=169)   458     62        1 (1609/ 450-1662/ 506) -56.3 W56-W62 S19-S16
     10 1513  303 (S30W65 L=175)   603     42        0 (1488/ 282-1535/ 330) -63.4 W61-W69 S32-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