im.makeImageFromSf(sfx=xrpix, sf=sfpix)
im.showImageI()
pylab.savefig('clouds_1.png', format='png')
im.image = im.imageI.real
#im.hanningFilter()
im.calcAll(min_npix=2, min_dr=1)
im.plotMore()
pylab.savefig('clouds_1_dat.png', format='png')
# Rescale sfx to be in physical (degrees)
im.sfx = im.sfx * pixscale
# make another image, as should see difference due to different random phases
im2 = PImagePlots(shift=True, nx= imsize, ny=imsize)
im2.invertSf(sfx=xrpix, sf=sfpix)
im2.invertAcovf1d()
im2.invertAcovf2d(useI=True)
im2.invertPsd2d(useI=True)
im2.invertFft(useI=True)
im2.showImageI()
pylab.savefig('clouds_2.png', format='png')
# Compare structure functions, but scaled (due to loss of amplitude info with random phases)
im2.image = im2.imageI.real
im2.calcAll(min_npix=2, min_dr=1)
im2.plotMore()
pylab.savefig('clouds_2_dat.png', format='png')
im2.sfx = im2.sfx * pixscale
im2.sf = im2.sf/im2.sf.max()
im.sf = im.sf/im.sf.max()
im.showSf(linear=True, comparison=im2, legendlabels=['Clouds 1 (scaled)', 'Clouds 2 (scaled)'])
