Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  708 1790 (N33E48 L=215)  2634    346        7 ( 597/1718- 854/1863)  49.4 E53-E39 N25-N40
      2 1007 1811 (N43E27 L=235)   461     72        1 ( 982/1779-1036/1842) -63.4 E30-E25 N39-N45
      3  790 1689 (N27E34 L=229)   962    184        3 ( 728/1658- 865/1727)  35.5 E39-E29 N24-N30
      4 1540 1440 (N31W22 L=285)  1895    486        0 (1394/1268-1715/1612) -23.4 W08-W38 N25-N39
      5  403 1570 (N12E57 L=205)   479    104        1 ( 366/1554- 455/1592)   6.3 E61-E52 N10-N13
      6  474 1535 (N11E49 L=214)   480     90        2 ( 435/1513- 509/1560)  53.1 E51-E46 N08-N12
      7 1782 1337 (N33W46 L=309)   603     99        0 (1729/1306-1815/1365)   0.0 W40-W50 N31-N34
      8 1184 1166 (N04W05 L=268)   372    108        1 (1127/1139-1227/1200)  50.2 W02-W07 N01-N07
      9  500 1131 (S12E34 L=229)   595    146        0 ( 448/1082- 552/1175)  49.4 E37-E31 S16-S09
     10 1912 1115 (N24W64 L=327)   540    109        0 (1881/1061-1942/1159) -15.3 W58-W70 N21-N25
     11  149  998 (S26E68 L=195)  2296    152        4 ( 119/ 920- 187/1072) -66.0 E73-E61 S30-S21
     12  322  927 (S28E46 L=217)   781    254        2 ( 224/ 884- 412/ 988)  22.8 E57-E38 S31-S22
     13 1175  887 (S11W12 L=275)   413    100        0 (1108/ 867-1222/ 903)  26.6 W08-W14 S13-S09
     14  541  790 (S32E26 L=237)  2296    486        6 ( 357/ 662- 746/ 866)  31.5 E43-E11 S43-S24
     15 1593  518 (S19W53 L=316)   749    166        0 (1525/ 475-1671/ 566)  52.1 W49-W57 S23-S14
     16 1299  530 (S28W32 L=295)   514    162        1 (1233/ 514-1380/ 547)  16.7 W27-W37 S30-S25
     17 1486  361 (S31W57 L=320)  2873    290        4 (1337/ 331-1610/ 400)  23.5 W44-W68 S36-S24
  

    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