Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1158 1798 (N50W13 L=122)   386     71        0 (1134/1762-1188/1830)  49.4 W10-W16 N46-N53
      2 1302 1808 (N52W28 L=137)   365     27        0 (1289/1781-1317/1829)  39.8 W25-W31 N48-N54
      3  654 1709 (N43E32 L= 78)   571    104        1 ( 611/1683- 698/1728) -16.7 E36-E26 N40-N44
      4 1650 1504 (N28W47 L=156)   506     62        1 (1621/1468-1671/1531)  45.0 W44-W50 N25-N29
      5 1753  718 (S21W54 L=163)   419     40        0 (1721/ 698-1780/ 739)  21.8 W50-W56 S21-S19
      6  234  639 (S25E65 L= 44)   455     58        0 ( 202/ 627- 267/ 653)   0.0 E71-E60 S25-S24
      7 1765  523 (S33W66 L=176)   818     77        1 (1737/ 480-1783/ 568) -19.7 W61-W75 S35-S29
      8  519  352 (S46E50 L= 59)   364     40        0 ( 503/ 326- 536/ 376)  26.6 E54-E45 S48-S44
  

    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