%Speckle Interferometry code for Averaging Power Spectrum
%copyright 2010 Arun Surya
%This program is free software: you can redistribute it and/or modify it %under the terms of the GNU General Public License as published by the Free %Software Foundation, either version 1 of the License. 
%This program is distributed in the hope that it will be useful, but WITHOUT %ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or %FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License for %more details.
%You should have received a copy of the GNU General Public License along
%with this program.  If not, see <http://www.gnu.org/licenses/>.



;This program usesSpeckle Interferometry technique  to reconstruct Object Auto correlation
;

;*****************************************************************
m = 200.0
n = 200.0
image=fltarr(m,n)
number = 100.;the no of frames
amp = dblarr(m,n)
ampi = dblarr(m,n)
fimagel = dcomplexarr(m,n)
bispm = dcomplexarr(3e7)
;*****************************************************************
;Calculations of  Powerspectrum

;*****************************************************************
print ,systime(0);
for jack=1 ,number,1 do begin
image(23,45) = 100
image(27,50) = 70
image(32,40) = 70
jacky=STRTRIM(jack, 2)
image=readfits('C:\imagerot'+jacky+'.fits')
;image=readfits('C:\image.fits')
image3 = image ;+ randomn(s,140,140)/10
fimage = shift(fft(shift(image3,m/2,n/2),1,double=3),m/2,n/2)
amp = amp+(abs(fimage)*abs(fimage));-float(fimage(m/2,n/2)))
print,jack
endfor


amp = amp/(number)

;reim = sqrt(amp)*ephase
reim = sqrt(amp)

imag = shift(fft(shift(reim,m/2,n/2),-1,double=3),m/2,n/2)
;writefits,'C:\bst\amp.fits',amp
writefits,'C:\bst\object.fits',image3
writefits,'C:\bst\image.fits',abs(imag)
print ,systime(0);
end