Анализ трехмерных алгоритмов реконструкции в рентгеновской компьютерной томографии

УДК 519.6

DOI: 10.14529/ctcr170202

АНАЛИЗ ТРЕХМЕРНЫХ АЛГОРИТМОВ РЕКОНСТРУКЦИИ В РЕНТГЕНОВСКОЙ КОМПЬЮТЕРНОЙ ТОМОГРАФИИ

Е.Н. Симонов, М.В. Аврамов, Д.В. Аврамов

Южно-Уральский государственный университет, филиал в г. Кыштыме

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

Ключевые слова: компьютерная томография, трехмерная реконструкция.

Введение

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

Так как общая формулировка задачи осталась неизменной - восстановление функции по ее интегралам вдоль прямых - подходы к решению трехмерной задачи схожи с теми, что применяются в классической двумерной КТ [1, 2].

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

Один из первых трехмерных алгоритмов реконструкции был предложен в 1984 году Фельд-кампом с соавторами. В нем полная геометрия конусной схемы учитывается только на этапе вычисления обратной проекции. Алгоритм дает приближенное решение для круговой траектории, обеспечивает хорошее качество восстановления, но не учитывает угол конуса (с увеличением угла растут ошибки реконструкции), хотя и превосходит по скорости классический подход и некоторые итерационные методы, требует больших затрат времени компьютерного счета. В связи с этим было предложено несколько модификаций алгоритма Фельдкампа. Основное внимание уделено созданию быстродействующих алгоритмов и исследованию возможных обобщений, связанных с усложнением способа сбора проекционных данных. Так в [3] предложен алгоритм, получивший название обобщенной параллельной обратной проекции (EPBP - extended parallel back-projection). От классического алгоритма Фельдкампа этот метод отличается применением пере-

паковки к параллельным данным и тем, что его использование не ограничено плоскими траекториями движения.

Альтернативой алгоритмам, основанным на методе Фельдкампа, может служить модернизированный алгоритм однослойной перепаковки (ASSR - advanced single-slice rebinning). Основная идея данного алгоритма заключается в том, что изображения не обязательно реконструировать в плоскостях, строго параллельных координатным плоскостям прямоугольной декартовой системы координат. Вместо этого реконструкция осуществляется в наклонных плоскостях, оптимально вписанных в спиральную траекторию фокусного пятна рентгеновской трубки [4].

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

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

1. Алгоритм обобщенной параллельной обратной проекции

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

Проекционные данные p0 (а, р,b) , где а - ракурс проекции, в - угол в пределах веера лучей,

b - продольная координата детектора, измеряются в области веера лучей.

Для реконструкции изображения удобно перевести данные в плоскость с параллельной геометрией. Параллельные лучи определяются их расстоянием до центра вращения и их углом наклона 0 к противоположному направлению оси OY. Соответствие координат векторов для веерного и параллельного пучка определяется простыми соотношениями:

S = a + P, а = S + arcsin

'О

V rf J

