Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1712 1464 (N27W56 L= 20)   674    227        0 (1642/1448-1770/1479)  -8.7 W49-W62 N25-N28
      2 1305 1331 (N22W20 L=344)  2113    646        6 ( 992/1213-1504/1434) -26.6 W00-W32 N13-N30
      3 1756 1280 (N15W55 L= 19)   490     95        1 (1733/1237-1777/1331)  80.5 W53-W56 N12-N18
      4  706 1183 (N15E20 L=304)   519    200        2 ( 666/1087- 747/1280) -71.6 E22-E17 N09-N21
      5 1793  892 (S09W55 L= 19)  1852    399        4 (1702/ 787-1867/ 963) -35.8 W46-W64 S16-S03
      6 1452  820 (S11W27 L=350)   358     87        2 (1425/ 774-1480/ 873)  63.4 W24-W28 S13-S07
      7 1558  588 (S27W36 L=  0)  2359    419        3 (1392/ 479-1685/ 753)  40.4 W24-W45 S33-S16
      8 1252  626 (S22W13 L=337)   462    138        1 (1185/ 579-1326/ 675)  29.1 W08-W18 S24-S18
  

    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