СПОСОБ РАЗДЕЛЕНИЯ ИЗОБРАЖЕНИЙ ВОДЫ И ЖИРА В МАГНИТНО-РЕЗОНАНСНОЙ ТОМОГРАФИИ Российский патент 2015 года по МПК G01N24/08 

Описание патента на изобретение RU2544387C1

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

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

Так, например, для томографа с магнитом 0.3 Т Δω=45 Гц, для томографа на сверхпроводящем магните 1,5 Т Δω=225 Гц.

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

где Ii - i-e изображение,

W - сигнал от воды,

F - сигнал от жира,

φi=Δω* ti - фаза сигнала от жира относительно сигнала от воды, накопленная за время ti.

Из (1) видно, что при известных φi теоретически для разделения жира и воды достаточно двух измерений I1 и I2. Соответствующий метод получил название двухточечного метода Диксона (Two point Dixon) и является наиболее эффективным с точки зрения временных затрат. В простейшем варианте по методу Диксона собираются 2 изображения, в одном из которых жир и вода находятся в фазе (φ1=0°), а во втором - в противофазе (φ2=180°) (фиг.1).

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

На фиг.2:

1 - возбуждающий радиоимпульс,

2 - градиент выделения слоя,

3 - градиент чтения,

4 - эхо-сигнал с временем эха ТЕ1, формирующий изображение I1 (жир и вода в фазе),

5 - эхо-сигнал с временем эха ТЕ2 формирующий изображение I2 (жир и вода в противофазе).

Подставляя в (1) значения φi соответственно 0° и 180°, имеем:

Из (2) жир и вода легко разделяются:

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

В этом случае уравнения (2) принимают вид:

где φ0 - начальная фаза векторов W и F, накопленная водой и жиром из-за неоднородности поля за время между возбуждением системы радиоимпульсом и первым сбором, φ - набег фазы, накопленный из-за неоднородности поля за время между двумя сборами

Получается 2 уравнения с 4 неизвестными W, F, φ0, φ.

Однако, если учесть, что I1 и I2 - комплексные числа, то число уравнений расширяется до 4-х и система в принципе разрешима.

Для решения проблемы предложен ряд способов.

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

Известен также способ разделения изображений воды и жира [2], основанный на априорном знании карты неоднородности магнитного поля (В0 map). Карта поля строится заранее, например, с использованием однородного объекта, а затем набег фаз, связанный с неоднородностью поля, корректируется с учетом полученной карты неоднородностей [2].

Наиболее близким к предложенному является способ разделения изображений воды и жира двухточечным методом Диксона [3], включающий получение двух комплексных изображений I1 и I2 с различными временами эха, в которых сигналы от воды и жира находятся соответственно в фазе и в противофазе, вычисление значений фазы 2φ комплексного вектора I2=(I2I*1/|I1|)2 для каждого пиксела матрицы изображений, построение матрицы "развернутой" фазы 2φu в диапазоне главных значений - 180°…180°, например, методом слияния регионов, определение знака комплексного вектора Ie-iφu в каждом пикселе матрицы и формирование изображения по воде как полусуммы абсолютного значения изображения в фазе и изображения в противофазе, умноженного на знак Ie-iφu и изображения жира как полуразности абсолютного значения изображения в фазе и изображения в противофазе, умноженного на знак Ie-iφu [3].

Недостатком этого способа является высокая вероятность "перепутывания" результирующих пикселов, содержащих воду и жир, в условиях неоднородности поля магнита и наличия шумов измерений при определении фазы комплексных чисел I2=(I2I*1/|I1|)2.

Целью предлагаемого способа является повышение надежности правильной классификации изображений по воде и жиру. Для достижении этой цели в способе разделения изображений воды и жира в магнитно-резонансной томографии, включающем получение двух комплексных изображений I1 и I2 с различными временами эха, в которых сигналы от воды и жира находятся соответственно в фазе и в противофазе, вычисление значений фазы 2φ комплексного вектора I2=(I2I*1/|I1|)2 для каждого пиксела матриц изображений, построение матрицы "развернутой" фазы 2φu в диапазоне главных значений - 180°…180°, например, методом слияния регионов, определение знака комплексного вектора Ie-iφu в каждом пикселе матрицы и формирование изображения по воде как полусуммы абсолютного значения изображения в фазе и изображения в противофазе, умноженного на знак Ie-iφu и изображения жира как полуразности абсолютного значения изображения в фазе и изображения в противофазе, умноженного на знак Ie-iφu, оценивают усредненные градиенты изменения фазы полученных изображений жира и воды по формулам:

