Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1233 1574 (N35W29 L=171)   405    115        1 (1181/1547-1289/1591) -29.7 W25-W32 N32-N36
      2 1023 1536 (N38W13 L=155)   223     71        1 (1005/1508-1046/1565)  36.9 W11-W15 N35-N39
      3  578 1384 (N38E24 L=119)   795    145        2 ( 525/1329- 638/1431) -45.0 E28-E19 N33-N41
      4  342 1151 (N27E45 L= 98)  1692    495        4 ( 287/ 953- 410/1307) -70.8 E50-E40 N14-N37
      5 1780 1180 (S02W54 L=197)   430     66        1 (1750/1158-1806/1202) -45.0 W51-W56 S04-N00
      6  867 1158 (N18E07 L=136)   461    175        0 ( 787/1129- 959/1201) -37.9 E11-E01 N14-N21
      7 1108  940 (N01W03 L=145)   932    257        3 (1023/ 825-1166/1036) -81.3 W00-W06 S06-N07
      8  160  979 (N18E65 L= 77)   510     76        1 ( 144/ 936- 177/1015) -78.7 E66-E63 N15-N20
      9 1721  937 (S14W43 L=186)   681    110        0 (1703/ 873-1745/ 990)  66.8 W41-W45 S18-S11
     10 1074  645 (S16E05 L=137)   738    149        1 (1003/ 589-1116/ 683) -54.5 E08-E02 S20-S12
     11  862  317 (S35E30 L=113)   980    174        1 ( 818/ 260- 898/ 374)  30.3 E36-E24 S38-S30
  

    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