Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1273 1653 (N39E05 L=340)   355     67        1 (1247/1627-1295/1685) -45.0 E06-E02 N37-N41
      2  360 1625 (N12E68 L=277)   653     83        0 ( 334/1596- 389/1659)  68.2 E69-E65 N10-N15
      3  776 1628 (N22E34 L=311)   257     45        1 ( 748/1609- 803/1643)  11.3 E36-E31 N20-N23
      4  416 1132 (S15E40 L=305)   344    114        2 ( 398/1070- 430/1187) -56.3 E42-E37 S18-S12
      5 1330 1083 (N06W15 L=360)   413    162        0 (1263/1044-1370/1153)  71.6 W12-W16 N01-N10
      6 1787 1083 (N21W47 L= 32)   377    110        0 (1750/1037-1811/1109)  -9.5 W43-W49 N18-N21
      7 1888 1045 (N22W58 L= 43)   886    167        1 (1851/ 997-1925/1095) -11.3 W53-W64 N20-N24
      8  389 1011 (S23E40 L=305)   645    177        1 ( 333/ 991- 464/1053)   0.0 E44-E33 S24-S21
      9 1148  499 (S31W25 L= 10)  1345    389        3 ( 981/ 423-1315/ 581)   6.3 W11-W39 S33-S28
  

    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