Q-сдвиговые фильтры Kingsbury
[ возвращает фильтры Кингсбери Q-сдвига для Q-сдвига комплексного преобразования двойного дерева. Целое число LoDa,LoDb,HiDa,HiDb,LoRa,LoRb,HiRa,HiRb] = qorthwavf(num)num относится к числу ненулевых коэффициентов (отводов) в фильтре. Допустимые параметры для num 6, 10, 14, 16 и 18. Все фильтры имеют четную длину, а фильтры дерева B представляют собой обратное время фильтров дерева A.
Получите фильтры Q-сдвига для случая с 10 ненулевыми коэффициентами.
[LoDa,LoDb,HiDa,HiDb,LoRa,LoRb,HiRa,HiRb] = qorthwavf(10);
Используйте dwtfilterbank создание двух дискретных блоков фильтров вейвлет-преобразования. Используйте фильтры анализа дерева A в первом блоке фильтров, а фильтры анализа дерева B во втором блоке фильтров.
fbTreeA = dwtfilterbank('Wavelet','Custom',... 'CustomScalingFilter',LoDa,... 'CustomWaveletFilter',HiDa); fbTreeB = dwtfilterbank('Wavelet','Custom',... 'CustomScalingFilter',LoDb,... 'CustomWaveletFilter',HiDb);
Постройте график наиболее масштабных вейвлетов каждого набора фильтров.
[psiA,t] = wavelets(fbTreeA); [psiB,~] = wavelets(fbTreeB); plot(t,psiA(end,:)) hold on plot(t,psiB(end,:)) grid on hold off legend('Tree A','Tree B')
Убедитесь, что оба блока фильтров ортогональны.
isOrthogonal(fbTreeA)
ans = logical
1
isOrthogonal(fbTreeB)
ans = logical
1
[1] Антонини, М., М. Барло, П. Матье и И. Даубехиес. «Кодирование изображения с помощью вейвлет-преобразования». Транзакции IEEE по обработке изображений 1, № 2 (апрель 1992 года): 205-20. https://doi.org/10.1109/83.136597.
[2] Кингсбери, Ник. «Сложные вейвлеты для инвариантного анализа сдвига и фильтрации сигналов». Прикладный и вычислительный гармонический анализ 10, № 3 (май 2001 года): 234-53. https://doi.org/10.1006/acha.2000.0343.
[3] Ле Галль, Д. и А. Табатабай. «Поддиапазонное кодирование цифровых изображений с использованием симметричных фильтров короткого ядра и методов арифметического кодирования». В ICASSP-88., Международная конференция по акустике, речи и обработке сигналов, 761-64. Нью-Йорк, Нью-Йорк, США: IEEE, 1988. https://doi.org/10.1109/ICASSP.1988.196696.