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



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

 


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

КОНИНКЛЕЙКЕ ФИЛИПС Н.В. (NL)

Группа изобретений относится к способам и устройствам для формирования дифференциальных фазовых контрастных изображений. Техническим результатом является обеспечение возможности корректировки количества артефактов в данных изображения. Для корректировки данных 52 дифференциального фазового изображения принимают данные 52 дифференциального фазового изображения, полученные с помощью излучения на разных энергетических уровнях, причем данные 52 дифференциального фазового изображения содержат пиксели 60, каждый пиксель 60 имеет значения 62a, 62b, 62c градиента фазы для каждого энергетического уровня. После этого определяют зависящее от энергии поведение значений 62a, 62b, 62c градиента фазы пикселя 60 и определяют скорректированное значение 68 градиента фазы для пикселя 60 на основании значений 62a, 62b, 62c градиента фазы пикселя 60 и модели для зависимости от энергии значений 62a, 62b, 62c градиента фазы. 5 н. и 10 з.п. ф-лы, 5 ил.

 

ОБЛАСТЬ ТЕХНИКИ, К КОТОРОЙ ОТНОСИТСЯ ИЗОБРЕТЕНИЕ

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

УРОВЕНЬ ТЕХНИКИ

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

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

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

В дифференциальной фазовой контрастной рентгеновской радиографии и томографии может возникнуть свертывание фазы. Если градиент фазового фронта находится вне диапазона [-π; π] на каждый период решетки, градиент свертывается в этот интервал. Эта ситуация может возникнуть, в частности, на краю объекта, например, вследствие большого перехода показателя преломления между воздухом и тканью. Свертывание фазы может возникнуть для умеренных шагов решетки, составляющих 2 мкм для фазовой решетки, уже для объектов с размером меньше 1 мм. В частности, в дифференциальной фазовой контрастной рентгеновской томографии это может привести к сильным артефактам перекрытия в воссозданном изображении.

РАСКРЫТИЕ ИЗОБРЕТЕНИЯ

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

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

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

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

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

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

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

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

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

Дополнительный аспект изобретения имеет отношение к системе формирования дифференциальных фазовых изображений, например к системе компьютерной томографии (CT).

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

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

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

Эти и другие аспекты изобретения будут очевидны и разобраны со ссылкой на описанные далее варианты осуществления.

КРАТКОЕ ОПИСАНИЕ ЧЕРТЕЖЕЙ

Ниже более подробно описаны варианты осуществления настоящего изобретения со ссылкой на приложенные чертежи.

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

Фиг. 2 показывает схему с абсолютными значениями градиента фазы в соответствии с вариантом осуществления изобретения.

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

Фиг. 4 схематично показывает данные изображения в соответствии с вариантом осуществления изобретения.

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

В основном на фигурах идентичным частям даны одинаковые обозначения.

ОСУЩЕСТВЛЕНИЕ ИЗОБРЕТЕНИЯ

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

Источник 10 излучения может содержать некогерентный источник 16 рентгеновского излучения, например рентгеновскую трубку 18, и решетку 20 источника для достижения пространственной когерентности пучка. Источник 10 излучения может быть выполнен с возможностью формировать пространственный когерентный пучок излучения.

Детектор 12 может содержать фазовую решетку 22, абсорбционную решетку 24 и элемент 26 детектора рентгеновского излучения, выполненный с возможностью обнаруживать данные изображения на основании рентгеновского излучения, излученного от источника излучения через интересующий объект 28.

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

Фазовая решетка 22 служит в качестве сдвигающего фазу разделителя пучка, который преобразовывает изменение фазового фронта пучка рентгеновского излучения в модуляцию интенсивности, то есть интерференционную картину с характерным масштабом длин в половину шага фазовой решетки 22 в местоположении абсорбционной решетки 24. Абсорбционная решетка 24 формирует муаровую интерференционную картину на сновании рентгеновского излучения, выходящего из фазовой решетки 22. Муаровая интерференционная картина на элементе 26 детектора содержит информацию о фазовом сдвиге отклоненного и смещенного по фазе рентгеновского излучения после прохождения и объекта 28 и фазовой решетки 22.

