Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  942 1323 (N18E09 L=171)   659    289        2 ( 907/1198- 976/1453) -76.0 E11-E04 N10-N26
      2 1798 1390 (N33W57 L=237)   701    237        1 (1742/1367-1841/1413)   0.0 W50-W62 N31-N34
      3 1866 1156 (N18W62 L=242)   454    101        1 (1845/1115-1886/1196) -68.2 W59-W65 N15-N20
      4 1756  932 (N03W52 L=232)   401     81        1 (1726/ 905-1778/ 954) -21.8 W49-W54 N01-N04
      5  491  811 (S18E33 L=147)   362    181        1 ( 446/ 751- 534/ 870)  53.1 E36-E29 S22-S14
      6 1846  664 (S12W71 L=251)  2560    249        4 (1812/ 586-1890/ 755)  70.0 W66-W76 S17-S06
      7 1756  707 (S11W57 L=237)   676    160        2 (1695/ 665-1787/ 752)  50.2 W52-W60 S13-S07
      8  202  695 (S30E62 L=118)  1170    248        2 ( 184/ 639- 229/ 732)  45.0 E65-E57 S34-S27
      9 1755  587 (S18W62 L=242)  1438    189        2 (1702/ 520-1783/ 651)  84.3 W58-W66 S23-S13
     10 1523  492 (S27W44 L=224)   382     50        1 (1500/ 468-1551/ 513)  71.6 W41-W45 S28-S25
     11 1430  274 (S45W50 L=229)   696     88        1 (1383/ 262-1478/ 288)  23.2 W45-W53 S46-S43
  

    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