Разработка вычислительного алгоритма реконструкции объемных изображений в рентгеновской компьютерной томографии с применением голографических методов

Abstract

Рассматривается вычислительный алгоритм решения обратной задачи рентгеновской томографии по реконструкции внутренней структуры микрообъектов в ближнем зоне Френеля с применением голографических методов визуализации объемных изображений. Голографические методы не дают прямого решения задачи по реконструкции внутренней структуры объекта. Они могут только решить задачу объемного отображения некоторой поверхности объекта. Однако, используя данные по поглощению рентгеновского излучения объекта и фазоконтрастных голографических сигналов в ближней зоне Френеля, показана возможность получения объемного голографического изображения внутренних слоев объекта. Решение этой сложной задачи потребовало использование трeхмерного (3D) преобразования Радона внутренней функции объекта и двухмерного (2D) преобразования Радона фазоконтрастной голографической проекции. Был получен алгоритм реконструкции фазоконтастных томографических изображений внутренней структуры объекта и на основании его разработан вычислительный алгоритм для практической реконструкции объемных томографических изображений внутренней структуры микрообъектов. Результаты исследований были подтверждены математическим моделированием алгоритма реконструкции объемных изображений. Для чего, была разработана математическая модель тестового фантома и для него смоделированы фазоконтрастные проекции с последующей реконструкцией по ним фазоконтрастных томографических изображений томографическими методами с использованием разработанного авторами программного реконструктора. Разработан компактный алгоритм реконструкции фазоконтрастных томографических изображений для экспериментальных установок по исследованию малоразмерных биологических объектов, когда для получения фазоконтрастных проекций осуществляется два массива измерений: измерения яркости от объекта на детекторе в ближней зоне Френеля и измерения от объекта на детекторе в зоне поглощения. Практическое применение разработанного подхода в получении фазоконтрастных томографических изображений томографическими методами может найти для исследования микрообъектов, где требуется высокое разрешение на изображении, чем при применении классических методов КТ. A computational algorithm for solving the inverse problem of x-ray tomography for reconstruction of the internal structure of micro objects in the near Fresnel zone using holographic methods of visualization of three-dimensional images is considered. Holographic methods do not provide a direct solution to the problem of reconstruction of the internal structure of the object. They can only solve the problem of volumetric display of some object surface. However, using data on the absorption of x-ray radiation of the object and phase-contrast holographic signals in the near Fresnel zone, the possibility of obtaining a three-dimensional holographic image of the inner layers of the object is shown. The solution of this complex problem required the use of three-dimensional (3D) Radon transformation of the internal function of the object and two-dimensional (2D) Radon transformation of phase-contrast holographic projection. An algorithm for reconstruction of phase-contrast tomographic images of the internal structure of the object was obtained and based on it a computational algorithm for the practical reconstruction of three-dimensional tomographic images of the internal structure of micro objects was developed. The results of the research were confirmed by mathematical modeling of the algorithm for reconstruction of three-dimensional images. For this purpose, a mathematical model of the test phantom was developed and phase-contrast projections were simulated for it with subsequent reconstruction of phase-contrast tomographic images by tomographic methods using the program reenactor developed by the authors. A compact algorithm for reconstruction of phase-contrast tomographic images for experimental facilities for the study of small biological objects, when two arrays of measurements are carried out to obtain phase-contrast projections: measurements of brightness from the object on the detector in the near Fresnel zone and measurements from the object on the detector in the absorption zone. The practical application of the developed approach in obtaining phase-contrast tomographic images by tomographic methods can be found for the study of micro objects where high resolution is required in the image than in the application of classical CT methods.