Decompose image into subbands (undecimated wavelet), denoise, and recompose again. fh = decomp_reconst_wavelet(im,Nsc,daub_order,block,noise,parent,covariance,optim,sig) im : image Nsc: Number of scales daub_order: Order of the daubec