Контроллер 16 содержит процессор 30 для регистрации или сбора данных изображения от элемента 26 датчика и для управления мотором 32 для перемещения абсорбционной решетки 24 в направлении, ортогональном по отношению к протяженности полос решеток 22, 24. Вследствие перемещения абсорбционной решетки 26 разные на элементе 26 датчика формируются муаровые рисунки. Эти разные муаровые картины могут быть зарегистрированы процессором 30 и преобразованы в данные дифференциального фазового изображения, которые могут быть сохранены в памяти 34 контроллера 14.

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

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

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

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

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

Может быть возможно, что система 10 представляет собой рентгеновскую систему 10 дифференциальной фазовой контрастной радиографии и/или рентгеновскую систему 10 дифференциальной фазой контрастной томографии. В последнем случае система 10 может вычислять сечения или трехмерные представления объекта 18 на основании данных изображения, которые были собраны от различных направлений, посредством перемещения устройства 12 источника излучения, детектора 14 и объекта 28 друг относительно друга.

Сформированные радиографические или томографические изображения могут быть отображены на дисплее 36 системы 10.

Вследствие пространственного изменения показателя преломления δ(x, y, z) объекта 28 два разных пучка рентгеновских лучей, проходящих через объект 28, могут подвергнуться разным фазовым изменениям, которые могут быть обнаружены с помощью детектора 12, и данные дифференциального фазового изображения могут быть вычислены на основании обнаруженных необработанных данных с помощью процессора 30 и сохранены в памяти 34.

Однако с помощью детектора 14 и последующего процесса преобразования необработанных данных не может быть однозначно определен градиент фазы. Другими словами, когда истинный градиент фазы находится вне интервала [-π; π], определенный градиент фазы свертывается в этот интервал, то есть может быть определен только истинный градиент фазы по модулю 2π.

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

Зависимость действительной части показателя преломления δ может быть описана с использованием потенциального закона

Поскольку фаза волнового фронта φ пучка излучения вдоль направления z составляет приблизительно

абсолютный градиент фазы g волнового фронта имеет следующую зависимость от энергии, то есть длины волны

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

Фиг. 2 показывает схему с истинными значениями градиента фазы для разных энергий и показывает зависимость градиента фазы от энергии. Ось X схемы изображает энергию в кэВ, ось Y - абсолютный градиент фазы в рад. Фиг. 2 иллюстрирует зависимость градиента фазы для разных значений. Проиллюстрированы пять кривых 40a, 40b, 40c, 40d и 40d, которые выбраны таким образом, что градиента g(E0) при расчетной энергии E0=20 кэВ составляет

g(E0)=(1+nπ)/(шаг фазовой решетки 22)

где n=0, 2, 4, 6, 8. Кривая 40a относится к n=0, кривая 40b к n=2 и так далее. Вследствие выбора приведенной выше формулы для кривых 40a, 40b, 40c, 40d, 40e значения градиента фазы оценивает g(E0) при E0=20 кэВ отличаются на 2π.

Однако обнаруженные (измеренные) и определенные данные дифференциального фазового изображения содержат не абсолютные (истинные) значения градиента фазы, а значения, которые свернуты в интервал [-π; π].

Это изображено на Фиг. 3, которая показывает схему со свернутыми значениями градиента фазы. Ось X схемы изображает энергию в кэВ, ось Y - свернутый градиент фазы в рад. Свернутые значения градиента фазы, относящиеся к одному и тому же значению для n, соединены линиями с целью иллюстрации.

На Фиг. 3 истинные градиенты фазы 40a, 40b, 40c, 40d, 40e, показанные на Фиг. 2, свернуты в интервал [-π; π] на каждый период решетки в относительные градиенты фазы 44a, 44b, 44c, 44d, 44e, 44f, 44g. Кривые 44a, 44b, 44c, 44d, 44e, 44f, 44g относятся к числам n=0, 2, 4, 6, 8, 10, 12, соответственно.

