Способ реконструкции изображений продольных срезов объекта



Способ реконструкции изображений продольных срезов объекта
Способ реконструкции изображений продольных срезов объекта
Способ реконструкции изображений продольных срезов объекта
Способ реконструкции изображений продольных срезов объекта
Способ реконструкции изображений продольных срезов объекта
Способ реконструкции изображений продольных срезов объекта
Способ реконструкции изображений продольных срезов объекта
Способ реконструкции изображений продольных срезов объекта
Способ реконструкции изображений продольных срезов объекта
Способ реконструкции изображений продольных срезов объекта
Способ реконструкции изображений продольных срезов объекта
Способ реконструкции изображений продольных срезов объекта

 

A61B6/00 - Приборы для радиодиагностики, например комбинированные с оборудованием для радиотерапии (рентгеноконтрастные препараты A61K 49/04; препараты, содержащие радиоактивные вещества A61K 51/00; радиотерапия как таковая A61N 5/00; приборы для измерения интенсивности излучения, применяемые в ядерной медицине, например измерение радиоактивности живого организма G01T 1/161; аппараты для получения рентгеновских снимков G03B 42/02; способы фотографирования в рентгеновских лучах G03C 5/16; облучающие приборы G21K; рентгеновские приборы и их схемы H05G 1/00)

Владельцы патента RU 2559688:

Закрытое акционерное общество "ИМПУЛЬС" (RU)

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

 

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

Описание предшествующего уровня техники

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

Из уровня техники известны способы получения томографических срезов объекта путем преобразования рентгенографических снимков, сделанных под разными углами.

Одним из самых первых и самых простых способов получения изображения томографических срезов является способ смещения и суммирования Shift-And-Add (SAA), который заключается в суммировании сдвинутых на определенное расстояние проекций объекта (James Т Dobbins III and Devon J Godfrey, Digital x-ray tomosynthesis: current state of the art and clinical potential, 2003 Phys. Med. Biol. 48 R72, section 3.2; doi: 10.1088/0031-9155/48/19/R01,).

Иллюстрацией способа служит фиг.1, где на фиг.1а изображены три положения 1, 2, 3 рентгеновской трубки и проекции круга из плоскости A и треугольника из плоскости B на детекторе; на фиг.1б приведены изображения структур плоскости A или B: четкие изображения структур плоскости A или B могут быть получены с помощью сдвига проекций и их суммирования, структуры вне интересующей плоскости будут размазаны по всему изображению (т.е. размыты).

На фиг.1 видно, как посредством сдвига и суммирования можно получить четкое изображение, например, треугольника из плоскости B: первую проекцию сдвигают вправо относительно центра детектора, центральную проекцию оставляют без изменений, а третью проекцию сдвигают влево. Суммирование трех сдвинутых проекций приводит к наложению треугольников друг на друга, тогда как круги оказываются «размазанными» по всему изображению.

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

Среди способов реконструкции изображений известен также способ обратных проекций (матриц) Matrix Inversion Tomosynthesis (MITS) (патент US №490320, опубл. 20.02.1990), который можно описать следующим образом.

Положим, что объект состоит из n томографических срезов. Обозначим через ti изображение i-го томографического среза, полученное методом SAA, где i=1, …, n. С другой стороны каждый томографический срез может быть описан в виде:

где sj - изображение действительной структуры объекта в томографическом срезе j, где j=1, …, n, а fij - функция размытия для изображения действительных структур среза j в срезе i.

Следует отметить, что ti и sj в системе уравнений (1) являются полноразмерными изображениями с размерами (X, Y) пикселей. Однако систему уравнений (1) следует решать непосредственно для каждого пиксела с координатами (x, y).

Переход в Фурье-область обеспечивает простую замену оператора свертки на оператор умножения. Для F, S, T - Фурье-спектров f, s и t соответственно получаем

Ввиду сделанного замечания о решении системы уравнений непосредственно для каждого пиксела с координатами (x, y) полученную систему уравнений (2) следует решать относительно пространственных частот kx и ky.

В матричном виде систему уравнений (2) можно записать так:

или

T=F·S.

Применив матричную алгебру, получаем

S=F-1·Т,

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

s=FT-1(S)=FT-1(F-1·Т).

Матрица F в данном случае является плохо обусловленной на низких пространственных частотах kx и ky. На нулевых частотах она является вырожденной. Поэтому для качественного восстановления низких и нулевых частот томографических срезов приходится задействовать дополнительные обработки, что является существенным недостатком данного способа. Также к недостаткам следует отнести и тот факт, что способ основан на реконструкциях SAA, что увеличивает не только продолжительность расчетов, но и количество возможных артефактов.

Раскрытие изобретения

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

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

