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

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

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

Пусть заданы функция f(x) = f(x₁, x₂, x₃) , точка S = (s₁, s₂, s₃) и вектор a = (a ₁, a ₂, a ₃). Лучевым преобразованием функции f(x) будем называть функцию

,

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

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

,

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

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

Пусть задана кривая, по которой движется источник, Ф(l ) = (Ф₁(l ), Ф₂(l ), Ф₃(l )) параметр l пробегает некоторый интервал Ù действительной прямой. Для любого a = (a ₁, a ₂, a ₃) и 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. Переход к функции, заданной во всем пространстве R³ при использовании преобразования Фурье приводит к обобщенным функциям. Преобразование Фурье в смысле обобщенных функций является линейным функционалом над соответствующим пространством. Подробнее об этом будет сказано в следующих параграфах. Здесь нам важно отметить, что не любой функционал задается с помощью регулярной функции. Для того, чтобы использовать формулы типа (2) для построения алгоритмов, необходимо показать, что задается с помощью регулярной функции и иметь для нее выражения через функцию . В работе [101] дается выражение, связывающее , при x отличном от нуля с помощью регулярных операций с искомой функций f(x), то есть фактически показано, что функционал задается с помощью регулярной функции. Однако для построения алгоритмов томографической реконструкции нужно выразить не через искомую функцию f(x), а через исходные данные .

Итак, перейдем к нахождению . Мы будем использовать то, что является однородной функцией по a фиксированном l . В [95] доказано следующее

Утверждение: Пусть есть преобразование Фурье в смысле обобщенных функций от однородной функции , тогда

. (2.1.3)

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

Представим в виде

,

(поскольку параметр l фиксирован, его на данном этапе можно опустить).

Как уже отмечалось выше, интеграл является расходящимся, тем не менее, переходя к сферическим координатам по обычным правилам, получаем:

,

где b = b (j ,q ) = (cosq cosj , sinq cosj , sinj ), j Î [-p /2, p /2], q Î [0, p ].

Учитывая, что , а также то, что интегрирование по углам j и q соответствует интегрированию по единичной сфере, приходим к выражению

.

Интеграл по r есть преобразование Фурье от r ₊+. Используя таблицы для преобразования Фурье обобщенных функций [19], приходим к выражению (2.1.3).

Для действительных функций f(x) в формуле (2) нужна мнимая часть :

.

Используя обобщенные функции, сосредоточенные на поверхности [19], получаем следующее следствие:

.

Здесь S(x ) = {g Î S^2½(x , g ) = 0), v производная по направлению x . Подставляя в (2.1.2) функции и , зависящие от параметра l , получаем формулу обращения, пригодную для построения численных алгоритмов:

(2.1.4)

Здесь S(x ) v окружность, являющаяся пересечением единичной сферы и плоскости P(b ). Плоскость P(b ) проходит через начало координат ортогональна вектору b . Символ W (x ) означает интегрирование по окружности. Оператор L(b , D) означает дифференцирование функции в направлении вектора b :

,

при этом l , зависящее от b и x, остается фиксированным.

Как и выше, b = b (q , j ) = (cosq cosj , cosq sinj , sinq ), l = l (q , j ) = l (x, b ) такое, что скалярное произведение (x, b ) равно (b , g (l )) и (b , g ^/(l )).

В формуле (4) используются регулярные функции, и она пригодна для построения численных алгоритмов.

Замечание. А.С. Денисюком независимо и другим методом, без явного использования преобразования Фурье обобщенных функций, получены формулы обращения функции g⁺ в Rⁿ . При n = 3 формулы А.С. Денисюка и формулы, получаемые изложенным способом из формулы Туя, совпадают.

Выше были получены формулы, позволяющие строить численные алгоритмы восстановления функции f(x) = f(x₁, x₂, x₃) по ее лучевому преобразованию

Далее мы будем опускать символ f и использовать обозначение .

При фиксированном S функция является функцией в трехмерном пространстве, но в силу ee однородности существуют поверхности, такие что полностью определяется своими значениями на них (поверхности расположения приемников излучения).

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

Плоский детектор.

Мы будем предполагать, что для источника, находящегося в точке S = (s₁, s₂, s₃), исходные данные регистрируются в плоскости P, определяемой уравнением xs₁ + ys₂ + zs₃ = -½ S½ . Плоскость P, определяется следующими условиями:

плоскость P перпендикулярна лучу, соединяющему источник с началом координат;

плоскость P проходит через точку S= (s₁, s₂, s₃.)

Расстояние D между плоскостью регистрации и источником равно удвоенному расстоянию от источника до начала координат. В плоскости регистрации будем использовать прямоугольную систему координат (p₁, p₂), начало которой находится в точке пересечения с лучем, соединяющим источник с точкой (0, 0, 0). Таким образом, если источник находится в точке S = (s₁, s₂, s₃), то начало системы координат (p₁, p₂) в плоскости наблюдения находится в точке с трехмерными координатами -s₁, -s₂, -s_{3 =}- S.

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

Траектория в виде двух окружностей.

