Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  904 1789 (N61W16 L=287)   617    113        2 ( 855/1760- 957/1813)   0.0 W09-W24 N58-N63
      2 1392 1719 (N40W52 L=324)   518     86        1 (1356/1689-1424/1761)   8.1 W46-W60 N38-N41
      3  452 1627 (N58E46 L=225)   428     45        1 ( 429/1609- 480/1643)  15.9 E50-E42 N56-N59
      4 1128 1593 (N39W22 L=294)  1041    331        1 (1041/1537-1285/1635) -11.9 W14-W33 N34-N43
      5 1736 1524 (N18W69 L=341)  1224    107        1 (1698/1488-1788/1548)  -2.7 W62-W83 N16-N19
      6 1631  685 (S27W34 L=306)   457    107        0 (1613/ 631-1644/ 739)  82.9 W32-W35 S31-S23
      7 1443  688 (S22W20 L=291)   247     76        0 (1422/ 655-1468/ 726)  36.9 W17-W21 S24-S20
  

    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