IDL Version 5.2 (linux x86). Research Systems, Inc.
Installation number: 82096-3.
Licensed for use by: University of Utrecht

For basic information, enter "IDLInfo" at the IDL> prompt.

% Compiled module: SKEY_LINUX.
use of Function-Keys:

F4  : Erase window        F8  : window 0 front & wset,0
F2  : Window front        F5  : -> print,
F3  : Iconify window      F6  : -> plot,
F9  : New free window     F7  : Toggle windows

IDL> .run spotmask
IDL> .run rubb
IDL> help,/rou
IDL> rubb, mov, cell=128, step=96, plot=4, /same
IDL> rubb, mov, cell=64, step=50, plot=4, /same
IDL> mov(*, *, 54) = (mov(*, *, 53)+mov(*, *, 55))/2
IDL> tvscl,mov(*,*,54)
IDL> tot=total(mov,3)/118.
IDL> tvscl,tot
IDL> tt=avg(tot)
IDL> print,tt
     0.841799
IDL> print,avg(mov(*,*,0))
     0.853297
IDL> print,avg(mov)
     0.721313
IDL> help,mov
MOV             FLOAT     = Array[600, 600, 118]
IDL> FOR i=0, num-1 DO mov(*, *, i) = mov(*, *, i) * (tt/avg(mov(*, *, i)))
IDL> for i=0,117,2 do tvscl,mov(*,*,i)
IDL> enhance, mov
IDL> mm1=alog(mov)        
IDL> mm1=scale_array(mm1)
IDL> write_mpeg,'spot.mpg',mm1
IDL> 
IDL> rubb, mm1, cell=32, step=25, plot=4, range=2, /same
IDL> rubb, mm1, cell=16, step=10, plot=4, range=2, /same
IDL> k=replicate(1b,5,5)
IDL> mask=tot lt .8
IDL> mask=dilate(erode(mask,k),k)
IDL> m1=erode(dilate(mask,k),k)
IDL> tvscl,tot
IDL> tmp=get_part(tot,/aspect,pos=pos)
IDL> print,pos
          76          66         552         542
IDL> help,tmp
TMP             FLOAT     = Array[477, 477]
IDL> print,477/2
     238
IDL> print,477/2+[66,76]
     304     314
IDL> kreis,314,304,238,238,/dev
IDL> tmp=polarint(shift(tot,314,304))
IDL> plot,tmp
IDL> plots,238,!y.crange
IDL> xt=findgen(300)
IDL> tvscl,tot
IDL> plots,xt+314,304+300*(tmp-tmp(0)),/dev
IDL> plots,314-xt,304+300*(tmp-tmp(0)),/dev
IDL> plots,314-xt,304+10*(tmp-tmp(0)),/dev
IDL> plots,314-xt,304+50*(tmp-tmp(0)),/dev
IDL> xx=shift(dist(628,608),314,304)
IDL> xx=xx(28:*,8:*)
IDL> 
IDL> plot,xx,tmp,psy=3
IDL> plots,238,!y.crange
IDL> oplot,xt,tmp,col=255
IDL> plot,xx(0:313,*),tot(0:313,*),psy=3
IDL> tmp1=hist_2d(xx(0:313,*),mov(0:313,*,0),bin2=0.01,max2=2.)
IDL> tvscl,tmp1
IDL> for i=1,117 do tmp1=tmp1+hist_2d(xx(0:313,*),mov(0:313,*,i),bin2=0.01,max2=2.)
IDL> tvscl,tmp1
IDL> tvscl,alog(tmp1+1)
IDL> tvscl,alog(tmp1+10)
IDL> tmp2=fltarr(417)
IDL> for i=0,416 do tmp2(i)=com(reform(tmp1(i,*)))
IDL> plots,238,!y.crange
IDL> plots,238,[0,200],/dev

IDL> tvscl,alog(tmp1+10)
IDL> plots,indgen(417),tmp2,/dev,col=255

IDL> tvscl,-alog(tmp1+10)
IDL> plots,indgen(417),tmp2,/dev,col=color('yellow')

IDL> tvimg,-alog(tmp1+10),findgen(417)/238,findgen(201)/100.
% Compiled module: TVIMG.
% Compiled module: RESCALE.
IDL> oplot,findgen(417)/238,tmp2/100.