( % \

= -Rp sin р,р = - arcsin

V rf J

где RF - расстояние от источника излучения до центра вращения.

Переход от веерной геометрии к параллельной называется операцией перепаковки и проиллюстрирован на рис. 1.

Рис. 1. Перепаковка веерной геометрии хода лучей в параллельную

Общее преобразование координат (а, в) в (£,, 0) дано выражением (1) и разбивается на два шага. Между каждым из шагов перепаковки EPBP выполняет продольную перепаковку.

На первом шаге осуществляется переход от переменной а к переменной 0, его называют азимутальной или угловой перепаковкой.

px (В,Р,b) = p0 (а,р,b)coss,

где а = а(0, в) согласно выражению (1); cos s - корректирующий длину луча коэффициент необходимый для учета наклона лучей относительно плоскости XY в конусном пучке, который опре-

2 2 1( 2 2 \

деляется из выражения cos s = RFD/ (b + RFD I; RFD - расстояние от источника до детектора.

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

Если определить новую продольную переменную l как b = l + с X = db/dтакое, что dl/d = 0 в направлении касательной к спирали, то продольную перепаковку можно записать в виде

p2 (В, р, l ) = p (В, р, b ), где b = l + Ц,.

Операция свертки, предшествующая обратному проецированию, требует доступ ко всем \ е [-RM,RM ], где RM - радиус области реконструкции. Для выполнения продольной перепаковки без потери части изначальной области детектора b £ [-5, B ] (B - граница детектора) необходимо использовать значения l е (B + |X| RM )[-1,1]. Согласно b = l + для этого необходимы значения b е B + 2 |X|RM [-1,1], которых, строго говоря, нет. Поэтому алгоритм EPBP для достижения 100 % использования дозы расширяет значения детектора за счет экстраполяции (рис. 2) [3].

/

Рис. 2. Продольная перепаковка с экстраполяцией данных детектора

Следует отметить, что экстраполированные данные используются исключительно для операции свертки. Обратное проецирование осуществляется по начальной области детектора. Далее данные переводятся к параллельной геометрии согласно выражению (1): Рз (В,I) = р2 (В,р,I). Этот шаг называют радиальной перепаковкой.

Для алгоритма ЕРВР могут быть использованы стандартные ядра свертки для параллельных данных, такие как ядро Шеппа - Логана. Операция свертки дает свернутую проекцию

Р(В,Ъ,I) = рз (В,Ъ,I)*kф.

Завершающим этапом алгоритма является взвешенная обратная проекция, которая представляет собой оценку

f (х, >>, 2 ) = { р (в, ъ I) У (в) а в,

где у (В) - нормализованная весовая функция.

Алгоритм ЕРВР использует воксельно-зависимую весовую функцию. Для удобства в [3] вводится диапазон видимости V(r) рассматриваемого воксела г = (х, у, z). V(r) представляет собой набор всех углов проекции, при которых измеряется г. Если необходимо выполнить стандартную реконструкцию, то можно использовать весовую функцию, которая постоянна в пределах V и равна нулю вне этого предела.

Результатом операции обратной проекции является значение воксела в (х, у, z). Эта процедура должна быть выполнена для каждого положения воксела г.

2. Алгоритм модифицированной односрезовой перепаковки

Алгоритм ASSR сформулирован как алгоритм перепаковки данных спирального сканирования с конусным пучком излучения в данные параллельной геометрии сканирования. При этом вместо плоскостей, перпендикулярных оси спирали, используются наклонные плоскости реконструкции. Перепакованные данные реконструируются, используя для получения наклонных изображений стандартный двумерный алгоритм фильтрованной обратной проекции ^ВР-алгоритм) для параллельной геометрии сканирования. Затем набор наклонных изображений переписывается с применением z-фильтрации для получения окончательного реконструированного объема данных в декартовой сетке.

Как и для алгоритма ЕРВР проекционные данные получены в веерной геометрии, описываемой ракурсом а и углом в пределах веера р. При перепаковке к параллельной геометрии осуществляется переход к параметрам 0 и £ соответственно.

Источник рентгеновского излучения движется по спирали «(а) (рис. 3). Область реконструкции представляет собой цилиндр радиусом Ям. Яр -расстояние от центра вращения до источника излучения, RD - расстояние от центра вращения до детектора, d - продвижение стола за один полный поворот системы источник - детектор на 360°.

Важным параметром является позиция реконструкции ая, который определяет положение фокуса излучения на спирали (ракурс). Набор данных каждой плоскости должен быть центрован относительно определенной угловой позиции ак. Также координата 0 в параллельной геометрии рассчитывается относительно ак: 0 + ак = = а + р. Для реконструкции в параллельных лучах необходима только часть полного 360° вращения/. Требованием полноты является/ > 1/2 (180°). Соответственно 0е [ак - /П, ак + /П].

Сперва необходимо задать плоскости реконструкции R, оптимально вписанные в спиральную траекторию. Пересечение R с цилиндром, описывающим спиральную траекторию движения фокуса, представляет собой эллипс е(а). Угол наклона плоскости реконструкции у относительно плоскости XY должен быть выбран таким, чтобы среднее расстояние от фокуса до эллипса минимизировалось:

( /* ^

d а

у = min | |e(а)-5(а) i

V" fл

Равенство (2) выполняется в случае, когда плоскость пересекает траекторию движения источника в трех точках при а = 0, а = ±а (рис. 4).

Рис. 3. Спиральная схема сканирования в конусе лучей с плоским детектором

Общая нормальная форма R-плоскости может быть записана как

R : х cos aR tan у + y sin аR tan у - z + d —R = 0.

2%

Далее необходимо определить, какие лучи из конусного пучка должны использоваться для получения проекционных данных виртуального сканера, который вращается в плоскости R вдоль e(a). Но, строго говоря, эти лучи не будут полностью лежать в R, за исключением случаев а = aR ± а и а = aR [4].

if 4 Rp tan 7

тг/2 -a* a* тг/2

— Ry tan -у

_ d

4

Рис. 4. Пересечение функции da/2n, соответствующей спиральной траектории (линейная кривая), с ÄFtanysina, соответствующей эллипсу реконструкции

Для аппроксимации можно использовать полные веерные пучки, которые исходят от источника, и все лучи в таком веере должны быть использованы для реконструкции. После выбора эти лучи перепаковываются в параллельные. При этом можно полагать, что все лучи проецируются на R только в продольном направлении (направление по оси

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

2. Для каждого S е [- f л, f л) и е

-i R 1r

2 m' 2 m

получить

in C0Sу í i \

Р (В с)= i 2 2 2 "C0S SP l + аr, uf, vf ),

•у/sin В + cos у cos В

с

где а 'L (В,С) = В + sin- -определяет позицию фокуса, которая используется для перепаковки

rf

луча (0, £,) в продольной аппроксимации; cos s - коррекция длинны луча для учета угла между измеренным (физическим) лучом и виртуальным лучом; uF и vF - координаты детектора для веерного пучка, используемые при аппроксимации.

3. Реконструировать перепакованный набор данных p(0, £) для получения соответствующей функции объекта fx, y, z) с z = z(x, y), заданное выражением (3):

z (x, y ) = x cos aR tan у + y sin aR tan у + d —R.

2%

После перепаковки данных в p(0, £) к ним могут применяться стандартные алгоритмы реконструкции, подробней о которых можно узнать в монографии [6].

Рис. 5. Интерполяция изображений, реконструированных в косой плоскости, в изображения, соответствующие стандартным плоскостям XV

Получив определенное количество изображений fx, y, z(x, y)), последним шагом алгоритма является z-интерполяция этих наклонных плоскостей в декартовую сетку (рис. 5).

Алгоритм ASSR получил коммерческую реализацию (алгоритм AMPR) и модификацию на случай неполного набора данных с углом охвата меньше 180° (алгоритм SMPR) [4, 7, 8].

Алгоритмы EPBP и ASSR относятся к приближенным алгоритмам реконструкции изображений в трехмерной КТ, обеспечивают высокое качество изображений при небольшом угле конуса и обладают рядом полезных свойств: легко расширяются на различные траектории движения источника, предоставляют гибкость в выборе питча в спиральной КТ, позволяют эффективно использовать дозу излучения (до 100 % площади детектора), легко расширяются на случай КТ исследования сердца, относительно просты в реализации. Однако приближенные алгоритмы могут уступать в качестве получаемых изображений точным алгоритмам реконструкции, что проявляется в значительном увеличении артефактов изображения с увеличением угла конуса.

3. Точное решение задачи трехмерной реконструкции изображений

по конусным проекциям

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

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

Для точного восстановления трехмерной функцииf r), r = (x, y, z) с Л3 по данным конусных проекций необходимо, чтобы они удовлетворяли условию полноты Кириллова - Смита - Туя. Для этого необходимо, чтобы каждая плоскость, проходящая через исследуемый объект, по крайней мере один раз должна пересекать траекторию движения источника излучения. Примеры схем сканирования с траекторией движения источника, удовлетворяющих этому условию: спиральная, две окружности, окружность и прямая, две окружности, лежащие во взаимно перпендикулярных плоскостях [5, 6, 7].

Математически задачу можно определить как поиск функции объекта f из уравнения

g (s, u ) = J f (s + tu ) dt,

которое описывает интегрирование по полулинии с началом в точке s и направлении и.

Одним из путей обращения интеграла выражения (4) является применение трехмерного преобразования Радона, согласно которому

р

где символ «•» обозначает скалярное произведение. Трехмерное преобразование Радона описывает интегрирование _Дг) по плоскости перпендикулярной к вектору п и лежащей на расстоянии 5 от начала координат.

Формула обращения уравнения (5) имеет вид

f (r ) = -

4 V

J p (n, r • n ) dn,

где V - Лапласиан г. Для использования (6) необходимо установить связь (4) с (5): у*р(п,8 • П) = |у(ы • П)g(8,и)а?и,

где у - функция, такая что y(ts) = t 2у(s) , * - свертка по второй величине преобразования Радона [6, 9].

Один из алгоритмов, вывод которых основан на связи преобразования в конусе лучей с преобразованием Радона в R3 и классическом трехмерном обратном преобразовании Радона, представлен в [5]. Представленная в нем формула обращения теоретически точна, дает решение для произвольной траектории, представлена двумерным интегралом, имеет структуру FBP-алгоритма (а при выборе весовой функции w = 1 дает сверточный FBP-алгоритм) и допускает реализацию фильтрации и обратного проецирования инвариантную к сдвигу.

Альтернативный подход к разработке аналитических алгоритмов трехмерной реконструкции представлен в [6]. В отличие от классических FBP-алгоритмов, предложенный там алгоритм сначала выполняет операцию обратного проецирования и затем фильтрацию. Алгоритм может быть описан двумя шагами:

1. Обратная проекция g (s,u ) ^ fs (r) .

2. Фильтрация fs (r) ^ f (r) .

На втором шаге алгоритма вводится уравнение визуализации, которое позволяет рассматривать задачу реконструкции как задачу восстановления изображения с инвариантной к сдвигу функцией рассеяния точки [6].

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

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

Выводы

Доминирующую роль в трехмерной КТ играют приближенные методы реконструкции объемных изображений. Они, хоть и уступают в точности аналитическим алгоритмам, обладают полезными свойствами такими, как предоставляют произвольный выбор питча спирального сканирования, позволяют эффективно использовать дозу излучения (вплоть до идеального 100%-ного значения), распространяются на случай синхронизированной реконструкции (КТ сердца).

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

Литература

1. Симонов, Е.Н. Физика визуализации изображений в рентгеновской компьютерной томографии /Е.Н. Симонов. - Челябинск: Издат. центр ЮУрГУ, 2013. - 479 с.

2. Пикалов, В.В. Сравнение алгоритмов спиральной томографии /В.В. Пикалов, А.В. Лихачев // Вычислительные методы и программирование. - 2004. - Т. 5. - С. 170-183.

3. Kachelrie, M. Extended parallel backprojection for standard three-dimensional and phase-correlated four-dimensional axial and spiral cone-beam CT with arbitrary pitch, arbitrary cone-angle, and 100% dose usage /M. Kachelrie, M. Knaup, W.A. Kalender //Medical Physics. - 2004. - Vol. 31, no. 6. - P. 1623-1641.

4. Kachelrie, M. Advanced single-slice rebinning in cone-beam spiral CT /M. Kachelrie, S. Schaller, W.A. Kalender //Medical Physics. - 2001. - Vol. 27, no. 4. - P. 1033-1041.

5. Katsevich, A. A general scheme for constructing inversion algorithms for cone beam CT / A. Kat-sevich // International Journal of Mathematics and Mathematical Sciences. - 2003. - No. 21. - P. 13051321.

6. Bronnikov, A.V. Cone-beam reconstruction by backprojection and filtering / A.V. Bronnikov // J. Opt. Soc. Am. A. - 2000. - Vol. 17, no. 11. - P. 1993-2000.

7. Календер, В. Компьютерная томография основы, техника, качество изображений и области клинического использования / В. Календер; пер. с англ. А.В. Кирюшина, А.Е. Соловченко; под ред. В.Е. Синицына. - М.: Техносфера, 2006. - 344 с.

8. Stienstorfer, K. Segmented multiple plane reconstruction - novel approximate reconstruction scheme for multislice spiral CT / K. Stienstorfer, T. Flohr, H. Bruder // Phys. Med. Biol. - 2002. -No. 47. - P. 2571-2851.

9. Наттерер, Ф. Математические аспекты компьютерной томографии / Ф. Наттерер; пер. с англ. И.В. Паламодова; под ред. В.П. Паламодова. - М. : Мир, 1990. - 286 с.

Симонов Евгений Николаевич, д-р техн. наук, профессор, заместитель директора по науке, Южно-Уральский государственный университет, филиал в г. Кыштыме; e.n.simonov@yandex.ru.

Аврамов Максим Витальевич, аспирант кафедры экономики, управления и информационных технологий, Южно-Уральский государственный университет, филиал в г. Кыштыме; avramov.maksim91@mail .ru.

Аврамов Денис Витальевич, аспирант кафедры экономики, управления и информационных технологий, Южно-Уральский государственный университет, филиал в г. Кыштыме; avramov.91@ mail.ru.

Поступила в редакцию 26 сентября 2016 г.

DOI: 10.14529/ctcr170202

COMPARISON OF 3D RECONSTRUCTION ALGORITHM IN X-RAY COMPUTED TOMOGRAPHY

E.N. Simonov, e.n.simonov@yandex.ru,

M.V. Avramov, avramov.maksim91@mail.ru,

D.V. Avramov, avramov.91@mail.ru

South Ural State University, Kyshtym, Russian Federation

The article analyzes the existing three-dimensional image reconstruction algorithms in cone-beam X-ray computed tomography with a spiral trajectory of the source. The main approaches to the design of the reconstruction algorithms are considered, their brief classification is given. The focus is on approximate Feldcamp-type algorithms and methods of rebinning cone data to parallel geometry. The extended parallel backprojection and advanced single-slice rebinning algorithms, their advantages and disadvantages are described. The authors have considered the exact methods for solving image reconstruction problem for the 3D data set in cone-beam CT. The author describes the additional conditions imposed on the projection data, the computational efficiency of the algorithms, and the quality of the images. The emphasis is on relation between the cone beam and Radon transform and the classical 3D Radon transform inversion. The general scheme for constructing inversion algorithms for cone beam CT is described. The algorithm of cone-beam reconstruction by backprojection and filtering is considered. The advantages and disadvantages of exact reconstruction algorithms are presented.

Keywords: computed tomography, three-dimensional image reconstruction.

References

1. Simonov E.N. Fizika vizualizatsii izobrazheniy v rentgenovskoy komp'yuternoy tomografii [Phisics Imaging in X-Ray Computed Tomography]. Chelyabinsk, South Ural St. Univ. Publ., 2013. 479 p.

2. Pikalov V.V., Lihachev A.V. [Comparison of Algorithms for Helical Tomography]. Computing Methods and Programming, 2004, vol. 5, pp. 170-183.

3. Kachelrie M., Knaup M., Kalender W.A. Extended Parallel Backprojection for Standard Three-Dimensional and Phase-Correlated Four-Dimensional Axial and Spiral Cone-Beam CT with Arbitrary Pitch, Arbitrary Cone-Angle, and 100 % Dose Usage. Medical Physics, 2004, vol. 31, no. 6, pp.1623-1641.

4. Kachelrie M., Schaller S., Kalender W.A. Advanced Single-Slice Rebinning in Cone-Beam Spiral CT. Medical Physics, 2001, vol. 27, no. 4, pp. 1033-1041.

5. Katsevich A. A General Scheme for Constructing Inversion Algorithms for Cone Beam CT. International Journal of Mathematics and Mathematical Sciences, 2003, no. 21, pp. 1305-1321.

6. Bronnikov A.V. Cone-Beam Reconstruction by Backprojection and Filtering. J. Opt. Soc. Am. A., 2000, vol. 17, no. 11, pp. 1993-2000.

7. Kalender V. Komp'yuternaya tomografiya: osnovy, tekhnika, kachestvo izobrazheniy i oblasti klinicheskogo ispol'zovaniya. [Computed Tomography: Fundamentals, System Technology, Image Quality, Applications]. Moscow, Technosfera Publ., 2006. 344 p.

8. Stienstorfer K., Flohr T., Bruder H. Segmented Multiple Plane Reconstruction - Novel Approximate Reconstruction Scheme for Multislice Spiral CT. Phys. Med. Biol, 2002, no. 47, pp. 2571-2851.

9. Natterer F. Matematicheskie aspekty komp'juternoj tomografii. [Mathematical Aspects of Computer Tomography]. Moscow, Mir, 1990. 288 p.

Received 26 September 2016

ОБРАЗЕЦ ЦИТИРОВАНИЯ

FOR CITATION

Симонов, Е.Н. Анализ трехмерных алгоритмов реконструкции в рентгеновской компьютерной томографии / Е.Н. Симонов, М.В. Аврамов, Д.В. Аврамов // Вестник ЮУрГУ. Серия «Компьютерные технологии, управление, радиоэлектроника». - 2017. - Т. 17, № 2. -

С. 24-32. DOI: 10.14529/ctcr170202

Simonov E.N., Avramov M.V., Avramov D.V. Comparison of 3D Reconstruction Algorithm in X-Ray Computed Tomography. Bulletin of the South Ural State University. Ser. Computer Technologies, Automatic Control, Radio Electronics, 2017, vol. 17, no. 2, pp. 24-32. (in Russ.) DOI: 10.14529/ctcr170202