На схеме точки данных или кривая 44a имеют отношение к случаю, когда свертывание не произошло. Кривая 44a показывает ожидаемое соотношение E0/E вокруг градиента фазы 1 рад на период решетки. Точки данных или кривая 44b имеют отношение к истинному градиенту фазы (1+2π) на период решетки при 20 кэВ, которые свернуты в 1 рад на период решетки. Поскольку разность градиента фазы между 20 кэВ и, например, 22 кэВ основана на истинном градиенте, свернутый градиент фазы изменяется более быстро с энергией, что делает возможным удаление свертывания.

Поскольку градиенты фазы g(E0) были выбраны таким образом, что они отличаются на 2π, свертывание приводит к тому, что при E0 все отображенные (то есть, измеренные) градиенты 42 фазы являются одинаковыми. Даже при том, что кривые 44a, 44b, 44c, 44d, 44e, 44f, 44g свертывания (которые представляют измеренные значения) испытывают свертывание фазы также и при других энергиях для n>0, ясно видно, что разные числа n свертывания дают в результате сильно отличающиеся распределения свернутых градиентов 46 фазы в пределах показанного диапазона энергий.

В частности, все кривые 44a, 44b, 44c, 44d, 44e, 44f, 44g имеют разные свернутые градиенты 46 фазы на разных энергетических уровнях (например, 16 кэВ, 18 кэВ, 22 кэВ, 24 кэВ) для нескольких кривых градиента 44a, 44b, 44c, 44d, 44e, 44f, 44g, которые отображены на одно и то же значение 42 при расчетной или опорной энергии (например, 20 кэВ).

Даже трех точек данных 42, 46 (и, таким образом, трех измерений на трех разных энергетических уровнях) может быть достаточно, чтобы различить все случаи (то есть кривые 44a, 44b, 44c, 44d, 44e, 44f, 44g), показанные на Фиг. 3, например точки данных энергетических уровней E-1, E0 и E1, соответствующие 18, 20 и 22 кэВ.

Фиг. 4 схематично показывает данные 50 изображения, которые могут быть приняты, обработаны и сохранены в памяти 34 контроллера 16 и, в частности, посредством 30.

Данные 50 изображения могут содержать данные 52 дифференциального фазового изображения, содержащие данные 54a, 54b, 54c дифференциального фазового изображения, собранные на разных энергетических уровнях E-1, E0 и E1 соответственно, и/или могут содержать данные 50 изображения затухания.

Данные 50 изображения составлены из пикселей 60, и каждый пиксель может иметь значение 62a, 62b, 62c градиента фазы для каждого энергетического уровня E-1, E0 и E1 и/или значение 64 затухания.

В соответствии с вариантом осуществления изобретения пиксель 60 имеет значения 62a, 62b, 62c градиента фазы, соответствующие по меньшей мере трем энергетическим уровням E-1, E0 и E1.

В соответствии с вариантом осуществления изобретения пиксель 60 имеет значение 64 затухания.

Данные 50 изображения могут содержать скорректированные данные 66 дифференциального фазового изображения, которые могут быть сформированы на основании данных 52, 58 изображения посредством способа, описанного выше и далее. В частности, каждый пиксель 60 может содержать скорректированное значение 68 градиента фазы, которое может быть определено на основании значений 62a, 62b, 62c, 64 посредством способа.

Следует отметить, что данные 50 изображения не обязательно должны представлять двухмерные изображения. Могут быть возможны другие представления, например линейные данные. Другими словами, элемент 26 датчика может представлять собой двухмерный детектором или линейный детектор.

Фиг. 5 показывает блок-схему последовательности операций для способа формирования и корректировки данных 52 дифференциального фазового изображения.

На этапе S10 источником 12 излучения управляет контроллер 16 для формирования излучения на первом энергетическом уровне E-1, например рентгеновского излучения с энергией 18 кэВ. Излучение проникает через объект 28 и падает на детектор 26.

В соответствии с вариантом осуществления изобретения излучение представляет собой электромагнитное излучение, например рентгеновское излучение.

