Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1638 1550 (N21W58 L=133)  2327    400        2 (1489/1515-1745/1600) -32.3 W47-W70 N15-N28
      2  327 1436 (N47E48 L= 28)   857    252        2 ( 268/1402- 372/1477) -28.3 E55-E42 N43-N50
      3 1324 1303 (N16W25 L=100)   491    161        0 (1249/1273-1362/1357)   6.3 W19-W28 N14-N18
      4  350 1139 (N27E43 L= 32)   598     96        0 ( 311/1106- 391/1180) -39.8 E46-E40 N24-N29
      5  954  646 (S13E13 L= 63)   758    139        2 ( 913/ 581- 994/ 716)  54.5 E16-E09 S16-S09
      6  344  559 (S08E59 L= 17)  1838    230        2 ( 273/ 500- 432/ 621) -19.7 E66-E51 S12-S04
      7 1745  547 (S44W45 L=121)   695    151        1 (1709/ 531-1785/ 567) -28.6 W40-W51 S47-S41
      8  676  303 (S31E45 L= 31)   470    102        1 ( 641/ 279- 698/ 336)  15.9 E48-E41 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