Современная медицинская диагностика включает в себя множество томографических методов: метод компьютерной томографии, магнитно-резонансная томография, позитронно-эмиссионная томография и однофотонная эмиссионная компьютерная томография. Все эти методы позволяют получать изображения внутренней анатомии пациентов. Такие методы как КТ и МРТ позволяют получать информацию о структуре внутренних органов человека, их плотности, коэффициентах линейного ослабления излучения и насыщенности атомами водорода. Также изображения, полученные при помощи этих методов, обладают высоким разрешением, большим, чем разрешение изображений, полученных при помощи ОФЭКТ и ПЭТ. Однако для точного диагноза необходима и информация о функциональных особенностях внутренних органов, которую позволяют получать такие методы как однофотонная эмиссионная компьютерная томография и позитронно-эмиссионная томография. Таким образом, методы радионуклидной диагностики являются неотъемлемой частью комплексного обследования и лечения пациентов.
Метод ПЭТ широко применяется при исследовании метаболизма и транспорта веществ в организме, экспрессии генов и других процессов. Исследование этих процессов было затруднено или невозможно до изобретения метода ПЭТ. Данный метод позволяет выявить болезнь на ее раннем этапе, что актуально при диагностике злокачественных опухолей и таких заболеваний как болезнь Альцгеймера и другие виды деменции. Очень часто ПЭТ-диагностика проводится совместно с другими видами томографических исследований, например компьютерной или магнитно-резонансной томографией.
Метод позитронно-эмиссионной томографии основан на регистрации пары гамма-квантов с энергией 511 кэВ, возникающих в результате позитронного бета-распада нуклида, входящего в состав вводимого пациенту радиофармпрепарата. При проведении ПЭТ-диагностики используют различные радиофармпрепараты, которые подбираются в зависимости от специфики исследования, что делает данный метод универсальным.
Детектирующая аппаратура в ПЭТ обычно представляет собой набор сцинтилляционных детекторов, расположенных в виде кольца (см. рис. 1). Сцинтилляторами называется группа веществ, испускающих гамма-кванты в видимом диапазоне длин волн при взаимодействии с излучением. В качестве детекторов используются различные сцинтилляторы, наиболее популярные в последнее время — кристаллы LSO (оксиортосиликат лютеция) и BGO (германат висмута).
Поток фотонов, испущенный сцинтилляторами, при помощи фотоэлектронных умножителей усиливается и преобразуется в электронный сигнал, который далее обрабатывается и, если сигнал соответствует определенному диапазону энергий, то компьютер сбора данных регистрирует единичное событие. Каждое единичное событие характеризуется индексом детектора, зарегистрировавшего его, и временной отметкой регистрации. На основе потока единичных событий регистрации строится список событий совпадения. Совпадением называется регистрация пары гамма-квантов, временные отметки которых лежат в определенном временном интервале, называемом временным окном совпадения. Величина окна выбирается в зависимости от параметров детектирующей аппаратуры и обычно не превышает десяти наносекунд.
Пара детекторов, зарегистрировавших совпадение, определяет линию отклика. В процессе сканирования вычисляется количество событий совпадения вдоль каждой линии отклика, полученные данные называются проекционными данными позитронно-эмиссионной томографии.
Параметры детектирующей аппаратуры, такие как количество составляющих ее детекторов и их геометрия, влияют на качество получаемых изображений, на разрешающую способность сканера. Поэтому актуальной является задача разработки методов моделирования сбора данных позитронно-эмиссионной томографии для различных конфигураций детекторов.
Различают двумерный и трехмерный режимы сбора данных позитронно-эмиссионной томографии (см. рис. 2, 3). При сборе данных в двумерном режиме кольца детекторов разделены септой коллиматора, поглощающей гамма-кванты, попадающие на детектор под углом, значительно отличающимся от прямого. В трехмерном же режиме септа коллиматора отсутствует, таким образом, возрастает чувствительность сканера, но в то же время увеличивается количество зарегистрированных случайных и рассеянных совпадений.
На качество изображений ПЭТ влияет множество факторов. Разрешающая способность метода ПЭТ ограничена в связи с тем, что до аннигиляции позитроны обладают небольшим (порядка нескольких миллиметров) свободным пробегом, который зависит от электронной плотности тканей. От качества изображений, полученных в ходе радионуклидной диагностики, зависит точность и корректность поставленного врачами диагноза, поэтому актуальной является задача разработки математических алгоритмов коррекции данных ПЭТ. Для тестирования алгоритмов коррекции необходимы проекционные данные, поэтому актуальной является также задача моделирования сбора проекционных данных ПЭТ с учетом различных факторов, влияющих на качество изображения.
Также модельные проекционных данные позитронно-эмиссионной томографии широко используются при разработке и проверке алгоритмов реконструкции томографических изображений.
К основным факторам, связанным с взаимодействием ионизирующего излучения и влияющим на качество изображения в ПЭТ, относятся ослабление излучения за счет поглощения и комптоновского рассеяния гамма-квантов. Также значительный шум в изображение вносят случайные и множественные совпадения, возникающие при регистрации гамма-квантов. Энергия гамма- квантов в ПЭТ составляет 511 кэВ, поэтому комптоновское рассеяние является основным эффектом при взаимодействии излучения и тканей органов пациента [19]. В результате регистрации рассеянных гамма-квантов ухудшается контрастность получаемого изображения. Случайные совпадения также отрицательно влияют на контрастность. Кроме того, помимо визуальной оценки изображений ПЭТ, зачастую врачи используют количественные методы обработки изображений, что требует предварительной коррекции шумов в изображениях, связанных с рассеянием, ослаблением и случайными совпадениями.
Данная работа посвящена моделированию сбора проекционных данных позитронно-эмиссионной томографии, а также разработке алгоритма коррекции одного из факторов, влияющих на качество изображений.
Результаты работы были представлены на различных конференциях: научно-практическая конференция «Радиационные технологии: достижения и перспективы. Ядерная медицина» (Ялта, 2014), международная научная конференция «Процессы управления и устойчивость» (Санкт-Петербург, 2014,
2015) , международная конференция «Устойчивость и процессы управления» в память В. И. Зубова (Санкт-Петербург, 2015), международная научная конференция «XXV Russian Particle Accelerator Conference» (Санкт-Петербург,
2016) , международная научная конференция «III International conference on Laser and Plasma Researches and Technologies» (Москва, 2017).
Результаты работы опубликованы в следующих работах:
1. Бажанов П. В. Моделирование сбора проекционных данных ПЭТ // Процессы управления и устойчивость, 2014, Т. 1 № 1. — С. 247-252.
2. Бажанов П. В. Исследование модели процесса сбора проекционных данных ПЭТ // Процессы управления и устойчивость, 2015, Т. 2. № 1. —
С. 276-281.
3. Bazhanov P.V. PET Projection Data Modeling and Scatter Correction // "Stability and Control Processes"in Memory of V.I. Zubov (SCP), 2015 International Conference, 2015. Article number 7342201, P. 516-517.
4. Бажанов П. В. Моделирование сбора проекционных данных ПЭТ и коррекция рассеяния // Устойчивость и процессы управления Материалы III международной конференции. 2015. С. 461-462.
5. Bazhanov P.V., Kotina E.D. Development of PET projection data correction algorithm // Journal of Physics: Conference Series, 2017, Volume 941 012097
6. Бажанов П.В, Котина Е. Д. Разработка алгоритма коррекции проекционных данных ПЭТ // Лазерные, плазменные исследования и технологии ЛаПлаз-2017 Сборник научных трудов III международной конференции. 2017. С. 91.
7. Бажанов П. В. Метод коррекции рассеяния и случайных совпадений в данных ПЭТ // Вестник Санкт-Петербургского государственного университета технологий и дизайна. Серия 1. Естественные и технические науки. 2017. С. 36-49.
8. Котина Е.Д., Овсянников Д.А., Плоских В.А., Латыпов В.Н., Бабин А.В., Широколобов А.Ю., Пасечная Г.А., Бажанов П.В. Программное обеспечение для обработки и визуализации данных однофотонной эмиссионной компьютерной томографии и позитронно-эмиссионной томографии //Вопросы атомной науки и техники. Серия: Техническая физика и автоматизация. 2015. — № 70. — С. 12-27.
9. E.D. Kotina, A.V. Babin, P.V. Bazhanov, D.A. Ovsyannikov, V.A. Ploskikh, A.Yu. Shirokolobov Mathematical and Computer Methods of Data Processing in Nuclear Medicine Studies // RuPAC2016 - Proceedings, Pages 480-482.
По теме исследования получено свидетельство о регистрации программы для ЭВМ № 2018611898 «Программа для моделирования сбора проекционных данных позитронно-эмиссионной томографии (AcqModeller)» от 08.02.2018.
[1] Tsoumpas C. Implementation and evaluation of scatter estimation algorithms in positron emission tomography. London: 2004. 118 pp.
[2] D. Strul et al. GATE (Geant4 Application for Tomographic Emission): a PET/SPECT general-purpose simulation platform // Nuclear Physics B (Proc. Suppl.). 2003. vol. 125. p. 75-79.
[3] Buvat I. and Castiglioni I. Monte Carlo simulations in SPET and PET // Q. J. Nucl. Med. 2002. vol. 46. p. 48-61.
[4] Harrison R. L., Vannoy S. D., Haynor D. R., Gillipsie S. B., Kaplan M.
S. and Lewellen T. K. Preliminary experience with the photon generator module of a public-domain simulation system for emission tomography // Conf. Rec. IEEE Nucl. Sci. Symp. and Med. Imag. Conf. (San Francisco). 1993. vol. 2. pp. 1154-1158.
[5] Bielajew A. F., Hirayama H., Nelson W. R. and Rogers D. W. O. History, overview and recent improvements of EGS4 NRCC Report PIRS-0436 // National Research Council. 1994. Ottawa. Canada. 26 p.
[6] Karine Assie et. al. Monte Carlo simulation in PET and SPECT instrumentation using GATE // Nuclear Instruments and Methods in Physics Research Section A: Accelerators, Spectrometers, Detectors and Associated Equipment. 2004. Volume 527, Issues 1-2, p. 180-189.
[7] S. Jan et. al. GATE: a simulation toolkit for PET and SPECT //Physics in medicine and biology. 2004. vol. 49. no. 19. pp. 4543-4561.
[8] Lewellen, T. K., R. L. Harrison. The SimSET program, in Monte Carlo Calculations in Nuclear Medicine. Applications in Diagnostic Imaging. S.-
E. S. M. Ljungberg, M.A. King. Bristol, UK. 1998. Institute of Physics Publication.
[9] Stute S. et. al. Monte Carlo simulations of clinical PET and SPECT scans: impact of the input data on the simulated images // Phys. Med. Biol. 2011. vol. 56(19). pp. 6441-6457.
[10] S.R. Meike, R.D. Badawy. Quantitative techniques in PET // Positron Emission Tomography: basic science and clinical practice. London: Springer. 2003. p. 115-146
[11] Hoffman E.J., Huang S., Phelps M.E., Kuhl D.E. Quantitation in positron¬emission computed tomography: 4. Effect of accidental coincidences. J. Comput Assist Tomography. 1981. vol. 5(3). p. 391-400.
[12] B.E. Cooke, A.C. Evans, E.O. Fanthome, R. Alarie and A.M. Sendyk. Performance figure and images from the Therascan 3128 positron emission tomograph // IEEE Trans. Nucl. Sci. 1984. vol. 31(1). p. 640-644.
[13] M.E. Casey and E.J. Hoffman. Quantitation in Positron Emission Computed Tomography: 7. A technique to reduce noise in accidental coincidence measurements and coincidence efficiency calibration // J. Comput. Assist. Tomogr. 1986. 10. p. 845-850.
[14] Badawi R.D., Miller M.P., Bailey D.L., Marsden P.K. Randoms variance reduction in 3D PET // Phys. Med. Biol. 1999. vol. 44(4). p. 941-954.
[15] David Brasse et.al. Correction Methods for Random Coincidences in Fully 3D Whole-Body PET: Impact on Data and Image Quality // J. Nucl. Med.
2005. vol. 46(5). p. 859-67.
[16] D. Brasse Correction methods for random coincidences in 3D wholebody PET imaging // Nuclear Science Symposium Conference Record. 2001. p. 923-930.
[17] J.F. Oliver, M. Rafecas. Modelling random coincidences in positron emission tomography by using singles and prompts: a comparison study // PLoS ONE, vol. 11(9). art. no. e0162096.
[18] Politte D.G., Snyder D.L. Corrections for Accidental Coincidences and Attenuation in Maximum-Likelihood Image Reconstruction for Positron¬Emission Tomography // IEEE Trans. Med. Imag. 1991. vol. 10. p.82-90.
[19] Гребенщиков В. В., Котина Е. Д. Физико-технические основы ядерной медицины. Спб.: СПбГУ, 2007. 172 с.
[20] Zaidi, H., Koral, K.F. Scatter modelling and compensation in emission tomography // Eur. J. Nucl. Med. Mol. Imaging. 2004. vol. 31. pp. 761¬782.
[21] J. M. Ollinger, Model-based scatter correction for fully 3D PET // Phys. Med. Biol. 1996. no. 41. pp. 153—176.
[22] Watson, C.C. New, faster, image-based scatter correction for 3D PET // IEEE Trans. Nucl. Sci. 2000. vol. 47 pp.1587-1594.
[23] Wollenweber S.D. Parameterization of a model-based 3-D PET scatter correction // IEEE Trans. Nucl. Sci. 2002. vol. 49. pp. 722-727.
[24] Accorsi R., Adam L.E., Werner M.E. et. al., Optimization of a fully 3D single scatter simulation algorithm for 3D PET // Phys. Med. Biol. 2004. vol. 49. pp. 2577-2598.
[25] Chen H.T., Kao C. M., Chen C.T. A fast, energy-dependent scatter reduction method for 3D PET imaging // IEEE Nuclear Science Symposium Conference Record. 2003. vol. 4. pp. 2630-2634.
[26] Adam L.E., Karp J.S., Freifelder R. Energy-based scatter correction for 3-D PET scanners using NaI(T1) detectors // IEEE Trans Med Imaging. 2000. vol. 19. pp. 513-521.
[27] Popescu L.M., Lewitt R.M., Matej S. et al. PET energy-based scatter estimation and image reconstruction with energy-dependent corrections // Phys. Med. Biol. 2006. vol. 51. pp. 2919-2937.
[28] Bazhanov P.V. PET Projection Data Modeling and Scatter Correction // Stability and Control Processes in Memory of V.I. Zubov (SCP 2015). Proceedings. 2015. pp. 516-517.
[29] Bazhanov P.V., Kotina E.D. Development of PET projection data correction algorithm // J. Phys.: Conf. Ser. vol. 941. 012097.
[30] Bailey D.L. Data acquisition and performance characterization in PET. In: Valk PE, Bailey DL,Townsend DW, et al, editors. Positron emission tomography: basic science and clinical practice. 1st edition. London: Springer; 2003. pp. 69-90.
[31] Floyd C.E., Jaszczak R.J., Greer K.L., et al. Inverse Monte Carlo as a unified reconstruction algorithm for ECT //J. Nucl. Med. 1986. vol. 27. pp. 1577-1585.
[32] Zaidi H., Herrmann Scheurer A., Morel C. An object-oriented Monte Carlo simulator for 3D positron tomographs // Comput. Meth. Prog. Biomed. 1999. vol. 58. pp. 133-145.
[33] Manavaki R., Reader A.J., Keller C., et. al. Scatter modeling for 3¬D PET list-mode EM reconstruction. IEEE Nuclear Science Symposium Conference Record 2002. vol.3. pp. 1863-1868.
[34] Кнут Д.Э. Искусство программирования. Том 1. Основные алгоритмы. Третье издание. Перевод с английского. Под общей редакцией Ю.В. Козаченко. Москва: Издательский дом «Вильямс». 2002. 720 с.
[35] Кнут Д. Э. Исскуство программирования. Том 2. Получисленные алгоритмы. Москва: Издательский дом «Вильямс», 2001. Т. 2. 832 с.
[36] Эддоус М., Стэнсфилд Р. Методы принятия решений. Пер. с английского под ред. член-корр. РАН И.И. Елисеевой. — М.: Аудит, ЮНИТИ, 1997. 590 с.
[37] Сивухин Д.В. "Атомная и ядерная физика". В 2-х ч. Ч.2. Ядерная физика. - М.: Наука. Гл. ред. физ. - мат. лит.1989. 416 с. - (Общий курс физики; Т.5).
[38] О.И.Василенко, Н.Г.Гончарова, В.К.Гришин, Ф.А.Живописцев,
Б.С.Ишханов, И.М.Капитанов, Э.И.Кэбин, Ж.М.Селиверстова, Н.А.Сотникова, В.Г.Сухаревский, Н.И.Тулинова, А.В.Шумаков “Субатомная физика. Вопросы. Задачи. Факты”: Учеб. Пособие/под ред. Б.С.Ишханова, - М: Изд-во МГУ,1994 - 224 с.