filt_btwl - Butterworth low-pass Fourier filter
output = filt_btwl(image, freql, freqh, pad)
- input image (can be either real or Fourier)
- low - pass-band frequency
- high - stop-band frequency
- logical flag specifying whether before filtering the image should be padded with zeroes in real space to twice the size (this helps avoiding aliasing artifacts). (Default pad = False).
All frequencies are in absolute frequency units `f_a` and their valid range is [0:0.5].
- filtered image. Output image is real when input image is real or Fourier when input image is Fourier
Fourier transform of the input image is multiplied by a radially symmetric Butterworth filter:
`B(f) = 1.0/sqrt(1.0+(f/(RAD))^(ORDER))`
Value of ORDER determines the filter falloff and RAD corresponds to the cut-off frequency. RAD and ORDER are calculated from the parameters specified by the user using following equations:
`ORDER = [2*log((eps)/sqrt(a^2-1))]/[log((low)/(high))]`
`RAD = (low)/((eps)^(2/(ORDER)))`
where low and high are the pass-band and stop-band frequencies, respectively, and parameters eps and a are set to 0.882 and 10.624, respectively. Note values of the filter B(f) deviate from 1.0 at frequency low by about 0.2 and from 0.0 at frequency high by about 0.09 (for a low-pass filter.)
Gonzalez, R. F., Woods, R. E., 2002. Digital Image Processing. Prentice Hall, Upper Saddle River, NJ.
Author / Maintainer
Pawel A. Penczek
- category 1
- category 2
- works for most people, has been tested; test cases/examples available.
None. It is perfect.