2-D FIR-фильтр методом окна 2-D
fwind2 конструирование функций 2-D фильтров FIR с использованием метода окна. fwind2 использует спецификацию окна 2-D для проектирования фильтра 2-D FIR на основе требуемой частотной характеристики. fwind2 работает только с окнами 2-D. Использовать fwind1 создать 2-й фильтр ЕЛИ из 1-D окна.
Вы можете применить 2-й фильтр ЕЛИ к изображениям при помощи filter2 функция.
В этом примере показано, как сконструировать приблизительно круговой симметричный двухмерный полосовой фильтр с использованием метода 2-D окна.
Создание векторов частотного диапазона f1 и f2 использование freqspace. Эти векторы имеют длину 21.
[f1,f2] = freqspace(21,'meshgrid');Вычислите расстояние каждой позиции от центральной частоты.
r = sqrt(f1.^2 + f2.^2);
Создание матрицы Hd который содержит требуемый полосовой отклик. В этом примере желаемая полоса пропускания находится в диапазоне от 0,1 до 0,5 (нормированная частота, где 1,0 соответствует половине частоты дискретизации, или радиан).
Hd = ones(21); Hd((r<0.1)|(r>0.5)) = 0;
Отображение идеального полосового отклика.
colormap(parula(64)) mesh(f1,f2,Hd)
Создайте 2-е Гауссовское использование окна fspecial. Нормализуйте окно.
win = fspecial('gaussian',21,2);
win = win ./ max(win(:));Откройте окно.
mesh(win)
С помощью окна 2-D спроектируйте фильтр, который наилучшим образом обеспечивает требуемую частотную характеристику.
h = fwind2(Hd,win);
Отображение фактической частотной характеристики этого фильтра.
freqz2(h)
fwind2 вычисляет h использование обратного преобразования Фурье и умножение на двумерное окно win.
ej
(n1, n2)
[1] Lim, Jae S., двумерная обработка сигналов и изображений, Englewood Cliffs, NJ, Prentice Hall, 1990, стр. 202-213.