| Image Processing Toolbox User's Guide | ![]() |
Design two-dimensional FIR filter using frequency transformation
Syntax
Description
h = ftrans2(b,t) produces the two-dimensional FIR filter h that corresponds to the one-dimensional FIR filter b using the transform t. (ftrans2 returns h as a computational molecule, which is the appropriate form to use with filter2.) b must be a one-dimensional, odd-length (Type I) FIR filter such
as can be returned by fir1, fir2, or remez in the Signal Processing Toolbox. The transform matrix t contains coefficients that define the frequency
transformation to use. If t is m-by-n and b has length Q, then h is size ((m-1)*(Q-1)/2+1)-by-((n-1)*(Q-1)/2+1).
h = ftrans2(b) uses the McClellan transform matrix t.
Remarks
The transformation below defines the frequency response of the two-dimensional filter returned by ftrans2,
where B(
) is the Fourier transform of the one-dimensional filter b,
and T(
1,
2) is the Fourier transform of the transformation matrix t.
The returned filter h is the inverse Fourier transform of H(
1,
2).
Example
Use ftrans2 to design an approximately circularly symmetric two-dimensional bandpass filter with passband between 0.1 and 0.6 (normalized frequency, where 1.0 corresponds to half the sampling frequency, or
radians):
ftrans2 transforms a one-dimensional FIR filter to create a two-dimensional filter, first design a one-dimensional FIR bandpass filter using the Signal Processing Toolbox function remez.
colormap(jet(64)) b = remez(10,[0 0.05 0.15 0.55 0.65 1],[0 0 1 1 0 0]); [H,w] = freqz(b,1,128,'whole'); plot(w/pi-1,fftshift(abs(H)))
ftrans2 with the default McClellan transformation to create the desired approximately circularly symmetric filter.
See Also
conv2, filter2, fsamp2, fwind1, fwind2
Reference
[1] Lim, Jae S., Two-Dimensional Signal and Image Processing, Englewood Cliffs, NJ, Prentice Hall, 1990, pp. 218-237.
| fspecial | fwind1 | ![]() |
© 1994-2005 The MathWorks, Inc.