Рассмотрим окружность, лежащую в плоскости z =0.

Направление оси p₂ в плоскости регистрации будет совпадать с направлением оси z.

Ось p₁ системы координат возьмем на линии пересечения плоскости регистрации с плоскостью, содержащей окружность, по которой движется источник. Для окончательного определения системы координат необходимо выбрать одно из двух возможных направлений оси p₁. Если s₃ = 0, s₁ = r cosl , s₂ = r sinl (источник движется в плоскости z =0), то положительный единичный вектор на оси p₁ выберем так, чтобы он совпадал с вектором (cos(l +p /2), sin(l +p /2), 0) = (-sinl , cosl , 0) = (-s₂/½ S½ , s₁/½ S½ , 0).

Точка, имеющая в плоскости регистрации координаты (p₁, p₂), имеет следующие пространственные координаты:

x = -p₁ sinl - r cosl = -p₁ s₂ /½ S½ - s₁ ,

y = p₁ cos l - r sinl = p₁ s₁ /½ S½ - s₂ , z = p₂.

В случае плоского детектора, исходными данными являются интегралы по лучам, соединяющим точки (p₁, p₂) в плоскости регистрации с источником S.

Регистрируемая функция g_r(p₁, p₂, l ) есть интеграл от искомой функции f(x) = f(x₁, x₂, x₃) вдоль луча исходящего из точки S = (s₁, s₂, s₃) = (rcosl , r sinl , 0) в направлении точки

P = (-p₁ sin l - rcosl , p₁ cosl - r sinl , p₂ ) = (-p₁ s₂/½ S½ v s₁, p₁ s₁/½ S½ v s₂, p₂).

Интегральная форма регистрируемой функции имеет вид:

При t = 0 луч проходит через точку S = (rcosl , rsinl , 0), при t = 1 v через точку P = (p_1, p₂) = (-p₁ sin l - rcosl , p₁ cosl - r sinl , p₂).

Итак, мы имеем соотношение между функциями g_r(p₁, p₂, l ) и :

,

.

Наряду с обозначением g_r(p₁, p₂, l ), мы будем использовать обозначения g_r(p₁, p₂, S(l )), g_r(p₁, p₂, S) и g_r(P, S) , здесь S(l ) точка на траектории источника, соответствующая параметру l , P = (p₁, p₂). Мы выразили функцию g_r(p₁, p₂, l ) через функцию = g⁺ (x , l ).

В формуле обращения лучевого преобразования используется функция g⁺ (x , l ) = для того, чтобы использовать g_r(p₁, p₂, l ), регистрируемую в случае плоского детектора, нужно выразить g⁺ (x , l ) используя g_r(p₁, p₂, l ).

Для дальнейшего нам потребуются координаты (p₁, p₂) (в системе координат плоскости регистрации) точки пересечения плоскости регистрации данных с лучем (S +tx ) = (s₁ + tx ₁, s₂ + tx ₂, s₃ + tx ₃). Эти координаты имеют вид:

.

.

Теперь мы можем выразить используя gr(p1, p2, l ):

= g+ (x , l ) = gr(2 ½ S(l )½ (s2(l )x 1 v s1(l )x 2) /, -2½ S(l )½ 2x 3 /,l ),

если < 0, = 0, если ³ 0.

Итак, мы имеем следующее соотношение между функциями:

g+ (P, l ) и = g+ (x , l ); P = (p1, p2), x = (x 1, x 2, x 3,);

= g+ (x , l ) =

= gr(2 ½ S(l )½ (s2(l )x 1 v s1(l )x 2) /, - 2½ S(l )½ 2x 3 /,l ),

если < 0,

= 0, если ³ 0.

При переходе от функции g+ (x , l ) = к функции gr (P, S) интегрирование по окружности S(l ) в трехмерном пространстве заменяется на интегрирование по прямым линиям в плоскости регистрации. Отметим, что формулы обращения лучевого преобразования, использующие интегрирование вдоль прямых в плоскости регистрации.

Похожие работы

Цифровые методы рентгенодиагностики

Скачать

50528

0

0

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

Цифровая рентгенография

Скачать

33118

0

7

... обычная система электронно-оптического преобразования для просвечивания, из ЭОП, телевизионного тракта с высоким разрешением, рентгеновского высоковольтного генератора и рентгеновского излучателя Рис.5 Цифровая рентгенография с экрана ЭОП 1-генератор; 2-рентгеновская трубка; 3-пациент; 4-ЭОП; 5-видеокамера; 6-аналого-цифровой преобразователь; 7-накопитель изображений; 8-видеопроцессор; 9-сеть; ...

Основные платформы ЭВМ

Скачать

127846

0

0

... банковских систем. Наиболее популярны сегодня смешанные решения, при которых часть модулей банковской системы разрабатывается компьютерным отделом банка, а часть покупается у независимых производителей. Основными платформами для банковских систем в настоящее время считаются: 1. ЛВС на базе сервера PC (10,7%); 2. Различные модели специализированных бизнес-компьютеров фирмы IBM типа AS/400 ...