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

В компьютерной рентгеновской томографии трехмерный объект представляется обычно в виде набора тонких срезов. Для восстановления плотности среза решается задача обращения двумерного преобразования Радона. Для исследования ряда объектов более естественной является другая схема, когда источник излучения движется по некоторой пространственной кривой. Каждой точке кривой соответствует конус лучей, проходящих через эту точку. Исходными данными являются данные об ослаблении излучения при прохождении через объект. Математически задача ставится как задача восстановления функции трех переменных по интегралам вдоль прямых, проходящих через заданную кривую. Была получена формула обращения для функций, имеющих финитный носитель, и для кривых, удовлетворяющих определенным условиям. Главным в этих условиях является то, что любая плоскость, пересекающая объект пересекает кривую, по которой движется источник. Примером кривой, удовлетворяющей условиям, является совокупность двух единичных окружностей, лежащих во взаимно перпендикулярных плоскостях. Однако построение численных алгоритмов непосредственно на основании этой формулы, затруднительно. Дело, в частности, в том, что формула обращения основана на преобразовании Фурье от однородной функции, получаемой из исходных данных. Причем преобразование Фурье понимается в смысле обобщенных функций, а преобразование Фурье в обычном смысле может не существовать. В настоящем параграфе приводятся выражения для используемого преобразования Фурье, позволяющие при построении численных алгоритмов использовать метод, изложенный в предыдущем параграфе. Устанавливаются также некоторые соотношения между результатами

Пусть заданы функция f(x) = f(x1, x2, x3) , точка S = (s1, s2, s3) и вектор a = (a 1, a 2, a 3). Лучевым преобразованием функции f(x) будем называть функцию

,

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

Наряду с функцией в некоторых ситуациях рассматривается функция

,

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

Множество точек S, для которых известно лучевое преобразование обычно является множеством точек, принадлежащих некоторой кривой, являющейся траекторией движения источника излучения.

Пусть задана кривая, по которой движется источник, Ф(l ) = (Ф1(l ), Ф2(l ), Ф3(l )) параметр l пробегает некоторый интервал Ù действительной прямой. Для любого a = (a 1, a 2, a 3) и l Î Ù определим функцию

.

Функция g(a ,l ) есть интеграл от функции f(x) вдоль проходящего через точку Ф(l ) в направлении вектора a . Отметим, что при любом фиксированном l функция является l однородной функцией a степени -1:

. (2.1.1)

Для функций, имеющих финитный носитель, в [101] получена формула:

. (2.1.2)

При фиксированном l функция G+(b ,l ) есть преобразование Фурье от функции по переменной a , b = (cosq cosf , sinq cosf , sinf ). В формуле (2.1.2) l зависит от x и b и выбирается из условий: скалярное произведение (b , x) равно (b ,f (l )), но (b ,Ф(l )) не равно нулю. Значение функции f(x) может быть восстановлено в точке x, если такое l существует для любого b . Геометрически это означает, что любая плоскость, пересекающая точку x носителя функции, пересекает кривую Ф(l ) так, что знаменатель в (2.1.2) не обращается в нуль. Примером кривой, удовлетворяющей условиям Кириллова-Туя, является совокупность двух единичных окружностей, лежащих во взаимно перпендикулярных плоскостях, если носитель лежит в единичном шаре. Для цилиндрических объектов можно использовать винтовую линию.

В формулу (2.1.2) входит G+(b ,l ) - преобразование Фурье от функции , однако преобразование Фурье, понимаемое в обычном смысле:

,

в данном случае не существует, так как является однородной и имеет на бесконечности порядок 1/ê a ê . Преоразование Фурье здесь понимается в смысле обобщенных функций. Поскольку однородная функция, то при любом фиксированном l исходные данные, полностью определяются своими значениями на поверхности ê a ê =1. Переход к функции, заданной во всем пространстве R3 при использовании преобразования Фурье приводит к обобщенным функциям. Преобразование Фурье в смысле обобщенных функций является линейным функционалом над соответствующим пространством. Подробнее об этом будет сказано в следующих параграфах. Здесь нам важно отметить, что не любой функционал задается с помощью регулярной функции. Для того, чтобы использовать формулы типа (2) для построения алгоритмов, необходимо показать, что задается с помощью регулярной функции и иметь для нее выражения через функцию . В работе [101] дается выражение, связывающее , при x отличном от нуля с помощью регулярных операций с искомой функций f(x), то есть фактически показано, что функционал задается с помощью регулярной функции. Однако для построения алгоритмов томографической реконструкции нужно выразить не через искомую функцию f(x), а через исходные данные .

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