Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  596 1166 (N15E27 L=200)   423    146        4 ( 524/1137- 658/1198) -18.4 E32-E23 N13-N17
      2  903  906 (S02E08 L=219)   269    132        0 ( 850/ 863- 961/ 932) -29.7 E11-E04 S04-N00
      3  237  650 (S17E67 L=161)   572     72        0 ( 206/ 617- 259/ 677) -26.6 E71-E64 S19-S15
      4 1714  628 (S26W51 L=279)   507     50        1 (1697/ 600-1736/ 661) -68.2 W48-W54 S28-S23
      5  478  634 (S18E41 L=187)   445     66        0 ( 448/ 607- 517/ 656)  18.4 E43-E37 S19-S16
      6  449  575 (S22E45 L=183)   350     60        1 ( 428/ 550- 478/ 599)  11.3 E47-E42 S23-S20
      7  413  518 (S26E51 L=176)   700     79        1 ( 385/ 490- 456/ 547) -31.0 E55-E47 S28-S24
      8 1478  459 (S36W33 L=260)   348     62        1 (1444/ 442-1511/ 479)   9.5 W29-W35 S37-S34
      9  424  425 (S33E56 L=171)   664    101        1 ( 407/ 392- 463/ 467) -53.1 E61-E52 S36-S30
  

    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