RU2672151C2 - Коррекция посторонних эхосигналов epi - Google Patents

Коррекция посторонних эхосигналов epi Download PDF

Info

Publication number
RU2672151C2
RU2672151C2 RU2016142190A RU2016142190A RU2672151C2 RU 2672151 C2 RU2672151 C2 RU 2672151C2 RU 2016142190 A RU2016142190 A RU 2016142190A RU 2016142190 A RU2016142190 A RU 2016142190A RU 2672151 C2 RU2672151 C2 RU 2672151C2
Authority
RU
Russia
Prior art keywords
image
ltr
rtl
data
phase
Prior art date
Application number
RU2016142190A
Other languages
English (en)
Other versions
RU2016142190A (ru
RU2016142190A3 (ru
Inventor
Михел Паул Юрриан ЮРРИССЕН
Йоханнес Петрус ГРУН
Миха ФЮДЕРЕР
Original Assignee
Конинклейке Филипс Н.В.
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Конинклейке Филипс Н.В. filed Critical Конинклейке Филипс Н.В.
Publication of RU2016142190A publication Critical patent/RU2016142190A/ru
Publication of RU2016142190A3 publication Critical patent/RU2016142190A3/ru
Application granted granted Critical
Publication of RU2672151C2 publication Critical patent/RU2672151C2/ru

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00—Arrangements or instruments for measuring magnetic variables
    • G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48—NMR imaging systems
    • G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • G01R33/561—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution by reduction of the scanning time, i.e. fast acquiring systems, e.g. using echo-planar pulse sequences
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00—Arrangements or instruments for measuring magnetic variables
    • G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48—NMR imaging systems
    • G01R33/4818—MR characterised by data acquisition along a specific k-space trajectory or by the temporal order of k-space coverage, e.g. centric or segmented coverage of k-space
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00—Arrangements or instruments for measuring magnetic variables
    • G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48—NMR imaging systems
    • G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • G01R33/561—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution by reduction of the scanning time, i.e. fast acquiring systems, e.g. using echo-planar pulse sequences
    • G01R33/5611—Parallel magnetic resonance imaging, e.g. sensitivity encoding [SENSE], simultaneous acquisition of spatial harmonics [SMASH], unaliasing by Fourier encoding of the overlaps using the temporal dimension [UNFOLD], k-t-broad-use linear acquisition speed-up technique [k-t-BLAST], k-t-SENSE
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00—Arrangements or instruments for measuring magnetic variables
    • G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48—NMR imaging systems
    • G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • G01R33/561—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution by reduction of the scanning time, i.e. fast acquiring systems, e.g. using echo-planar pulse sequences
    • G01R33/5615—Echo train techniques involving acquiring plural, differently encoded, echo signals after one RF excitation, e.g. using gradient refocusing in echo planar imaging [EPI], RF refocusing in rapid acquisition with relaxation enhancement [RARE] or using both RF and gradient refocusing in gradient and spin echo imaging [GRASE]
    • G01R33/5616—Echo train techniques involving acquiring plural, differently encoded, echo signals after one RF excitation, e.g. using gradient refocusing in echo planar imaging [EPI], RF refocusing in rapid acquisition with relaxation enhancement [RARE] or using both RF and gradient refocusing in gradient and spin echo imaging [GRASE] using gradient refocusing, e.g. EPI
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00—Arrangements or instruments for measuring magnetic variables
    • G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48—NMR imaging systems
    • G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • G01R33/565—Correction of image distortions, e.g. due to magnetic field inhomogeneities
    • G01R33/56545—Correction of image distortions, e.g. due to magnetic field inhomogeneities caused by finite or discrete sampling, e.g. Gibbs ringing, truncation artefacts, phase aliasing artefacts
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00—Arrangements or instruments for measuring magnetic variables
    • G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48—NMR imaging systems
    • G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • G01R33/565—Correction of image distortions, e.g. due to magnetic field inhomogeneities
    • G01R33/56554—Correction of image distortions, e.g. due to magnetic field inhomogeneities caused by acquiring plural, differently encoded echo signals after one RF excitation, e.g. correction for readout gradients of alternating polarity in EPI
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00—Arrangements or instruments for measuring magnetic variables
    • G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48—NMR imaging systems
    • G01R33/58—Calibration of imaging systems, e.g. using test probes, Phantoms; Calibration objects or fiducial markers such as active or passive RF coils surrounding an MR active material
    • G01R33/583—Calibration of signal excitation or detection systems, e.g. for optimal RF excitation power or frequency
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B2217/00—General characteristics of surgical instruments
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61K—PREPARATIONS FOR MEDICAL, DENTAL OR TOILETRY PURPOSES
    • A61K6/00—Preparations for dentistry
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R1/00—Details of instruments or arrangements of the types included in groups G01R5/00 - G01R13/00 and G01R31/00

Landscapes

  • Physics & Mathematics (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Condensed Matter Physics & Semiconductors (AREA)
  • General Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Radiology & Medical Imaging (AREA)
  • Engineering & Computer Science (AREA)
  • Signal Processing (AREA)
  • Magnetic Resonance Imaging Apparatus (AREA)

Abstract

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

Description

ОБЛАСТЬ ИЗОБРЕТЕНИЯ
Изобретение относится к способу магнитно-резонансной томографии, который содержит последовательность сбора данных эхопланарной томографии (EPI).
ПРЕДПОСЫЛКИ ИЗОБРЕТЕНИЯ
Сбор данных EPI из магнитно-резонансных сигналов известен из статьи "Image-based ghost correction for interleaved EPI" in MRM 45(2001)96-108 by M.H. Buoncore and D.C. Zhu.
Известный способ магнитно-резонансной томографии направлен на решение проблемы, связанной с тем, что фазовые искажения в EPI последовательностях вызывают посторонние эхосигналы в виде функции общего вида направлений (x) считывания и направления (y) фазового кодирования. Изображения реконструируют с использованием только линий k-пространства, проходящих слева направо или справа налево (при этом недостающие линии заменяют нулями). Оптимальное фазовое искажение в каждом местоположении (x,y) находят из итерационного решения.
Кроме того, в международной заявке WO2012/047771 раскрыт способ магнитно-резонансной томографии, в котором сглаженные изображения формируют из данных, полученных в многочисленных EPI-снимках.
СУЩНОСТЬ ИЗОБРЕТЕНИЯ
Задача настоящего изобретения состоит в том, чтобы обеспечить способ магнитно-резонансной томографии со сбором данных EPI с коррекцией посторонних эхосигналов, которая требует менее громоздких вычислений.
Данная задача решена с помощью способа магнитно-резонансной томографии по изобретению, который содержит:
- последовательность сбора данных эхопланарной томографии (EPI), которая включает в себя дискретизацию k-пространства для магнитно-резонансных сигналов для сбора:
- набора (mltr) ltr-данных, реконструированного из линий в k-пространстве, просканированных вдоль положительного направления обхода в k-пространстве, и
- набора (mrtl) rtl-данных, реконструированного из линий в k-пространстве, просканированных вдоль отрицательного направления обхода в k-пространстве,
- причем магнитно-резонансные сигналы из набора ltr-данных и набора rtl-данных получают с помощью нескольких РЧ приемных антенн, имеющих профили пространственной чувствительности;
- получение доступа к пространственному распределению фазовых ошибок для набора ltr-данных и для набора rtl-данных;
- формирование матрицы (Sltr) ltr-кодирования из (i) фазового кодирования линий в k-пространстве, просканированных вдоль положительного направления обхода в k-пространстве, (ii) пространственного распределения фазовых ошибок для набора ltr-данных и (iii) профилей пространственной чувствительности катушек;
- формирование матрицы (Srtl) rtl-кодирования из (i) фазового кодирования линий в k-пространстве, просканированных вдоль отрицательного направления обхода в k-пространстве, (ii) пространственного распределения фазовых ошибок для набора rtl-данных и (iii) профилей пространственной чувствительности катушек;
- объединение матрицы ltr-кодирования и матрицы rtl-кодирования в общую матрицу S кодирования;
- реконструкцию диагностического магнитно-резонансного изображения (p) путем решения соотношения кодирования между наборами ltr- и rtl-данных для пиксельных значений (pj(r)) магнитно-резонансного изображения:
Figure 00000001
.
Идея изобретения состоит в том, что из-за фазовых ошибок матрицы кодирования для магнитно-резонансных изображений, реконструированных, соответственно, из набора ltr-данных и набора rtl-данных, которые измеряются (дискретизируются) из линий k-пространства, которые пересекаются вдоль положительного и отрицательного направления распространения вдоль направления считывания в k-пространстве, являются различными. Более того, фазовые ошибки зависят как от позиции вдоль направления (x) считывания, так и от направления (y) фазового кодирования. Магнитно-резонансные сигналы получают с помощью многочисленных приемных антенн (катушек или катушечных элементов), имеющих профиль пространственной чувствительности. Соответственно, существует соотношение кодирования, которое устанавливается по аналогии с подходом параллельного формирования изображения SENSE между набором ltr-данных и окончательным магнитно-резонансным изображением, а также между набором rtl-данных и диагностическим магнитно-резонансным изображением. Это соотношение кодирования объединяет градиентное (фазовое) кодирование, фазовые ошибки, которые добавляются к фазовому кодированию и профилям пространственной чувствительности, которые определяют линейные комбинации пиксельных значений из-за субдискретизации в k-пространстве, которые вызывают свертку в пространстве изображений. Как из набора rtl-данных, так и из набора ltr-данных могут быть реконструированы отдельные магнитно-резонансные изображения и происходит развертка субдискретизации события, что выполняется на базе профилей пространственной чувствительности. Это обеспечивает достаточную избыточность для решения распределения фазовых ошибок. Распределение фазовых ошибок может быть первоначально оценено, измерено на этапе калибровки или получено из предыдущей итерации или более ранней динамики. На основании профилей пространственной чувствительности и пространственного распределения фазовых ошибок устанавливается соотношение глобального кодирования между объединенными набором ltr-данных и набором rtl-данных и окончательным магнитно-резонансным изображением (пиксельными значениями окончательного магнитно-резонансного изображения). Другая идея настоящего изобретения состоит в том, что это соотношение глобального кодирования является, как правило, избыточным в зависимости от плотности дискретизации в направлении фазового кодирования. Эта избыточность обеспечивает стабильное разрешение пиксельных значений (pj(r)) магнитно-резонансного изображения из соотношения кодирования в наборах ltr- и rtl-данных. Это выполняется с помощью обобщенных методов инверсии матриц для объединенной матрицы
Figure 00000002
 глобального кодирования. Альтернативно, подход со среднеквадратической минимизацией может быть использован для получения окончательного магнитно-резонансного изображения в виде наилучшего приближения к соотношению кодирования. То есть отыскивается решение соотношения глобального кодирования, которое имеет минимальное среднеквадратическое отклонение от точного решения. Такие методы разрешения этих соотношений кодирования с большой избыточностью между измеренными МР данными и пикселями окончательного магнитно-резонансного изображения обычно известны сами по себе из области технологий параллельного формирования изображений типа SENSE.
Эти и другие аспекты изобретения будут в дальнейшем представлены со ссылкой на варианты осуществления, определенные в зависимых пунктах формулы изобретения.
Различные фазовые ошибки, которые возникают в магнитно-резонансных сигналах, полученных из k-пространства вдоль положительного и отрицательного направлений обхода, соответственно, могут быть получены в виде субтракционного фазового изображения из развернутого ltr-изображения и развернутого rtl-изображения. Субтракционное фазовое изображение образует пространственное распределение фазовых ошибок. Развернутое ltr-изображение получают из набора ltr-данных, а развернутое rtl-изображение получают из набора rtl-данных. На практике эти развернутые ltr-изображения и rtl-изображения реконструируют из нечетных и четных линий в k-пространстве траектории k-пространства EPI. Реконструированные набор ltr-данных и набор rtl-данных могут иметь относительно низкое отношение сигнал/шум из-за более низкого переопределения или даже недоопределения их соотношений кодирования в окончательно реконструированном магнитно-резонансном изображении. Это вызвано относительно большим интервалом дискретизации в k-пространстве дискретизированных линий фазового кодирования. Из-за относительно низкой пространственной частоты фазовых ошибок пространственное распределение фазовых ошибок можно точно получить для субтракционного фазового изображения, например, путем фильтрации субтракционного фазового изображения для выравнивания вариаций, вызванных шумом.
В дополнительном воплощении изобретения более точное субтракционное фазовое изображение получают итерационным способом. На этапе калибровки измеряют начальное пространственное распределение фазовых ошибок. В простом подходе этап калибровки включает в себя сбор данных EPI путем сканирования вдоль линий k-пространства, в котором каждая линия сканируется вдоль чередующихся направлений распространения вдоль направления считывания. Затем начальное распределение фазовых ошибок используют для коррекции фазовых ошибок в развернутом ltr-изображении и развернутом rtl-изображении. Результирующие развернутое ltr-изображение с фазовой коррекцией и развернутое rtl-изображение с фазовой коррекцией затем используют для получения более точного субтракционного фазового изображения, которое обеспечивает пространственное распределение фазовых ошибок текущей итерации, которое, в свою очередь, может быть использовано в следующей итерации опять же для коррекции фазовых ошибок в развернутом ltr-изображении и развернутом rtl-изображении. Пространственное распределение фазовых ошибок может быть использовано для дополнительной коррекции фазовых ошибок в развернутом ltr-изображении и развернутом rtl-изображении в следующей итерации. Вместо коррекции фазовых ошибок в ltr-изображении и rtl-изображении фазовая коррекция может быть непосредственно применена к диагностическому магнитно-резонансному изображению. На практике достаточно нескольких итераций, чтобы достичь достаточно точной коррекции фазовых ошибок, которые пространственно распределены как в направлении считывания, так и в направлении фазового кодирования. На практике часто возникает ситуация, когда уже после первой итерации получается точное распределение фазовых ошибок. Таким образом, когда производится динамическая серия сбора данных EPI, то для каждой динамики достаточно только одной итерации для того, чтобы достичь точной динамической серии распределений фазовых ошибок для конкретной динамической серии. Такая динамическая серия обычно получается при функциональном обследовании МРТ (fMRI), в котором исследуется временное поведение BOLD-сигналов в кровеносных сосудах мозга пациента.
В другом воплощении изобретения на этапе калибровки сбор данных EPI выполняется при наличии одного набора линий в k-пространстве, каждая из которых сканируется сначала вдоль положительного, а затем вдоль отрицательного направления распространения, и при наличии другого набора линий в k-пространстве, каждая из которых сканируется сначала вдоль отрицательного, а затем вдоль положительного направления распространения. То есть на этапе калибровки траектории k-пространства аналогичны тем, которые приводят к набору rtl-данных и наборам ltr-данных. Этот подход частично корректирует модуляцию сигналов из-за различных моментов времени сбора данных вдоль линий k-пространства, проходящих вдоль положительного и отрицательного направления распространения в k-пространстве. Следует отметить, что этот подход уравнивает средние периоды времени сбора данных из набора ltr-данных и набора rtl-данных. Предложенное сканирование с калибровкой имеет своей целью обеспечить отдельные 2D фазовые карты для двух направлений считывания данных. С этой целью используется меньший набор фазовых кодирований, и все линии получаются для обоих направлений считывания данных. Дополнительный обход линии ky=0 выполняется для того, чтобы достичь одинакового среднего таймирования для обоих наборов. Это позволяет избежать повреждения разности фаз, например, за счет сдвигов фаз с наведенным смещением B0. Два ltr-сигнала линии ky=0 необходимо объединить соответствующим образом.
В другой уточненной реализации линия вдоль нулевого фазового кодирования сканируется в k-пространстве сначала вдоль положительного обхода, затем вдоль отрицательного обхода и, наконец, снова вдоль положительного обхода. Это позволяет иметь средние периоды времени сбора данных набора ltr-данных и набора rtl-данных в однокадровом сборе данных EPI.
В данном изобретении в последовательности сбора данных EPI для магнитно-резонансных сигналов k-пространство сканируют вдоль наборов линий в k-пространстве вдоль противоположных направлений распространения, например, нечетных и четных линий в k-пространстве. Фазовые ошибки, которые возникают из-за противоположных направлений распространения, корректируют при параллельной реконструкции изображения типа SENSE. Распределение фазовых ошибок в пространстве изображений может быть сначала оценено, а затем вычислено из разности фаз между изображениями, реконструированными из магнитно-резонансных сигналов, полученных из соответствующих наборов линий k-пространства или из более ранней динамики.
Изобретение дополнительно относится к компьютерной программе, включающей в себя инструкции для управления системой магнитно-резонансного исследования для выполнения способа по изобретению. Компьютерная программа по изобретению может быть предоставлена на носителе информации, таком как диск CD-ROM или карта памяти USB, или компьютерную программу по изобретению можно загрузить из сети передачи данных, такой как Интернет. При установке на компьютере, который включен в систему магнитно-резонансной томографии, система магнитно-резонансной томографии позволяет работать согласно изобретению и позволяет обеспечить сбор данных EPI с коррекцией посторонних эхосигналов.
КРАТКОЕ ОПИСАНИЕ ЧЕРТЕЖЕЙ
Эти и другие аспекты изобретения будут объяснены со ссылкой на варианты осуществления, описанные в дальнейшем, и со ссылкой на сопроводительные чертежи, на которых:
на фиг.1 показано схематичное представление выполнения способа магнитно-резонансной томографии по изобретению;
на фиг.2 показано представление примера двумерной EPI-траектории k-пространства для этапа калибровки для измерения начального пространственного распределения фазовых ошибок, которое необходимо измерить;
на фиг.3 показано схематичное представление системы магнитно-резонансного исследования, в которую включено настоящее изобретение.
ПОДРОБНОЕ ОПИСАНИЕ ИЗОБРЕТЕНИЯ
На фиг.1 показано схематичное представление выполнения способа магнитно-резонансной томографии по настоящему изобретению. При однокадровом сборе 101 данных EPI двумерное k-пространство сканируется вдоль линий в положительном направлении распространения (слева направо) и в отрицательном направлении распространения (справа налево). Данные, полученные таким образом, образуют набор ltr-данных и набор rtl-данных. Как правило, нечетные линии или четные линии сканируются вдоль противоположных направлений. Из данных, полученных путем сканирования вдоль положительных направлений распространения (например, нечетные линии), то есть из набора ltr-данных, формируют ltr-изображение 102, как правило, путем быстрого преобразования Фурье набора ltr-данных. Из данных, полученных путем сканирования вдоль положительных направлений распространения (например, нечетных линий), то есть из набора rtl-данных, реконструируется свернутое rtl-изображение 103, как правило, путем быстрого преобразования Фурье набора ltr-данных. Обычно, набор ltr-данных (rtl-данных), образованный только нечетными (четными) линиями, будет субдискретизироваться в направлении фазового кодирования в k-пространстве, приводя к возникновению зубцеобразного дефекта, который проявляется в виде сложных артефактов. Осуществляют развертку свернутого ltr-изображения и свернутого rtl-изображения 104, 105 посредством реконструкции SENSE, которая позволяет использовать профиль 111 чувствительности катушки. Эта реконструкция SENSE сама по себе является известной и обычно используется для реконструкции данных, которые субдискретизируются в k-пространстве для того, чтобы уменьшить время сбора данных. Из развернутого ltr-изображения 104 и развернутого rtl-изображения 105 получают фазоразностное изображение 112. Это фазоразностное изображение содержит разности фаз набора ltr-данных относительно набора rtl-данных. Эта разности фаз создаются из различий между фазовыми кодированиями и различия между фазовыми ошибками в наборе ltr-данных и в наборе rtl-данных, которые вызваны эффектами, такими как задержки при градиентном переключении и токе Фуко. Так как фазовое кодирование в наборе ltr-данных и в наборе rtl-данных применяется управляемым образом, пространственное распределение 113 фазовых ошибок может быть получено из фазоразностного изображения 112. Для того чтобы инициировать процесс, на этапе 114 калибровки производится первое измерение пространственного распределения фазовых ошибок путем сканирования двумерного k-пространства вдоль набора линий фазового кодирования, где каждая линия сканируется вдоль противоположных направлений. Разность между фазой магнитно-резонансных сигналов, полученных в противоположных направлениях распространения, представляет собой фазовые ошибки. Так как измерение производится при различных значениях фазового кодирования, фазовые ошибки измеряются как в зависимости от направления считывания, так и направления частотного кодирования. Это измеренное распределение фазовых ошибок образует начальную оценку пространственного распределения фазовых ошибок. Распределение фазовых ошибок может быть получено различными способами из фазоразностного изображения, измерения на этапе калибровки и/или пространственного распределения 113 фазовых ошибок, доступного из предыдущей итерации или из предыдущей динамики. Распределение фазовых ошибок могло бы быть просто заменено на доступное в текущий момент времени фазоразностное изображение, но текущее фазоразностное изображение может быть также объединено с результатом, полученным на этапе калибровки из предыдущей динамики или сохраненной информации. Таким образом, принимается в расчет возможно более высокий уровень шума в фазоразностных изображениях, а также результат, полученный на этапе калибровки, который становится относительно устаревшим.
Диагностическое магнитно-резонансное изображение реконструируют 106 на основании профилей пространственной чувствительности, пространственного распределения фазовых ошибок, а соотношение глобального кодирования устанавливают между расположенными рядом набором ltr-данных и набором rtl-данных и окончательным магнитно-резонансным изображением (пиксельными значениями окончательного магнитно-резонансного изображения). Пиксельные значения диагностического магнитно-резонансного изображения могут быть использованы для подхода регуляризации при развертке (104, 105) rtl-свернутого изображения и ltr-свернутых изображений. Эта регуляризация повышает численную устойчивость решения развертки rtl-свернутого изображения и ltr-свернутого изображения и уменьшает уровень шума в развернутых изображениях. Следовательно, в распределении фазовых ошибок уменьшается уровень шума.
В итерационном подходе диагностическое магнитно-резонансное изображение может периодически обновляться в виде последовательных итераций. Таким образом, на текущей итерации регуляризация может быть выполнена над относительно точными пиксельными значениями имеющейся в текущий момент версии диагностического магнитно-резонансного изображения, которое может периодически обновляться с тем, чтобы учесть распределение фазовых ошибок, которое становится более точным с увеличением количества итераций.
Альтернативно, в динамическом подходе (схематично показанном пунктирными линиями) производятся однокадровые сборы данных EPI. Таким образом, сбор данных в k-пространстве производится в следующие один за другим моменты времени, которые в дальнейшем представляют динамическое изменение, которое может происходить в объекте из-за перемещения, вызванного сердцебиением и дыханием пациента. В таком динамическом подходе в качестве оценки распределения фазовых ошибок можно применять результат для распределения фазовых ошибок, полученный из более ранней, предпочтительно предыдущей динамики. Оказывается, что распределение фазовых ошибок только медленно варьируется в зависимости от динамики, и в итерационном подходе также оказывается, что часто единственной итерации достаточно для точного определения распределения фазовых ошибок.
Более конкретно, с математической точки зрения, реконструкция диагностического магнитно-резонансного изображения из однокадрового сбора (однокадровых сборов) данных EPI выглядит следующим образом. Для нечетных и четных эхосигналов (т.е. для набора ltr-данных и набора rtl-данных) уравнения SENSE можно записать в виде:
,
где S - матрица чувствительности катушки, m - результаты измерений, и p - пиксели окончательного развернутого ltr-изображения 104 и развернутого rtl-изображения 105.
Φ - диагональная матрица, содержащая дельта фазовое кодирование ΔΦenc и 2D EPI фазовые ошибки Φepi:
Figure 00000003
.
При этом N=2*R, где R - коэффициент уменьшения SENSE, и 2 - дополнительный "коэффициент SENSE" при делении на нечетные и четные эхосигналы. Кодирование чувствительности, оставшееся фазовое кодирование и фазовые ошибки можно объединить в одну матрицу Seven и Sodd "чувствительности".
.
Предположим, что четные эхосигналы содержат профиль k=0 (поэтому отсутствует дельта фазовое кодирование), и только дельта 2D EPI фазовая коррекция применяется к нечетным эхосигналам, таким образом, Φ для четных и нечетных эхосигналов имеет вид:
Figure 00000004
Figure 00000005
 при
Figure 00000006
.
Член ΔΦenc намеренно введен для того, чтобы иметь лучшее разделение между наборами уравнений; это производится с помощью расстояния k-пространства между нечетными и четными линиями k-пространства. Теперь уравнение SENSE можно записать в виде:
Figure 00000007
,
().
Объединение нечетных и четных эхосигналов в одном так называемом ядре реконструкции SENSE-IRIS имеет вид:
Figure 00000008
,
где pall - диагностическое магнитно-резонансное изображение, возникающее в результате из нечетных-четных эхосигналов, в сочетании с разверткой реконструкции, включая 2D EPI фазовую коррекцию.
Теперь, чтобы определить EPI фазовые изменения в зависимости от динамики, нечетные и четные эхосигналы могут быть дополнительно реконструированы по отдельности. Например, для четных эхосигналов имеет место:
Figure 00000009
Здесь pall используется для регуляризации (R) реконструкции. Функция f должна быть настроена и "инструктирует" реконструкцию SENSE относительно того, насколько близким должно быть решение по отношению к предыдущему полному решению (pall). Вклад шума представляет собой nn в измеренных данных. Реконструкция развернутых пикселей может включать в себя неопределенности и ошибки реконструкции, которые учитываются членом np. В подходе регуляризации при развертке предполагается, что решение развертки близко к контрольной точке pall. Эту контрольную точку можно получить из предыдущего решения, такого как (1) решение предыдущей итерации или (2) решение предыдущей динамики.
Матрица R регуляризации, как правило, зависит от модуля pall:
Figure 00000010
 .
Размер матрицы эквивалентен числу свернутых пикселей (Sf) SENSE. Если предполагается, что решение будет отличаться на 10% от pall, то F=0,1. Если для регуляризации используется предыдущая динамика, и ожидаемые изменения (за счет нагревания) являются маленькими (например, 1%), то можно выбрать меньшее значение F (например, 0,01). При итерационной реконструкции f можно уменьшить (более сильная регуляризация) для более высоких итераций.
Решение SENSE для четных и нечетных эхосигналов имеет вид:
Figure 00000011
,
Figure 00000012
.
Таким образом, фазоразностное изображение вычисляется путем простого вычитания:
Figure 00000013
.
Отдельная реконструкция SENSE нечетных и четных эхосигналов приводит к более низкому отношению сигнал/шум из-за более высоких геометрических факторов SENSE. Поэтому эти изображения не являются оптимальными для клинического использования, но достаточными для определения глобального 2D EPI фазового изменения, например, для соответствия ΔΦEPI 2D линейной фазовой ошибке или сильного сглаживания карты ΔΦEPI.
Оценочное значение ΔΦEPI добавляется к 2D EPI фазовым ошибкам и используется в реконструкции SENSE-IRIS следующей динамики, что приводит к динамически обновленной 2D EPI фазовой коррекции, интегрированной в реконструкции SENSE.
На фиг.3 схематично показана система магнитно-резонансной томографии, в которой применимо настоящее изобретение. Система магнитно-резонансной томографии включает в себя набор основных катушек 10, посредством которых генерируется постоянное, однородное магнитное поле. Например, основные катушки сконструированы так, чтобы они обхватывали пространство для обследования в форме туннеля. Подлежащий обследованию пациент размещается на подвижном столе для пациента, который плавно перемещается в этом пространстве для обследования в форме туннеля. Система магнитно-резонансной томографии также включает в себя ряд градиентных катушек 11, 12, посредством которых генерируются магнитные поля, проявляющие пространственные вариации особенно в форме временных градиентов в конкретных направлениях так, чтобы накладываться на однородное магнитное поле. Градиентные катушки 11, 12 соединены с управляемым блоком 21 источника питания. Градиентные катушки 11, 12 возбуждаются за счет подачи электрического тока посредством блока 21 источника питания; с этой целью блок источника питания оснащен электронной схемой градиентного усиления, которая подает электрический ток в градиентные катушки для того, чтобы сгенерировать градиентные импульсы (которые также упоминаются как "градиентные сигналы") подходящей временной формы. Управление интенсивностью направлением и длительностью градиентов осуществляется путем управления от блока источника питания. Система магнитно-резонансной томографии также включает в себя передающей и приемной катушки 13, 16 для генерации РЧ импульсов возбуждения и для захвата магнитно-резонансных сигналов, соответственно. Передающая катушка 13 предпочтительно сконструирована в виде катушки 13 для тела, посредством которой может быть окружен подлежащий обследованию объект (часть обследуемого объекта). Катушка для тела обычно размещается в системе магнитно-резонансной томографии так, чтобы подлежащий обследованию пациент 30 был окружен катушкой 13 для тела, когда он или она размещается в системе магнитно-резонансной томографии. Катушка 13 для тела действует как передающая антенна для передачи РЧ импульсов возбуждения и РЧ импульсов повторной фокусировки. Предпочтительно, катушка 13 для тела включает в себя пространственно однородное распределение интенсивности передаваемых РЧ импульсов (RFS). Такая же катушка или антенна обычно используется альтернативно в качестве передающей катушки и приемной катушки. Кроме того, передающая и приемная катушка обычно имеет форму катушки, но также возможны и другие геометрические формы, при этом передающая и приемная катушка действует как передающая и приемная антенна для РЧ электромагнитных сигналов. Передающая и приемная катушка 13 соединена с электронной передающей и приемной схемой 15.
Следует отметить, что альтернативно можно использовать отдельные приемные и/или передающие катушки 16. Например, поверхностные катушки 16 можно использовать в качестве приемных и/или передающих катушек. Такие поверхностные катушки имеют высокую чувствительность при сравнительно маленьком объеме. Приемные катушки, такие как поверхностные катушки, соединены с демодулятором 24, и принимаемые магнитно-резонансные сигналы (MS) демодулируются посредством демодулятора 24. Демодулированные магнитно-резонансные сигналы (DMS) подаются в блок реконструкции. Приемная катушка соединена с предусилителем 23. Предусилитель 23 усиливает РЧ резонансный сигнал (MS), принятый приемной катушкой 16, и усиленный РЧ резонансный сигнал подается в демодулятор 24. Демодулятор 24 демодулирует усиленный РЧ резонансный сигнал. Демодулированный резонансный сигнал содержит фактическую информацию, относящуюся к локальной плотности спинов в части изображаемого объекта. Кроме того, передающая и приемная схема 15 соединена с модулятором 22. Модулятор 22 и передающая и приемная схема 15 возбуждают передающую катушку 13 для того, чтобы передать РЧ импульсы возбуждения повторной фокусировки. Блок реконструкции получает один или более сигналов изображения из демодулированных магнитно-резонансных сигналов (DMS), сигналы изображения которых представляют информацию в виде изображений изображенной части обследуемого объекта. На практике блок 25 реконструкции сконструирован предпочтительно в виде блока 25 обработки цифровых изображений, который запрограммирован таким, чтобы получать из демодулированных магнитно-резонансных сигналов сигналы изображения, которые представляют информацию в виде изображений части объекта, подлежащей отображению. Сигнал на выходе монитора 26 реконструкции, такого как монитор, может отображать магнитно-резонансное изображение. Альтернативно, можно сохранить сигнал, подаваемый из блока 25 реконструкции в буферный блок 27 во время ожидания дальнейшей обработки.
Система магнитно-резонансной томографии согласно изобретению также содержит блок 20 управления, например, в виде компьютера, который включает в себя (микро)процессор. Блок 20 управления управляет выполнением РЧ возбуждения и прикладывания временных градиентных полей. С этой целью компьютерная программа согласно изобретению загружается, например, в блок 20 управления и блок 25 реконструкции.

Claims (42)

1. Способ магнитно-резонансной томографии, содержащий:
- последовательность сбора данных эхопланарной томографии (EPI), которая включает в себя дискретизацию k-пространства для магнитно-резонансных сигналов для сбора:
- набора (mltr) ltr-данных, реконструированного из линий в k-пространстве, просканированных вдоль положительного направления обхода в k-пространстве, и
- набора (mrtl) rtl-данных, реконструированного из линий в k-пространстве, просканированных вдоль отрицательного направления обхода в k-пространстве;
- причем магнитно-резонансные сигналы из набора ltr-данных и набора rtl-данных собирают с помощью нескольких РЧ приемных антенн, имеющих профили пространственной чувствительности;
- получение доступа к пространственному распределению фазовых ошибок для набора ltr-данных и для набора rtl-данных;
- формирование матрицы (Sltr) ltr-кодирования из (i) фазового кодирования линий в k-пространстве, просканированных вдоль положительного направления обхода в k-пространстве, (ii) пространственного распределения фазовых ошибок для набора ltr-данных и (iii) профилей пространственной чувствительности катушек;
- формирование матрицы (Srtl) rtl-кодирования из (i) фазового кодирования линий в k-пространстве, просканированных вдоль отрицательного направления обхода в k-пространстве, (ii) пространственного распределения фазовых ошибок для набора rtl-данных и (iii) профилей пространственной чувствительности катушек;
- объединение матрицы ltr-кодирования и матрицы rtl-кодирования в общую матрицу S кодирования;
- реконструкцию диагностического магнитно-резонансного изображения (p) путем решения соотношения кодирования между наборами ltr- и rtl-данных для пиксельных значений (pj(r)) магнитно-резонансного изображения:
Figure 00000014
, и причем набор ltr-данных и/или набор rtl-данных субдескретизируют в k-пространстве; и
- развернутое ltr-изображение получают из набора ltr-данных в виде решения соотношения кодирования между набором ltr-данных и пиксельными значениями развернутого ltr-изображения:
Figure 00000015
,
- развернутое rtl-изображение получают из набора rtl-данных в виде решения соотношения кодирования между набором rtl-данных и пиксельными значениями развернутого rtl-изображения:
Figure 00000016
, и
- субтракционное фазовое изображение формируют из развернутого ltr-изображения и развернутого rtl-изображения и используют субтракционное фазовое изображение в качестве пространственного распределения фазовых ошибок, и причем
- начальное пространственное распределение фазовых ошибок измеряют на этапе калибровки;
- развернутое ltr-изображение с фазовой коррекцией формируют тем, что развернутое ltr-изображение или набор ltr-данных корректируют по начальному пространственному распределению фазовых ошибок;
- развернутое rtl-изображение с фазовой коррекцией формируют тем, что развернутое ltr-изображение или набор rtl-данных корректируют по начальному пространственному распределению фазовых ошибок; и
- субтракционное фазовое изображение формируют из развернутого ltr-изображения с фазовой коррекцией и развернутого rtl-изображения с фазовой коррекцией, и причем ltr-изображение с фазовой коррекцией и rtl-изображение с фазовой коррекцией вырабатывают итерационным образом.
2. Способ магнитно-резонансной томографии по п.1, причем результирующие развернутое ltr-изображение с фазовой коррекцией и развернутое rtl-изображение с фазовой коррекцией используют для получения более точного субтракционного фазового изображения, которое обеспечивает пространственное распределение фазовых ошибок текущей итерации, которое, в частности, в свою очередь, используют в следующей итерации для повторной коррекции фазовых ошибок в развернутом ltr-изображении и развернутом rtl-изображении.
3. Способ магнитно-резонансной томографии по п.1, причем начальное пространственное распределение фазовых ошибок измеряют на этапе калибровки, на котором производят сбор данных EPI, имеющих один набор линий в k-пространстве, каждую из которых сканируют сначала вдоль положительного, а затем вдоль отрицательного направления распространения, и имеющих другой набор линий в k-пространстве, каждую из которых сканируют сначала вдоль отрицательного, а затем вдоль положительного направления распространения.
4. Способ магнитно-резонансной томографии по п.3, причем линию вдоль нулевого фазового кодирования сканируют в k-пространстве сначала вдоль положительного обхода, затем вдоль отрицательного обхода и, наконец, снова вдоль положительного обхода.
5. Носитель информации, содержащий компьютерную программу для управления системой магнитно-резонансного исследования, содержащую инструкции для:
- применения последовательности сбора данных эхопланарной томографии (EPI), которая включает в себя дискретизацию k-пространства для магнитно-резонансных сигналов для сбора:
- набора (mltr) ltr-данных, реконструированного из линий в k-пространстве, просканированных вдоль положительного направления обхода в k-пространстве, и
- набора (mrtl) rtl-данных, реконструированного из линий в k-пространстве, просканированных вдоль отрицательного направления обхода в k-пространстве;
- сбора магнитно-резонансных сигналов из набора ltr-данных и набора rtl-данных с помощью нескольких РЧ приемных антенн, имеющих профили пространственной чувствительности;
- получения доступа к пространственным распределениям фазовых ошибок для набора ltr-данных и для набора rtl-данных;
- формирования матрицы (Sltr) ltr-кодирования из (i) фазового кодирования линий в k-пространстве, просканированных вдоль положительного направления обхода в k-пространстве, (ii) пространственного распределения фазовых ошибок для набора ltr-данных и (iii) профилей пространственной чувствительности катушек;
- формирования матрицы (Srtl) rtl-кодирования из (i) фазового кодирования линий в k-пространстве, просканированных вдоль отрицательного направления обхода в k-пространстве, (ii) пространственного распределения фазовых ошибок для набора rtl-данных и (iii) профилей пространственной чувствительности катушек;
- объединения матрицы ltr-кодирования и матрицы rtl-кодирования в общую матрицу S кодирования;
- реконструкции магнитно-резонансного изображения (p) путем решения соотношения кодирования между наборами ltr- и rtl-данных и пиксельными значениями (pj(r)) магнитно-резонансного изображения:
Figure 00000017
и
- для субдискретизации набора ltr-данных и/или набора rtl-данных в k-пространстве; и
- для получения развернутого ltr-изображения из набора ltr-данных в качестве решения соотношения кодирования между набором ltr-данных и пиксельными значениями развернутого ltr-изображения:
Figure 00000015
- для получения развернутого rtl-изображения, которое получается из набора rtl-данных в качестве решения соотношения кодирования между набором rtl-данных и пиксельными значениями развернутого rtl-изображения:
Figure 00000016
и
- для формирования субтракционного фазового изображения из развернутого ltr-изображения и развернутого rtl-изображения и использования субтракционного фазового изображения в качестве пространственного распределения фазовых ошибок; и
- для измерения начального пространственного распределения фазовых ошибок на этапе калибровки,
- для формирования развернутого ltr-изображения с фазовой коррекцией с тем, чтобы развернутое ltr-изображение или набор ltr-данных корректировались с учетом начального пространственного распределения фазовых ошибок;
- для формирования развернутого rtl-изображения с фазовой коррекцией с тем, чтобы развернутое ltr-изображение или набор rtl-данных корректировались с учетом начального пространственного распределения фазовых ошибок; и
- для формирования субтракционного фазового изображения из развернутого ltr-изображения с фазовой коррекцией и развернутого rtl-изображения с фазовой коррекцией, в котором ltr-изображение с фазовой коррекцией и rtl-изображения с фазовой коррекцией вырабатываются итерационным образом.
RU2016142190A 2014-03-28 2015-03-20 Коррекция посторонних эхосигналов epi RU2672151C2 (ru)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
EP14162169.8 2014-03-28
EP14162169 2014-03-28
PCT/EP2015/055862 WO2015144569A1 (en) 2014-03-28 2015-03-20 Epi ghost correction involving sense

Publications (3)

Publication Number Publication Date
RU2016142190A RU2016142190A (ru) 2018-05-03
RU2016142190A3 RU2016142190A3 (ru) 2018-06-20
RU2672151C2 true RU2672151C2 (ru) 2018-11-12

Family

ID=50382355

Family Applications (1)

Application Number Title Priority Date Filing Date
RU2016142190A RU2672151C2 (ru) 2014-03-28 2015-03-20 Коррекция посторонних эхосигналов epi

Country Status (6)

Country Link
US (1) US10241184B2 (ru)
EP (1) EP3123192B1 (ru)
JP (1) JP2017508556A (ru)
CN (1) CN106133546B (ru)
RU (1) RU2672151C2 (ru)
WO (1) WO2015144569A1 (ru)

Families Citing this family (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP3364207A1 (en) 2017-02-17 2018-08-22 Koninklijke Philips N.V. Phase error correction for bipolar read out gradients
EP3413074A1 (en) 2017-06-09 2018-12-12 Koninklijke Philips N.V. Parallel multi-slice mr imaging
EP3447520A1 (en) * 2017-08-22 2019-02-27 Koninklijke Philips N.V. Data-driven correction of phase depending artefacts in a magnetic resonance imaging system
CN108378851B (zh) * 2018-01-03 2021-08-06 上海东软医疗科技有限公司 磁共振弥散加权成像方法和装置
CN108872904B (zh) * 2018-06-06 2021-01-05 上海东软医疗科技有限公司 磁共振成像方法和装置
CN109342981B (zh) * 2018-09-20 2021-10-12 上海东软医疗科技有限公司 一种基于epi序列的动态成像方法和装置
US11002815B2 (en) 2019-03-21 2021-05-11 University Of Cincinnati System and method for reducing artifacts in echo planar magnetic resonance imaging
CN110570487B (zh) * 2019-08-26 2022-11-25 上海联影智能医疗科技有限公司 欠采样模型生成方法、图像重建方法、装置和计算机设备
US11568584B2 (en) 2019-08-26 2023-01-31 Shanghai United Imaging Intelligence Co., Ltd. Systems and methods for magnetic resonance imaging
CN111352055B (zh) * 2020-03-17 2022-05-06 无锡鸣石峻致医疗科技有限公司 一种前瞻性相位校正平面回波成像方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2002253526A (ja) * 2001-03-02 2002-09-10 Hitachi Medical Corp 磁気共鳴イメージング装置
JP2006255046A (ja) * 2005-03-16 2006-09-28 Hamano Life Science Research Foundation 磁気共鳴映像法および画像処理装置
US20080157767A1 (en) * 2007-01-02 2008-07-03 The Board Of Trustees Of The Leland Stanford Junior University Mri data acquisition using propeller k-space data acquisition
WO2012047771A2 (en) * 2010-10-07 2012-04-12 Duke University Multi-dimensional iterative phase-cycled reconstruction for mri images
WO2012145754A1 (en) * 2011-04-21 2012-10-26 The Board Of Trustees Of The University Of Illinois Correcting constant, linear and oblique phase errors in mr images acquired using propeller echo planar imaging
RU2011130555A (ru) * 2008-12-22 2013-01-27 Конинклейке Филипс Электроникс Н.В. Мр-томография с усилением контраста при переносе насыщения в химическом обмене (cest)

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20050033154A1 (en) * 2003-06-03 2005-02-10 Decharms Richard Christopher Methods for measurement of magnetic resonance signal perturbations
JP4619674B2 (ja) * 2004-03-24 2011-01-26 株式会社東芝 磁気共鳴イメージング装置
CN100409807C (zh) * 2005-11-29 2008-08-13 东南大学 基于均匀标记物校正回波平面成像技术中幽灵伪影的方法
CN101484823B (zh) 2006-07-06 2013-03-27 皇家飞利浦电子股份有限公司 磁共振成像设备和方法
US8649846B2 (en) * 2007-07-11 2014-02-11 Kabushiki Kaisha Toshiba Magnetic resonance imaging apparatus and magnetic resonance imaging method
US9606209B2 (en) * 2011-08-26 2017-03-28 Kineticor, Inc. Methods, systems, and devices for intra-scan motion correction
US9081074B2 (en) * 2012-08-21 2015-07-14 Beth Israel Deaconess Medical Center, Inc. (Bidmc, Inc.) Method and apparatus for accelerated phase contrast magnetic resonance angiography and blood flow imaging

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2002253526A (ja) * 2001-03-02 2002-09-10 Hitachi Medical Corp 磁気共鳴イメージング装置
JP2006255046A (ja) * 2005-03-16 2006-09-28 Hamano Life Science Research Foundation 磁気共鳴映像法および画像処理装置
US20080157767A1 (en) * 2007-01-02 2008-07-03 The Board Of Trustees Of The Leland Stanford Junior University Mri data acquisition using propeller k-space data acquisition
RU2011130555A (ru) * 2008-12-22 2013-01-27 Конинклейке Филипс Электроникс Н.В. Мр-томография с усилением контраста при переносе насыщения в химическом обмене (cest)
WO2012047771A2 (en) * 2010-10-07 2012-04-12 Duke University Multi-dimensional iterative phase-cycled reconstruction for mri images
WO2012145754A1 (en) * 2011-04-21 2012-10-26 The Board Of Trustees Of The University Of Illinois Correcting constant, linear and oblique phase errors in mr images acquired using propeller echo planar imaging

Also Published As

Publication number Publication date
US10241184B2 (en) 2019-03-26
CN106133546A (zh) 2016-11-16
RU2016142190A (ru) 2018-05-03
JP2017508556A (ja) 2017-03-30
RU2016142190A3 (ru) 2018-06-20
EP3123192B1 (en) 2020-10-14
US20170108571A1 (en) 2017-04-20
CN106133546B (zh) 2019-04-30
EP3123192A1 (en) 2017-02-01
WO2015144569A1 (en) 2015-10-01

Similar Documents

Publication Publication Date Title
RU2672151C2 (ru) Коррекция посторонних эхосигналов epi
US9797974B2 (en) Nonrigid motion correction in 3D using autofocusing with localized linear translations
JP6018401B2 (ja) 拡散強調エコープラナー撮像法において高次渦電流に誘発された歪みを予測補正するためのシステムおよび方法
CN109239633B (zh) 磁共振成像方法和***、校正扩散加权磁共振数据的方法
US10444315B2 (en) MRI with motion correction using navigators acquired using a dixon technique
US8648599B2 (en) Magnetic resonance imaging apparatus and magnetic resonance imaging method
US9050018B2 (en) Means and methods for providing high resolution MRI
CN103238082B (zh) 使用多点Dixon技术和低分辨率校准的MR成像
US9482732B2 (en) MRI reconstruction with motion-dependent regularization
JP6243522B2 (ja) 正則化された検出再構成を使用するマルチエコーディクソン水−脂肪分離及びｂ０歪補正による並列ｍｒｉ
JP6333293B2 (ja) 金属耐性ｍｒ画像化
EP3041418B1 (en) Method for correcting errors associated with asynchronous timing offsets between transmit and receive clocks in mri wireless radiofrequency coils
US20120002859A1 (en) Magnetic resonance partially parallel imaging (ppi) with motion corrected coil sensitivities
CN106574954A (zh) 针对epi的具有奈奎斯特伪影校正的并行mr成像
JP4679158B2 (ja) 磁気共鳴イメージング装置
US10690741B2 (en) Method and systems for reducing artifacts in magnetic resonance imaging
US20220057467A1 (en) Epi mr imaging with distortion correction
CN111164444B (zh) 具有经改进的脂肪位移校正的Dixon型水/脂肪分离MR成像
JP6348355B2 (ja) 磁気共鳴信号処理装置及び磁気共鳴装置並びにプログラム
JP2005288026A (ja) 磁気共鳴イメージング装置、渦磁場分布推定方法、及び静磁場補正方法
RU2773486C2 (ru) Магнитно-резонансная визуализация с разделением вода/жир методом диксона
JP2006522638A (ja) ヒストグラムベース位相補正を用いた磁気共鳴画像化

Legal Events

Date Code Title Description
MM4A The patent is invalid due to non-payment of fees

Effective date: 20210321