GF=(|I1|-|I2|)2/NF при Ie-iφu<0

Gw=(|I1|-|I2|)2/NW при Ie-iφu>0,

сравнивают значения GF и GW и, в случае, если GF<GW, пиксели изображений жира и воды обменивают местами.

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

I1=(W+F)eiφ0 жир и вода в фазе

и I2=(W-F)ei(φ0+φ) жир и вода в противофазе

где φ0 - начальная фаза векторов W и F, накопленная водой и жиром из-за неоднородности поля за время между возбуждением системы радиоимпульсом и первым сбором, φ - набег фазы, накопленный из-за неоднородности поля за время между двумя сборами

Для решения системы уравнений относительно W и F вводится вспомогательная комплексная величина

где * - знак комплексного сопряжения

Определить однозначно искомую фазу φ как фазу комплексного числа I нельзя из-за неопределенности знака W-F, так как заранее не известно, содержит ли данный пиксел больше воды или жира. Эта неоднозначность устраняется возведением обеих частей (6) в квадрат:

В принципе теперь фаза определяется как

Однако выражение (8) справедливо только при условии, что все значения угла 2φ лежат в диапазоне -180°<2φ<180°, что в реальности маловероятно при наличии неоднородности поля и шумов измерений эхо-сигналов.

Эту ситуацию для одномерного случая иллюстрирует фиг.3.

Как видно из фиг.3, при переходе через 180° фаза испытывает скачок, который надлежит сгладить, прежде чем определять фазу φ и решать уравнение (6).

Существует ряд методов двумерного "разворачивания" фазы, среди которых можно использовать, например, метод слияния регионов [4].

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

На фиг.4 слева показано изображение фазы 2φ, рассчитанной по формуле (8), а справа - изображение "развернутой" фазы φu, полученной методом слияния регионов.

Определив оценку фазы φu, из (6) можно получить оценку разности W-F

Так как левая часть (9) является действительным числом, то фаза комплексного числа Ie-iφu теоретически может принимать только 2 значения 0° или 180°.

Поэтому (9) можно переписать в виде:

где g принимает значения ±1 в зависимости от знака Ie-iφu

Сумма воды и жира также является действительной и положительной:

Из (10) и (11) находим искомые изображения W и F:

В реальности из-за наличия шумов коэффициент g не получается строго равным 1 или -1. Поэтому более гладкое решение дает аппроксимация:

Как видно из (12), при ошибке в определении знака коэффициента g происходит "перепутывание" воды и жира, что с вероятностью, близкой к 50%, происходит в условиях неоднородного поля и шумов измерений.

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

Усредненные градиенты изменения фазы оцениваются по исходным матрицам изображений "в фазе" I1 и в "противофазе" I2 и матрице коэффициентов g:

Иллюстрация к (14) представлена на фиг.5.

Порядок выполнение действий предлагаемым способом схематически представлен на фиг.6.

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

На фиг.7:

а - суммарное изображение воды и жира в противофазе,

b - суммарное изображение воды и жира в фазе,

с - изображение воды,

d - изображение жира.

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

Источники информации

[1] Патент ЕР 2414860 B1, "TWO-POINT DIXON TECHNIQUE WITH FLEXIBLE CHOICE OF ECHO TIMES".

[2] Патент ЕР 2610632 A1, "MRI with Dixon-type water/fat separation and prior knowledge about inhomogeneity of the main magnetic field".

[3] Bernstein, Zhou, King "Handbook-of-MRI-Pulse-Sequence". Part 17.3.1 Two point Dixon, ELSEVIER ACADEMIC PRESS, 2004.

