Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  879 1658 (N50W02 L=243)   287     47        1 ( 863/1628- 897/1682)  53.1 W00-W04 N48-N52
      2 1475 1476 (N25W41 L=281)   951    330        0 (1321/1423-1601/1527) -31.4 W30-W48 N19-N30
      3  335 1094 (N20E47 L=193)   543     89        1 ( 297/1070- 383/1119) -18.4 E51-E42 N17-N21
      4 1553  589 (S31W29 L=270)  3712    575        6 (1337/ 447-1728/ 763)  20.0 W11-W45 S37-S21
      5  451  637 (S09E45 L=196)   625    106        2 ( 417/ 596- 488/ 675)  38.7 E47-E41 S11-S06
      6  290  588 (S11E64 L=177)   451     66        2 ( 275/ 557- 306/ 620)  63.4 E65-E61 S13-S09
      7  435  446 (S22E56 L=185)   500     68        1 ( 415/ 419- 462/ 466) -56.3 E60-E53 S23-S20
      8  702  346 (S32E38 L=203)   685    103        0 ( 675/ 301- 725/ 395)  35.0 E43-E33 S35-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