Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1263 1809 (N52E07 L= 83)   621     68        2 (1231/1787-1302/1830)  51.3 E10-E03 N48-N55
      2 1364 1770 (N52W03 L= 94)   393     56        0 (1349/1735-1387/1793)  82.9 W01-W05 N47-N55
      3  828 1741 (N32E34 L= 56)   579    115        0 ( 787/1689- 864/1787) -26.6 E40-E28 N29-N35
      4 1432 1730 (N51W11 L=101)   368     38        0 (1413/1709-1453/1752)  59.0 W09-W12 N48-N53
      5  454 1644 (N18E58 L= 33)  1532    113        3 ( 409/1608- 520/1693)  18.4 E64-E50 N15-N21
      6  858 1656 (N27E27 L= 63)   248     37        1 ( 841/1636- 883/1674)  63.4 E29-E25 N25-N28
      7  258 1076 (S19E54 L= 37)   390     49        0 ( 241/1044- 279/1108)  63.4 E55-E51 S20-S16
      8 1569  364 (S29W64 L=155)   832    137        1 (1517/ 308-1620/ 404)  90.0 W62-W68 S33-S24
      9  990  222 (S59W33 L=124)  1411    201        1 ( 926/ 146-1053/ 302) -13.2 W19-W53 S62-S54
  

    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