Фильтр реконструкции
В системе смешанных сигналов ( аналоговых и цифровых ) фильтр реконструкции , иногда называемый антивизуальным фильтром , используется для формирования сглаженного аналогового сигнала с цифрового входа, как в случае цифро-аналогового преобразователя ( ЦАП ). или другое устройство вывода выборочных данных.
Фильтры восстановления выборочных данных
[ редактировать ]Теорема выборки описывает, почему на входе АЦП требуется нижних частот аналоговый электронный фильтр , называемый фильтром сглаживания : дискретизированный входной сигнал должен быть ограничен по полосе пропускания , чтобы предотвратить наложение спектров (здесь это означает, что волны более высокой частоты записываются как более низкие частоты). .
По той же причине на выходе ЦАП требуется аналоговый фильтр нижних частот, называемый фильтром реконструкции, поскольку выходной сигнал должен быть ограничен по полосе пропускания, чтобы предотвратить формирование изображения (это означает, что коэффициенты Фурье восстанавливаются как ложные высокочастотные «зеркала»). Это реализация интерполяционной формулы Уиттекера-Шеннона .
В идеале оба фильтра должны быть фильтрами кирпичной стены , постоянной задержкой фазы в полосе пропускания с постоянной ровной частотной характеристикой и нулевой характеристикой частоты Найквиста . Этого можно достичь с помощью фильтра с импульсной характеристикой sinc .
Выполнение
[ редактировать ]Хотя теоретически ЦАП выдает серию дискретных импульсов Дирака , на практике реальный ЦАП выдает импульсы с конечной полосой пропускания и шириной. И идеализированные импульсы Дирака, и удерживаемые ступени нулевого порядка , и другие выходные импульсы, если их не фильтровать, будут содержать ложные высокочастотные копии « или изображения » исходного сигнала с ограниченной полосой пропускания. Таким образом, фильтр реконструкции сглаживает форму сигнала, чтобы удалить частоты изображения (копии) выше предела Найквиста . При этом он восстанавливает непрерывный временной сигнал (независимо от того, был ли он первоначально дискретизирован или смоделирован цифровой логикой), соответствующий цифровой временной последовательности.
Практические фильтры имеют неровную частотную или фазовую характеристику в полосе пропускания и неполное подавление сигнала в других местах. Идеальная форма сигнала sinc имеет бесконечную реакцию на сигнал как в положительном, так и в отрицательном направлении времени, что невозможно реализовать в реальном времени, поскольку для этого потребуется бесконечная задержка. Следовательно, реальные фильтры восстановления обычно либо допускают некоторую энергию выше скорости Найквиста, либо ослабляют некоторые внутриполосные частоты, либо и то, и другое. По этой причине можно использовать передискретизацию , чтобы гарантировать точное воспроизведение интересующих частот без выброса избыточной энергии за пределы полосы пропускания.
В системах, которые имеют оба, фильтр сглаживания и фильтр реконструкции могут иметь идентичную конструкцию. Например, как вход, так и выход аудиооборудования могут дискретизироваться с частотой 44,1 кГц. В этом случае оба аудиофильтра максимально блокируют частоты выше 22 кГц и максимально пропускают частоты ниже 20 кГц.
Альтернативно, система может не иметь восстанавливающего фильтра и просто допускать трату некоторой энергии на воспроизведение более высокочастотных изображений спектра первичного сигнала.
Обработка изображений
[ редактировать ]При обработке изображений фильтры цифровой реконструкции используются как для воссоздания изображений из образцов, так и в медицинской визуализации. [1] и для повторной выборки . [2] Был проведен ряд сравнений по различным критериям; [1] [2] [3] [4] одно наблюдение состоит в том, что реконструкцию можно улучшить, если производная сигнала, помимо амплитуды также известна [3] и наоборот, выполнение реконструкции по производной может улучшить методы реконструкции сигнала. [1]
Повторную выборку можно назвать прореживанием или интерполяцией , соответственно, когда частота дискретизации уменьшается или увеличивается - как и при выборке и реконструкции, в обоих случаях обычно применяются одни и те же критерии, и, следовательно, можно использовать один и тот же фильтр.
В принципе, при повторной выборке аналоговое изображение восстанавливается, а затем дискретизируется, и это необходимо для общего изменения разрешения. Для целочисленных отношений частоты дискретизации можно упростить выборку импульсной характеристики фильтра непрерывной реконструкции для создания дискретного фильтра передискретизации, а затем использовать дискретный фильтр передискретизации для прямой передискретизации изображения. Для прореживания на целое число необходим только один фильтр выборки; для интерполяции целым числом для разных фаз необходимы разные выборки – например, если выполняется повышающая дискретизация в 4 раза, то один дискретный фильтр используется для средней точки, а другой дискретный фильтр используется для точка 1/4 пути от одной точки до другой.
Тонкость обработки изображений заключается в том, что (линейная) обработка сигнала предполагает линейную яркость – удвоение значения пикселя удваивает яркость выходного сигнала. Однако изображения часто кодируются гамма-гаммой , особенно в цветовом пространстве sRGB , поэтому яркость не является линейной.Таким образом, чтобы применить линейный фильтр, необходимо сначала гамма-декодировать значения, а в случае повторной выборки необходимо гамма-декодировать, повторно дискретизировать, а затем гамма-кодировать.
Общие фильтры
[ редактировать ]Наиболее распространенными повседневными фильтрами являются: [5]
- интерполяция ближайшего соседа , с ядром и коробочным фильтром – для понижения дискретизации, что соответствует усреднению;
- билинейная интерполяция с ядром палаточного фильтра;
- бикубическая интерполяция с ядром в виде кубического сплайна – последний имеет свободный параметр, причем каждое значение параметра дает другой фильтр интерполяции.
Они расположены в порядке возрастания подавления полосы задерживания (сглаживания) и уменьшения скорости.
Для целей реконструкции используются различные ядра, многие из которых можно интерпретировать как аппроксимацию функции sinc: [4] либо с помощью окон, либо путем аппроксимации сплайнами, либо кубами, либо сплайнами более высокого порядка. В случае оконных sinc-фильтров частотную характеристику восстанавливающего фильтра можно понимать как частотную характеристику окна, поскольку частотная характеристика оконного фильтра представляет собой свертку исходной характеристики (для sinc — кирпичного фильтра). стены) с частотной характеристикой окна. Среди них окно Ланцоша и окно Кайзера часто хвалят .
Другой класс фильтров реконструкции включает фильтры Гаусса для различной ширины: [2] или кардинальные B-сплайны более высокого порядка - коробчатый фильтр и палаточный фильтр представляют собой кардинальные B-сплайны 0-го и 1-го порядка. Эти фильтры не являются интерполирующими фильтрами, поскольку их импульсная характеристика не исчезает во всех ненулевых исходных точках выборки - при повторной выборке 1:1 они не тождественны, а скорее размыты. С другой стороны, будучи неотрицательными, они не вносят каких-либо артефактов перерегулирования или звона , а будучи шире во временной области, они могут быть уже в частотной области (по принципу неопределенности Фурье ), хотя и за счет размытия, которое отражается в спаде полосы пропускания («зубце»).
В фотографии существует множество интерполяционных фильтров. [6] некоторые проприетарные, мнения по которым неоднозначны. Оценка часто бывает субъективной, реакции различаются, и некоторые утверждают, что при реалистичных коэффициентах повторной выборки разница между ними незначительна по сравнению с бикубической. [7] хотя для более высоких коэффициентов передискретизации поведение более разнообразно. [8]
Фильтры реконструкции вейвлета
[ редактировать ]Фильтры реконструкции также используются при «реконструкции» формы сигнала или изображения из набора вейвлет- коэффициентов.В медицинской визуализации обычным методом является использование нескольких двумерных рентгеновских фотографий или МРТ-сканирований для «реконструкции» трехмерного изображения.
См. также
[ редактировать ]Ссылки
[ редактировать ]- ^ Jump up to: а б с Тойсль, Томас; Хаузер, Хельвиг; Греллер, Мейстер Эдуард (октябрь 2000 г.). Освоение Windows: улучшение реконструкции (PDF) . Симпозиум IEEE/ACM SIGGRAPH по объемной визуализации . Солт-Лейк-Сити, Юта, США. стр. 101–108. дои : 10.1109/ВВ.2000.10002 . ISBN 1-58113-308-1 . ( Веб-страница проекта )
- ^ Jump up to: а б с Турковски, Кен (1990). «Фильтры для распространенных задач повторной выборки» (PDF) .
- ^ Jump up to: а б Митчелл, Дон П.; Нетравали, Арун Н. (август 1988 г.). Реконструкция фильтров в компьютерной графике (PDF) . Международная конференция ACM SIGGRAPH по компьютерной графике и интерактивным технологиям . Том. 22. С. 221–228. дои : 10.1145/54852.378514 . ISBN 0-89791-275-6 .
- ^ Jump up to: а б Мейеринг, Эрик Х.В.; Ниссен; Плуим; Виргевер. Количественное сравнение Sinc-аппроксимирующих ядер для интерполяции медицинских изображений . Вычисление медицинских изображений и компьютерное вмешательство - MICCAI '99: вторая международная конференция, Кембридж, Великобритания, 19–22 сентября 1999 г., протоколы .
- ^ dpreview: Интерполяция , Винсент Бокарт
- ^ Обзор интерполяции цифровых фотографий
- ^ Интерполяция - Часть I , Рон Бигелоу
- ^ Фильтр изображений - Сепия