Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  871 1728 (N31E32 L=140)  1513    360        2 ( 787/1618- 960/1787)  90.0 E37-E23 N21-N36
      2 1279 1674 (N40E03 L=169)   639    156        2 (1218/1621-1322/1725) -21.8 E09-W01 N36-N42
      3  466 1480 (N06E48 L=124)   829    169        3 ( 388/1452- 523/1532)  41.2 E53-E43 N03-N10
      4 1792 1415 (N41W48 L=220)   579     98        1 (1753/1376-1818/1447) -23.2 W43-W51 N38-N42
      5  300  892 (S31E48 L=124)  1064    332        1 ( 222/ 782- 364/1036)  52.1 E58-E42 S39-S21
      6 1695  964 (N08W42 L=214)   285     75        1 (1674/ 928-1720/ 996) -21.8 W40-W45 N06-N09
      7  820  885 (S20E08 L=164)   458    177        2 ( 747/ 815- 872/ 928) -16.7 E13-E03 S22-S18
      8  443  767 (S36E33 L=139)   689    225        1 ( 405/ 682- 502/ 858) -51.3 E37-E27 S40-S30
      9 1575  736 (S08W41 L=213)   411     75        3 (1556/ 699-1596/ 773) -76.0 W38-W43 S09-S05
     10 1138  667 (S24W17 L=189)   198     44        1 (1119/ 644-1157/ 695) -26.6 W14-W18 S25-S22
     11 1416  459 (S27W44 L=216)   994    283        0 (1287/ 426-1532/ 493)  31.0 W34-W51 S30-S21
     12 1098  344 (S45W31 L=203)   414     97        0 (1075/ 304-1129/ 383) -60.3 W29-W34 S47-S40
     13  661  313 (S64E05 L=167)  1497    160        5 ( 593/ 272- 732/ 365)   0.0 E15-W05 S68-S60
     14 1564  339 (S28W67 L=239)  1448    177        4 (1480/ 314-1643/ 363)  14.6 W57-W81 S31-S24
     15 1334  288 (S39W53 L=225)   340     67        1 (1314/ 268-1362/ 317)  -9.5 W49-W56 S40-S37
     16 1260  245 (S45W53 L=225)  1818    182        3 (1191/ 209-1350/ 290)  12.5 W43-W62 S47-S41
  

    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