|
ru.algorithms
- RU.ALGORITHMS ----------------------------------------------------------------
From : Valentin Davydov 2:5020/400 24 Jan 2002 08:25:34
To : Evgenij Masherow
Subject : Re: Ищется алгоритм чистки радиокарты
--------------------------------------------------------------------------------
> From: "Evgenij Masherow" <EMasherow@nsi.ru>
> Date: Tue, 15 Jan 2002 20:28:52 +0300
>
>Задача распадается на две части:
>1. Оценка "вслепую" апертурной функции.
>2. Деконволюция.
Hеобязательно. Можно придумать некий итеративный алгоритм, постепенно
"улучшающий" изображение и одновременно выдающий всё более точные приближения
для аппаратной функции (типа, разложить её по ортогональным функциям, а потом
деконволюировать карту с каждым членом разложения до "упора в ноль").
>2. Hаиболее эффективна в вычислительном отношении деконволюция через Фурье.
>Она основана на теореме "Фурье от свертки равно произведению Фурье исходных
>последовательностей". При деконволюции Фурье от свернутой (наблюдаемой) карты
>делится на Фурье от апертурной функции.
>Примечания:
>а. В теореме фигурирует циклическая свертка. Поэтому нужно либо использовать
>окна, убывающие к краям, либо дополнять нулями.
>б. Делить на ноль интересно, но не очень полезно. Если верить астроному
>Стругацкому - от этого на ушах шерсть растет... Поэтому если в Фурье-образе
>нули - надо регуляризовать...
Дык, Тихонов именно это и делает. Переписывает выражение F(A)/F(B) в
виде F(A)*cnjg(F(B))/abs(F(B))^2, после чего имеет в знаменателе
вещественную неотрицательную функцию, каковую можно сделать заведомо
ненулевой посредством добавки регуляризующей константы. Это Тихонов.
А по-правильному там должна быть не константа, а положительная функция,
отражающая статистические свойства шума, присутствующего в B и/или A.
>в. Вычитать ли предварительно среднее - надо подумать.
IMHO, не надо, по причине того, что радиоисточники некогерентные, то есть как
"истинная" карта, так и аппаратная функция суть существенно неотрицательные.
Вал. Дав.
--- ifmail v.2.15dev5
* Origin: Demos online service (2:5020/400)
Вернуться к списку тем, сортированных по: возрастание даты уменьшение даты тема автор
|
Re: Ищется алгоритм чистки радиокарты 