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

Алгоритм обратного проецирования (англ. backprojection algorithm) — это метод цифровой обработки сигналов и реконструкции изображений, используемый для восстановления двумерного или трёхмерного распределения некоторой физической величины (например, плотности, радиоизлучения или отражающей способности) по набору её проекций, полученных под различными углами. Алгоритм является основой для ряда методов в компьютерной томографии, радиолокации (в частности, в радиолокации с синтезированной апертурой — РСА) и сейсморазведке. В отличие от более распространённого метода фильтрованных обратных проекций, классический алгоритм обратного проецирования не включает этап фильтрации, что приводит к характерному размытию (артефакту в виде «звёздного» шума) в восстановленном изображении, но компенсируется простотой реализации и возможностью работы с неполными или нерегулярными данными.

Принцип работы

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

Математическая формулировка

Пусть \( f(x, y) \) — функция, описывающая неизвестное распределение (например, коэффициент ослабления рентгеновского излучения в томографии). Проекция \( p(\theta, s) \) — это интеграл от \( f(x, y) \) вдоль прямой линии, заданной углом \( \theta \) и расстоянием \( s \) от начала координат:

\[ p(\theta, s) = \int_{-\infty}^{\infty} f(x, y) \, \delta(x \cos\theta + y \sin\theta - s) \, dx \, dy \]

где \( \delta \) — дельта-функция Дирака.

Алгоритм обратного проецирования восстанавливает изображение \( \hat{f}(x, y) \) следующим образом:

\[ \hat{f}(x, y) = \int_{0}^{\pi} p(\theta, x \cos\theta + y \sin\theta) \, d\theta \]

То есть для каждой точки \( (x, y) \) вычисляется сумма (интеграл) значений проекции \( p \) для всех углов \( \theta \) при \( s = x \cos\theta + y \sin\theta \). В дискретном случае интеграл заменяется суммой по конечному числу углов.

История

Идея обратного проецирования восходит к работам Иоганна Радона, который в 1917 году опубликовал математическую основу преобразования, названного впоследствии преобразованием Радона. Радон показал, что функцию двух переменных можно однозначно восстановить по её интегралам по всем прямым линиям. Однако практическая реализация стала возможной только с развитием вычислительной техники в середине XX века.

В 1960-х годах Годфри Хаунсфилд и Аллан Кормак независимо друг от друга разработали первые компьютерные томографы, в которых использовались алгоритмы, основанные на обратном проецировании. Хаунсфилд в 1972 году представил первый коммерческий КТ-сканер, который использовал метод последовательных приближений, близкий к обратному проецированию. Позднее, для устранения размытия, был разработан метод фильтрованных обратных проекций, который стал стандартом в медицинской томографии.

В радиолокации алгоритм обратного проецирования начал применяться в 1970-х годах для обработки данных радиолокаторов с синтезированной апертурой (РСА). Он позволяет формировать изображение земной поверхности с высоким разрешением, компенсируя искажения, вызванные движением носителя.

Виды и модификации

Классическое обратное проецирование (без фильтрации)

Этот метод даёт размытое изображение, так как каждая точка проекции «размазывается» по всей линии обратного проецирования. В результате возникает эффект «звёздного» шума, когда яркие точки объекта создают крестообразные артефакты. Для устранения этого эффекта применяется фильтрация.

Фильтрованное обратное проецирование (Filtered Back Projection, FBP)

Является наиболее распространённой модификацией. Перед обратным проецированием проекции подвергаются фильтрации высокочастотным фильтром (например, фильтром Рамачандрана — Лакшминараяна или фильтром Шеппа — Логана). Это позволяет скомпенсировать размытие и получить чёткое изображение. FBP лежит в основе большинства медицинских КТ-сканеров.

Полярное обратное проецирование

Применяется в радиолокации, когда данные получены в полярной системе координат (угол — дальность). Алгоритм работает непосредственно в полярных координатах, что упрощает интерполяцию и уменьшает вычислительные затраты.

Взвешенное обратное проецирование

Используется для учёта неравномерности плотности проекций или для подавления шума. Каждой проекции присваивается весовой коэффициент, зависящий от угла или расстояния.

Применение

