Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  982 1903 (N63W49 L=199)   519     64        0 ( 962/1872-1003/1935)   8.5 W40-W61 N60-N65
      2  421 1457 (N47E35 L=116)   410     23        1 ( 399/1430- 447/1476)  26.6 E37-E31 N44-N48
      3  997 1316 (N22W06 L=157)   232     88        0 ( 957/1290-1043/1332) -11.3 W03-W08 N20-N22
      4  330 1113 (N26E44 L=106)   642    196        0 ( 291/1037- 403/1156)  16.7 E48-E37 N22-N28
      5 1685  905 (S20W38 L=189)   423    127        0 (1644/ 856-1733/ 939)  14.0 W34-W42 S21-S18
      6  279  839 (N11E54 L= 97)   334     72        1 ( 262/ 806- 293/ 886)  53.1 E56-E52 N09-N13
      7 1617  682 (S32W29 L=180)   528    153        0 (1556/ 642-1677/ 739)  10.3 W23-W35 S32-S29
      8 1223  582 (S26E00 L=151)   234     69        1 (1194/ 540-1251/ 610)  45.0 E02-W02 S27-S24
      9  826  561 (S17E24 L=127)   273    102        0 ( 784/ 533- 867/ 597)   8.1 E27-E20 S17-S15
     10  431  513 (S10E54 L= 97)   860    162        0 ( 408/ 434- 464/ 582)  52.1 E57-E48 S15-S06
     11 1401  485 (S38W09 L=160)   687    152        0 (1335/ 446-1473/ 531)   4.8 W03-W15 S39-S36
     12  571  284 (S28E59 L= 91)  1828    210        4 ( 494/ 238- 649/ 332) -11.3 E67-E48 S32-S23
  

    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