Научная статья на тему 'Математический подход к описанию апертурной функции сканирующей системы в комптоновской томографии'

Математический подход к описанию апертурной функции сканирующей системы в комптоновской томографии Текст научной статьи по специальности «Математика»

CC BY
164
20
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
труды учёных ТПУ / электронный ресурс / математический подход / описание / апертурная функция / сканирующие системы / комптоновская томография / математические зависимости / рост / рассеивающий объем / аналитические соотношения / графические зависимости / сканирование

Аннотация научной статьи по математике, автор научной работы — Короткова И. А., Капранов Борис Иванович

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

i Надоели баннеры? Вы всегда можете отключить рекламу.

Похожие темы научных работ по математике , автор научной работы — Короткова И. А., Капранов Борис Иванович

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

Текст научной работы на тему «Математический подход к описанию апертурной функции сканирующей системы в комптоновской томографии»

УДК: 620.179.15

И.А. КОРОТКОВА, Б.И. КАПРАНОВ

МАТЕМАТИЧЕСКИЙ ПОДХОД К ОПИСАНИЮ АПЕРТУРНОЙ ФУНКЦИИ СКАНИРУЮЩЕЙ СИСТЕМЫ В КОМПТОНОВСКОЙ ТОМОГРАФИИ

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

Системы пространственного отбора комптоновской томографии конструируются с применением двух или более коллимационных каналов. Пересечение первичного пучка излучения и зоны чувствительности детекторов дает возможность получения рассеивающего объема (РО), которым проводится сканирование объекта контроля. Размер и форма рассеивающего объема описываются апертурной функцией S(x,y,z), участвующей в дальнейшей обработке изображения. Исследования показали, что экспериментальные измерения S(x,y,z) сопряжены с большими трудностями. Так, при снятии самой характеристики S(x,y,z) очень сложно реализовать "точечный" рассеиватель. При вычислении S(x,y,z), как производной от переходной характеристики V(x,y,z), сложно использовать "полубесконечную" геометрию по одному направлению, не нарушая геометрий по двум другим осям. Наибольшего соответствия, а следовательно, и большей точности реконструкции, можно добиться с помощью детального математического описания апертурной функции рассеивающего объема (АФРО) по конструктивным данным коллимационной системы.

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

Общие подходы для описания АФРО сформулированы в статьях [1] для установки, используемой в НПО "АЛТАЙ" г. Бийска и [2] для томографа "TomScan-200", г. Москвы.

В данной статье объектом исследования является система CBS LBD китайского производства [3]. Для построения АФРО необходимо детальное описание геометрической фигуры, которая образуется при взаимном пересечении первичного пучка излучения, фокусируемого узким цилиндрическим коллиматором и зоной чувствительности детекторов (4шт.), которая формируется коллимационной системой, состоящей из набора стальных пластин (рис. 1). Пластины собраны в секции по 11шт. с каждой стороны от центральной оси. Строение секции показано на рис. 2.

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

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

y0(z) = -1.1667-z +177.184 0 (1) y,(z) = -1.1667 z +175.866

Рис. 2. Секция фокусирующего коллиматора томографа CBS LBD

Рис. 1. Сканирующая система томографа CBS LBD

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

у2 (z) = 0.217 - z-36.97, у3 (z) = 0.099-z-10.835.

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

(2)

Сканирование вдоль оси г

При сканировании в данном направлении АФРО будет состоять из трех участков: 1. Зона входа (между точками С и Э)

АФРО определяется площадью сегмента, расположенного в сечении рассеивающего объема при пошаговом перемещении системы сканирования вдоль оси г. Параметры сканирующей системы: а = 1 мм - диаметр первичного пучка излучения;

а = 49,4° - угол, образованный пересечением центральных осей пучков коллиматора источника и детектора;

Н = 203мм - расстояние от фокуса рентгеновского аппарата до преколлиматора; И = \2мм - длинна преколлиматора; ^ = 151, Ъмм - фокус постколлиматора; Л = 80мм - длина постколлиматора;

р - 0,5мм - малая ось эллипса, находящегося в сечении первичного пучка излучения;

г = 0,659лш - большая ось эллипса, находящегося в сечении первичного пучка излучения;

у0 = 9,399 - протяженность рассеивающего объема вдоль оси у;

/0 = 4,467мм - расстояние от оси z до правой границы РО;

/, = 4,932лш - расстояние от оси z до левой границы РО;

р = -12,279мм - протяженность РО вдоль оси z. sin а

Исходя из параметров системы можно записать уравнение эллипса, находящегося в сечении РО:

x(z) = -yjr2-y{zf

(3)

где 7(2) = г-((-1,1667-2 + 177,184) + /0) - текущая координата г при сканировании. Тогда апертурная функция будет определяться по формуле

р ■ г - arceos

y(z)

-y(z)-x(z)

(4)

Для построения функции роста рассеивающего объема достаточно проинтегрировать полученную апертурную функцию на промежутке от (го,г\), который определяется координатами точек С и О и соответственно равен го=155,696, г|=154,567.

F,(z)= )s¿z)dz.

(5)

2. Средняя зона (между точками Э и В)

Искомая апертурная функция системы будет определяться площадью эллипса, расположенного в сечении РО:

8г{г) = л-р-г, (6)

где И(х) = г(х) ■ бш а - приращение функции при сканировании вдоль оси г, определяемое шагом сканирования

