Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  403 1621 (N22E64 L=199)   745     48        0 ( 371/1594- 452/1655)  -3.8 E72-E57 N20-N23
      2  358 1478 (N12E58 L=205)  1505    174        7 ( 321/1417- 414/1543) -66.0 E64-E52 N07-N16
      3  388 1380 (N07E51 L=212)   450    114        0 ( 359/1332- 417/1430)  81.9 E51-E48 N03-N10
      4 1609 1318 (N29W34 L=297)   979    135        3 (1561/1259-1647/1383) -35.5 W29-W37 N25-N33
      5 1868  925 (N11W64 L=327)   573     70        3 (1842/ 894-1889/ 958)  45.0 W61-W67 N08-N12
      6  533  841 (S23E28 L=235)   427    159        0 ( 484/ 796- 569/ 912)  58.0 E31-E25 S26-S17
      7  263  774 (S32E54 L=209)   394     67        0 ( 245/ 750- 286/ 810) -45.0 E56-E50 S33-S30
      8 1027  703 (S21W07 L=270)   521    128        2 ( 970/ 672-1085/ 741)  39.8 W03-W09 S23-S17
      9  724  581 (S35E10 L=253)   984    302        2 ( 617/ 545- 839/ 627)   6.7 E19-E02 S37-S31
     10  278  576 (S47E60 L=203)   955    115        2 ( 223/ 554- 333/ 600)   7.1 E82-E50 S49-S44
     11 1645  444 (S22W64 L=327)  2134    254        8 (1565/ 370-1718/ 540)  85.6 W59-W71 S28-S15
     12 1549  479 (S23W51 L=314)   619    109        2 (1518/ 433-1580/ 521)  90.0 W49-W53 S25-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