pro contrast,init,numim,x,y,arr ;logarr ;logarr - array with seq. nos and time (made by readlog.pro) ;init - initial image no. ;numim - number of images to be processed ;x,y = lower left corner coordinates of the selected box. ;subimages 256x256 will be extracted and contrast computed ;arr - resulting array with image numbers and contrasts ;parametros de entrada a modificar con editor n=700 ;dimension en x m=500 ;dimension en y size=256 ;size of subimages st1='/home1/msobotka/lap93/grand/spot' st2='.dat' im1=intarr(n,m,/nozero) ;numim=n_elements(logarr(0,*)) arr=fltarr(2,numim) FOR i = 0,numim-1 DO BEGIN ;seq=strtrim(string(fix(logarr(0,i))),1) seq=strtrim(i+init,1) print,'Now reading image:' print,st1+seq+st2 openr,1,st1+seq+st2 readu,1,im1 close,1 dev=stdev(im1(x:x+size-1,y:y+size-1),mea) arr(0,i)=i+init arr(1,i)=dev/mea print,seq,arr(1,i) ENDFOR ;plot,arr,yra=[min(arr)*0.95,max(arr)*1.05],/yst end