comp.lang.idl-pvwave archive: archive

Home » Public Forums » archive » Wiener filter

Show: Today's Messages :: Show Polls :: Message Navigator
E-mail to friend

Re: Wiener filter [message #28534 is a reply to message #28424]

Wed, 19 December 2001 12:10

Richard Tyc
Messages: 69
Registered: June 1999

Member

> Optimal Wiener filtering of a one-dimensional data set is described in
> section 12.6 of "Numerical Recipes in C", by Preuss et.al. It cites
> three books on signal processing as references. The basic result is that
> if you have a corrupted signal with the fourier spectrum S(f),
> containing noise with a fourier specturm N(f), it can be shown
> rigorously that the optimal (in the sense of a least-squares fit)
> frequency filter for removing the noise is:
>
> |S(f)|^2
> phi(f) = --------------------
> |S(f)|^2 + |N(f)|^2
>

What exactly is |S(f)|^2

If I have a 2D corrupted image, say I(x,y)

Is it ABS( FFT(I) ) ^ 2 or the magnitude of the complex FFT result
squared (Power Spectrum) squared ?

> The procedure is straightforward. Estimate the fourier spectrum of the
> noise. Calculate the fourier spectrum of the corrupted signal. Calculate
> the corresponding filter function. Multiply the fourier spectrum of the
> corrupted signal by the filter function. Do an inverse fourier transform
> on the resulting function, to get an optimum estimate to the uncorrupted
> signal.

It seems some knowledge of the noise is required. What if it was modeled as
'white noise' where it would be constant at all spatial frequencies.

A paper I am using that discusses this in the context of my problem points
out, "....Assuming that noise power spectrum is white, the mean spectral
density at high spatial frequencies was calculated and subtracted from P(f)
(the power spectral density of the corrupted image) to estimate S(f) (power
spectral density of uncorrupted image). Can you shed any light on this in
terms of IDL code ??

Thanks for the help
Rich

Report message to a moderator

[Message index]

		Wiener filter By: Richard Tyc on Mon, 17 December 2001 12:41
		Re: Wiener filter By: James Kuyper Jr. on Wed, 19 December 2001 13:48
		Re: Wiener filter By: James Kuyper Jr. on Wed, 19 December 2001 13:50
		Re: Wiener filter By: Richard Tyc on Fri, 21 December 2001 10:40
		Re: Wiener filter By: Richard Tyc on Wed, 19 December 2001 12:10
		Re: Wiener filter By: James Kuyper Jr. on Wed, 19 December 2001 07:03
		Re: Wiener filter By: jeyadev on Wed, 19 December 2001 13:34
		Re: Wiener filter By: jeyadev on Tue, 18 December 2001 12:30
		Re: Wiener filter By: James Kuyper Jr. on Fri, 21 December 2001 12:26

Previous Topic:	Re: FindFile for more than one filetype
Next Topic:	looking for some idl + seawifs (seadas) help

-=] Back to Top [=-

[ Syndicate this forum (XML) ] [

] [

]

Current Time: Wed Oct 08 14:56:05 PDT 2025

Total time taken to generate the page: 0.00486 seconds