Томографическая реконструкция

Продолжительность: 25 секунд. 0:25

Томографическая реконструкция: проекция, обратная проекция и обратная проекция с фильтром.

Томографическая реконструкция — это тип многомерной обратной задачи , задача которой состоит в том, чтобы получить оценку конкретной системы на основе конечного числа проекций . Математическая основа томографической визуализации была заложена Иоганном Радоном . Ярким примером применения является реконструкция компьютерной томографии (КТ), при которой изображения поперечного сечения пациентов получаются неинвазивным способом. Недавние разработки показали, что преобразование Радона и его обратное значение используются для задач, связанных с реалистичной вставкой объектов, необходимых для тестирования и оценки использования компьютерной томографии в службах безопасности аэропортов .

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

Представляем формулу

Рисунок 1: Геометрия параллельного луча, используемая в томографии и томографической реконструкции. Каждая проекция, полученная в результате томографии под определенным углом, состоит из набора интегралов по линиям, проходящим через объект.

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

Проекция объекта, полученная в результате томографического измерения под заданным углом. $\theta$, состоит из набора линейных интегралов (см. рис. 1). Набор многих таких проекций под разными углами, организованных в 2D, называется синограммой (см. рис. 3). В рентгеновской КТ линейный интеграл представляет собой общее ослабление луча рентгеновских лучей при его прохождении через объект по прямой линии. Как уже говорилось выше, полученное изображение представляет собой 2D (или 3D) модель коэффициента затухания . То есть мы хотим найти изображение $\mu (x,y)$. Самый простой и легкий способ визуализации метода сканирования — это система параллельного проецирования , использовавшаяся в первых сканерах. Для этого обсуждения мы рассматриваем данные, которые должны быть собраны как серия параллельных лучей в позиции $r$, поперек проекции под углом $\theta$. Это повторяется для разных углов. Затухание происходит экспоненциально в тканях :

$I=I_{0}\exp \left({-\int \mu (x,y)\,ds}\right)$

где $\mu (x,y)$ — коэффициент затухания как функция положения. Поэтому, как правило, общее затухание $p$ луча в положении $r$, на проекции под углом $\theta$, определяется линейным интегралом:

$p_{\theta }(r)=\ln \left({\frac {I}{I_{0}}}\right)=-\int \mu (x,y)\,ds$

Используя систему координат рисунка 1, значение $r$ на что точка $(x,y)$ будет проецироваться под углом $\theta$ дается:

$x\cos \theta +y\sin \theta =r\$

Таким образом, приведенное выше уравнение можно переписать как

$p_{\theta }(r)=\int _{-\infty }^{\infty }\int _{-\infty }^{\infty }f(x,y)\delta (x\cos \theta +y\sin \theta -r)\,dx\,dy$

где $f(x,y)$ представляет $\mu (x,y)$ и $\delta ()$ – дельта-функция Дирака . Эта функция известна как преобразование Радона (или синограмма ) 2D-объекта.

Преобразование Фурье проекции можно записать как

$P_{\theta }(\omega )=\int _{-\infty }^{\infty }\int _{-\infty }^{\infty }f(x,y)\exp[-j\omega (x\cos \theta +y\sin \theta )]\,dx\,dy=F(\Omega _{1},\Omega _{2})$ где $\Omega _{1}=\omega \cos \theta ,\Omega _{2}=\omega \sin \theta$

$P_{\theta }(\omega )$ представляет собой фрагмент двумерного преобразования Фурье $f(x,y)$ под углом $\theta$. Используя обратное преобразование Фурье , можно легко вывести формулу обратного преобразования Радона.

$f(x,y)={\frac {1}{2\pi }}\int \limits _{0}^{\pi }g_{\theta }(x\cos \theta +y\sin \theta )d\theta$

где $g_{\theta }(x\cos \theta +y\sin \theta )$ является производной Гильберта преобразования $p_{\theta }(r)$

