Анализ алгоритма построения имитационного рентгеновского изображения для инспекционно-досмотровых комплексов. Особенности настройки основных параметров, деталей моделирования движущегося рентгеновского источника и сборки результирующего изображения.
Аннотация к работе
Необходимо задать характеристики источника излучения: ориентацию в пространстве, расположение, форму и размер, мощность излучения, определить характеристики излучения, а также характеристики детектора: количество пикселей, размер и расположение. Также требуется учесть свойства объектов контроля: материалы, из которых они изготовлены, форму и расположение в пространстве. Построение имитационных рентгеновских изображений осуществляется в рамках одной из наиболее распространенных моделей ИДК HCV-Mobile, тактико-технические характеристики (ТТХ) которого представлены в Таблице 1. Между последовательными парами точек (вход в объект - выход из объекта) вычисляется длина траектории внутри объекта. Имитационное рентгеновское изображение строится как результат моделирования воздействия рентгеновского излучения на детекторную сетку набора источников, ослабленного взаимодействием с объектами контроля.Отметим, что данный и подобные ему алгоритмы имеют широкий спектр применения и, в частности, могут быть использованы для моделирования рентгеновских изображений. Реализация, описанная в статье, нашла свое применение в обучающей системе Симурен, одним из разработчиков которой является автор [5-7]. Сильной стороной продукта является получение исходной информации о геометрии и материалах сканируемых объектов из STL-файлов, которые создаются с использованием встроенного в приложение трехмерного редактора.
План
Содержание
Аннотация
Введение
Построение имитационного изображения
Алгоритм проекции
Заключение
Литература
Аннотация
Введение
Одним из актуальных направлений деятельности таможенных органов является внедрение технологий, содействующих ускорению и упрощению таможенных процедур и одновременно способствующих пресечению и профилактике правонарушений. Для этого в деятельность таможенных органов активно внедряется «Концепция создания системы таможенного контроля крупногабаритных грузов и транспортных средств» (Приказ ФТС России от 24.01.2005 № 52 // СПС «Консультант Плюс»).
В настоящее время таможенные органы Российской Федерации оснащены современными инспекционно-досмотровыми комплексами (ИДК) - техническим средством, использующим ионизирующее излучение для получения рентгеновского изображения крупногабаритных грузов и транспортных средств в целях таможенного контроля. ИДК позволяет без вскрытия и разгрузки грузового транспортного средства получить рентгеновское изображение транспортного средства и перевозимых в нем товаров. При этом имеется возможность идентифицировать перевозимые товары, обнаруживать в них предметы, запрещенные к перевозке.
ИДК должны эксплуатироваться только специально подготовленными сотрудниками, обладающими навыками распознавания теневых рентгеновских изображений объектов контроля. Обучение соответствующим компетенциям проводится в ростовском филиале Российской таможенной академии (РТА) с помощью специально разработанного учебного пособия [1].
Существующая в настоящее время методика подготовки операторов опирается на программные продукты компании «Smiths Heimann» [2,3]. Следует отметить, что указанные программные продукты функционируют исключительно под операционной системой Windows XP, уже не поддерживаемой разработчиком Microsoft. Более того, в настоящее время не существует компьютерных симуляторов, имитирующих процесс получения рентгеновского изображения с помощью технологии ИДК. Мировой опыт показывает высокую эффективность программ подготовки, обучения и контроля знаний на основе применения компьютерных моделей, объектов трехмерной графики, методов и алгоритмов имитации условий работы оператора [4]. Поэтому виртуальный симулятор рентгеновских изображений мог бы существенно способствовать развитию и совершенствованию базовых и продвинутых навыков анализа рентгеновских изображений.
Целью этой статьи является анализ алгоритма построения имитационного рентгеновского изображения для ИДК. Алгоритм был разработан и применен в ходе опытно-конструкторской работы по созданию программного продукта Симурен в Ростовском филиале РТА [5-7]. Представлено математическое и алгоритмическое обеспечение процесса построения имитационных рентгеновских изображений, приведены примеры имитационных изображений, построенных на основе реалистичных трехмерных изображений в стереолитографическом формате (STL).
Построение имитационного изображения
Для имитации рентген-изображения требуется определить значения различных параметров. Необходимо задать характеристики источника излучения: ориентацию в пространстве, расположение, форму и размер, мощность излучения, определить характеристики излучения, а также характеристики детектора: количество пикселей, размер и расположение. Также требуется учесть свойства объектов контроля: материалы, из которых они изготовлены, форму и расположение в пространстве.
Характеристики источника и излучения. Построение имитационных рентгеновских изображений осуществляется в рамках одной из наиболее распространенных моделей ИДК HCV-Mobile, тактико-технические характеристики (ТТХ) которого представлены в Таблице 1. Согласно характеристикам, для определенности выбрано верхнее положение источника. Геометрические параметры системы HCV-Mobile представлены на Рис. 1.
Таблица №1. Тактико-технические характеристики системы HCV-Mobile
ПараметрыХарактеристики
Энергетический уровень ускорителя3.8 МЭВ
Проникающая способностьМинимум 280 мм стали
Подсистема обнаружения рентгеновского излучения
Число датчиков детекторной линейки1216 датчиков на детекторной линейке Г - образной формы.
Высота сканированияВ нижнем положении источника: максимальная высота 4,45 м при расстоянии от земли 0,2 м; в стандартном: 4,50 м, 0,26 м; в верхнем: 4,75 м, 0,46 м.
Рис. 1 - Геометрические параметры системы HCV-Mobile.
Так как считаем, что в модели отсутствует геометрическая нерезкость изображения, обусловленная соотношением размеров и расположением источника и объектов контроля, источник излучения полагаем точечным [8]. Модель движущегося источника излучения представляет собой дискретное множество точечных равномощных источников, расположенных вдоль прямой линии. Энергия источника излучения принята за 3,8МЭВ.
Моделирование объектов контроля. В качестве объекта контроля выбран 20-футовый стальной контейнер, внешние длина, ширина и высота которого составляют 6096 мм, 2370 мм, 2591 мм. Толщина стенок контейнера - 5 мм. Для определенности контейнер расположен на расстоянии 500мм от вертикального детектора и 795 мм от источника излучения.
Объекты контроля, в том числе моделируемые грузы, являются замкнутыми трехмерными непересекающимися поверхностями, созданными в системах автоматизированного проектирования (CAD-системах, англ. CAD - computer aided design) [9] и хранящимися в STL-формате. Стереолитографический файл является полигональным представлением трехмерной поверхности. Поверхность представляется в виде серии треугольников. Каждая грань описывается вектором нормали, указывающим внешнее направление и тремя точками - вершинами треугольника. Эти данные используются в алгоритмах трассировки для построения центральных проекций. Указанные алгоритмы позволяют определить количество точек пересечения модели рентгеновского луча с поверхностями объектов контроля.
Каждый объект состоит из одной или нескольких частей из однородных материалов. Характеристики используемых материалов: плотность и массовые коэффициенты ослабления рентгеновского излучения для уровней энергии 3МЭВ и 4МЭВ - хранятся в специальной базе данных. Исходные данные взяты из справочника [10].
Моделирование рентгеновского излучения. Из каждого точечного источника испускается набор лучей к каждому центру пикселя соответствующей полосы детектора; количество лучей в пучке для каждого точечного источника равно количеству пикселей по вертикали детекторной сетки. Луч может пересекать элементы полигональной сетки поверхности некоторого объекта контроля. Каждый объект имеет четное число точек пересечения с лучом. Между последовательными парами точек (вход в объект - выход из объекта) вычисляется длина траектории внутри объекта.
Моделирование детектора. Все возможные положения детекторной линейки моделируются как детекторная сетка; один пиксель соответствует одному датчику. Детектор расположен параллельно движению источника излучения.
Построение имитационного изображения
Имитационное рентгеновское изображение строится как результат моделирования воздействия рентгеновского излучения на детекторную сетку набора источников, ослабленного взаимодействием с объектами контроля.
Предварительный этап
Настройки материалов. Производится загрузка информации из базы данных материалов, с помощью линейной интерполяции вычисляется линейный коэффициент ослабления для уровня энергии соответствующего ТТХ ИДК.
IMG_58af6834-5a2e-4d73-b1f5-4e03800b1cd8
IMG_e6019442-3de2-4fe9-84fa-af2dec0549ce
, где m - индекс материала, E - энергия излучения,
IMG_f013c8b8-0afa-4c38-96a4-08e04bbc83c5
IMG_097ff0c7-bc84-411e-9e11-4a937a356441 - коэффициент массового ослабления, ? - плотность материала.
Настройки объектов контроля. Определяются размеры ограниченной трехмерной области, в которой располагаются объекты контроля; объекты и их расположение загружаются из STL-файла, происходит сопоставление объектов с материалами.
Настройки источника. Определяется количество точечных источников (движущийся источник моделируется как набор точечных). Задаются расположения относительно контролируемой области (расстояние до контейнера, координаты начальной и конечной позиции источника - расстояние между ними совпадает с длиной контейнера).
Настройки детектора. Определяется расположение относительно источника и контейнера и размеры (один детектор - один пиксель).
Настройки грузового контейнера. Описываются в математической системе координат расположение контейнера между источником и детектором, внешние и внутренние размеры, толщина и материал стенок.
Также задаются Jmax (соответствует белому цвету на мониторе) и Jmin (соответствует черному цвету) приведенные интенсивности рентгеновского излучения, регистрируемые детектором для энергии излучения 3,8МЭВ на основе ТТХ ИДК из Таблицы 1.
Рентгеновские изображения можно рассматривать как специального вида проекции объектов контроля, совмещающие в себе принципы параллельного и центрального проектирования. После построения координаты проекции переводятся в экранные координаты. Далее описан алгоритм построения проекции объекта контроля на экранную плоскость.
Алгоритм проекции
Каждый точечный источник рассматривается как центр проекции. Проекционная плоскость - плоскость, параллельная источнику и проходящая через вертикальную детекторную линейку.
IMG_46043acd-5182-4b4f-9328-4652386bebcf
Рис. 2. - Источник излучения и объект контроля (вид сбоку).
Задаем параметры экранной системы координат (максимальное число пикселей в строке Xmax и максимальное число строк Ymax): начало координат находится в левом верхнем углу экрана, координата Хэ возрастает слева направо, а координата Уэ - сверху вниз.
Вводим мировую ортогональную трехмерную систему координат как показано на Рис.4
IMG_d5c1c48b-3acd-4589-b19c-7e203a5032fe
Рис. 4. - Расположение и размеры контейнера, детекторной сетки, источника излучения и проецируемой точки A (x, y, z) объекта контроля в заданной системе координат.
Описываем источник в выбранной мировой системе координат: начальное и конечное положения источника: (0,0,0), (xmax,0,0).
Определяем значения координат контейнера, размеры детектора (xmin, ymin, zmax) - (xmax, ymax, zmax), где ymax - сумма длин вертикальной и горизонтальной детекторных линеек.
Задаем ymid - максимально возможную координату проекции по y на вертикальную детекторную линейку.
Найдем проекцию точки A (x,y,z) объекта контроля на экран. Направляем луч из источника с координатами (x, y0, 0) в точку A и находим координаты проекции Аґ (хґ, уґ, zґ) на плоскость Z=zmax по формулам:
IMG_6e893d55-6fc6-41ee-b9aa-63b84506701f
Если луч проходит через вертикальную детекторную линейку, т.е. уґ?ymid, то искомая проекция Аэ (Хэ,Уэ) имеет следующие координаты:
IMG_ab7e145d-cb2a-41b3-9fd1-4029d4288e8d
В противном случае (уґ>y<mid):>
IMG_5bc1bf5b-369e-4cef-b66e-0e189311ce49
Таким образом, экран разбивается на две части: верхней соответствуют проекции на горизонтальный детектор, а нижней - на вертикальный.
Определение интенсивности серого цвета проекции
Основой алгоритма является закон ослабления рентгеновского излучения в веществе:
IMG_36a30d91-5f91-4ee7-9416-d0738703ef76
где I0 - изначальная интенсивность излучения, ?i - линейные коэффициенты ослабления материалов, di - расстояния, которые прошло излучение, в каждом материале, I - итоговая интенсивность излучения.
Алгоритм построения имитационного рентгеновского изображения, связанный с определением длины пути рентгеновского луча в веществе, выглядит следующим образом. Пусть интенсивность пикселя серого цвета меняется в диапазоне от 0 до Imax.
Для каждого пикселя детектора с координатами (Хэ,Уэ): Переводим экранные координаты в трехмерную математическую систему координат:
IMG_df662872-a3ee-4a9d-bddf-f65eb2691c6f
IMG_ff6f9298-61a5-4caa-b3d9-a3b9fab36eab
если Y < ymid, то Z=zmax; если Y ? ymid, то Z=zmax ymid-Y
Выпускаем виртуальный луч, идущий в обратном направлении из точки Аґ (X,Y,Z) к источнику излучения O (X,y0,0).
Для каждого объекта контроля i (i=1,…,N) выполняем запрос на пересечение луча с ограничивающим параллелепипедом соответствующего изображаемого объекта (CAD-модели). В случае успеха выполняем запрос на пересечение луч - треугольник с использованием отсортированного списка треугольников, полученных от i-й CAD-модели. Сортируем все точки пересечения по расстоянию от источника излучения. Так как существует четное число пересечений (в силу замкнутости полигональной сетки, моделирующей объект), вычислим расстояния между последовательными точками пересечения (входа и выхода из объекта) и найдем их сумму - длину траектории di рентгеновского луча внутри объекта контроля i.
Определяем значение ?i - коэффициента линейного ослабления материала i-го объекта контроля.
Вычисляем значение интенсивности пикселя I (Хэ,Уэ) приведенное к диапазону 0…Imax с помощью формулы:
IMG_2d5f9d69-f699-4683-86f4-d67a843589a9
Результирующее изображение строится с использованием полинейной сборки изображений для каждого положения источника излучения.
Примеры исходного трехмерного изображения и построенного с его использованием имитационного рентгеновского приведены на рисунках 5,6.
IMG_013b278f-8015-4cc3-868b-e05e6dad27a1
Рис.5. - Трехмерная сцена с объектами контроля: «20-футовый контейнер» и «ящики».
IMG_07141b9d-5fe8-492c-9bc8-a0cf1bd253ca
Рис. 6. - Результат имитационного рентгеновского сканирования. Видны скрытые объекты: «автомат» и «ящик с алкоголем».
Вывод
Отметим, что данный и подобные ему алгоритмы имеют широкий спектр применения и, в частности, могут быть использованы для моделирования рентгеновских изображений. Реализация, описанная в статье, нашла свое применение в обучающей системе Симурен, одним из разработчиков которой является автор [5-7].
Данная система применяется в ростовском филиале РТА. Сильной стороной продукта является получение исходной информации о геометрии и материалах сканируемых объектов из STL-файлов, которые создаются с использованием встроенного в приложение трехмерного редактора. Такой подход позволяет, изменяя исходные данные, видеть влияние этих изменений на результирующее изображение, что предоставляет широкие возможности по обучению работе с ИДК.
Кроме того, разработанный алгоритм и подобные ему могут быть использованы для проведения имитационных экспериментов с целью анализа возможностей ИДК. Другим преимуществом моделирования является возможность одновременного использования симулятора несколькими членами команды, что невозможно при работе с реальным оборудованием, которое, как правило, доступно в единичном экземпляре [11].
Другим направлением развития является построение трехмерных рентгеновских изображений. Это потребует разработки более сложной модели, включающей несколько источников рентгеновского излучения, и позволяющей восстанавливать трехмерное изображение по многоракурсному набору двумерных.
Таким образом применение разработанного алгоритма и его возможное развитие позволит повысить эффективность применения ИДК как средства таможенной проверки грузов, а также будет способствовать дальнейшему совершенствованию технологии.
Список литературы
1. Вербов В.Ф., Гамидуллаев С.Н., Кириллов Д.В., Шевцов А.В. Мобильный инспекционно-досмотровый комплекс HCV-MOBILE: основы эксплуатации и применения: учебное пособие / Ростов-на-Дону: Ростовский филиал Российской таможенной академии, 2011.178 с.
2. Heimann CARGOVISION mobile. Руководство по эксплуатации. Smiths Heimann, 2004.63 с.
3. Heimann CARGOVISION mobile. Руководство оператора по использованию программного обеспечения станции RIW. Smiths Detection, 2007.56 с.
5. Заключительный отчет о научно-исследовательской работе по теме: «Изучение алгоритмических особенностей проектирования компетентностного тренажера для инспекционно-досмотровых комплексов»/ Рук.А.С. Гречко. ?inoia-ia-Aiio: III IIO «EIAAEC Nenoain», 2013.195 n.
6. Заключительный отчет об опытно-конструкторской работе по теме: «Исследование возможностей программной реализации виртуального симулятора рентгеновских изображений в контексте технологий мобильных инспекционно-досмотровых комплексов / Рук.А.С. Гречко. ?inoia-ia-Aiio: III IIO «EIAAEC Nenoain», 2013.151 n.
7. Кудрявцев О.Е., Трофимчук А.М. Сетевая версия виртуального стимулятора, имитирующего рентгеновские изображения, полученные с использованием инспекционно-досмотровых комплексов // О проблемах импортозамещения в таможенных органах и новых разработках в сфере информационно-коммуникационных технологий: сборник материалов Межведомственной научной конференции. E?aa?ou: ?inneeneay oaii?aiiay aeaaaiey, 2016. N.43-47.
8. Duvauchelle, P., N. Freud, V. Kaftandjian and D. Babot, 2000. A computer code to simulate X-ray imaging techniques. Nuclear Instruments and Methods in Physics Research, 170 (B): pp.245-258.
9. Петров М.Н., Молочков В. П Компьютерная графика. Учебник. Серия: Учебник для вузов. NIA.: Ieoa?, 2004.812 n.
10. Creagh, D. C. X-Ray Absorption (or Attenuation) Coefficients / D. C. Creagh, J. H. Hubbell // Sec.4.2.4 in International Tables for Crystallography,Vol. C, A. J. C. Wilson, ed. Kluwer Academic Publishers, Dordrecht, 1992. pp.189-206.
11. Гуренко Б.В., Федоренко Р.В., Береснев М.А., Сапрыкин Р.В., Переверзев В.А. Разработка симулятора автономного необитаемого подводного аппарата // Инженерный вестник Дона. 2014. ?3. URL: ivdon.ru/ru/magazine/archive/n3y2014/2504