Привет, Хабр! На связи Smart Engines. Вы же помните, что кроме распознавания документов мы занимаемся и томографией.
Сегодня мы продолжаем нашу серию статей про фантомы для томографии. Мы уже рассмотрели медицинские тестово‑калибровочные фантомы, антропоморфные фантомы и ответили на вопрос: "Что такое пространственное разрешение в томографии и как его измерить?"
Дальнейшее продолжение предполагает рассказ о фантомах, полученных с помощью математических формул. Сегодня мы расскажем о самых первых математических фантомах, которые моделируют голову человека и его мозг, в частности. Такие фантомы нужны для различных расчетов, а для каких именно - рассказываем под катом.
Стилизованные антропоморфные фантомы
Первые математические фантомы появились для точных дозиметрических расчетов, обеспечивающих выработку безопасных протоколов рентгеноскопии, компьютерной рентгеновской томографии и ядерной медицины. В первых попытках математического моделирования использовали геометрические описания в виде сфер и эллипсоидов различного размера, но вскоре стало понятно, что эта информация имеет ограниченное применение при прогнозировании, поскольку фактическая геометрия тканей и органов не может быть хорошо смоделирована с помощью этого упрощенного предположения. Одним из первых (1969 г.) стилизованных (на основе уравнений) математических фантомов всего тела, в котором для моделирования внутренних органов использовался более широкий диапазон геометрических форм стал фантом, известный как MIRD-5 (Medical Internal Radiation Dosimetry (MIRD) Committee of the Society of Nuclear Medicine, Pamphlet No.5).
![Рис. схема дозиметрического фантома MIRD-5, [DeWerd-2014] Рис. схема дозиметрического фантома MIRD-5, [DeWerd-2014]](https://habrastorage.org/getpro/habr/upload_files/5a0/cb7/deb/5a0cb7debf9f4389b2db8fc5b2a783ee.png)
Этот фантом простыми уравнениями грубо моделировал скелет, пары легких и оставшиеся части мягких тканей для представления «среднего» здорового взрослого мужчины (эталонного человека), как это было определено Международной комиссией по радиологической защите (ICRP) на основе обширного обзора медицинской и другой научной литературы о популяциях Европы и Северной Америки. Новая серия стилизованных фантомов была разработана в 1980-1987 гг. и представляла собой целую «семью», которая состояла из взрослого мужчины и женщины, новорожденного и детей в возрасте 1, 5, 10 и 15 лет.
![Рис. визуализация математических фантомов семейства MIRD, адаптировано из [Xu-2000], [DeWerd-2014] Рис. визуализация математических фантомов семейства MIRD, адаптировано из [Xu-2000], [DeWerd-2014]](https://habrastorage.org/getpro/habr/upload_files/691/6cb/578/6916cb5788ed6a97e7bf9726c9f8c484.png)
Каждый фантом состоит из трех типов тканей различной плотности: кости, мягкие ткани и легкие. Основные части тела определены как эллиптические цилиндры (руки, туловище и бедра), усеченный эллиптический конус (ноги и ступни) и эллиптический цилиндр (голову и шею).
Фантомы Шеппа и Логана (Shepp-Logan Phantom, 2d и 3d)
Первый компьютерный рентгеновский томограф был представлен в 1971 г., но ещё много лет вычислительных возможностей компьютеров не хватало для детального моделирования процесса томографического сканирования тела человека. Тем не менее, развитие методов реконструкции требовало референтных математических моделей. Первые КТ сканеры использовались для пространственного представления внутренней структуры головы, поэтому именно модель головы стала первым специализированным математическим КТ фантомом. Разумеется, речь идет о широко известном фантоме Шеппа и Логана, предложенном в 1974 г. в классической работе [Shepp-1974].
![Рис. фантом Шеппа и Логана - пример визуализации “уровней серого”, схема и таблица параметров [Shepp-1974] Рис. фантом Шеппа и Логана - пример визуализации “уровней серого”, схема и таблица параметров [Shepp-1974]](https://habrastorage.org/getpro/habr/upload_files/f06/449/aa4/f06449aa4b51e06cdfc5d2d2604207db.png)
Фантомы Шеппа и Логана формируются из эллипсоидов, они удобны для аналитического моделирования проекций (особенно для методов, опирающихся на фурье-преобразование). Оригинальный 2D фантом моделирует один поперечный срез головы и состоит из 10 эллипсов, которые за счет разных размеров, ориентации и “уровней серого” моделируют основные анатомические структуры на срезе. Возможность аналитического моделирования “исходного сигнала” оказалась крайне важной, и для исследования численного метода реконструкции МР-изображения на основе набора измерений в параллельных плоскостях. Уже в 1980 г. в работе [Shepp-1980] предложена обобщение фантома головы Шеппа-Логана на 3D.
![Рис. схема 3d фантома Шеппа и Логана для вида спереди, сбоку и сверху, [Shepp-1980] Рис. схема 3d фантома Шеппа и Логана для вида спереди, сбоку и сверху, [Shepp-1980]](https://habrastorage.org/getpro/habr/upload_files/2bb/70a/c81/2bb70ac81b72ad4256a3c19bddabe37d.png)
Этот фантом содержит уже 17 эллипсоидов, которые соответствуют геометрии и характеристикам различных геометрических структур (например, нос, глаза, сгустки крови, желудочки, опухоли и многие другие). Позже упрощенную версию из 10 эллипсоидов для трехмерного (или усложненную версию оригинального двумерного) фантома головы использовали, например [Kak-1988], при тестировании алгоритмов КТ реконструкции для конусной схемы получения проекций.
![Рис. иллюстрация вертикального среза реконструированного объема фантома и схемы горизонтальных срезов в плоскостях, [Kak-1988] Рис. иллюстрация вертикального среза реконструированного объема фантома и схемы горизонтальных срезов в плоскостях, [Kak-1988]](https://habrastorage.org/getpro/habr/upload_files/535/4e1/01b/5354e101be2c1f851160da0e7aee4d07.png)
Наиболее полно аналитические свойства 3D фантома Шеппа и Логана раскрыты в работе [Koay-2007], которая придерживается классической геометрии и вносит небольшие вариации в уровни серого.
![Рис. визуализации срезов и рендеринг 3d фантома Шеппа и Логана по версии [Koay-2007] Рис. визуализации срезов и рендеринг 3d фантома Шеппа и Логана по версии [Koay-2007]](https://habrastorage.org/getpro/habr/upload_files/a17/2f7/186/a172f7186e8d9bc3d0a2ae5943c5562e.png)
Фантом мозга (Hoffman’s 3D brain phantom)
Голова человека много лет остается предметом интереса исследователей. Для изучения процессов диффузии растворов в головном мозге для развития методов позитронно-эмиссионной томографии (ПЭТ) в 1990 г. [Hoffman-1990] был предложен физический 3D фантом, который стал “золотым стандартом” в этой области.
![Рис. схема пластин для моделирования поперечных срезов мозга [Hoffman-1990] Рис. схема пластин для моделирования поперечных срезов мозга [Hoffman-1990]](https://habrastorage.org/getpro/habr/upload_files/d55/0f9/454/d550f94545e3b0905f1b6d0d52526289.png)
Фантом состоит из пластин с вырезами сложной формы, которые чередуются в определенном порядке и физически моделируют систему кровотока головного мозга человека образуя систему полостей и каналов.
![Рис. визуализация эталонных срезов и результатов регистрации фантомы, [Hoffman-1990] Рис. визуализация эталонных срезов и результатов регистрации фантомы, [Hoffman-1990]](https://habrastorage.org/getpro/habr/upload_files/bfb/6d0/c10/bfb6d0c1077edd3dc4b6130cb0669321.png)
Дополнительные вставки позволяют имитировать “горячие” и “холодные” сферические поражения, делать количественную оценку пространственного разрешения и размеров исследуемого объекта, оценивать эффекты ослабления и рассеяния.
В 2020 г. для физического фантома был создан цифровой двойник [Harrison-2020], при этом авторы провели серьезную работу по уточнению параметров цифрового фантома и согласования математической и физических моделей.
![Рис. сопоставление схем физического и математического фантомов Hoffman’s 3D brain, [Harrison-2020] Рис. сопоставление схем физического и математического фантомов Hoffman’s 3D brain, [Harrison-2020]](https://habrastorage.org/getpro/habr/upload_files/1e5/d38/0fe/1e5d380fe340e93a4f8cdd49b337212f.png)
Фантом является воксельным, доступен в DICOM и RAW форматах, может применяться для тестирования или проверки программного обеспечения анализа изображений в качестве достоверного эталона, для которого в настоящее время не существует общепринятых стандартов.
Многослойный (дисковый) фантом (Defrise phantom)
Для развития спиральной схемы сканирования с широким конусным пучком, которую и сейчас активно используют в КТ, необходимы свои тестовые объекты. Дисковый фантом является стандартным “простым объектом” для выявления проблем и анализа качества реконструкции при использовании широкого конусного пучка. Он упоминается под разными названиями и приводится в качестве иллюстрации в большом количестве работ, но мы затрудняемся указать оригинальную публикацию, поскольку ссылки на первоисточник отсутствуют. Отметим, что в некоторых публикациях диски имеют вид сильно сплюснутых эллипсоидов для упрощения численного моделирования.
![Рис. моделирование результатов реконструкции фантома для разных вариантов ширины конуса [Zeng-2010] Рис. моделирование результатов реконструкции фантома для разных вариантов ширины конуса [Zeng-2010]](https://habrastorage.org/getpro/habr/upload_files/fc9/d70/a38/fc9d70a38ecb4e7b87cf6389df70bda4.png)
Рис. моделирование результатов реконструкции фантома для разных вариантов ширины конуса [Zeng-2010]В других работах многослойность достигается при помощи обычных плоских цилиндрических дисков, которые легко могут быть воспроизведены в физическом воплощении фантома, или даже при помощи цилиндрических слоев для других направлений.
![Рис. схема усиройства дискового фантома [Zeng-2010] Рис. схема усиройства дискового фантома [Zeng-2010]](https://habrastorage.org/getpro/habr/upload_files/054/eac/f28/054eacf2879af907c88b5895aa46e381.png)
Фантомы для радиомики (Stanford DRO Toolkit)
Радиомика возникла на стыке онкологии, радиологии, обработки и анализа изображений, однако этот метод можно применить к любому медицинскому исследованию, при котором можно визуализировать патологический процесс. В основе подхода подтвержденная исследованиями гипотеза, что количественный анализ характеристик изображений позволяет улучшить компьютерную диагностику, улучшить выбор стратегий лечения и сделать более точным прогнозирование ответа на терапию. В радиомике основные характеристики изображений относятся к 5 классам дескрипторов: размер, форма, интенсивность, текстура и резкость границ, при этом часть значимых признаков может быть не видна человеческому глазу. Применение метода сильно усложняется тем, что при разных протоколах сканирования, реконструкции и постобработки для одних и тех же объектов изображения могут сильно отличаться, то есть под большим вопросом оказывается воспроизводимость и повторяемость вычисления признаков, а также возможность построения достоверных моделей. Для возможности анализа влияния всех многочисленных факторов на извлекаемые признаки и стандартизации между различными производителями оборудования и программного обеспечения были разработаны [McNitt-Gray-2020] цифровые референсные объекты (Digital Reference Objects, DRO).
![](https://habrastorage.org/getpro/habr/upload_files/90d/556/c94/90d556c94a10a05915e6f6e6e5a049bd.png)
Строгие математические описания позволяют управлять основными признаками, а сопровождающее программное обеспечение позволяет сгенерировать воксельные представления для заданных параметров квантования, а также эталонную сегментацию.
Заключение
Фантомы, о которых мы рассказали в этой статье, служат нам для количественной оценки разрабатываемых нами алгоритмов реконструкции, алгоритмов подавления артефактов или алгоритмов уточнения экспериментальных параметров. После тщательной проверки работоспособности методов они попадают в продуктизацию в наш продукт для томографической реконструкции Smart Tomo Engine. Видео с результатами томографической реконструкцией можно посмотреть в нашем ютуб-канале. Почитать об одном из прикладных применений томографии можно в нашей статье.
Список литературы
[DeWerd-2014] DeWerd L. A. The phantoms of medical and health physics. – Berlin : Springer, 2014.
[Harrison-2020] Harrison R. L. et al. Technical Note: A digital reference object representing Hoffman’s 3D brain phantom for PET scanner simulations // Medical physics. – 2020. – V. 47. – №. 3. – P. 1174-1180.
[Hoffman-1990] Hoffman E. J. et al. 3-D phantom to simulate cerebral blood flow and metabolic images for PET //IEEE Transactions on Nuclear Science. – 1990. – V. 37. – №. 2. – P. 616-620.
[Kak-1988] Kak A.C., Slaney M. Principles of computerized tomographic imaging. IEEE Press, New York, 1988.
[Koay-2007] Koay CG, Sarlls JE, Ozarslan E. Three-dimensional analytical magnetic resonance imaging phantom in the Fourier domain. Magn ResonMed 2007;58:430–436
[McNitt-Gray-2020] McNitt-Gray M. et al. Standardization in quantitative imaging: a multicenter comparison of radiomic features from different software packages on digital reference objects and patient data sets //Tomography. – 2020. – V. 6. – №. 2. – P. 118-128.
[Shepp-1974] Shepp LA, Logan BF. The Fourier reconstruction of a head section.IEEE Trans Nucl Sci 1974;21:21–43.
[Shepp-1980] Shepp L. A. et al. Computerized tomography and nuclear magnetic resonance //J. Comput. Assist. Tomogr. – 1980. – V. 4. – №. 1. – P. 94-107.
[Xu-2000] Xu X. G., Chao T. C., Bozkurt A. VIP-Man: an image-based whole-body adult male model constructed from color photographs of the Visible Human Project for multi-particle Monte Carlo calculations //Health physics. – 2000. – V. 78. – №. 5. – P. 476-486.
[Zeng-2010] Zeng G. L. Medical image reconstruction. – Heidelberg : Springer, 2010. – V. 530.