Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1125 1718 (N38E11 L= 84)   997    155        2 (1060/1667-1188/1750) -10.3 E16-E05 N35-N40
      2  371 1645 (N19E65 L= 29)  1154     91        3 ( 333/1611- 411/1687) -71.6 E70-E60 N16-N22
      3  449 1622 (N18E54 L= 40)   343     69        1 ( 419/1593- 478/1649)  53.1 E56-E52 N15-N20
      4 1675 1452 (N35W36 L=130)  1221    269        3 (1549/1397-1771/1486)  12.8 W24-W46 N29-N38
      5  151 1129 (S15E67 L= 28)  4730    427       13 ( 107/ 953- 233/1263)  64.7 E76-E59 S25-S06
      6 1232 1188 (N06W08 L=103)   191     41        1 (1209/1170-1262/1202)  14.0 W06-W10 N05-N07
      7  463 1106 (S13E36 L= 58)   472    140        1 ( 395/1069- 520/1134)   0.0 E41-E31 S14-S11
      8 1899 1049 (N16W61 L=156)   947    152        1 (1843/1003-1940/1083)   0.0 W54-W69 N13-N18
      9 1075  909 (S12W05 L=100)   984    526        1 ( 845/ 829-1275/ 987)  34.3 E06-W16 S21-S05
     10  469  930 (S23E33 L= 61)   445    137        0 ( 399/ 910- 533/ 959)  21.8 E38-E28 S25-S20
     11 1933  844 (N07W73 L=168)  1238    186        2 (1906/ 739-1950/ 918) -73.3 W70-W76 N00-N10
     12 1324  531 (S28W32 L=127)   367     66        1 (1282/ 506-1353/ 553)  11.3 W28-W34 S29-S26
  

    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