stringtranslate.com

Слепая деконволюция

В электротехнике и прикладной математике слепая деконволюция — это деконволюция без явного знания функции импульсного отклика , используемой в свертке . Обычно это достигается путем принятия соответствующих предположений о входных данных для оценки импульсного отклика путем анализа выходных данных. Слепая деконволюция неразрешима без принятия предположений о входных данных и импульсном отклике. Большинство алгоритмов для решения этой задачи основаны на предположении, что и входные данные, и импульсный отклик находятся в соответствующих известных подпространствах. Однако слепая деконволюция остается очень сложной невыпуклой задачей оптимизации даже с этим предположением.

Иллюстрация слепой деконволюции
Верхнее левое изображение: NGC224 с космического телескопа Хаббл . Верхний правый контур: наилучшее соответствие функции рассеяния точки (PSF) (априори). [1] Среднее левое изображение: деконволюция по максимальной апостериорной оценке (MAP), 2-я итерация. Средний правый контур: оценка PSF по MAP, 2-я итерация. Нижнее левое изображение: деконволюция по MAP, окончательный результат. Нижний правый контур: оценка PSF по MAP, окончательный результат.

В обработке изображений

В обработке изображений слепая деконволюция — это метод деконволюции, который позволяет восстановить целевую сцену из одного или набора «размытых» изображений при наличии плохо определенной или неизвестной функции рассеяния точки (ФРТ). [2] Обычные линейные и нелинейные методы деконволюции используют известную ФРТ. Для слепой деконволюции ФРТ оценивается по изображению или набору изображений, что позволяет выполнить деконволюцию. Исследователи изучали методы слепой деконволюции в течение нескольких десятилетий и подходили к проблеме с разных сторон.

Большая часть работ по слепой деконволюции началась в начале 1970-х годов. Слепая деконволюция используется в астрономической и медицинской визуализации.

Слепая деконволюция может быть выполнена итеративно, при этом каждая итерация улучшает оценку PSF и сцены, или неитеративно, когда одно применение алгоритма, основанное на внешней информации, извлекает PSF. Итеративные методы включают максимальную апостериорную оценку и алгоритмы максимизации ожидания . Хорошая оценка PSF полезна для более быстрой сходимости, но не обязательна.

Примерами неитеративных методов являются SeDDaRA, [3] преобразование кепстра и APEX. Методы преобразования кепстра и APEX предполагают, что PSF имеет определенную форму, и необходимо оценить ширину формы. Для SeDDaRA информация о сцене предоставляется в виде опорного изображения. Алгоритм оценивает PSF, сравнивая информацию о пространственной частоте в размытом изображении с информацией о целевом изображении.

Примеры

Любое размытое изображение может быть передано в качестве входных данных для алгоритма слепой деконволюции, он может устранить размытость изображения, но основное условие для работы этого алгоритма не должно быть нарушено, как обсуждалось выше. В первом примере (изображение фигур) восстановленное изображение было очень четким, в точности похожим на исходное изображение, потому что L > K + N. Во втором примере (изображение девушки) L < K + N, поэтому основное условие нарушено, следовательно, восстановленное изображение сильно отличается от исходного изображения.

Размытое изображение, полученное сверткой исходного изображения с ядром размытия. Входное изображение лежит в фиксированном подпространстве вейвлет-преобразования, а ядро ​​размытия лежит в случайном подпространстве.

В обработке сигналов

Сейсмические данные

В случае деконволюции сейсмических данных исходный неизвестный сигнал состоит из пиков, поэтому его можно охарактеризовать с помощью ограничений разреженности [4] или регуляризаций, таких как отношения l 1 norm / l 2 norm norm [5] , предложенные WC Gray в 1978 году. [6]

Аудио деконволюция

Аудио деконволюция (часто называемая дереверберацией ) — это уменьшение реверберации в аудиосмесях. Это часть аудиообработки записей в некорректных случаях, таких как эффект коктейльной вечеринки . Одной из возможностей является использование ICA . [7]

В общем

Предположим, что у нас есть сигнал, передаваемый по каналу. Канал обычно можно смоделировать как линейную инвариантную к сдвигу систему , поэтому рецептор получает свертку исходного сигнала с импульсным откликом канала. Если мы хотим обратить эффект канала, чтобы получить исходный сигнал, мы должны обработать полученный сигнал второй линейной системой, инвертируя отклик канала. Эта система называется эквалайзером .

