Wavelet Transformation in IDL

PRO WAVELETTRANSFORM


READ_GIF,'d:\Image\bird.gif',f
sz=size(f,/dimensions)
DEVICE,DECOMPOSED = 0
LOADCT, 0
WINDOW,0,XSIZE=SZ[0],YSIZE=SZ[1],TITLE='Original Image'
tvscl, congrid(f,sz[0],sz[1])
f1 = wtn(f,20)
f2=abs(f1)^2
f2 = alog10(f2)
WINDOW,1,XSIZE=SZ[0],YSIZE=SZ[1],TITLE='Transformed Image'
tvscl, congrid(f2,sz[0],sz[1])
f3 = wtn(f1,20,/inverse)
WINDOW,2,XSIZE=SZ[0],YSIZE=SZ[1],TITLE='inverse Transformed Image'
tvscl, congrid(f3,sz[0],sz[1])

END