Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  403 1545 (N30E57 L=144)  2038    214        3 ( 342/1487- 480/1622) -13.6 E74-E45 N26-N34
      2  258  894 (S12E55 L=146)   568    114        0 ( 230/ 850- 285/ 948)  49.4 E58-E51 S15-S08
      3  199  833 (S17E64 L=137)   488     95        1 ( 181/ 794- 214/ 874) -68.2 E65-E61 S19-S14
      4  444  758 (S19E38 L=163)   371    138        0 ( 384/ 733- 480/ 800)  20.6 E43-E35 S20-S16
      5  737  774 (S16E17 L=185)   240     42        1 ( 716/ 755- 763/ 793)  33.7 E17-E14 S16-S14
      6 1223  724 (S15W15 L=216)   504    188        1 (1146/ 669-1291/ 752) -11.3 W09-W19 S18-S13
      7 1652  625 (S19W50 L=251)  1846    319        5 (1571/ 541-1726/ 685)  39.3 W43-W57 S25-S15
      8  760  383 (S44E16 L=185)   348     62        0 ( 732/ 360- 785/ 416)  50.2 E19-E13 S46-S40
      9  519  349 (S50E45 L=156)   511     75        1 ( 486/ 328- 554/ 376) -15.9 E49-E41 S51-S47
  

    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