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

0 y (расстояние)

−20 −Hv

Рис. 2.25. График функции «тонкой структуры» модели μ0 (0, y) и функции

μg (0, y) , полученной с помощью регуляризированного решения

Из рис. 2.25 видно, что регуляризированное решение достаточно хорошо описывает «тонкую структуру» решения μ0 (x, y) .

176

2.4.Анализ технических задач проектирования рентгеновских компьютерных томографов

2.4.1. Общие требования к рентгено-оптическому тракту

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

Теория взаимодействия рентгеновского излучения с объектом, а также экспериментальные данные для основных компонентов живого организма (Н2О, Са) позволяют определить эффекты и явления, которые необходимо учитывать при создании томографических систем и задании требований к ним.

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

Решение уравнения переноса излучения достаточно точно определяет геометрию их распространения от источника излучения до блока детектирования. В зависимости от диапазона энергии фотонов, испускаемых источником излучения, составляющая рассеянного излучения на выходе из объекта исследования может составлять различные значения. Для диапазона энергий 17–150 кэВ, который выбирается исходя из принципа «полупрозрачности», высокого контраста мягких (Н2О) и костных (Са) тканей человека, минимальной дозы облучения, на исследуемом фантоме человека (Н2О) диаметром 30 см составляет около 3 % от прямолинейно распространяемых фотонов для неколлимированного пучка. Эта доля рассеянного излучения достаточно высока, чтобы реконструировать

177

изображение с «шумом» не более 0,5 %, что требуется для томографических систем.

Для уменьшения доли рассеянного излучения необходимо использовать точечный источник излучения с коллимацией луча как в плоскости сканирования, так и в плоскости, перпендикулярной плоскости сканирования. Более того, эти условия создания геометрии рентгеновского луча определяют также требуемое физическое приближение решения уравнения переноса с учетом дополнительного условия – моноэнергетического источника излучения. В сумме эти условия дают возможность найти приближенные решения уравнения переноса и с помощью этого приближения решить уравнение КТ при практической реализации. Практическая реализация заключается в установлении «лучевых сумм» вдоль прямолинейного распространения рентгеновского луча, представлении этих «сумм» в виде проекции, равной логарифму отношения измеренной интенсивности на входе и измеренной интенсивности на выходе исследуемого объекта.

Необходимо отметить, что если бы преобладала криволинейная (рассеянная) составляющая излучения, то вопрос создания рентгеновского компьютерного томографа до сих пор был бы открытым, так как решение уравнения переноса (например, для лазерного излучения или СВЧ-излучения, где преобладает рассеянная (дифракционная) составляющая), для создания этих видов томографов в настоящее время не найдено. Исходя из физического приближения решения уравнения переноса для рентгеновского луча идеальной схемой сканирования является использование точечного моноэнергетического коллимированного в двух плоскостях источника излучения с единичным детектором. Эта схема исключает влияние рассеянного излучения, но при этом осуществляется параллельновращательное движение системы источник–детектор: параллельное движение – для набора дискретных значений единичной проекции, вращательное – для набора дискретных ракурсов.

при фиксированном угле θj выполняется идеальным образом для единичной проекции, так как измерение J0i (входной интенсивности из-

178

лучения перед объектом исследования) и Jxi (выходной интенсивности после объекта исследования) для единичной проекции проводится одним и тем же детектором. Это исключает влияние аддитивных и мультипликативных ошибок измерения детекторной системы и всего рентгено-оптического тракта преобразования измеряемого сигнала.

Учитывая, что параллельная схема сканирования является мало производительной, возможно использование веерной схемы с вращательным движением системы источник–детектор. Здесь используется множество единичных детекторов, и поэтому влиянием рассеянного излучения в объекте исследования здесь пренебрегать нельзя. Уменьшение влияния рассеянного излучения на качество изображения в этой схеме сканирования можно достичь конструктивным и алгоритмическим путем. В первом случае – это использование коллимированных ионизационных детекторных камер с достаточно большими в длину по направлению излучения катодными стенками из тяжелых (с большим Z) металлов (тантал, вольфрам, свинец). Если используется полупроводниковые детекторы (они имеют, как правило, достаточно малую длину в направлении излучения), то на окно детекторной линейки можно установить сетку-коллиматор с входным окном, равным входному окну единичного детектора.

Во втором случае используется корректировка соотношения ln J0i Jxi путем введения дополнительного соотношения в лога-

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

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

179

ется моноэнергетический, например, изотопный источник С060, то первичного коллиматора ставить не требуется.

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

Если учитывать, что при практическом проектировании томографических систем по критерию вес – размеры – мощность, как правило, используются источники излучения, основанные на принципе облучения тяжелой мишени электронами (рентгеновские трубки, бетатроны, линейные ускорители), т.е. полиэнергетические источники рентгеновского излучения, то при переносе излучения в исследуемом объекте на выходе происходит ужесточение рентгеновского пучка, т.е. смещение спектра в сторону высоких энергий по отношению к входному спектру. Этот эффект приводит к нелинейности в определении томографического параметра μ – линейного коэффи-

циента ослабления, который зависит от E – энергии фотонов. Эта нелинейность выражается в виде нелинейного уменьшения значения μ с увеличением толщины (расстояния от периферии до цен-

тра) исследуемого объекта. И здесь требуется корректировка проекции с целью уменьшения влияния полихроматичности. Корректировка проекции проводится в виде изменения соотношения J0i/Jxi в логарифме или введением дополнительного соотношения в логарифм на основании расчетных моноэнергетических «лучевых сумм» для сплошного водяного фантома человека и экспериментальных калибровочных данных на фантоме реальных полиэнергетических «лучевых сумм».

Суть корректировки полихроматичности должна заключаться в таком видоизменении рабочей проекции ln J0i Jxi , которое бы приводило к μ = const , на всем изображении томограммы (и в центре, и на

периферии) и значение которого бы соответствовало некоторой эффективной энергии фотонов, которая бы определяла соотношение

∫μ(Eэф)dl =∫μ(E)dl , где левый интеграл – лучевая моноэнергети-

L L

180