Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  410 1542 (N30E56 L=145)  1618    141        4 ( 347/1489- 480/1589) -12.8 E68-E46 N26-N32
      2  259  885 (S13E55 L=145)   528     87        0 ( 232/ 849- 287/ 924)  39.8 E58-E51 S15-S10
      3  451  759 (S19E38 L=163)   412    123        0 ( 396/ 733- 488/ 801)  26.6 E42-E34 S20-S16
      4  744  773 (S16E16 L=184)   253     53        1 ( 716/ 755- 773/ 788)  33.7 E18-E14 S16-S14
      5 1214  734 (S15W14 L=215)   355    123        1 (1153/ 710-1265/ 751)  -8.1 W10-W17 S15-S13
      6 1699  643 (S18W54 L=254)   977    126        2 (1650/ 599-1731/ 681)  23.2 W48-W57 S20-S15
      7 1620  617 (S20W47 L=247)   513    143        1 (1590/ 554-1653/ 680) -83.7 W45-W49 S24-S15
      8  768  381 (S44E16 L=185)   321     74        0 ( 737/ 358- 790/ 415)  50.2 E18-E13 S46-S40
      9  524  346 (S50E45 L=156)   492     79        0 ( 487/ 326- 559/ 368) -23.2 E48-E40 S51-S48
  

    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