Deblur image using Wiener filter

`J = deconvwnr(I,psf,nsr)`

`J = deconvwnr(I,psf,ncorr,icorr)`

deconvolves image `J`

= deconvwnr(`I`

,`psf`

,`nsr`

)`I`

using the Wiener filter algorithm,
returning deblurred image `J`

. `psf`

is the
point-spread function (PSF) with which `I`

was convolved.
`nsr`

is the noise-to-signal power ratio of the additive
noise. The algorithm is optimal in a sense of least mean square error between the
estimated and the true images.

The output image

`J`

could exhibit ringing introduced by the discrete Fourier transform used in the algorithm. To reduce the ringing, use`I = edgetaper(I,psf)`

before calling`deconvwnr`

.

[1] "Digital Image Processing", R. C. Gonzalez & R. E. Woods, Addison-Wesley Publishing Company, Inc., 1992.

`deconvblind`

| `deconvlucy`

| `deconvreg`

| `edgetaper`

| `otf2psf`

| `padarray`

| `psf2otf`