На этапе S12 контроллер 16 управляет детектором 14 и, в частности, решеткой 24, в результате чего элементом 26 детектора обнаруживаются разные муаровые интерференционные картины.

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

На этапе S14 муаровые интерференционные картины преобразовываются в данные 52a дифференциального фазового изображения и данные 58 изображения затухания процессором 30. Эти данные 52a, 58 изображения могут быть сохранены в памяти 34.

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

Этапы S10-S14 повторяются для каждого дополнительного энергетического уровня E0, E1. Данные 58 изображения затухания могут быть получены только для одного энергетического уровня, например для основного энергетического уровня E0.

В соответствии с вариантом осуществления изобретения способ содержит этап, на котором формируют излучение на разных энергетических уровнях E-1, E0, E1.

В соответствии с вариантом осуществления изобретения энергетические уровни E-1, E0 и E1 излучения содержат опорный энергетический уровень E0 и два соседних энергетических уровня E-1, E1, отличающиеся приблизительно на 8-12%, например, на 10% от опорного энергетического уровня E0, например уровни составляют 18, 20 и 22 кэВ.

Может быть возможно, что излучение на одном энергетическом уровне формируется в разные моменты времени и излучение обнаруживается в эти разные моменты времени, как объяснено выше. Однако также может быть возможно, что излучение с разными энергетическими уровнями формируется и обнаруживается одновременно, например, с помощью мультихромного источника 12 излучения и спектроскопического детектора 14.

На этапе S16 данные 50 изображения принимаются процессором 30 из памяти 34 (например, на основе пикселей).

В соответствии с вариантом осуществления изобретения способ содержит этап, на котором принимают данные 52 дифференциального фазового изображения, полученные на основании излучения на разных энергетических уровнях E-1, E0 и E1, причем данные 52 дифференциального фазового изображения содержат пиксели 60, каждый пиксель 60 имеет значения 62a, 62b, 62c градиента фазы для каждого энергетического уровня E-1, E0 и E1 и дополнительно значение 64 затухания.

Затем следующие этапы S18-S22 могут быть повторены для каждого пикселя 60 из данных 50 изображения.

На этапе S18 зависящее от энергии поведение значений 62a, 62b, 62c градиента фазы выбранного пикселя 60 определяется при помощи кривых 44a-44g свертывания, показанных на Фиг. 3.

В соответствии с вариантом осуществления изобретения способ содержит этап, на котором определяют зависящее от энергии поведение 44a-44g значений 62a, 62b, 62c градиента фазы пикселя 60.

Кривые 44a-44g свертывания могут быть представлены в памяти 30 с помощью точек 41, 46 данных. Например, для каждой кривой 44a-44g могут быть сохранены по меньшей мере три предварительно вычисленных значения 42, 46 данных.

В соответствии с вариантом осуществления изобретения кривая 44a-44g свертывания представлена значениями 42, 46 градиента фазы, соответствующими энергетическим уровням E-1, E0 и E1.

В соответствии с вариантом осуществления изобретения значения 42, 46 градиента фазы кривой 44a-44g свертывания предварительно вычислены и/или сохранены в контроллере 16, например, в памяти 34.

В соответствии с вариантом осуществления изобретения значения 42, 46 градиента фазы кривых 44a-44g свертывания определяются таким образом, что значения 42 градиента фазы на опорном энергетическом уровне E0 равны. Затем посредством процессора 30 значения 62a, 62b, 62c градиента фазы могут быть отображены на каждую кривую свертывания, сохраненную в контроллере 16. Например, для каждой кривой 44a-44g свертывания вычисляется значение наилучшего приближения, которое указывает, насколько значения 62a, 62b, 62c градиента фазы пикселя соответствуют кривой 44a-44g свертывания. В частности, значение 62a, 62b, 62c градиента фазы пикселя 60 может быть отображено на соответствующие значения 42, 46 градиента фазы кривой свертывания. После этого из отображенных значений 62a, 62b, 62c, 42, 46 может быть вычислено численное приближение.

В соответствии с вариантом осуществления изобретения зависящее от энергии поведение определяется посредством приближения значений 62a, 62b, 62c градиента фазы пикселя 60 к множеству кривых 44a-44g свертывания фазы.

