Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1060 1553 (N29W01 L=301)   445    152        1 ( 978/1533-1143/1566)  -5.7 E04-W06 N28-N30
      2 1807 1497 (N30W61 L=  1)   646    125        0 (1777/1467-1827/1520) -33.7 W57-W63 N27-N31
      3  745 1425 (N20E19 L=281)   342     51        1 ( 710/1403- 775/1449)  36.9 E21-E16 N18-N21
      4 1213 1365 (N17W11 L=311)   539    134        2 (1177/1289-1241/1430) -77.5 W08-W12 N12-N21
      5  994 1308 (N13E03 L=297)   997    236        0 ( 875/1244-1109/1388) -29.7 E10-W04 N09-N17
      6  462 1274 (N10E37 L=263)   241     73        0 ( 442/1245- 485/1298)  90.0 E38-E35 N08-N12
      7  783 1215 (N07E15 L=285)   189     78        0 ( 742/1201- 813/1246)  26.6 E17-E13 N06-N08
      8 1717 1069 (N02W46 L=346)   263     59        0 (1688/1052-1742/1090)  26.6 W43-W47 N00-N02
      9 1870  920 (S06W62 L=  2)  1817    302        4 (1827/ 847-1919/ 981)  69.4 W57-W69 S10-S02
     10  192  859 (S14E61 L=239)  1118    181        4 ( 161/ 790- 225/ 922)  41.2 E65-E56 S17-S09
     11 1683  874 (S10W44 L=344)   368     59        0 (1652/ 843-1708/ 912)  68.2 W41-W46 S12-S07
     12 1316  825 (S15W19 L=319)   354    149        0 (1286/ 759-1364/ 884)  74.1 W16-W21 S18-S11
     13  220  763 (S20E59 L=241)  2179    274        4 ( 176/ 679- 290/ 851)  42.3 E64-E50 S24-S14
     14  273  769 (S20E53 L=247)   360     72        1 ( 254/ 747- 295/ 803) -71.6 E55-E51 S20-S17
     15  990  659 (S26E01 L=299)   876    161        3 ( 914/ 598-1040/ 760)  65.6 E06-W01 S30-S19
     16 1843  628 (S24W66 L=  6)   775     56        0 (1819/ 599-1868/ 654)  63.4 W63-W68 S25-S21
     17  275  548 (S33E59 L=241)  1187    239        1 ( 258/ 471- 302/ 594)  48.8 E63-E55 S38-S30
  

    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