Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  853 1812 (N49E38 L=254)   478     72        0 ( 825/1782- 888/1829)  78.7 E41-E35 N46-N51
      2  646 1718 (N37E48 L=245)   763    131        2 ( 618/1661- 688/1776) -22.4 E57-E40 N33-N41
      3  957 1690 (N41E20 L=272)   283     28        0 ( 934/1667- 978/1710) -26.6 E22-E18 N39-N42
      4  440 1641 (N28E62 L=231)   481     44        2 ( 415/1624- 464/1661)   0.0 E65-E57 N26-N29
      5 1559 1559 (N44W34 L=327)   330     39        2 (1536/1537-1580/1576)  38.7 W31-W36 N42-N46
      6 1384 1444 (N32W17 L=310)   243     66        0 (1362/1421-1414/1467) -26.6 W15-W19 N30-N33
      7 1060 1370 (N21E04 L=289)   303    101        0 (1022/1344-1109/1400) -18.4 E06-W00 N19-N22
      8  267 1221 (S01E57 L=236)   989    160        2 ( 216/1172- 317/1295) -19.7 E65-E50 S03-N02
      9 1772  652 (S10W63 L=356)   575     85        0 (1726/ 631-1810/ 679)  26.6 W58-W68 S11-S07
     10  916  569 (S31W01 L=294)   296     86        1 ( 873/ 551- 957/ 588)  36.9 E01-W04 S32-S28
     11  970  379 (S44W10 L=303)   396    106        0 ( 933/ 340-1001/ 414)  74.1 W08-W12 S47-S40
     12 1561  345 (S33W64 L=357)   606     54        2 (1540/ 328-1587/ 368) -71.6 W60-W69 S34-S31
  

    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