Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  825 1894 (N69E00 L=327)   881     82        1 ( 802/1852- 853/1929)  63.4 E04-W06 N62-N74
      2  681 1441 (N30E17 L=310)   766    140        4 ( 636/1398- 726/1484)  45.0 E20-E13 N27-N32
      3 1173 1378 (N19W14 L=342)   834    324        1 (1052/1254-1253/1460) -61.9 W08-W17 N10-N25
      4 1806 1417 (N13W65 L= 33)  1183    156        3 (1740/1368-1860/1453) -10.6 W57-W73 N09-N14
      5  313 1250 (N23E47 L=280)   557    111        1 ( 294/1191- 334/1304) -74.1 E48-E44 N19-N26
      6 1187 1085 (N01W10 L=338)   321     71        1 (1163/1051-1212/1122)  36.9 W08-W12 N00-N03
      7 1542  562 (S36W31 L=359)  2867    463        7 (1348/ 439-1718/ 697)  17.9 W15-W47 S42-S28
      8  435  548 (S20E49 L=278)  2098    315        2 ( 383/ 426- 527/ 681)  86.8 E52-E45 S29-S11
      9  576  368 (S34E47 L=280)  1038    172        2 ( 537/ 285- 611/ 424)  37.6 E56-E43 S40-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