Wavelet denoising for tomographically reconstructed image

We have developed a wavelet denoising (thresholding) method for a tomographically reconstructed image to which the conventional wavelet methods are not necessarily applicable because of their limitation of applicable noise models. The basic idea of our new method is that noise variance is. in general. spatially varying and the threshold must be adapted to it. Specifically, our algorithm includes two key steps: The first is to estimate local variances in image space to produce a "sigma-map". The second is to calculate the standard deviations of individual wavelet coefficients from the a-map by a formula of "covariance propagation". Spatially adaptive thresholds are then Liven as those proportional to the standard deviations. Our method is applicable to a wider range of noise models. and numerical experiments have shown that it can yield a denoised image with 1070 less residual error than that in the boxcar smoothing or the median filtering. (c) 2006 The Optical Society of Japan