Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1260 1576 (N40W15 L=338)   589    136        1 (1209/1533-1314/1599)  -6.3 W10-W19 N37-N41
      2 1585 1419 (N30W41 L=  4)   278     43        1 (1559/1403-1610/1439) -14.0 W38-W42 N29-N31
      3 1149 1283 (N19W07 L=330)   574    241        2 (1050/1215-1231/1333) -26.6 W00-W12 N15-N21
      4  701 1236 (N14E22 L=301)   635    179        2 ( 632/1199- 765/1303) -33.7 E27-E17 N11-N17
      5 1644 1042 (N06W42 L=  5)   445     70        1 (1619/1006-1668/1070)  53.1 W39-W43 N04-N08
      6  676  725 (S18E21 L=302)   494    142        1 ( 635/ 661- 712/ 789)  63.4 E24-E18 S22-S14
      7 1608  705 (S15W43 L=  6)   350     72        0 (1582/ 672-1629/ 737)  63.4 W40-W44 S16-S12
      8 1706  649 (S18W54 L= 17)   371     47        0 (1688/ 624-1723/ 673) -45.0 W51-W55 S19-S16
      9  873  641 (S23E08 L=315)   315    139        0 ( 812/ 624- 932/ 662)   0.0 E11-E03 S23-S21
     10  337  524 (S36E58 L=265)  1049    129        1 ( 282/ 505- 399/ 552)   0.0 E66-E49 S37-S33
  

    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