Восстановленное изображение после применения алгоритма слепой деконволюции. Этот алгоритм в основном решает задачу оптимизации с использованием минимизации ядерной нормы. L=65536, K=65 и N=44838,

Если нам дан исходный сигнал, мы можем использовать метод наблюдения, например, нахождение фильтра Винера , но без него мы все равно можем исследовать то, что мы знаем о нем, чтобы попытаться восстановить его. Например, мы можем отфильтровать полученный сигнал, чтобы получить желаемую спектральную плотность мощности . Это то, что происходит, например, когда известно, что исходный сигнал не имеет автокорреляции , и мы « отбеливаем » полученный сигнал.

Отбеливание обычно оставляет некоторые фазовые искажения в результатах. Большинство методов слепой деконволюции используют статистику более высокого порядка сигналов и позволяют корректировать такие фазовые искажения. Мы можем оптимизировать эквалайзер, чтобы получить сигнал с PSF, приближающимся к тому, что мы знаем об исходном PSF.

Исходное изображение
Размытое изображение: получено после свертки исходного изображения с ядром размытия. Исходное изображение лежит в фиксированном подпространстве вейвлет-преобразования, а размытие — в случайном подпространстве. L=65536, K=200, N=65400
Восстановленное изображение: сильно отличается от исходного изображения, поскольку нарушено основное условие для алгоритма слепой деконволюции с использованием минимизации ядерной нормы. L=65536, K=200, N=65400

Статистика высшего порядка

Алгоритмы слепой деконволюции часто используют статистику высокого порядка с моментами выше двух. Это может быть неявным или явным. [8]

Смотрите также

Внешние ссылки

Ссылки

  1. ^ Бармби, Полин; Маклафлин, Дин Э.; Харрис, Уильям Э.; Харрис, Гретхен Л. Х.; Форбс, Дункан А. (2007). «Структурные параметры шаровых скоплений в M31 и обобщения для фундаментальной плоскости» (PDF) . The Astronomical Journal . 133 (6): 2764–2786. arXiv : 0704.2057 . Bibcode :2007AJ....133.2764B. doi :10.1086/516777. S2CID  58913061.
  2. ^ Лам, Эдмунд Й.; Гудман, Джозеф В. (2000). «Итеративный статистический подход к слепой деконволюции изображений». Журнал оптического общества Америки A. 17 ( 7): 1177–1184. Bibcode : 2000JOSAA..17.1177L. doi : 10.1364/JOSAA.17.001177. PMID  10883969.
  3. ^ Caron, James N.; Namazi, Nader M.; Rollins, Chris J. (2002). «Безытеративное слепое восстановление данных с использованием извлеченной функции фильтра». Applied Optics . 41 (32): 6884–9. Bibcode : 2002ApOpt..41.6884C. doi : 10.1364/AO.41.006884. PMID  12440543.
  4. ^ L. Tonellot, T.; K. Broadhead, M. (2010). «Разреженная сейсмическая деконволюция методом ортогонального сопоставления». 72-я конференция и выставка EAGE, включающая SPE EUROPEC 2010. Европейская ассоциация геологов и инженеров. doi :10.3997/2214-4609.201401250.
  5. ^ Barmby, P.; McLaughlin, DE; Harris, WE; Harris, GLH; Forbes, DA (2015). «Евклид в такси: разреженная слепая деконволюция со сглаженной регуляризацией l1/l2». IEEE Signal Processing Letters . 22 (5): 539–543. arXiv : 1407.5465 . Bibcode : 2015ISPL...22..539R. doi : 10.1109/LSP.2014.2362861. S2CID  9605797.
  6. ^ Gray, WC (1978). Деконволюция переменной нормы (PDF) (диссертация доктора философии). Стэнфордский университет. Архивировано из оригинала (PDF) 2015-04-09.
  7. ^ Колдовский, Збынек; Тихавский, Петр (2007). «Слепое разделение источников звука во временной области с использованием передовых методов ICA». Труды 8-й ежегодной конференции Международной ассоциации речевой коммуникации (Interspeech 2007) . стр. 846–849.
  8. ^ Кардосо, Ж.-Ф. (1991). «Суперсимметричное разложение тензора кумулянта четвертого порядка. Слепая идентификация большего количества источников, чем датчиков». [Труды] ICASSP 91: 1991 Международная конференция по акустике, речи и обработке сигналов . Том 5. С. 3109–3112. CiteSeerX 10.1.1.8.9380 . doi :10.1109/ICASSP.1991.150113. ISBN  978-0-7803-0003-3. S2CID  7972548.