Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1119 1770 (N47W39 L=117)   299     64        0 (1091/1755-1153/1784) -31.0 W36-W41 N45-N49
      2  919 1704 (N50W19 L= 97)   319     89        1 ( 896/1674- 959/1728)   7.1 W15-W23 N48-N51
      3 1671 1636 (N19W71 L=150)  1358    144        4 (1635/1572-1708/1677)  63.4 W67-W80 N14-N22
      4  446 1603 (N60E35 L= 43)   489     47        1 ( 420/1576- 476/1629) -20.0 E41-E30 N57-N63
      5  504 1358 (N39E27 L= 52)   459    138        0 ( 477/1298- 534/1423) -73.3 E28-E24 N34-N44
      6  824 1278 (N26E05 L= 74)   473    186        1 ( 738/1258- 906/1296) -37.9 E09-W00 N22-N29
      7 1901 1175 (S12W68 L=146)  1034     88        2 (1859/1150-1945/1192) -12.5 W61-W79 S14-S10
      8  597 1040 (N18E25 L= 54)   890    343        3 ( 528/ 909- 685/1130) -72.6 E28-E22 N08-N24
      9  507 1065 (N21E31 L= 48)   393    128        1 ( 482/1012- 533/1116) -90.0 E32-E28 N18-N24
     10  326  997 (N21E47 L= 32)   417    128        0 ( 286/ 956- 374/1036) -40.6 E50-E43 N18-N24
     11  521  872 (N10E34 L= 45)   663    128        2 ( 497/ 802- 553/ 927)  60.3 E36-E30 N06-N13
     12 1564  778 (S23W26 L=104)   886    324        0 (1442/ 685-1660/ 893)  14.0 W15-W35 S25-S18
     13 1811  765 (S33W49 L=127)   808    111        0 (1764/ 724-1844/ 808)  12.5 W43-W53 S34-S31
     14 1251  630 (S22W02 L= 81)   371    200        1 (1194/ 571-1316/ 720)  31.0 E02-W08 S24-S18
     15 1707  647 (S37W37 L=115)   339     40        2 (1690/ 624-1727/ 668)  33.7 W35-W39 S39-S35
     16  777  437 (S21E33 L= 46)   577    114        0 ( 738/ 397- 827/ 465)   7.1 E36-E28 S23-S19
     17  683  278 (S30E50 L= 28)   540    133        0 ( 635/ 251- 716/ 331)  15.9 E57-E43 S31-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