Теоретически обратное преобразование Радона даст исходное изображение. Теорема о срезах проекций говорит нам, что если бы у нас было бесконечное количество одномерных проекций объекта, снятых под бесконечным числом углов, мы могли бы идеально восстановить исходный объект. $f(x,y)$. Однако на практике будет доступно лишь ограниченное число прогнозов.

Предполагая $f(x,y)$ имеет эффективный диаметр $d$ и желаемое разрешение $R_{s}$, эмпирическое правило для количества проекций, необходимых для реконструкции, равно $N>\pi d/R_{s}$

Алгоритмы реконструкции

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

Алгоритм восстановления Фурье-области

Реконструкцию можно произвести с помощью интерполяции. Предполагать $N$ прогнозы $f(x,y)$ генерируются под одинаково расположенными углами, каждый из которых производится с одинаковой частотой. Дискретное преобразование Фурье (ДПФ) в каждой проекции дает выборку в частотной области. Объединение всех проекций с частотной дискретизацией создает полярный растр в частотной области. Полярный растр разрежен, поэтому для заполнения неизвестных точек ДПФ используется интерполяция, а реконструкцию можно выполнить с помощью обратного дискретного преобразования Фурье . Производительность реконструкции можно улучшить за счет разработки методов изменения разреженности полярного растра, повышая эффективность интерполяции.

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

$\theta '={\frac {R_{0}}{\max\{|\cos \theta |,|\sin \theta |\}}}$

где $R_{0}$ самая высокая частота, подлежащая оценке.

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

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

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

В дискретной системе обратное преобразование Радона имеет вид

$f(x,y)={\frac {1}{2\pi }}\sum _{i=0}^{N-1}\Delta \theta _{i}g_{\theta _{i}}(x\cos \theta _{i}+y\sin \theta _{i})$

$g_{\theta }(t)=p_{\theta }(t)\cdot k(t)$

где $\Delta \theta$ - угловое расстояние между выступами и $k(t)$ представляет собой ядро ​​Радона с частотной характеристикой $|\omega |$.

Название «обратная проекция» происходит от того факта, что одномерная проекция должна быть отфильтрована с помощью одномерного ядра Радона (обратная проекция), чтобы получить двумерный сигнал. Используемый фильтр не содержит усиления по постоянному току, поэтому добавление смещения по постоянному току может быть желательным . Реконструкция с использованием обратной проекции обеспечивает лучшее разрешение, чем метод интерполяции, описанный выше. Однако он вызывает больший шум, поскольку фильтр склонен усиливать высокочастотный контент.

Алгоритм итерационной реконструкции

Основная статья: Итеративная реконструкция

Итерационный алгоритм требует больших вычислительных ресурсов, но позволяет включать априорную информацию о системе. $f(x,y)$.

Позволять $N$ быть числом проекций и $D_{i}$ быть оператором искажения для $i$проекция, снятая под углом $\theta _{i}$. $\{\lambda _{i}\}$ представляют собой набор параметров для оптимизации преобразования итераций.

$f_{0}(x,y)=\sum _{i=1}^{N}\lambda _{i}p_{\theta _{i}}(r)$

Веерная реконструкция Шеппа-Логана Фантома с другим расстоянием между датчиками. Меньшее расстояние между датчиками обеспечивает более точную реконструкцию. Рисунок был создан с помощью MATLAB.

$f_{k}(x,y)=f_{k-1}(x,y)+\sum _{i=1}^{N}\lambda _{i}[p_{\theta _{i}}(r)-D_{i}f_{k-1}(x,y)]$

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

Реконструкция веерного луча

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

Реконструкция глубокого обучения

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

