Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1688 1407 (N08W53 L=314)   983    189        2 (1653/1334-1714/1481)  45.0 W49-W58 N03-N11
      2 1622 1190 (S02W40 L=301)   254     80        0 (1598/1168-1658/1225) -45.0 W38-W42 S03-N00
      3  506  931 (N10E33 L=228)   243     77        0 ( 476/ 899- 533/ 950) -53.1 E35-E31 N07-N11
      4  658  818 (N00E26 L=235)   539    121        0 ( 584/ 794- 719/ 841) -29.7 E29-E22 S01-N03
      5  278  774 (N06E55 L=206)  1018    156        2 ( 233/ 738- 334/ 799)  -5.7 E59-E49 N03-N07
      6 1340  737 (S20W11 L=273)   208     81        0 (1307/ 713-1372/ 778)  21.8 W09-W14 S21-S18
      7 1890  744 (S37W58 L=319)   815    196        0 (1867/ 729-1926/ 770) -29.1 W54-W63 S39-S34
      8  484  455 (S17E51 L=210)  1554    283        7 ( 416/ 372- 577/ 526) -69.0 E55-E47 S24-S11
      9  396  438 (S16E60 L=201)  2401    318        7 ( 342/ 364- 465/ 504) -79.7 E65-E55 S22-S11
     10  721  313 (S32E42 L=219)   300     78        0 ( 707/ 284- 737/ 340)  36.9 E43-E38 S33-S29
  

    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