|
ru.algorithms
- RU.ALGORITHMS ----------------------------------------------------------------
From : Evgenij Masherov 2:5020/175.2 24 Jan 2002 11:01:46
To : Valentin Davydov
Subject : Re: Ищется алгоритм чистки радиокарты
--------------------------------------------------------------------------------
Thu Jan 24 2002 07:25, Valentin Davydov wrote to Evgenij Masherow:
>> Задача распадается на две части:
>> 1. Оценка "вслепую" апертурной функции.
>> 2. Деконволюция.
VD> Hеобязательно. Можно придумать некий итеративный алгоритм, постепенно
VD> "улучшающий" изображение и одновременно выдающий всё более точные
VD> приближения для аппаратной функции (типа, разложить её по ортогональным
VD> функциям, а потом деконволюировать карту с каждым членом разложения до
VD> "упора в ноль").
Можно. И даже сделать некую интерактивность. И останавливаться, когда получено
самое лучшее решение ("Господин коллекционер! Расчищая абстрактную картину,
которую Вы приказали, дабы обмануть таможню, нарисовать поверх купленной Вами
картины Леонардо да Винчи, нашел под последней портрет Муссолини! Hа каком
слое остановимся?"). Вопрос остановки и здесь самый интересный...
>> 2. Hаиболее эффективна в вычислительном отношении деконволюция через
>> Фурье. Она основана на теореме "Фурье от свертки равно произведению Фурье
>> исходных последовательностей". При деконволюции Фурье от свернутой
>> (наблюдаемой) карты делится на Фурье от апертурной функции.
>> Примечания:
>> а. В теореме фигурирует циклическая свертка. Поэтому нужно либо
>> использовать окна, убывающие к краям, либо дополнять нулями.
>> б. Делить на ноль интересно, но не очень полезно. Если верить астроному
>> Стругацкому - от этого на ушах шерсть растет... Поэтому если в
>> Фурье-образе нули - надо регуляризовать...
VD> Дык, Тихонов именно это и делает. Переписывает выражение F(A)/F(B) в
VD> виде F(A)*cnjg(F(B))/abs(F(B))^2, после чего имеет в знаменателе
VD> вещественную неотрицательную функцию, каковую можно сделать заведомо
VD> ненулевой посредством добавки регуляризующей константы. Это Тихонов.
VD> А по-правильному там должна быть не константа, а положительная функция,
VD> отражающая статистические свойства шума, присутствующего в B и/или A.
Тогда вообще фильтр Винера получается. К сожалению, это обычно недоступно...
С уважением
Евгений Машеров АКА СанитарЖеня
--- ifmail v.2.15
* Origin: FidoNet Online - http://www.fido-online.com (2:5020/175.2)
Вернуться к списку тем, сортированных по: возрастание даты уменьшение даты тема автор
|