В настоящее время методы глубокого обучения широко применяются для реконструкции изображений и достигли впечатляющих результатов в различных задачах реконструкции изображений, включая низкодозное шумоподавление, реконструкцию с разреженным изображением, томографию с ограниченным углом и уменьшение металлических артефактов. Отличный обзор можно найти в специальном выпуске. транзакции IEEE по медицинской визуализации. Одна группа алгоритмов реконструкции глубокого обучения применяет нейронные сети постобработки для достижения реконструкции изображений, где входные изображения реконструируются обычными методами реконструкции. Одним из примеров применения является уменьшение артефактов с использованием U-Net в томографии с ограниченным углом. Однако в изображении, восстановленном таким полностью управляемым данными методом, могут возникнуть неправильные структуры. как показано на рисунке. Таким образом, интеграция известных операторов в архитектуру нейронных сетей представляется полезной, как описано в концепции точного обучения. Например, прямую реконструкцию изображения по данным проекции можно изучить с помощью фильтрованной обратной проекции. Другой пример — создание нейронных сетей путем развертывания алгоритмов итеративной реконструкции. За исключением точного обучения, с использованием традиционных методов реконструкции с предварительной реконструкцией глубокого обучения. Это также альтернативный подход к улучшению качества изображения при реконструкции глубокого обучения.

Программное обеспечение для томографической реконструкции

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

Большинство коммерческих пакетов программного обеспечения, доступных для приобретения, ориентированы на обработку данных для настольных систем конусно-лучевой компьютерной томографии. Некоторые из этих программных пакетов включают Volume Graphics , InstaRecon , iTomography , Livermore Tomography Tools (LTT) и Cone Beam Software Tools (CST) .

Некоторые заслуживающие внимания примеры программного обеспечения для реконструкции с открытым исходным кодом включают в себя: Reconstruction Toolkit (RTK), КОНРАД, ТомоПи, набор инструментов АСТРА, ПИРО-НН, ОДЛ, ТИГР, и ЛЕАП.

Галерея

В галерее показан полный процесс томографии простого объекта и последующей томографической реконструкции на основе ART.

• 
  Рис. 2: Объект-фантом , два квадрата с кошачьими углами.
• 
  Рис. 3: Синограмма фантомного объекта (рис. 2), полученная по данным томографии. Было взято 50 проекционных срезов под углом 180 градусов, отобранных на равном расстоянии (только по совпадению ось X отмечает смещение на -50/50 единиц).
• 
  Рис.4: ART Томографическая реконструкция синограммы на рис.3 на основе , представленная в виде анимации процесса итеративной реконструкции. Исходный объект можно аппроксимативно реконструировать, поскольку полученное изображение имеет некоторые визуальные артефакты .

См. также

• Операция компьютерной томографии#Томографическая реконструкция
• Реконструкция конусной балки
• Промышленная компьютерная томография
• Промышленные томографические системы» ООО «

Дальнейшее чтение

• Авинаш Как и Малкольм Слейни (1988), Принципы компьютерной томографической визуализации, IEEE Press, ISBN   0-87942-198-3 .
• Брюянт, П. П. «Алгоритмы аналитической и итеративной реконструкции в ОФЭКТ» Журнал ядерной медицины 43 (10): 1343-1358, 2002 г.

Внешние ссылки

• Слейни, AC Как и Малкольм. «Принципы компьютерной томографии» . Slaney.org . Проверено 7 сентября 2018 г.
• Инструментарий Insight; программное обеспечение для томографической поддержки с открытым исходным кодом
• «TomoPy — документация TomoPy 1.1.3» . Tomopy.readthedocs.org . Проверено 7 сентября 2018 г.
• Набор инструментов ASTRA (Всемасштабная томографическая реконструкция Антверпена); очень гибкое и быстрое программное обеспечение с открытым исходным кодом для компьютерной томографической реконструкции
• НифтиРек; комплексное программное обеспечение для томографической реконструкции с открытым исходным кодом; Возможность написания сценариев Matlab и Python
• Инструмент томографической реконструкции и визуализации с открытым исходным кодом
• «ITS plc — Томография электрических процессов для промышленной визуализации» . Itoms.com . Проверено 7 сентября 2018 г.