При этом для нахождения функции роста рассеивающего объема необходимо проинтегрировать полученную апертурную функцию на промежутке (гьгг), который определяется координатами точек О и В и соответственно равен г)=154,567, г1=147,64.

2\

3. Зона выхода (между точками В и А) Определяется аналогично зоне входа

г 4

где у(г) = -г + (1{ -(-1,1667-2 + 175,866)) - текущая координата г при сканировании. Тогда апертурная функция будет определяться по формуле

(7)

(8)

p-r- arceos

К*

\ \ -y(z)-x(z)

(9)

При интегрировании апертурной функции на промежутке (г2,г3), определяемом координатами точек А и В, соответственно 22=147,64, 23=146,54, получим формулу для определения функции роста РО.

(Ю)

22

Объединив все три рассмотренных участка, получим зависимости АФРО и функции роста рассеивающего объема при сканировании вдоль оси г, которые представлены на рис. 3 (а, Ь).

«56Д4 151-54 1-Й

1ЯМ 1ЯЛ и 54

0 0.5

мы

1.5 мм' 1 / \ 10 Ум. мм3 I а^ 15 мм' 10 1 Г\

0 / , \ 0 / 1 0 1 , \

0 0.5

мм

е Г

Рис. 3. Функции роста РО и АФРО при сканировании: а, Ь - вдоль оси ¿\ с, А - вдоль оси х; е, Г - вдоль оси у

Сканирование вдоль оси х

АФРО будет определяться площадью параллелограмма, расположенного в сечении РО при сканировании

8(х) = 2рИ(х), (11)

где И(х)■

• 5ш ос - высота параллелограмма.

При интегрировании полученной апертурной функции на промежутке, определяемом малой осью эллипса, получим функцию роста РО:

У(х)= ¡$(х)с£с. (12)

Полученные зависимости приведены на рис. 3 (с, с!).

Сканирование вдоль оси у

Апертурная функция будет определяться площадью эллипса, расположенного в се-

чении РО при сканировании:

Путем интегрирования 5(_у), получим функцию роста РО:

У(у)= ^Б(у)с1у. Полученные зависимости приведены на рис. 3 (е, 1).

(13)

АФРО представляет собой сложную объемную фигуру, образованную пересечением слаборасходящегося цилиндра и Х-образной протяженной фигуры. Получены координаты всех особых точек, математические выражения и массивы чисел, описывающие переходные функции V(x), V(y), V(z) и апертурные функции по всем направлениям

ад, s(y), ад.

Исходя из приведенных зависимостей, геометрические размеры РО имеют порядок 1 мм по оси х, 9,4мм по оси у, 9,2мм по оси z. Анализ графиков показывает, что наиболее оптимально сканирование объекта вдоль оси х, так как в этом случае апертурная функция имеет выраженный максимум и отсутствуют прямолинейные участки, которые вносят значительный вклад в потерю информации при восстановлении внутренней структуры контролируемого объекта. Также статистика квантов, определяемых функцией роста рассеивающего объема в данном направлении, максимальна по сравнению с остальными.

Полученные результаты будут использованы при реконструкции изображений в томографе CBS LBD, а примененный математический подход позволяет исследовать АФРО для произвольных геометрий, используемых в комптоновской томографии.

СПИСОК ЛИТЕРАТУРЫ

1. Kapranov B.I., Varga V.V. Spatial resolution compton of systems, Сб. докладов Russian-Korean International Symposium on Science and Technology, June 22-25, 1999, Novosibirsk, Vol 2, p. 672-676.

2. Kapranov B.I., Korotkova I. A . Increase of spatial resolution the compton tomographtom-scan-200 due to optimization of aperture function of disseminating volume. Сб. Докладов 5th Korea - Russia International Symposium on Science and Technology, June 26 - July 3, 2001, Tomsk, Vol 1, p. 360 - 362.

3. В. Фирстов, А Бронников, С. Пронин, А Кузнецов, Gu Benli, Huang Yan Wang, Zhou Wenming Xiaochun. Метрологические аспекты радиационных интроско-пических систем на основе регистрации обратнорассеянного излучения. Контроль. Диагностика, № 11 1999, с. 18-25.

УДК 620.179

/О. В. АЛХИМОВ, В. К. КУЛЕШОВ

МОЗАИЧНЫЙ ДЕТЕКТОР РЕНТГЕНОВСКОГО ИЗЛУЧЕНИЯ НА ОСНОВЕ ГРП ДЛЯ КОНТРОЛЯ ОБЪЕКТОВ БОЛЬШОГО РАЗМЕРА

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

Работа выполнена в рамках программы сотрудничества Министерства образования Российской Федерации и Министерства Российской Федерации по атомной энергии по направлению "Научно-инновационное сотрудничество" 2002 г.

Для рентгеновского контроля объектов большого размера необходимо использовать преобразователи рентгеновского изображения в видимое большой площади. Предложено использовать для увеличения площади контроля мозаичный газоразрядный экран. Мозаичный экран составлен из нескольких газоразрядных преобразователей (ГРП) рентгеновского излучения в видимое изображение. ГРП представляет собой отпаянную плоскую камеру, которая наполнена ксеноном при атмосферном давлении. На двух противоположных стенках камеры расположены электроды. При взаимодействии рент-

i Надоели баннеры? Вы всегда можете отключить рекламу.