[4] Salah Karout "Two-Dimensional Phase Unwrapping", General Engineering Research Institute (GERI), Liverpool John Moores University, 2007.

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

название год авторы номер документа
ДВУХТОЧЕЧНЫЙ МЕТОД ДИКСОНА С ГИБКИМ ВЫБОРОМ ВРЕМЕНИ ЭХА 2010
  • Брендель Бернард
  • Эггерс Хольгер
  • Дейндам Адри Й.
RU2521751C2
МАГНИТНО-РЕЗОНАНСНАЯ (MR) ВИЗУАЛИЗАЦИЯ С В1-ОТОБРАЖЕНИЕМ 2013
  • Нерке Кай
  • Бурнерт Петер
  • Граесслин Ингмар
  • Качер Ульрих
  • Лойсслер Кристоф
  • Эггерс Хольгер
RU2611082C2
СПОСОБ ДИАГНОСТИКИ САКРОИЛЕИТА НА ВЫСОКОПОЛЬНОМ МАГНИТНО-РЕЗОНАНСНОМ ТОМОГРАФЕ 2018
  • Завылова Ксения Александровна
  • Шахов Борис Евгеньевич
  • Сафонов Дмитрий Владимирович
  • Моровов Сергей Викторович
  • Завылов Дмитрий Михайлович
RU2712310C1
МАГНИТНО-РЕЗОНАНСНАЯ ВИЗУАЛИЗАЦИЯ С РАЗДЕЛЕНИЕМ ВОДА/ЖИР МЕТОДОМ ДИКСОНА 2018
  • Эггерс, Хольгер
RU2773486C2
МАГНИТНО-РЕЗОНАНСНАЯ ВИЗУАЛИЗАЦИЯ ХИМИЧЕСКИХ СОЕДИНЕНИЙ С ПОМОЩЬЮ СПЕКТРАЛЬНОЙ МОДЕЛИ 2011
  • Эггерс Хольгер
RU2552408C2
МАГНИТНО-РЕЗОНАНСНАЯ ТОМОГРАФИЯ С РАЗДЕЛЕНИЕМ ВОДЫ И ЖИРА ПО МЕТОДУ ДИКСОНА 2017
  • Эггерс Хольгер
RU2739479C2
СПОСОБ ДИФФЕРЕНЦИАЛЬНОЙ ДИАГНОСТИКИ РАКА ПОДЖЕЛУДОЧНОЙ ЖЕЛЕЗЫ 2019
  • Кит Олег Иванович
  • Франциянц Елена Михайловна
  • Иозефи Дмитрий Ярославович
  • Винидченко Михаил Александрович
RU2742744C1
ПАРАЛЛЕЛЬНАЯ МР-ТОМОГРАФИЯ С КАРТИРОВАНИЕМ ЧУВСТВИТЕЛЬНОСТИ С ПОМОЩЬЮ РЧ-КАТУШКИ 2016
  • Бернерт Петер
  • Нерке Кай
RU2683605C1
СПОСОБ ОЦЕНКИ ОБЪЕМА ЖИРОВОЙ ТКАНИ В ТЕЛЕ ЧЕЛОВЕКА 2008
  • Анисимов Николай Викторович
  • Пирогов Юрий Андреевич
RU2373840C1
СПОСОБ И ОПТИЧЕСКАЯ СИСТЕМА ДЛЯ ПОЛУЧЕНИЯ ТОМОГРАФИЧЕСКОГО РАСПРЕДЕЛЕНИЯ ВОЛНОВЫХ ФРОНТОВ ЭЛЕКТРОМАГНИТНЫХ ПОЛЕЙ 2018
  • Фернандес Вальдивия, Хуан Хосе
  • Трухильо Севилья, Хуан Мануэль
  • Гомес Карденес, Оскар
RU2733822C1

Иллюстрации к изобретению RU 2 544 387 C1

Реферат патента 2015 года СПОСОБ РАЗДЕЛЕНИЯ ИЗОБРАЖЕНИЙ ВОДЫ И ЖИРА В МАГНИТНО-РЕЗОНАНСНОЙ ТОМОГРАФИИ