В соответствии с вариантом осуществления изобретения значение 42, 46 градиента фазы энергетического уровня E-1, E0 и E1 кривой 44a-44g свертывания приближается с помощью значения 62a, 62b, 62c градиента фазы пикселя 60 на этом энергетическом уровне.

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

На этапе S18 процессор 30 может определить число свертывания n для пикселя 60 при помощи наилучшего приближения значений 62a, 62b, 62c градиента фазы к кривым 44a-44g свертывания. Поскольку каждая кривая свертывания может соответствовать числу свертывания, число свертывания n для пикселя может быть выбрано в качестве числа свертывания кривой свертывания с наилучшим приближением.

В соответствии с вариантом осуществления изобретения способ содержит этап, на котором определяют число свертывания n пикселя 60 на основании зависящего от энергии поведения.

В соответствии с вариантом осуществления изобретения число свертывания n пикселя 60 определяется посредством определения кривой 44a-44g свертывания с наилучшим приближением приближенных значений 62a, 62b, 62c градиента фазы пикселя 60 и посредством выбора числа свертывания, соответствующего определенной кривой 44a-44g свертывания.

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

В соответствии с вариантом осуществления изобретения значение 42 градиента фазы при энергии E0 обозначается как g0 и вычисляется посредством приближения методом наименьших квадратов

где σ2 - дисперсии измеренных свернутых значений градиента фазы gi при энергиях Ei, w обозначает операцию свертывания, и || ||π обозначает специальную операцию расстояния, а именно

В итоге для каждого значимого числа свертывания n (например, 0, 2, 4,…, 12) может быть сделано численное приближение для получения наилучшего приближения g(E0) в предположении, что фаза свернута n/2 раз при энергии E0. Наконец, среди этих результатов может быть выбран тот, который имеет наилучшее общее приближение (например, с использованием среднеквадратичного показателя).

На этапе S20 процессор 30 определяет абсолютное или скорректированное значение 68 градиента фазы на основании определенного числа свертывания n. Например, к значению 62b градиента фазы на опорном энергетическом уровне E0 может быть добавлено nπ для вычисления значения 68.

В соответствии с вариантом осуществления изобретения способ содержит этап, на котором определяют скорректированное значение 68 градиента фазы для пикселя 60 от значений 62a, 62b, 62c градиента фазы пикселя 60 и число свертывания пикселя 60.

В соответствии с вариантом осуществления изобретения скорректированное значение 68 градиента фазы определяется посредством смещения значения 62b градиента фазы пикселя 60, соответствующего опорному энергетическому уровню E0, на определенное число свертывания пикселя 60.

Когда множество кривых 44a-44g свертывания предварительно не выбрано, способ может дать возможность различать только достаточно малые значения для возможных значений обертывания n. Однако если n становится чрезмерно большим, результаты снова могут стать неоднозначными. В этой ситуации изображение 58 затухания может обеспечить грубую оценку числа свертывания n.

Изображение 58 затухания может затем использоваться для разрешения оставшейся неопределенности.

Факультативно на этапе S18 с помощью процессора 30 на основании данных изображения 58 затухания может быть определено значение градиента в пикселе 60. На основании значения градиента могут быть определены значимые числа свертывания (например, n= от 20 до 30). Например, в памяти 34 может храниться функция или таблица, с помощью которых значение градиента может быть отображено на предполагаемую область чисел свертывания.

Затем множество кривых свертывания ограничивается кривыми свертывания, соответствующими числам свертывания в этой области.

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

Где функция штрафа P используется для определения допустимого диапазона от -gmax до gmax для результата, например, посредством использования дискретного штрафа

Или

для большого значения m для приближения дискретного штрафа с помощью гладкой функции. В приведенном выше случае, когда число свертывания оценивается на основании изображения затухания как, например, имеющее порядок n0, что эквивалентно тому, что g имеет порядок nπ, штраф может быть выбран как

Или

В соответствии с вариантом осуществления изобретения способ содержит этап, на котором определяют градиент значения 64 затухания в пикселе 60.

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