Технический результат в заявляемом способе реконструкции изображений продольных срезов объекта, который заключается в том, что с помощью томографической системы заданной геометрии выполняют множество проекций объекта, снятых под разными углами, решают систему линейных уравнений с помощью обратных матриц, неизвестными в которой выступают Фурье-образы S изображений реконструируемых срезов, и затем восстанавливают по ним изображения продольных срезов объекта, достигается тем, что в качестве коэффициентов системы линейных уравнений выступают Фурье-образы D дельта-функций от величин смещений для каждой проекции и каждого среза, при этом величину смещения для каждой проекции и каждого среза рассчитывают, исходя из геометрии системы; в качестве свободных членов выступают Фурье-образы P проекций объекта, а Фурье-образы S изображений реконструируемых срезов определяют методом наименьших квадратов с использованием регуляризации по Тихонову по формуле

S=(DTD+λ·I)-1·DTP,

где λ - параметр регуляризации, I - единичная матрица, DT - транспонированная матрица D.

Осуществление изобретения

Указанные преимущества, а также особенности изобретения поясняются лучшим вариантом его выполнения со ссылками на прилагаемые фигуры. На фигурах изображены:

фиг.1 - схематичное изображение способа смещения и суммирования Shift-And-Add (SAA) - предшествующий уровень техники;

фиг.2 - схематичное изображение способа по заявляемому техническому решению;

фиг.3 - блок-схема алгоритма реконструкции продольных срезов объекта по изобретению.

Томографическая система, реализующая способ, представляет собой установку для линейной томографии, основные принципы работы которой проиллюстрированы на фиг.2. Объект 1 размещают на столе 2. Над столом располагают рентгеновскую трубку 3, которая может двигаться вдоль стола и поворачиваться вокруг оси вращения. Детектор 4, расположенный под столом 2, перемещается согласно геометрии в сторону, противоположную направлению движения трубки так, чтобы центральный луч L от трубки проходил через ось вращения в центр детектора. Тогда за один проход трубки слева направо получают n проекций объекта 1, снятых под разными углами (на фиг.2 показана геометрия съемки для получения трех проекций объекта). В начале работы (шаг 1) трубка 3 и детектор 4 расположены симметрично относительно оси вращения. На данном шаге получают Проекцию 1. Далее (шаг 2) трубка и детектор смещаются в противоположные стороны, получают Проекцию 2 (центральную). На следующем шаге (шаг 3) трубка и детектор снова смещаются, в результате чего получают Проекцию 3. Каждая проекция представляет собой рентгенографический снимок размером dx на dy пикселей, т.е. зависимость яркости от координат (x, y).

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

pi - проекция объекта, где i=1, …, n, n - число проекций;

sd - расстояние от рентгеновской трубки 3 до детектора 4 в центральной проекции;

sc - расстояние от трубки 3 до оси вращения в центральной проекции;

a i - угол наклона трубки 3 к вертикали в точке оси вращения для каждой проекции pi;

nx - ширина детектора 4, мм;

dx - ширина изображения, пиксель;

sj - томографический срез, где j=1, …, m, m - число срезов;

hj - высота над поверхностью детектора, на которой требуется реконструировать томографический срез sj.

Расстояния sd от трубки до детектора и sc от трубки до оси вращения в центральной проекции, угол наклона a i трубки к вертикали в точке оси вращения для каждой проекции pi, ширина детектора nx в мм, ширина изображения dx в пикселях и сами проекции pi являются входными параметрами алгоритма.

Пусть требуется реконструировать m томографических срезов. Каждый томографический срез sj объекта задается высотой hj над поверхностью детектора и представляет собой, так же, как и любая проекция, рентгенографический снимок размером (dx, dy) пикселей.

Каждый пиксел pi(x,y) изображения проекции pi может быть представлен в виде суммы пикселей sj(x,y) срезов sj, смещенных относительно друг друга на величину Δi,j:

p1(x,y)=δ(Δ1,1)⊗s1(x,y)+δ(Δ1,2)⊗s2(x,y)+…+δ(Δ1,m)⊗sm(x,y),

где δ - дельта-функция, аргументом которой является величина смещения Δi,j, ⊗ - оператор свертки.

Величина смещения Δi,j рассчитывается для каждого среза sj (на высоте hj) и для каждой проекции pi (снятой под углом a i) по формуле

Для работы непосредственно с изображением необходимо перевести смещение в пиксели:

.

Таким образом, получаем систему уравнений (координаты (x,y) опущены)

Переход в Фурье-область обеспечивает простую замену оператора свертки на оператор умножения. Для Di,j(kx,ky), Pi(kx,ky), Sj(kx,ky) - Фурье-спектров δ(Δi,j(x,y)), pi(x,y), sj(x,y) соответственно, где kx и ky - пространственные частоты, - система уравнений (4) имеет вид (частоты (kx,ky) опущены)

Перепишем полученную систему уравнений (5) в матричном виде:

или

P=D·S.

Применив матричную алгебру, получаем

а применив обратное преобразование Фурье к обеим частям выражения, получаем искомые томографические срезы объекта:

s=FT-1(S)=FT-1(D-1·P).

