Зафиксируем какое-либо значение энергии 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
эффициент ослабления равен нулю во всей области, кроме части области А, которая заштрихована (рис. 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