Затем на этапе S20 факультативно для приближения значений 62a, 62b, 62c градиента фазы пикселя 60 используются только выбранные кривые свертывания.

В соответствии с вариантом осуществления изобретения способ содержит этап, на котором выполняют приближение значений 62a, 62b, 62c градиента фазы пикселя 60 к множеству кривых 44a-44g свертывания с помощью соответствующих чисел свертывания в предполагаемой области.

Если система 10 является системой 10 томографии, система может сформировать данные изображения, показывающие сечения или трехмерные представления объекта 28. Данные изображения могут быть отображены на дисплее 36 контроллера 16.

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

Таким образом, на этапе S24 контроллер 16 может изменить оптическую ось A относительно объекта, например, посредством поворота размещения источника 12 излучения и детектора 14 вокруг объекта 18 и может повторить сбор данных 50 изображения относительно измененного направления.

После того как собрано и скорректировано достаточно данных 50 изображения, контроллер 16 и, в частности, процессор 30 могут сформировать данные томографического изображения на этапе S26.

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

Следует понимать, что этапы на Фиг. 5 не обязательно должны выполняться в порядке, описанном относительно Фиг. 5.

Также следует отметить, что различие между разными значениями 62a, 62b, 62c, 42, 46 может потребовать, чтобы отношение сигнал-шум было ориентировочно достаточно большим, чтобы правильно различать разные точки данных. Более определенно, для измерения при расчетной энергии E0 и последующего измерения при энергии E0+ΔE градиенты фазы отличаются следующим образом:

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

В зависимости от объекта 28 интервал между энергетическими уровнями E-1, E0 и E1, используемыми в способе, может быть настроен для достижения такой возможности.

Описанный выше и далее способ может представлять собой компьютерную программу, исполняемую в процессоре 30 и хранящуюся в памяти 34. Компьютерная программа может быть сохранена на машиночитаемом носителе, таком как гибкий диск, жесткий диск, устройство хранения на универсальной последовательной шине (USB), оперативное запоминающее устройство (RAM; ОЗУ), постоянное запоминающее устройство (ROM; ПЗУ) и стираемое программируемое постоянное запоминающее устройство (EPROM; СППЗУ). Машиночитаемый носитель также может представлять собой сеть передачи данных, например Интернет, которая позволяет загружать программный код.

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

1. Способ корректировки данных (52) дифференциального фазового изображения, причем способ содержит этапы, на которых:
принимают данные (52) дифференциального фазового изображения, полученные с помощью излучения на разных энергетических уровнях (E-1, E0, E1), причем данные (52) дифференциального фазового изображения содержат пиксели (60), причем каждый пиксель (60) имеет значение (62a, 62b, 62c) градиента фазы для каждого энергетического уровня (E-1, E0, E1);
определяют зависящее от энергии поведение (44a-44g) значений (62a, 62b, 62c) градиента фазы пикселя (60);
определяют скорректированное значение (68) градиента фазы для пикселя (60) на основании значений (62a, 62b, 62c) градиента фазы пикселя (60) и модели для зависимости от энергии значений (62a, 62b, 62c) градиента фазы.

2. Способ по п. 1,
в котором зависящее от энергии поведение определяется посредством приближения значений (62a, 62b, 62c) градиента фазы пикселя (60) к множеству кривых (44a-44g) свертывания фазы, причем каждая кривая свертывания соответствует числу свертывания;
при этом число свертывания пикселя (60) определяется посредством определения кривой (44a-44g) свертывания наилучшим приближением приближенных значений (62a, 62b, 62c) градиента фазы пикселя (60) и посредством выбора числа свертывания, соответствующего определенной кривой (44a-44g) свертывания.

3. Способ по п. 2, в котором кривая (44a-44g) свертывания представлена значениями (42, 46) градиента фазы, соответствующими энергетическим уровням (E-1, E0, E1);
при этом значение (42, 46) градиента фазы энергетического уровня (E-1, E0, E1) кривой (44a-44g) свертывания приближается с помощью значения (62a, 62b, 62c) градиента фазы пикселя (60) на этом энергетическом уровне.