Строго говоря, решение данной системы уравнения возможно только при n=m, то есть для квадратной матрицы D. Для того чтобы обойти данное ограничение, т.е. находить решение, когда число проекций не равно числу реконструируемых срезов n≠m, применяется метод наименьших квадратов, который заменяет выражение (6) на:

где DT - транспонированная матрица D.

Матрица D является плохо обусловленной на низких пространственных частотах kx и ky. На нулевых частотах она является вырожденной. Поэтому для качественного восстановления низких и нулевых частот томографических срезов применяется метод регуляризации по Тихонову (Кирьянов Д.В., Кирьянова Е.Н. Вычислительная физика. - М.: Полибук Мультимедиа, 2006. - 352 с.: ил), который заменяет выражение (7) на:

S=(DTD+λ·I)-1·DTP,

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

Алгоритм описывается блок-схемой, изображенной на фиг.3. После получения массива проекций p1..n объекта и параметров sd, sc, a 1…n, nx, dx, h1…m, λ рассчитывают матрицу смещений Δ1…n, 1…m для каждой проекции из массива p1…n и каждого томографического среза из массива s1…m на определенной высоте из массива h1…m по приведенной выше формуле (3). Здесь индекс вида 1…n обозначает одномерный массив (вектор) размерностью n, индекс вида (1…n, 1…m) обозначает двумерный массив размерностями (n, m.)

Далее выполняют преобразование Фурье для каждой проекции из массива p1…n только по горизонтальным частотам (т.е. в направлении сдвига проекций) и получают вектор Фурье-образов P1..n, где каждый Фурье-образ есть набор пространственных частот.

Затем выполняют два цикла: i - внешний цикл по строкам (в блок-схеме nf - частота Найквиста), и j - внутренний цикл (в блок-схеме ny - количество столбцов в проекциях p1..n):

1. в цикле по i в зависимости от пространственной частоты kxi формируется матрица A1…n, 1…m, которая соответствует Фурье-преобразованию дельта-функции от смещения Δ1…n, 1…m:

A1…n, 1…m=exp(Δ1…n, 1…m·kxi).

Исходя из метода наименьших квадратов и регуляризации по Тихонову, рассчитывают матрицу Ar

,

где - транспонированная матрица A, I1…n, 1…m - единичная матрица.

2. в цикле по j выбирается вектор B1…n=Pij, 2…n, соответствующий вектору частот с индексами i и j всех проекций от 1 до n и производится расчет вектора частот с индексами i и j всех томографических срезов от 1 до m:

.

Таким образом получают все частоты пространственного спектра томографического среза S1…m. Выполняя обратное преобразование Фурье над S1…m (только по горизонтальным частотам, т.е. в направлении сдвига проекций) получают непосредственно изображение томографического среза sm.

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

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

Способ реконструкции изображений продольных срезов объекта, заключающийся в том, что с помощью томографической системы заданной геометрии получают множество проекций объекта, снятых под разными углами, решают систему линейных уравнений с помощью обратных матриц, неизвестными в которой являются Фурье-образы S изображений реконструируемых срезов, и затем восстанавливают по ним изображения продольных срезов объекта, отличающийся тем, что в качестве коэффициентов системы линейных уравнений выступают Фурье-образы D дельта-функций от величин смещений для каждой проекции и каждого среза, при этом величину смещения для каждой проекции и каждого среза рассчитывают, исходя из геометрии системы; в качестве свободных членов выступают Фурье-образы P проекций объекта, а Фурье-образы S изображений реконструируемых срезов определяют методом наименьших квадратов с использованием регуляризации по Тихонову по формуле
S=(DTD+λ·I)-1DТP,
где λ - параметр регуляризации, I - единичная матрица, DТ - транспонированная матрица D.



 

Похожие патенты:

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

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

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

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

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

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

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

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

Изобретение относится к компьютерной технике, а именно к системам генерирования изображения с иллюзией. Техническим результатом является автоматическое генерирование иллюзии с изображением из произвольного изображения.

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

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

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

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

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

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

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

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

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

Изобретение относится к медицинской технике, а именно к средствам формирования функциональных изображений. Способ содержит получение первого изображения накопления первого контрастного вещества в ткани пациента, не являющейся объектом исследования, при этом первое изображение генерируется на основе первых данных от первого средства формирования изображений, получение второго изображения накопления второго контрастного вещества в исследуемой ткани пациента и ткани пациента, не являющейся объектом исследования, при этом второе изображение генерируется на основе вторых данных от другого второго средства формирования изображений, генерирование первой маски изображения на основе первого изображения, генерирование первого изображения особенности на основе второго изображения и первой маски изображения и отображение первого изображения особенности, которое не включает в себя накопление контрастного вещества в исследуемой ткани, не накапливающей контрастное вещество. Устройство формирования изображений содержит генератор масок и блок выделения особенностей, генерирующий изображение особенности. Использование изобретения позволяет снизить ошибки при оценке результатов терапии. 2 н. и 12 з.п. ф-лы, 6 ил.
Наверх