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

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

Иллюстрации

Показать все

Реферат

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

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

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

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

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

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

В дифференциальной фазовой контрастной рентгеновской радиографии и томографии может возникнуть свертывание фазы. Если градиент фазового фронта находится вне диапазона [-π; π] на каждый период решетки, градиент свертывается в этот интервал. Эта ситуация может возникнуть, в частности, на краю объекта, например, вследствие большого перехода показателя преломления между воздухом и тканью. Свертывание фазы может возникнуть для умеренных шагов решетки, составляющих 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(E₀) при расчетной энергии E₀=20 кэВ составляет

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

где n=0, 2, 4, 6, 8. Кривая 40a относится к n=0, кривая 40b к n=2 и так далее. Вследствие выбора приведенной выше формулы для кривых 40a, 40b, 40c, 40d, 40e значения градиента фазы оценивает g(E₀) при E₀=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 показывает ожидаемое соотношение E₀/E вокруг градиента фазы 1 рад на период решетки. Точки данных или кривая 44b имеют отношение к истинному градиенту фазы (1+2π) на период решетки при 20 кэВ, которые свернуты в 1 рад на период решетки. Поскольку разность градиента фазы между 20 кэВ и, например, 22 кэВ основана на истинном градиенте, свернутый градиент фазы изменяется более быстро с энергией, что делает возможным удаление свертывания.

Поскольку градиенты фазы g(E₀) были выбраны таким образом, что они отличаются на 2π, свертывание приводит к тому, что при E₀ все отображенные (то есть, измеренные) градиенты 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, E₀ и E₁, соответствующие 18, 20 и 22 кэВ.

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

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

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

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

В соответствии с вариантом осуществления изобретения пиксель 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 повторяются для каждого дополнительного энергетического уровня E₀, E₁. Данные 58 изображения затухания могут быть получены только для одного энергетического уровня, например для основного энергетического уровня E₀.

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

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

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

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

В соответствии с вариантом осуществления изобретения способ содержит этап, на котором принимают данные 52 дифференциального фазового изображения, полученные на основании излучения на разных энергетических уровнях E_-1, E₀ и E₁, причем данные 52 дифференциального фазового изображения содержат пиксели 60, каждый пиксель 60 имеет значения 62a, 62b, 62c градиента фазы для каждого энергетического уровня E_-1, E₀ и E₁ и дополнительно значение 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, E₀ и E₁.

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

В соответствии с вариантом осуществления изобретения значения 42, 46 градиента фазы кривых 44a-44g свертывания определяются таким образом, что значения 42 градиента фазы на опорном энергетическом уровне E₀ равны. Затем посредством процессора 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, E₀ и E₁ кривой 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 градиента фазы при энергии E₀ обозначается как g₀ и вычисляется посредством приближения методом наименьших квадратов

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

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

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

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

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

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

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

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

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

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

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

Или

для большого значения m для приближения дискретного штрафа с помощью гладкой функции. В приведенном выше случае, когда число свертывания оценивается на основании изображения затухания как, например, имеющее порядок n₀, что эквивалентно тому, что 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 может потребовать, чтобы отношение сигнал-шум было ориентировочно достаточно большим, чтобы правильно различать разные точки данных. Более определенно, для измерения при расчетной энергии E₀ и последующего измерения при энергии E₀+ΔE градиенты фазы отличаются следующим образом:

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

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

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

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

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

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