Приведение формул обращения томографической реконструкции в конусе лучей к виду, позволяющему строить численные алгоритмы
Научные материалы / Средства визуализации изображений в компьютерной томографии и цифровых рентгенографических системах / Математические основы компьютерной томографии / Приведение формул обращения томографической реконструкции в конусе лучей к виду, позволяющему строить численные алгоритмы
Страница 6

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

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

Напомним определение лучевого преобразования, которое было дано в предыдущих параграфах.

Лучевым преобразованием функции f(x) = f(x1, x2, x3) называется функция

, (2.2.8)

являющаяся интегралом от f(x) вдоль луча, исходящего из точки S = (s1, s2, s3) в направлении вектора a = (a 1, a 2, a 3).

Как уже отмечалось выше, в наряду с функцией рассматривается функция

,

являющаяся интегралом по всей прямой или, что тоже самое, суммой интегралов вдоль лучей из точки S в направлениях a и - a .

Обе функции являются однородными степени -1, то есть для них выполняются равенства

, .

Отметим также, что является четной, а функция таковой не является.

Понятие однородности степени l можно естественным образом расширить на обобщенные функции, если взять за основу равенство g(g x) = g l g(x). В терминах действия на основную функцию j равенство запишется в виде (g, j (x/g ) = g l +n (g, j (x)), здесь g v любое вещественное число большее нуля, n n - размерность пространства, в котором заданы основные функции. В интегральном представлении обобщенных функций показатель n возникает при соответствующей замене переменных в dx.

Известно, что преобразование Фурье однородной обобщенной функции, тоже является однородной обобщенной функцией.

Для интегрируемых, ограниченных и имеющих ограниченный носитель, функций f их лучевое преобразование является регулярной однородной функцией. Из результатов работ следует, что в трехмерном пространстве преобразование Фурье таких функций, понимаемое в обобщенном смысле, задается регулярной функцией. Регулярная однородная функция задается своими значениями на единичной сфере. Таким образом, в практических ситуациях при инвертировании лучевого преобразования нас интересует соотношение между двумя функциями. Одна из них является сужением на единичную сферу лучевого преобразования, а другая - сужением на единичную сферу преобразование Фурье лучевых данных, понимаемого в смысле обобщенных функций. Подобное преобразование между функциями, заданными на единичной сфере естественно назвать преобразованием Семянистого, поскольку в его работе впервые получены подобные соотношения для симметричных однородных функций в n-мерных пространствах. Как уже отмечалось выше, функция не является симметричной, для нее соответствующие соотношения для функций на единичных сферах в трехмерном пространстве были получены в предыдущих параграфах.

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

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

Лучевыми данными называется функция

,

Ф = (Ф1, Ф2, Ф3) Î R3, b Î S2 (S2 v единичная сфера). (Не трудно видеть, что в наших обозначениях это функция ).

В формулах обращения используются следующие функции:

(2.2.9)

, (2.2.10)

(S2/2 - половина единичной сферы), - скалярное произведение векторов и .

Страницы: 1 2 3 4 5 6 7