Медицинская компьютерная томография (КТ)

В КТ алгоритм обратного проецирования (в форме FBP) используется для реконструкции поперечных срезов тела пациента. Проекции получаются путём измерения ослабления рентгеновского излучения детекторами, расположенными вокруг пациента. Алгоритм позволяет восстановить распределение коэффициента ослабления, что даёт возможность визуализировать внутренние органы, кости и патологии.

Радиолокация с синтезированной апертурой (РСА)

В РСА алгоритм обратного проецирования применяется для формирования изображений земной поверхности. Радиолокатор, установленный на самолёте или спутнике, излучает импульсы и принимает отражённые сигналы. Для каждой точки земной поверхности вычисляется задержка сигнала, и соответствующие значения суммируются по всем положениям носителя. Это позволяет достичь высокого пространственного разрешения, независимо от дальности.

Сейсморазведка

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

Ультразвуковая томография

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

Оптическая когерентная томография (ОКТ)

В ОКТ, используемой в офтальмологии и дерматологии, алгоритм обратного проецирования применяется для восстановления трёхмерных изображений биологических тканей.

Достоинства и недостатки

Достоинства

• Простота реализации: алгоритм не требует сложных математических преобразований, легко реализуется на цифровых вычислителях.
• Универсальность: работает с любыми геометриями сбора данных (параллельными, веерными, коническими пучками).
• Робастность: устойчив к неполным данным (например, при ограниченном угле обзора).
• Параллелизуемость: вычисления для каждой точки изображения независимы, что позволяет эффективно использовать многопроцессорные системы и графические процессоры (GPU).

Недостатки

• Размытие изображения: без фильтрации даёт нечёткую картинку с артефактами.
• Высокая вычислительная сложность: для больших объёмов данных (например, в 3D-томографии) требует значительных вычислительных ресурсов, хотя современные GPU позволяют это компенсировать.
• Чувствительность к шуму: при низком отношении сигнал/шум качество восстановления ухудшается.

Сравнение с другими методами

В медицинской томографии алгоритм обратного проецирования (FBP) долгое время был стандартом, но в 2010-х годах всё большее распространение получают итеративные методы реконструкции (например, метод максимального правдоподобия — MLEM, или метод статистической реконструкции). Итеративные методы позволяют лучше подавлять шум и работать с меньшими дозами облучения, но требуют большего времени вычислений. В радиолокации алгоритм обратного проецирования остаётся одним из основных, особенно для систем с нерегулярной траекторией носителя.

Интересные факты

• Алгоритм обратного проецирования является дискретной реализацией обратного преобразования Радона.
• В ранних КТ-сканерах 1970-х годов использовался именно классический алгоритм без фильтрации, что приводило к значительному размытию изображений. Пациентам приходилось делать снимки с большим количеством проекций (до 180) для получения приемлемого качества.
• В радиолокации алгоритм обратного проецирования позволяет формировать изображения с разрешением, превышающим физическую апертуру антенны в десятки и сотни раз, что и дало название методу «синтезированная апертура».
• В 2009 году за разработку методов компьютерной томографии, включая алгоритмы обратного проецирования, Нобелевская премия по физиологии и медицине была присуждена Полу Лотербуру и Питеру Мэнсфилду (за магнитно-резонансную томографию), однако Хаунсфилд и Кормак получили Нобелевскую премию за КТ ещё в 1979 году.

Источники

• Радон, И. (1917). «Об определении функций по их интегралам вдоль некоторых многообразий». Сообщения Саксонской академии наук.
• Хаунсфилд, Г. Н. (1973). «Компьютерная томография: история и принципы». Британский журнал радиологии.
• Кормак, А. М. (1963). «Представление функции по её интегралам по прямым линиям с некоторыми приложениями в радиологии». Журнал прикладной физики.
• Каку, А. С., Слейни, М. (1988). «Принципы компьютерной томографии». IEEE Press.
• Каррара, У. и др. (1995). «Алгоритмы обратного проецирования для РСА». IEEE Transactions on Aerospace and Electronic Systems.
• Барретт, Х. Х., Суинделл, У. (1996). «Радиологическое изображение: теория формирования изображения». Academic Press.