4. Способ по п. 3,
в котором значения (42, 46) градиента фазы кривых (44a-44g) свертывания определены таким образом, что значения (42) градиента фазы на опорном энергетическом уровне (E0) равны.

5. Способ по одному из пп. 1-4,
в котором пиксель (60) имеет значение (64) затухания;
причем способ содержит этапы, на которых:
определяют градиент значения (64) затухания в пикселе (60);
выбирают предполагаемую область чисел свертывания;
выполняют приближение значений (62a, 62b, 62c) градиента фазы пикселя (60) к множеству кривых (44a-44g) свертывания, соответствующих числам свертывания в предполагаемой области.

6. Способ по одному из пп. 1-4,
в котором скорректированное значение (68) градиента фазы определяется посредством смещения значения (62b) градиента фазы пикселя (60), соответствующего опорному энергетическому уровню (E0), с определенным числом свертывания пикселя (60).

7. Способ по одному из пп. 1-4,
в котором пиксель имеет значения (62a, 62b, 62c) градиента фазы, соответствующие по меньшей мере трем энергетическим уровням (E-1, E0, E1).

8. Способ по одному из пп. 1-4,
в котором излучение представляет собой электромагнитное излучение.

9. Способ по одному из пп. 1-4,
в котором энергетические уровни E-1, E0, E1 излучения содержат опорный энергетический уровень (E0) и два соседних энергетических уровня (E-1, E1), отличающиеся на 8-12% от опорного энергетического уровня (E0).

10. Способ формирования скорректированных данных (66) дифференциального фазового изображения, причем способ содержит этапы, на которых:
формируют излучение на разных энергетических уровнях (E-1, E0, E1);
обнаруживают сформированное излучение, проникающее через интересующий объект (28);
получают данные (52) дифференциального фазового изображения на основании обнаруженного излучения;
исполняют этапы способа по одному из пп. 1-9 со сформированными данными (52) дифференциального фазового изображения.

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

12. Машиночитаемый носитель (34), на котором сохранена компьютерная программа для корректировки данных (52) дифференциального фазового изображения, которая при ее исполнении процессором (30) выполнена с возможностью выполнять этапы способа по одному из пп. 1-9.

13. Машиночитаемый носитель (34), на котором сохранена компьютерная программа для корректировки данных (52) дифференциального фазового изображения, которая при ее исполнении процессором (30) выполнена с возможностью выполнять этапы способа по одному из пп. 10-11.

14. Система (10) формирования дифференциальных фазовых изображений, содержащая:
источник (12) излучения;
детектор (14); и
контроллер (16);
причем источник (12) излучения выполнен с возможностью формировать излучение разных энергетических уровней;
причем детектор (14) выполнен с возможностью обнаруживать данные (52) дифференциального фазового изображения интересующего объекта (28), через который проникает излучение;
причем контроллер (16) выполнен с возможностью выполнять способ по одному из пп. 1-11.

15. Система (10) формирования дифференциальных фазовых изображений по п. 14,
причем система формирования дифференциальных фазовых изображений представляет собой систему (10) рентгеновской компьютерной томографии (CT).



 

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

Использование: для реконструкции рентгеновской двухэнергетической компьютерной томографии. Сущность изобретения заключается в том, что способ реконструкции рентгеновской двухэнергетической CT согласно настоящему изобретению содержит: (a) оценку энергетического спектра и создание двухэнергетической таблицы поиска; (b) сбор данных высокой энергии и данных низкой энергии системы формирования изображений двухэнергетической CT с использованием детектора системы формирования изображений двухэнергетической CT; (c) получение изображений проекции и масштабированных изображений и согласно полученным данным высокой энергии и данным низкой энергии ; (d) реконструкцию масштабированного изображения с использованием первого условия ограничения кусочной гладкости и, тем самым, получение изображения электронной плотности; и (e) реконструкцию масштабированного изображения с использованием второго условия ограничения кусочной гладкости и, тем самым, получение изображения эквивалентного атомного номера.
Наверх