Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  422 1602 (N38E54 L=238)   498    119        2 ( 397/1578- 455/1625) -14.0 E57-E49 N36-N39
      2  199 1287 (N16E66 L=225)   656    152        1 ( 169/1266- 226/1311)  15.9 E71-E63 N14-N17
      3  663 1226 (N14E24 L=268)   818    243        3 ( 564/1149- 754/1302)  33.7 E30-E18 N09-N19
      4  519 1187 (N11E34 L=258)   285     77       -1 ( 498/1164- 539/1214) -71.6 E35-E32 N10-N13
      5 1772  777 (S12W56 L=347)  3465    522        4 (1645/ 559-1874/ 980)  63.4 W47-W66 S26-N00
      6  155  914 (S08E69 L=223)  2015    238        5 ( 121/ 845- 179/ 979)  71.6 E76-E65 S12-S03
      7  546  892 (S07E31 L=261)   302    102        0 ( 508/ 867- 595/ 910)  -9.5 E33-E27 S08-S05
      8  278  815 (S13E54 L=238)   652    155        5 ( 253/ 758- 304/ 863) -71.6 E56-E51 S16-S10
      9  170  809 (S14E68 L=224)   759    188        0 ( 142/ 782- 193/ 849)  32.0 E73-E63 S16-S11
     10  862  649 (S21E10 L=282)  1235    325        3 ( 747/ 559- 936/ 778)  51.3 E18-E04 S27-S12
     11  587  584 (S27E30 L=262)   253    112        0 ( 573/ 538- 603/ 631)  66.8 E31-E28 S30-S23
     12  396  420 (S41E54 L=238)  1330    302        1 ( 350/ 391- 441/ 455)  23.6 E62-E46 S44-S37
  

    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