Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1240 1814 (N54E15 L=244)   590     79        0 (1201/1799-1287/1829)  14.0 E19-E10 N51-N56
      2  866 1795 (N38E40 L=219)   822    126        1 ( 810/1772- 939/1818)  15.3 E45-E33 N36-N39
      3  504 1703 (N22E61 L=198)   485    113        0 ( 482/1668- 536/1748) -54.5 E63-E56 N19-N26
      4  292 1526 (N07E70 L=189)  1912    169        5 ( 237/1467- 355/1577)  45.0 E77-E65 N02-N10
      5 1314 1351 (N22W08 L=267)   488    146        0 (1274/1316-1372/1405) -23.2 W04-W11 N19-N24
      6  213 1293 (S08E66 L=193)   643     77        2 ( 176/1271- 249/1318)  26.6 E70-E61 S10-S06
      7 1846  881 (N11W61 L=319)  1097    163        2 (1805/ 816-1872/ 961) -77.5 W57-W63 N06-N16
      8 1508  595 (S15W43 L=302)   277     65        1 (1480/ 567-1531/ 629)  90.0 W42-W44 S16-S12
      9  738  559 (S40E05 L=254)   371    110        1 ( 710/ 513- 781/ 607) -29.7 E07-W00 S41-S37
     10 1400  515 (S22W40 L=298)   363    113        0 (1358/ 498-1452/ 540)  26.6 W36-W42 S23-S19
  

    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