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

Внимание! Если размещение файла нарушает Ваши авторские права, то обязательно сообщите нам

Зафиксируем какое-либо значение энергии E* (E1, E2 ) и пред-

ставим функцию μ(r, E) M0 в виде

μ(r, E ) = μ(r )+μ1 (r, E ),

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

μ

E

)

,

 

 

A ,

 

 

 

 

 

 

 

где

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r

,

Ai

= suppμi ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

μ(r ) = μ(r, E *), μ1 (r, E ) = i (

 

 

 

 

 

 

 

 

 

 

 

i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0, r

A ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A A, i =1,2,...n,

μ

(E ) > 0, E < E*,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i

 

 

 

 

 

 

 

 

 

i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

< 0, E* < E.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

), где

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

*

 

Пусть

 

μ(r, E) M1

 

и

 

μ(r, E )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= μ(r )+ μ(E, E

 

 

 

 

 

 

 

*

 

 

 

 

 

 

 

 

 

 

*

 

 

 

 

 

 

*

 

 

 

 

 

 

 

 

 

 

 

 

 

*

 

 

 

 

 

*

 

 

 

 

 

 

 

)= μ

(r, E

),

μ

(E, E

)= μ(E )

−μ(E

) ,

 

E

[E1, E2 ],

μ(r ) = μ(r, E

 

 

 

 

 

 

 

тогда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

*

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

max

 

 

 

μ1 (r, E )−μ(E, E

)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

max

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

μ −μ

= E [E , E

]

r A

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= μmax

−μmin ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

max

 

 

 

 

 

 

 

 

 

 

 

 

 

min μ(E )

 

 

 

 

где μmax = E [E1, E2 ]μmax (E ),

 

μmin = E [E1, E2 ].

 

 

 

 

Для мягких биотканей m ≈0,0108 см-1 и εm ≈ 0,175 при diamA =

= 20 см.

Заметим, что при наличии костных тканей m ≈ 0,8 см-1 и

εm ≈ 16.

Из этих оценок следует, что с удовлетворительной для

практики точностью

Jμ (ρ)

может быть аппроксимирован

Jμ (ρ)

только при диагностике мягких тканей. В этом случае имеем

 

E2

 

 

 

 

 

 

 

 

 

 

 

 

 

Jμ (ρ) = ∫ F0 (ρ, E)exp −μ(E, E* ) ∫

dL _ exp −

∫ μ(r , E* ) dL dE =

E1

 

L(ρ)

 

 

L(ρ)

 

 

 

 

 

 

 

 

 

 

 

∫

 

,

(2.44)

= Φ(ρ, E* ) exp −

μ(r , E* ) dL

 

L(ρ)

 

 

 

 

 

 

 

где

106

 

(

 

)

 

E2

0

(ρ, E)exp

{

 

(

 

)

 

0

}

Ф

ρ, E*

=

∫

−μ

E, E*

Z

 

 

 

F

 

 

 

 

(ρ) dE , (2.45)

 

 

 

 

 

E1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Z0 (ρ) =

 

∫

dL .

 

 

 

 

 

 

 

 

 

 

 

 

L(ρ)

 

 

 

 

 

 

Будем называть величину μ(r , E* ) эффективным линейным ко-

эффициентом ослабления рентгеновского излучения, а величину Ф(ρ, E* )– эффективной мощностью источника излучения. Введем

эффективный поток J* (ρ) рентгеновских фотонов с помощью со-

отношения

 

J* (ρ) = Ф(ρ, E* ) exp{− ∫ μ(r , E* ) dL .

(2.46)

L(ρ)

 

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

ным потоком J*(ρ) рентгеновских фотонов, ослабляющихся в ве-

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

описываются функцией Ф(ρ, E* ). На рис. 2.14 показаны графики

Zs = Zs (μl ) , где

E

 

 

Zs = ∫2

Z λ(E) (μe) dE

(2.46′)

E1

 

 

для однородного водного фантома радиуса R0

в случаях парал-

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

107

достаточной коллимации существенен эффект накопления вторичных фотонов (угол α – раствор пучка в направлении детекторов).

Zs , усл. ед

α = 44°

1,1

α = 0,1°

1

 

−μ

R0

 

 

μ

R0

μ

l

 

0

Рис. 2.14. График функции Zs = Zs (μl )

 

 

 

Ф(

 

l, E*) J0

 

 

 

 

 

 

 

 

 

μ

Zmax

= 1,25

 

 

 

 

 

 

 

 

 

 

 

 

 

1,03

 

 

 

 

 

 

1,2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1,01

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1,0

 

 

−μ

R

1,1

 

μ

R

μ

l

 

 

0

0

 

 

 

 

 

 

0,99

 

 

 

 

 

 

1,05

1,00

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0,92

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рис. 2.15. Зависимость нормированной эффективной мощности источника излучения от безразмерной величины μ l при фиксированной Z

Также представляет практический интерес проанализировать зависимость Ф(μl, E* ) от величины Zmax = Z . Графики этой зави-

симости показаны на рис. 2.15, анализируя которые можно сделать вывод, что если величина Zmax достаточно мала (~1,05–1,2), то

108

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

2.2.3. Эффект частично заполненного объема

Когда строят алгоритм реконструкции, то необходимо знать моноэнергетическую лучевую сумму вдоль линии L(l,θ) (см. рис.

2.1) для определенных l и θ, т. е. для определенных прямых L(l,θ) .

Моноэнергетическую лучевую сумму

P(M )

получают из (2.20)

 

 

JH (l,θ, E0 )

 

(l,θ)

 

 

(M )

 

 

∫

μ(x, y)dl .

 

P(l,θ)

= −ln

 

=

(2.46″)

C0

 

 

 

L(l,θ)

 

 

Статистика фотонов и изменение энергетического спектра рентгеновского излучения при прохождении через вещество являются причинами того, что физические измерения дают нам возможность

только приближенно оценить значения P((lM,θ)) .

Существуют другие погрешности, влияющие на получение P((l,Mθ )) . Одна из погрешностей связана с тем, чт о источник (фокаль-

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

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

Пусть имеется точечный источник моноэнергетического излучения И и детектор Д в виде отрезка линии, и пусть линейный ко-

109

результате рабочих
AE0 ≈ 568 .
Д
Рис. 2.16. Пояснение эффекта частично заполненного объема
500000 a
A
b
И

эффициент ослабления равен нулю во всей области, кроме части области А, которая заштрихована (рис. 2.16), где он равен 2. Примем, что длина заключенного между ad и bc отрезка любой из прямых, идущих от И к Д, равны 1.

Допустим, кроме того, что эталонный материал имеет линейный коэффициент ослабления, равный нулю (вакуум), и что число фотонов, регистрируемых эталонным детектором за время как калибровочного, так и рабочего измерения, равно 1000. Следовательно, число фотонов, вылетающих из источника в направлении детектора за время проведения калибровки, и соответствующее число фотонов за время проведения рабочего измерения приблизительно равны. Пусть это число составляет 1000000. Таким образом, в результате калибровочного измерения, получим на основании (2.25)

СЕ0 ≈ 1000 .

500000

d

1

c

Если разбить рентгеновский пучок на две равные половины, как это показано на рис. 2.16, то очевидно, в каждую из половин А войдут примерно по 500000 фотонов. Из левой части, где линейный коэффициент ослабления равен 0 и, следовательно, коэффициент пропускания равен 1, все 500000 фотонов дойдут до детектора. Из правой же части, где линейный коэффициент ослабления равен 2 и, следовательно, коэффициент пропускания составляет е–2 ≈ ≈ 0,135, число фотонов, которые дойдут до детектора, будет равно

примерно 500000 e−2 ≈ 68000 .

Следовательно, общее число зарегистрированных детектором фотонов составит около 568000, и в

измерений на основании (2.26) получим

110

Источник: https://studfile.net/preview/16708896/