Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  693 1734 (N27E46 L= 30)   579     82        0 ( 674/1699- 711/1789) -33.7 E51-E42 N23-N30
      2 1104 1614 (N30E13 L= 63)   676    278        0 (1022/1492-1184/1690) -26.6 E21-E03 N23-N32
      3  718 1631 (N20E37 L= 39)   472    119        2 ( 690/1601- 759/1659)  11.3 E39-E34 N18-N22
      4  819 1559 (N18E27 L= 48)   274     86        1 ( 793/1523- 843/1599) -33.7 E30-E24 N16-N20
      5  155 1050 (S25E67 L=  9)  1820    180       10 ( 108/ 997- 225/1121) -15.5 E76-E56 S27-S20
      6 1497  586 (S17W43 L=118)   275    100        1 (1455/ 561-1548/ 611)  45.0 W39-W45 S19-S13
      7 1506  442 (S24W52 L=128)   431    101        0 (1463/ 407-1541/ 485)  90.0 W51-W53 S27-S20
      8 1642  455 (S18W63 L=139)   425     46        1 (1617/ 435-1673/ 470)  36.9 W59-W66 S19-S16
      9 1628  366 (S23W71 L=147)  1280     99        1 (1598/ 321-1662/ 406) -40.6 W67-W79 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