Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1261 1742 (N46W09 L=260)  1352    203        2 (1172/1689-1351/1774) -11.3 W01-W16 N42-N48
      2  363 1531 (N20E56 L=196)  1263    164        3 ( 340/1449- 402/1609) -45.0 E62-E50 N14-N25
      3  830 1535 (N24E19 L=233)   231     73        0 ( 805/1512- 862/1558)  56.3 E20-E17 N22-N25
      4 1035 1431 (N19E04 L=248)  1056    302        2 ( 920/1328-1147/1531) -33.7 E12-W03 N13-N24
      5 1529 1395 (N23W29 L=281)   996    272        5 (1437/1333-1648/1471) -17.4 W22-W38 N20-N27
      6 1739 1324 (N23W47 L=299)  1504    297        3 (1699/1188-1775/1459)  86.8 W44-W51 N13-N31
      7 1842 1203 (N17W57 L=309)  1951    375        4 (1785/1036-1914/1349) -59.5 W51-W67 N07-N26
      8  361  730 (S29E46 L=206)   624    165        0 ( 322/ 667- 420/ 789)  33.7 E51-E39 S33-S25
      9 1822  665 (S15W67 L=318)  2950    250        9 (1713/ 583-1887/ 726)   2.0 W52-W81 S19-S10
     10 1193  608 (S29W17 L=268)  1281    398        2 (1044/ 488-1352/ 716)  51.7 W08-W24 S38-S19
     11  866  461 (S43E04 L=247)   232     50        0 ( 850/ 433- 882/ 493)  90.0 E05-E03 S45-S40
  

    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