Материал: Симонов Томографические измерителные информационные системы 2011

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

вой сумме P(l,θ), прикладывается ко всем точкам, которые обра-

зуют этот луч. После того, как это сделано для всех проекций, получается приближенная аппроксимация исследуемого объекта. Так как в этом методе проекции как бы растягиваются обратно через восстанавливаемое сечение, метод и назван обратной проекцией. Для каждой точки изображения восстановленная μ(х,у) является суммой всех лучевых проекций, которые проходят эту точку. Поэтому метод обратной проекции называется методом суммирования или методом линейной суперпозиции.

P (l, θ)

2

Y Y

P (l, θ)

1

Рис. 3.40. Восстановление объекта методом обратной проекции: 1 – объект; 2 – его проекции; справа – восстановленный объект

Математическое описание метода обратной проекции может быть представлено выражением

j=1

281

где x cosθj + y sin θj = l в соответствии с (3.122) и суммирование

ходят через точку (х,у), коэффициент Δθj представляет угловое рас-

стояние между соседними проекциями, М – количество проекций. Символ μ указывает на то, что величина, полученная с помо-

щью уравнения (3.126), не идентична истинной μ . Как видно из

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

Можно показать, в чем состоит причина невысокой точности прямого применения метода обратной проекции.

Предположим, что объект исследования состоит из одной точки. Тогда результат восстановления по проекциям методом суммирования (3.126) будет представлять собой не точку, а многолучевую звезду, центр которой находится в восстанавливаемой точке

(рис. 3.41).

1

2

Рис. 3.41. Восстановление точечного объекта методом обратной проекции: 1 – объект; 2 – его проекция; справа – восстановленный объект

Очевидно, что точка будет представлена наиболее ярко, но в то же время на окружающее пространство эта точка будет наклады-

282

γ=−∞

где Σ (ν,θm ) – результат суммирования совокупности функций, получаемых из функции (ν) путем сдвига ее аргумента.

Учитывая (3.124), что P(l, θm ) = 0 при l >T , т. е. функция P(l, θm ) не является функцией с ограниченной шириной спектра, следовательно выборочные значения P(n l, θm ) и Σ (ν,θm ) под-

вержены эффекту наложения, который был рассмотрен в гл. 2. Влияние этого эффекта можно уменьшить фильтрацией функции

P(l, θm ) , осуществив ее до дискретизации. В реальности усред-

няющее действие такого рода осуществляется за счет апертурных ограничений (конечной ширины окна единичного детектора) блока детекторов. Даже в том случае, когда наложение мало, из выборочного характера данных следует, что функцию (ν,θm ) можно оце-

нить лишь внутри конечной полосы спектра шириной ν <12 l .

Это является одной из главных причин восстановления точки, через прямое применение метода обратной проекции, в виде “лучевой звезды”.

Рассмотрим проекционные теоремы.

Обобщенная проекционная теорема. Если h(l) – произвольная функция одного переменного, для которой существуют нижеследующие интегралы, то при любыхуглахθ справедливо равенство

Т

∫ Р(l,θ)h(l )dl = ∫∫μ(x, y)h(x cosθ + уsin θ)dx dy . (3.132)

−Т Ω

Доказательство. Рассмотрим левую часть равенства (3.132). Подставим в нее выражение для P(l,θ) из равенства (3.120) и, воспользовавшись выражением (3.122), сделаем замену переменных, перейдя от координат (l, t ) к прямоугольным координатам (x, y) .

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

284

екцией при угле θ эквивалентна определенной операции над исходным объектом μ(x, y). Если операция на μ(x, y) обратима, как например, преобразование Фурье, то из этого следует способ нахождения μ(x, y) по заданной P(l,θ).

Теорема о проекционном слое или центральном сечении.

Пусть (ν,θ) одномерное преобразование Фурье функции P(l,θ) по первому аргументу, определяемое формулой (3.130). Тогда

где μˆ – одномерное преобразование фурье-функции μ(x, y) . Доказательство. В качестве функции h(l) в выражении (3.132)

возьмем функцию exp(−i 2πνl ) . Результат следует из определений

одномерного преобразования Фурье по (3.130).

В соответствии с физической интерпретацией этот результат называется теоремой о центральном сечении. Действительно, пусть θ – заданный фиксированный угол. Тогда в теореме утверждается, что при изменении ν величина (ν,θ) равна значению μ при ра-

диусе ν и угле θ в фурье-пространстве. Перейдем к выводу формул реконструкции.

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

π ∞

μ(х, у) = ∫ ∫ (ν,θ)exp(i 2πν(x cosθ + y sin θ)) ν dνdθ . (3.134)

0 −∞

Подстановка выражения (3.130) для (ν,θ) в (3.134) дает фор-

мулу реконструкции функции μ(х,у) по заданным P(l,θ).

Из приведенной теоремы о центральном сечении следует, что при реконструкции изображения по выборочным проекциям не следует ожидать, что метод (3.134) даст лучшую оценку функции

μ(х,у), чем ее аппроксимация μ(х, у) с ограниченной шириной спектра. Преобразование фурье-функции μ(х, у) можно записать в виде произведения фурье-функции μˆ на функцию окна W (ν) :