Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1371 1797 (N39W55 L=152)  1206    149        2 (1304/1761-1414/1824) -24.0 W50-W60 N36-N41
      2  746 1304 (N30E09 L= 87)   584    182        1 ( 675/1251- 799/1368) -53.1 E12-E05 N26-N35
      3 1750 1078 (S11W47 L=143)  3711    522        5 (1706/ 852-1797/1319)  58.0 W39-W54 S23-N01
      4 1167 1094 (N07W10 L=106)   598    228        2 (1051/1044-1273/1151) -42.3 W04-W15 N02-N12
      5  704 1041 (N16E18 L= 78)   237     83        0 ( 664/1024- 746/1059) -31.0 E20-E15 N14-N17
      6  226  999 (N23E56 L= 41)  1628    320        4 ( 143/ 956- 311/1033) -20.6 E64-E48 N19-N26
      7 1134  617 (S19E05 L= 92)   201     65        0 (1114/ 588-1167/ 638)  26.6 E06-E02 S20-S18
      8  748  509 (S16E31 L= 65)  1097    306        2 ( 713/ 384- 791/ 629)  60.6 E35-E24 S24-S08
      9  518  402 (S18E52 L= 44)   352     99        0 ( 484/ 373- 569/ 436) -56.3 E54-E49 S20-S14
  

    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