Использование: для разделения изображений воды и жира в магнитно-резонансной томографии. Сущность изобретения заключается в том, что осуществляют получение двух комплексных изображений I1 и I2 с различными временами эха, в которых сигналы от воды и жира находятся соответственно в фазе и в противофазе, вычисление значений фазы 2φ комплексного вектора I 2 = ( I 2 I 1 * / | I 1 | ) 2 для каждого пиксела матриц изображений, построение матрицы "развернутой" фазы 2φ и в диапазоне главных значений -180°…180° определение знака комплексного вектора Ie-iφu в каждом пикселе матрицы, формирование изображения по воде как полусуммы абсолютного значения изображения в фазе и изображения в противофазе, умноженного на знак Ie-iφu, изображения жира как полуразности абсолютного значения изображения в фазе и изображения в противофазе, умноженного на знак Ie-iφu, при этом оценивают усредненные градиенты изменения фазы полученных изображений жира и воды по формулам:

GF=(|I1|-|I2|)2/NF при Ie-iφu<0

GW=(|I1|-|I2|)2/NW при Ie-iφu<0, сравнивают значения GF и GW и, в случае, если GF<GW, пиксели изображений жира и воды обменивают местами. Технический результат: повышение надежности правильной классификации изображений по воде и жиру. 7 ил.

Формула изобретения RU 2 544 387 C1

Способ разделения изображений воды и жира в магнитно-резонансной томографии, включающий получение двух комплексных изображений I1 и I2 с различными временами эха, в которых сигналы от воды и жира находятся соответственно в фазе и в противофазе, вычисление значений фазы 2φ комплексного вектора I2=(I2I*1/|I1|)2 для каждого пиксела матриц изображений, построение матрицы "развернутой" фазы 2φ и в диапазоне главных значений -180°…180°, определение знака комплексного вектора I е-iφu в каждом пикселе матрицы и формирование изображения по воде как полусуммы абсолютного значения изображения в фазе и изображения в противофазе, умноженного на знак I e-iφu, и изображения жира как полуразности абсолютного значения изображения в фазе и изображения в противофазе, умноженного на знак I е-iφu, отличающийся тем, что оценивают усредненные градиенты изменения фазы полученных изображений жира и воды по формулам:
GF=(|I1|-|I2|)2/NF при I е-iφu<0
GW=(|I1|-|I2|))2W при I е-iφu>0, сравнивают значения GF и GW и, в случае, если GF<GW, пиксели изображений жира и воды обменивают местами.

Документы, цитированные в отчете о поиске Патент 2015 года RU2544387C1

Bernstein, Zhou, King "Handbook-of-MRI-Pulse-Sequence", Part 17.3.1 Two point Dixon, ELSEVIER ACADEMIC PRESS, 2004
RU 2011143800A, 10.05.2013
СПОСОБ ОЦЕНКИ ОБЪЕМА ЖИРОВОЙ ТКАНИ В ТЕЛЕ ЧЕЛОВЕКА 2008
  • Анисимов Николай Викторович
  • Пирогов Юрий Андреевич
RU2373840C1
Сигнатурный анализатор 1989
  • Терехов Владимир Георгиевич
  • Савченко Олег Владимирович
  • Козлов Игорь Владимирович
SU1635184A1
СПОСОБ ОБРАБОТКИ ЛОЖА ЖЕЛЧНОГО ПУЗЫРЯ ПОСЛЕ ХОЛЕЦИСТЭКТОМИИ ПРИ ДЕСТРУКТИВНЫХ ФОРМАХ ОСТРОГО КАЛЬКУЛЕЗНОГО ХОЛЕЦИСТИТА 2009
  • Журихин Александр Викторович
  • Одишелашвили Гиви Доментиевич
  • Курбанов Рауль Мурадханович
RU2414860C1
US 2003062900A1, 03.04.2003

RU 2 544 387 C1

Авторы

Мишкинис Александр Борисович

Пистрак Григорий Михайлович

Бочкарев Серафим Николаевич

Даты

2015-03-20Публикация

2013-12-24Подача