logo
Лекции - Восстановление и реконструкция изображ

Алгоритм обратного проецирования

Рассмотрим еще одно соотношение, позволяющее вычислять обратное преобразование Радона. Пусть на плоскости XOY задана функция f(x,y), преобразованием Радона которой является функция R(s,). Функция R(s,) представляет собой интеграл от функции f(x,y) вдоль прямой, заданной уравнением

.

Зададим на плоскости XOY точку (x,y). В общем случае прямых, проходящих через эту точку и удовлетворяющих соотношению бесконечно много. Для любого угла всегда можно определить значение s, при котором точка (x,y) будет принадлежать заданной прямой. Из всего множества значений функции R(s,) выберем те которые получены при интегрировании вдоль прямых проходящих через точку (x,y) и просуммируем эти значения. В результате получим функцию

,

которая называется суммарной обратной проекцией.

Выразим R(s,) через обратное преобразование Фурье:

Подставим последнее соотношение в выражение . В результате получим

.

В соответствии с теоремой о центральном сечении . После подстановки этого соотношения в выражение получим:

.

Поменяем пределы интегрирования в последнем выражении:

.

Объединяем пределы интегрирования:

.

С другой стороны, по определению спектра обратной проекции

.

Перейдем в последнем выражении к полярным координатам

.

.

Сравним подынтегральные выражения в и . Видно, что

,

или

.

Вычислив обратное преобразование Фурье, получим

.

Получено еще одно соотношение для вычисления преобразования Радона. Несмотря на то, что при выводе выражения и выражения использовались различные предпосылки, эти выражения полностью эквивалентны.

  1. Yandex.RTB R-A-252273-3
    Yandex.RTB R-A-252273-4