МАГНИТНО-РЕЗОНАНСНАЯ ВИЗУАЛИЗАЦИЯ С ПОДАВЛЕНИЕМ АРТЕФАКТОВ ПОТОКА
MAGNETIC RESONANCE IMAGING WITH SUPPRESSION OF FLOW ARTEFACTS

ОБЛАСТЬ ТЕХНИКИ, К КОТОРОЙ ОТНОСИТСЯ ИЗОБРЕТЕНИЕ
Изобретение относится к области магнитно-резонансной визуализации. Оно касается способа МР-визуализации участка тела, помещенного в объем исследования МР-устройства. Изобретение также относится к МР-устройству и компьютерной программе, работающей в МР-устройстве.
УРОВЕНЬ ТЕХНИКИ
Способы формирования МР-изображений, в которых используется взаимодействие между магнитными полями и ядерными спинами с целью построения двумерных или трехмерных изображений широко применяются в настоящее время, в частности в области медицинской диагностики, поскольку для визуализации мягких тканей они во многих отношениях превосходят другие способы визуализации, не требуют ионизирующего излучения и, как правило, являются неинвазивными.
Согласно МР-способу, в общем, тело пациента, которое требуется обследовать, располагают в сильном, однородном магнитном поле Bо, направление которого одновременно определяет ось (обычно ось z) координатной системы, на которой основывается измерение. Магнитное поле Bосоздает различные уровни энергии для индивидуальных ядерных спинов в зависимости от напряженности магнитного поля, которые далее могут возбуждаться (спиновый резонанс) путем наложения переменного электромагнитного поля (РЧ-поле) определенной частоты (т.н. Ларморовой частоты или МР-частоты). С макроскопической точки зрения распределение индивидуальных ядерных спинов создает общую намагниченность, которая может отклоняться от состояния равновесия путем наложения электромагнитного импульса соответствующей частоты (РЧ-импульса), когда силовые линии магнитного поля Bo проходят перпендикулярно z-оси, так что вектор намагниченности совершает прецессионное движение вокруг z-оси. Прецессионное движение описывает поверхность конуса, угол при вершине которого называют углом отклонения вектора намагниченности. Величина угла отклонения вектора намагниченности зависит от силы и продолжительности приложенного электромагнитного импульса. В случае т.н. 90°-импульса спины отклоняются от z-оси к трансверсальной плоскости (угол отклонения вектора намагниченности составляет 90°).
1
2
По окончании РЧ-импульса намагниченность возвращается в исходное состояние равновесия, при котором снова устанавливается намагниченность в направлении оси z согласно первой временной постоянной T (время спин-решеточной или продольной релаксации), а намагниченность в направлении, перпендикулярном z-направлению, релаксирует согласно второй временной постоянной T (время спин-спиновой или поперечной релаксации). Изменение намагниченности может детектироваться с помощью приемных РЧ катушек, расположенных и сориентированных в объеме исследования МР-устройства так, что изменение намагниченности измеряется в направлении, перпендикулярном z-оси. Затухание поперечной намагниченности сопровождается, после того как приложен, например, 90°-импульс, переходом ядерных спинов (индуцированным локальными неоднородностями магнитного поля) из упорядоченного состояния с одинаковой фазой в состояние, при котором все фазовые углы равномерно распределены (дефазировка). Дефазировка может быть скомпенсирована с помощью рефокусирующего импульса (например, 180°-импульса). Это порождает эхо-сигнал (спиновое эхо) в приемных катушках.
Для осуществления пространственного разрешения в теле на однородное магнитное поле Bонакладываются постоянные градиенты магнитного поля, направленные вдоль трех главных осей, что приводит к линейной пространственной зависимости частоты спинового резонанса. Сигнал, принятый в приемных катушках, в этом случае содержит компоненты различных частот, которые можно связать с различными местоположениями в теле. Данные сигнала, полученные посредством приемных катушек, соответствуют пространственно-частотной области и носят название данных k-пространства. В число данных k-пространства обычно входят множество линий, полученных с использованием различного фазового кодирования. Каждая линия оцифровывается путем сбора множества образцов. Набор данных k-пространства преобразуется в МР-изображение посредством преобразования Фурье.
В реферате ISMRM-2000 (стр.1829) публикации «Effects of gradient moment nulling in 3D half-Fourier FCE bright blood imaging», Y. Machida и другие, раскрыта технология визуализации артериального кровотока с использованием быстрого спин-эха (FSE). В известной FSE-последовательности обнуление момента градиента достигается путем приложения дополнительной градиентной составляющей, приводящей к нулевому моменту градиента (первого порядка).
РАСКРЫТИЕ ИЗОБРЕТЕНИЯ
В МР-визуализации часто требуется получить информацию об относительном вкладе различных химических веществ, таких как вода и жир, в общий сигнал, либо для подавления вклада некоторых из них, либо для отдельного или совместного анализа вклада каждого из них. Эти вклады можно рассчитать, если объединить информацию двух или более соответствующих эхо-сигналов, полученных в различное время появления эхо-сигналов. Это можно рассматривать как кодирование химического сдвига, при котором определяется и кодируется некоторая дополнительная величина, а именно величина химического сдвига, путем получения пары изображений в немного отличающиеся времена появления эхо-сигналов. В частности, для разделения вода-жир такого типа эксперименты часто называют измерениями по Диксону. Путем визуализация по Диксону или визуализации вода/жир по Диксону можно выполнить разделение вода-жир с помощью расчета вкладов от воды и жира из двух или более соответствующих эхо-сигналов, полученных в различное время появления эхо-сигналов. Вообще такое разделение возможно, поскольку существует известная разность частот прецессии водорода в жире и воде. В простейшей форме изображения воды и жира генерируются либо путем сложения, либо путем вычитания совпадающих по фазе и не совпадающих по фазе массивов данных. Однако такое разделение различных веществ на основе кодирования химического сдвига не ограничивается только веществами вода/жир. Могут также рассматриваться другие вещества с другими химическими сдвигами.
Известно, что может иметь место множество артефактов в МР-изображении, полученном из некоторой области обследуемого организма, в которой присутствуют явления потока, например артериального кровотока. Два наиболее поразительных эффекта, связанных с потоком, - снижение интенсивности МР-сигнала вследствие некогерентной добавки фаз индивидуальных спинов в заданный элемент изображения (пиксель/воксель), а также образование множества артефактов, связанных с двоением изображения, вследствие пульсации кровотока в ходе цикла сердечного сокращения. Поскольку оба упомянутых эффекта пропорциональны 1-му моменту времени используемых переключаемых градиентов магнитного поля, соответствующие артефакты можно уменьшить посредством компенсации потока путем управления переключением градиентов, так чтобы 1-й момент был по существу равен нулю. Такой вид компенсации потока иногда в литературе упоминается как рефазировка момента градиента (gradient moment rephasing), рефокусировка градиентного движения (gradient motion refocusing) или устранения момента градиента (gradient moment cancellation).
В то время как вышеупомянутые технологии позволяют преимущественно подавить отрицательное влияние потока в МР-изображениях, они влекут за собой другие проблемы. Главным образом, упомянутые технологии рефазировки момента градиента обычно требуют более продолжительного времени появления эхо-сигналов и времени повторения. Это, в частности, практически не позволяет их использовать в ангиографии «первого прохода» с применением контрастного вещества.
Из вышесказанного ясно, что существует потребность в усовершенствованной технологии МР-визуализации. Таким образом, задача изобретения заключается в создании способа, позволяющего эффективно компенсировать артефакты потока, в особенности для ангиографии «первого прохода» с применением контрастного вещества в сочетании с разделением вода/жир по Диксону.
Согласно изобретению раскрыт способ МР-визуализации тела, помещенного в объем исследования МР-устройства. Способ содержит этапы: a) генерирования, по меньшей мере, двух градиентных эхо-сигналов в два различных момента времени появления эха путем подвергания участка тела воздействию визуализирующей последовательности РЧ-импульсов и переключаемых градиентов магнитного поля, при этом 0-й момент считываемого градиента магнитного поля по существу исчезает во время появления первого градиентного эха, 1-й момент считываемого градиента отличен от нуля во время появления первого градиентного эха, в то время как и 0-й, и 1-й моменты считываемого градиента магнитного поля по существу исчезают во время появления второго градиентного эха;
b) сбора градиентных эхо-сигналов;
c) повторения этапов a) и b) для множества этапов фазового кодирования;
d) реконструкции первого МР-изображения из градиентных эхо-сигналов первого градиентного эха, а также второго МР-изображения из градиентных эхо-сигналов второго градиентного эха; а также
e) идентификации артефактов, связанных с двоением изображения, на первом и/или втором МР-изображении путем сравнения первого и второго МР-изображений.
Согласно изобретению осуществляется сбор, по меньшей мере, двух независимых наборов данных сигнала, при этом первый набор данных сигнала содержит градиентные эхо-сигналы первого градиентного эха, а второй набор данных сигнала содержит градиентные эхо-сигналы второго градиентного эха. В изобретении используется эффект рефазировки момента градиента, который по своей природе подавляет большинство артефактов потока во втором наборе данных сигнала. Таким образом, различие между первым и вторым наборами данных сигнала позволяет идентифицировать артефакты, связанные с двоением изображения, в первом наборе данных сигнала, обусловленные потоком, и соответственно ввести поправку на эти артефакты. Таким образом, изобретение позволяет устранить артефакты, связанные с двоением изображения, без продления времени появления эха и времени повторения при сборе МР-изображений.
Например, первое МР-изображение реконструируется из градиентных эхо-сигналов первого градиентного эха, в то время как второе МР-изображение реконструируется из градиентных эхо-сигналов второго градиентного эха. Артефакты, связанные с двоением изображения, далее можно идентифицировать на первом и/или втором МР-изображении путем сравнения первого и второго МР-изображений. Например, артефакты, связанные с двоением изображения, можно идентифицировать на основе потери и/или усиления локальной интенсивности на первом МР-изображении по сравнению со вторым МР-изображением. После того как артефакты, связанные с двоением изображения, идентифицированы на МР-изображении, их можно непосредственно устранить на первом МР-изображении, например восстановив потерянную локальную интенсивность или ограничив избыточную локальную интенсивность на первом МР-изображении.
После устранения артефактов, связанных с двоением изображения, на первом МР-изображении в предпочтительном варианте осуществления изобретения можно рассчитать, по меньшей мере, одно конечное МР-изображение, объединив первое и второе МР-изображения, например, согласно вышеупомянутой технологии Диксона.
Следует отметить, что термины "первое градиентное эхо" и "второе градиентное эхо" следует толковать в самом широком смысле в свете изобретения, при этом каждый термин относится к любой из последовательностей двух или более градиентных эхо-сигналов, сгенерированных визуализирующей последовательностью. В частности, термины "первое градиентное эхо" и "второе градиентное эхо" не всегда подразумевают, что два соответствующих градиентных эхо-сигнала генерируются в конкретной временной последовательности. При этом "первое градиентное эхо" не обязательно должно быть начальным эхо-сигналом в последовательности, состоящей из нескольких градиентных эхо-сигналов. Необязательно также, чтобы "первое градиентное эхо" и "второе градиентное эхо" были эхо-сигналами, следующими друг за другом.
Градиентные эхо-сигналы, собранные согласно изобретению, могут содержать вклады, по меньшей мере, от двух химических веществ, обладающих неодинаковыми МР-спектрами, при этом моменты времени появления эхо-сигналов выбираются так, что вклады сигналов, по меньшей мере, от двух химических веществ более совпадают по фазе во время второго градиентного эха, чем во время первого градиентного эха. Визуализация по Диксону может применяться для разделения вкладов сигналов, по меньшей мере, двух химических веществ (например, воды и жира) так, что, по меньшей мере, одно конечное МР-изображение содержит вклады только от одного из химических веществ. Подобные способы по Диксону уже использовались в МР-ангиографии для устранения вкладов сигналов от жира и, таким образом, для усиления контраста между сосудистой системой и фоновой тканью. Способы по Диксону - перспективные претенденты для замены способов вычитания, поскольку они уже используются для той же цели. Как пояснялось выше, способы по Диксону основываются на сборе, по меньшей мере, двух эхо-сигналов в различные времена появления эхо-сигналов, при этом вклады сигнала от воды и жира более совпадают по фазе во время одного из двух эхо-сигналов, чем во время другого эхо-сигнала. Способ по Диксону предпочтительно может использоваться в сочетании с подходом к компенсации потока по изобретению, который также основывается на сборе (по меньшей мере) двух эхо-сигналов в различные моменты времени появления эхо-сигналов. Способы по Диксону, в общем, не требуют большего совпадения по фазе вкладов сигнала от двух химических веществ во время второго градиентного эха, чем во время первого градиентного эха. Они способны разделить воду и жир независимо от этого условия. Однако подавление артефакта потока, предложенное в изобретении, получает выгоду от этого условия. Это объясняется тем, что происходит рефазировка сигналов от двух химических веществ в отношении времени второго градиента, т.е. величина суммарного сигнала либо остается постоянной (если только одно из химических веществ присутствует в соответствующем пикселе/вокселе), либо она увеличивается (в соответствующем пикселе/вокселе присутствуют оба химических вещества). При этом если величина суммарного сигнала падает, это можно интерпретировать как указание на артефакт потока согласно изобретению.
Следует отметить, однако, что способ по изобретению можно также использовать независимо от способов по Диксону, т.е. во всех случаях, когда артефакты потока следует устранить.
Способ по изобретению может быть также распространен на случай более двух эхо-сигналов. В этом случае МР-изображения, реконструированные из МР-сигналов четных градиентных эхо-сигналов, суммируются для избавления от артефактов, связанных с двоением изображения, благодаря рефазировке момента градиента. Артефакты, связанные с двоением изображения, будут присутствовать на МР-изображениях, реконструированных из МР-сигналов нечетных эхо-сигналов. Это различие может быть переведено в расчетные весовые значения для оценки совместимости в качестве основы для итеративного установления совместимости. Таким образом, артефакты, связанные с двоением изображения, могут подавляться особенно эффективно.
Помимо этого, известные процедуры для компенсации вихревых токов, устраняющие изменение фазы МР-сигнала между нечетными и четными эхо-сигналами, могут использоваться для способа по изобретению с целью устранения сдвига в амплитуде МР-сигнала между нечетными и четными градиентными эхо-сигналами.
Вышеописанный способ по изобретению еще не касался вопроса снижения общей интенсивности изображения в месте существования потока на первом МР-изображении. Это снижение интенсивности может приводить к потере сигнала в изображении жира при комбинировании первого и второго МР-изображений согласно подходу Диксона. Однако, если это необходимо, кровеносный сосуд, в котором наблюдается поток, может детектироваться в рамках изображения воды, при этом интенсивность сигнала в изображении жира может быть восстановлена в соответствии с изображением воды во всех затронутых пикселях или вокселях.
Следует отметить, что термин "химические вещества" следует понимать в широком смысле в контексте изобретения как химическое вещество любого рода или ядро любого рода, обладающее МР-свойствами. В простом примере осуществляется сбор МР-сигналов двух химических веществ, при этом химические вещества представляют собой протоны в "химических композициях" воды и жира. В более сложном примере спектральная модель с множеством пиков в действительности описывает ядра во множестве различных химических композиций, присутствующих в известных относительных количествах.
Описанный способ по изобретению может осуществляться с помощью МР-устройства, включающего в себя, по меньшей мере, одну главную магнитную катушку для генерирования однородного постоянного магнитного поля Bo в пределах объема исследования, множество градиентных катушек для генерирования переключаемых градиентов магнитного поля в различных пространственных направлениях в пределах объема исследования, по меньшей мере, одну РЧ катушку для тела, предназначенную для генерирования РЧ-импульсов в пределах объема исследования и/или для приема МР-сигналов от тела пациента, помещенного в объем исследования, блок управления для управления временной последовательностью РЧ-импульсов и переключаемыми градиентами магнитного поля, а также блок реконструкции изображения для реконструкции МР-изображения из полученных МР-сигналов. Способ по изобретению может быть реализован путем соответствующего программирования блока реконструкции изображения и/или блока управления МР-устройства.
Способ по изобретению предпочтительно может осуществляться на большинстве МР-устройств, применяемых в клинической практике в настоящее время. Для этого достаточно использовать компьютерную программу, с помощью которой МР-устройство управляется так, чтобы выполнить вышеуказанные этапы способа по изобретению. Компьютерная программа может находиться либо на носителе данных, либо присутствовать в сети передачи данных, чтобы загружаться для установки в блоке управления МР-устройства.
КРАТКОЕ ОПИСАНИЕ ЧЕРТЕЖЕЙ
Прилагаемые чертежи раскрывают предпочтительные варианты осуществления настоящего изобретения. Следует, однако, понимать, что чертежи представлены лишь в целях иллюстрации, но не для обозначения границ изобретения. На чертежах:
на Фиг. 1 показано МР-устройство для осуществления способа по изобретению;
на Фиг. 2 проиллюстрировано переключение считываемого градиента для генерирования градиентных эхо-сигналов согласно изобретению;
на Фиг. 3 показаны первое и второе МР-изображения, реконструированные согласно изобретению, совместно с разностным изображением первого и второго МР-изображений;
на Фиг. 4 показаны изображения воды и жира, полученные традиционным способом из первого и второго МР-изображений, представленных на Фиг. 3;
на Фиг. 5 показаны изображения воды и жира, полученные из первого и второго МР-изображений, представленных на Фиг. 3, с устраненными артефактами, связанными с двоением изображения, согласно изобретению.
ОСУЩЕСТВЛЕНИЕ ИЗОБРЕТЕНИЯ
На Фиг. 1 показано МР-устройство 1. Устройство содержит сверхпроводящие или резистивные главные магнитные катушки 2, так чтобы создавалось по существу однородное, постоянное во времени основное магнитное поле Bo вдоль z-оси по всему объему исследования. Устройство дополнительно содержит набор шиммирующих катушек 2′ (1-го, 2-го и в соответствующих случаях 3-го порядка), при этом ток, протекающий через индивидуальные шиммирующие катушки набора 2′, может контролироваться с целью минимизации отклонений Bo в пределах объема исследования.
Система для создания магнитного резонанса и управления им прикладывает серию РЧ-импульсов и переключаемые градиенты магнитного поля для опрокидывания или возбуждения ядерно-магнитных спинов, индуцирования магнитного резонанса, рефокусировки магнитного резонанса, управления магнитным резонансом, пространственного или иного кодирования магнитного резонанса, насыщения спинов и т.п. для проведения МР-визуализации.
В частности, усилитель 3 градиентных импульсов прикладывает импульсы тока к выбранным градиентным катушкам 4, 5 и 6 для всего тела вдоль x, y и z-осей объема исследования. Цифровой передатчик 7 радиочастот передает РЧ-импульсы или пакеты импульсов через переключатель 8 приема-передачи на РЧ катушку 9 для тела с целью передачи РЧ-импульсов в объем исследования. Обычная МР-визуализирующая последовательность состоит из пакета сегментов РЧ-импульсов малой продолжительности, которые совместно с любыми приложенными градиентами магнитного поля позволяют избирательно управлять ядерным магнитным резонансом. РЧ-импульсы используются для насыщения, возбуждения резонанса, инвертирования намагниченности, рефокусировки резонанса либо управления резонансом, а также выбора участка тела 10, располагаемого в объеме исследования. МР-сигналы также принимаются РЧ катушкой 9 для тела.
Для генерирования МР-изображений ограниченных областей тела 10 посредством параллельной визуализации набор локальных матричных РЧ катушек 11, 12, 13 располагают смежно с областью, выбранной для визуализации. Матричные катушки 11, 12, 13 могут использоваться для приема МР-сигналов, индуцированных передачей РЧ сигнала катушки для тела.
Получаемые в результате МР-сигналы принимаются РЧ катушкой 9 для тела и/или матричными РЧ катушками 11, 12, 13 и демодулируются приемником 14, который предпочтительно включает в себя предварительный усилитель (не показан). Приемник 14 соединен с РЧ катушками 9, 11, 12 и 13 посредством переключателя 8 прием-передача.
Главный компьютер 15 управляет шиммирующими катушками 2′, а также усилителем 3 градиентных импульсов и передатчиком 7 для генерирования любой из МР-визуализирующих последовательностей, например, эхо-планарной визуализации (EPI), объемной эхо-визуализации, визуализации с использованием градиентного и спинового эха, визуализации с использованием быстрого спин-эха и т.п. Для выбранной последовательности приемник 14 получает одну или множество МР-линий данных в быстрой последовательности вслед за каждым RF импульсом возбуждения. Система 16 сбора данных выполняет аналого-цифровое преобразование принятых сигналов и преобразует каждую МР-линию данных в цифровой формат, пригодный для дальнейшей обработки данных. В современных МР-устройствах система 16 сбора данных представляет собой отдельный компьютер, адаптированный для сбора исходных данных изображения.
И наконец, исходные данные изображения реконструируются в представление изображения с помощью процессора 17 реконструкции изображений, использующего преобразование Фурье или иные соответствующие алгоритмы реконструкции, например SENSE или SMASH. МР-изображение может представлять планарный срез тела пациента, набор параллельных планарных срезов, трехмерный объем и т.п. Изображение далее загружается в память для хранения изображений, где оно доступно для преобразования срезов, проекций или других участков представления изображения в соответствующий формат для визуализации, например с помощью видеомонитора 18, обеспечивающего воспринимаемое человеком отображение итогового МР-изображения.
readout
0
1
0
1
0
0
1
1
Согласно изобретению генерируется, по меньшей мере, два градиентных эхо-сигнала в два различных момента времени появления эхо-сигналов, при этом 0-й момент считываемого градиента магнитного поля по существу исчезает во время первого градиентного эха, первый момент считываемого градиента отличен от нуля во время первого градиентного эха, при этом как и 0-й, и 1-й момент считываемого градиента магнитного поля по существу исчезает во время второго градиентного эха. Это схематично показано на Фиг. 2. Показаны считываемый градиент G магнитного поля по схеме сбора с двойным эхо по изобретению, а также 0-й момент M и 1-й момент M считываемого градиента магнитного поля. Нулевой момент M исчезает в оба момента времени появления эхо-сигналов (обозначены пунктирными линиями на Фиг. 2), однако 1-й момент M исчезает только во второй момент времени появления эхо-сигнала. То же самое относится к нечетным и четным эхо-сигналам при множественном сборе эхо-сигналов согласно изобретению. Характерная компенсация потока второго градиентного эха достигается благодаря эффекту рефазировки момента градиента. 0-й момент M считываемого градиента магнитного поля определяется как интеграл от считываемого градиента магнитного поля по времени. Величина M равна нулю в моменты времени первого и второго градиентных эхо-сигналов. 1-ый момент M считываемого градиента магнитного поля определяется как интеграл от произведения считываемого градиента магнитного поля на время, взятый по времени. Величина M равна нулю только в момент времени второго градиентного эха. Это делает градиентные эхо-сигналы первого градиентного эха значительно более подверженными артефактам потока, чем градиентные эхо-сигналы второго градиентного эха.
На Фиг. 3 показано первое (левое) и второе (среднее) МР-изображения, реконструированные из градиентных эхо-сигналов первого градиентного эха и градиентных эхо-сигналов второго градиентного эха соответственно. Приведенные первое и второе МР-изображения представляют собой выбранные срезы, полученные при проведении трехмерной ангиографии с применением контрастного вещества. Существенное появление посторонних изображений (отмеченных белыми стрелками) присутствует на первом МР-изображении (левом), как и ожидалось из приведенного выше обсуждения способа по изобретению. На втором МР-изображении (среднем) артефакты, связанные с двоением изображения, не просматриваются. Правое изображение на Фиг. 3 представляет собой разность первого (левого) и второго (среднего) МР-изображений. Разностное изображение можно использовать для идентификации артефактов, связанных с двоением изображения, на первом МР-изображении, поскольку разностное изображение показывает потерю локальной интенсивности, а также усиление локальной интенсивности на первом МР-изображении вследствие появления посторонних изображений. Артефакты, связанные с двоением изображения, далее могут быть устранены на первом МР-изображении путем минимизации идентифицированных потерь или усилений интенсивности.
Представленные МР-изображения брюшной полости пациента получены с использованием модулированной в трехмерном пространстве двойной градиент-эхо последовательности при напряженности основного магнитного поля 1,5 Т, при этом время появления эхо-сигналов первого и второго градиентных эхо-сигналов составляло 1,8 мс и 3,0 мс соответственно. МР-сигналы, относящиеся к воде и жиру, чуть более сдвинуты по фазе в первое эхо-время, чем во второе эхо-время.
На Фиг. 4 показано МР-изображение воды (левое) и МР-изображение жира (правое), полученные из первого и второго МР-изображений, представленных на Фиг. 2, с использованием традиционного двухточечного способа разделения по Диксону, т.е. без устранения артефактов, связанных с двоением изображения. Как можно видеть на Фиг. 4, артефакты, связанные с двоением изображения и обусловленные потоком, распространяются в оба изображения (белые стрелки).
На Фиг. 5 показаны МР-изображение воды (левое) и МР-изображение жира (правое), полученные опять же из первого и второго МР-изображений, представленных на Фиг. 2, при этом использован способ по изобретению для устранения артефактов, связанных с двоением изображения. Путем детектирования пониженной локальной интенсивности изображения на втором МР-изображении по сравнению с первым МР-изображением и путем устранения локальных потерь интенсивности изображения на первом МР-изображении артефакты, связанные с двоением изображения, преимущественно подавляются. Поскольку время появления эхо-сигналов выбирается так, что вклады сигналов от воды и жира более совпадают по фазе во время второго градиентного эха, чем во время первого градиентного эха, осуществляется рефазировка воды и жира, а значит, они не вызывают потери сигнала по отношению ко второму времени эхо-задержки. Таким образом, потери интенсивности могут быть отнесены к артефактам, связанным с двоением изображения и присутствующим на первом МР-изображении, но подавленным на втором МР-изображении благодаря эффекту рефазировки момента градиента. Путем устранения идентифицированных артефактов, связанных с двоением изображения, достигается существенное улучшение качества изображения по сравнению с изображениями, показанными на Фиг. 4, без продления времени появления эха и времени повторения.