[go: nahoru, domu]

CN111714097B - Bimodal magnetic resonance temperature measurement method based on multi-gradient echo sequence - Google Patents

Bimodal magnetic resonance temperature measurement method based on multi-gradient echo sequence Download PDF

Info

Publication number
CN111714097B
CN111714097B CN202010621882.9A CN202010621882A CN111714097B CN 111714097 B CN111714097 B CN 111714097B CN 202010621882 A CN202010621882 A CN 202010621882A CN 111714097 B CN111714097 B CN 111714097B
Authority
CN
China
Prior art keywords
echo
magnetic resonance
heating
pixel point
temperature
Prior art date
Legal status (The legal status 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 status listed.)
Active
Application number
CN202010621882.9A
Other languages
Chinese (zh)
Other versions
CN111714097A (en
Inventor
曹鹏
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Hangzhou Jialiang Medical Technology Co ltd
Original Assignee
Hangzhou Jialiang Medical Technology Co ltd
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 Hangzhou Jialiang Medical Technology Co ltd filed Critical Hangzhou Jialiang Medical Technology Co ltd
Priority to CN202010621882.9A priority Critical patent/CN111714097B/en
Publication of CN111714097A publication Critical patent/CN111714097A/en
Priority to PCT/CN2021/102229 priority patent/WO2022001841A1/en
Application granted granted Critical
Publication of CN111714097B publication Critical patent/CN111714097B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/01Measuring temperature of body parts ; Diagnostic temperature sensing, e.g. for malignant or inflamed tissue
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/0033Features or image-related aspects of imaging apparatus classified in A61B5/00, e.g. for MRI, optical tomography or impedance tomography apparatus; arrangements of imaging apparatus in a room
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/05Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves 
    • A61B5/055Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves  involving electronic [EMR] or nuclear [NMR] magnetic resonance, e.g. magnetic resonance imaging

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Medical Informatics (AREA)
  • Surgery (AREA)
  • Engineering & Computer Science (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Biophysics (AREA)
  • Molecular Biology (AREA)
  • Pathology (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Radiology & Medical Imaging (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Magnetic Resonance Imaging Apparatus (AREA)

Abstract

The invention discloses a bimodal magnetic resonance temperature measurement method based on a multi-gradient echo sequence: acquiring magnetic resonance signals before and during heating through a multi-gradient echo magnetic resonance sequence, and simultaneously reserving a phase diagram and an amplitude diagram of each echo; according to the phase diagrams of different echoes before and during heating, calculating the temperature change before and after heating of each pixel point on the phase diagrams
Figure DDA0002563308520000011
Calculating the signal intensity of each pixel point TE when TE is 0 according to amplitude maps of different echoes in heating; calculating the temperature change delta T (M) of each pixel point based on the amplitude value according to the signal intensity of the amplitude diagram before heating and the signal intensity of each pixel point TE when the signal intensity is 0; according to
Figure DDA0002563308520000012
And Δ T (M) determines the final actual temperature T of each pixel point or issues an error warning. According to the temperature measurement method, the amplitude diagram and the phase diagram in multi-gradient-echo magnetic resonance temperature imaging are synchronously used, so that bimodal temperature imaging is realized, and the accuracy of magnetic resonance temperature imaging is further improved.

Description

Bimodal magnetic resonance temperature measurement method based on multi-gradient echo sequence
Technical Field
The invention relates to a magnetic resonance imaging technology, in particular to a bimodal magnetic resonance temperature measurement method based on a multi-gradient echo sequence.
Background
The real-time magnetic resonance temperature imaging can realize noninvasive, rapid and accurate temperature measurement on human tissues, is mainly applied to real-time temperature detection of minimally invasive surgeries such as tumor tissue thermal ablation and the like, and judges whether the temperature reaches an ideal value and temperature spatial distribution, thereby helping to judge the effect of the minimally invasive surgeries and early warning risks in the surgeries such as local over-temperature and the like.
Gradient echo magnetic resonance sequences (GRE, or other deformation sequences) are the mainstream technique of magnetic resonance temperature imaging methods, mainly used for proton resonance frequency shift (PRF) based magnetic resonance temperature imaging. Proton resonance frequency shift magnetic resonance temperature imaging is the mainstream imaging technology of magnetic resonance temperature imaging at present because of better space-time resolution, higher sensitivity and accuracy under high field and near linear sensitivity to higher temperature. The method mainly acquires space phase images through a gradient echo magnetic resonance sequence, acquires phase images before and during heating respectively, and calculates the difference (phase difference image for short) between the phase images to obtain a temperature image. For example, chinese patent publication No. CN102258362a discloses a method for reducing magnetic resonance temperature measurement errors, which is used in a high-intensity focused ultrasound apparatus for magnetic resonance imaging monitoring, and the method includes: before the high-intensity focused ultrasound equipment heats a heating area, acquiring a magnetic resonance phase image as a reference image; during or after the high-intensity focused ultrasound equipment is heated, acquiring another magnetic resonance phase image as a heating image; calculating the temperature change of a heating area according to the heating image and the reference image; the method further comprises the following steps: measuring the magnetic field change caused by the position change of an ultrasonic transducer of the high-intensity focused ultrasound equipment, and compensating the temperature change according to the magnetic field change. For example, chinese patent publication No. CN107468251a discloses a method for correcting low-field magnetic resonance temperature imaging phase drift, which is applied to monitoring the area of a heating region and an accurate temperature change value, and the method includes: before a microwave ablation instrument ablates a target region, a magnetic resonance image is acquired by using a GRE sequence or an SPGR sequence and is used as a reference image; acquiring another magnetic resonance image during or after MW ablation as a heating image; selecting a simulated heating area and an unheated area in the heating image; fitting out phase change caused by non-temperature change in the heating area by using a first-order polynomial model of a weighted least square method according to the phase shift of the unheated area; calculating a temperature difference value according to the phase difference value; and obtaining the area of the heating area according to the phase difference diagram and the temperature difference diagram.
In recent years, a multi-gradient-echo magnetic resonance sequence (the traditional sequence is a single echo) is introduced into magnetic resonance temperature imaging, and the phase diagrams of a plurality of echoes are collected, so that the accuracy of the magnetic resonance temperature imaging is improved to a certain extent. The main disadvantage of phase-based magnetic resonance temperature imaging is that the phase difference map is simultaneously influenced by non-temperature factors, such as the offset of the main magnetic field, and so on, with a certain risk of measurement deviation. Meanwhile, a phase diagram and an amplitude diagram can be obtained simultaneously by a multi-gradient echo magnetic resonance sequence, while the amplitude diagram is usually abandoned in the traditional magnetic resonance temperature imaging, only the phase diagram is used, and the temperature information in the amplitude diagram is not fully utilized.
Disclosure of Invention
The invention aims to disclose a bimodal magnetic resonance temperature measurement method based on a multi-gradient echo sequence, which realizes bimodal temperature imaging by synchronous use of an amplitude diagram and a phase diagram in multi-gradient echo magnetic resonance temperature imaging, and further improves the accuracy of the magnetic resonance temperature imaging.
The invention provides the following technical scheme:
a bimodal magnetic resonance temperature measurement method based on a multi-gradient echo sequence comprises the following steps:
(1) Acquiring magnetic resonance signals before and during heating through a multi-gradient echo magnetic resonance sequence, and simultaneously reserving a phase diagram and an amplitude diagram of each echo;
(2) According to the phase diagrams of different echoes before and during heating, calculating the temperature change before and after heating of each pixel point on the phase diagrams
Figure BDA0002563308500000031
(3) Calculating the signal intensity M when each pixel TE is 0 according to amplitude maps of different echoes in heating;
(4) Calculating the temperature change delta T (M) of each pixel point based on the amplitude value according to the signal intensity of the amplitude diagram before heating and the signal intensity M in the step (3);
(5) According to step (2)
Figure BDA0002563308500000032
And in step (4)
Figure BDA0002563308500000033
And judging the final temperature change delta T of each pixel point or sending out error warning.
In the step (1), before heating, scanning a multi-gradient echo magnetic resonance sequence for 5-10 frames, wherein the number of echoes is N, and N is more than or equal to 2; and averaging the phase diagram and the amplitude diagram of each echo obtained by the multi-gradient-echo magnetic resonance sequence scanning before heating to obtain the phase diagram and the amplitude diagram of each echo before heating.
Preferably, in step (1), the sequence repetition Time (TR) should be as long as possible, e.g. more than 200ms, and the flip angle should be as small as possible, e.g. less than 5 degrees, in order to eliminate the T1 effect.
Preferably, in step (1), if two echoes are used, the first echo time should be as short as possible, e.g. within 2ms, and the second echo time should be as close as possible to the time constant T2 * Thereby improving the estimation accuracy of the signal strength M of step (3).
In the step (2), the temperature change before and after each pixel point on the phase diagram is heated is calculated
Figure BDA0002563308500000036
The method comprises the following steps:
(2-1) according to the formula
Figure BDA0002563308500000034
Calculating the temperature change Delta T (i) of each echo, wherein gamma represents the gyromagnetic ratio of hydrogen proton, B 0 Showing the static magnetic field intensity, TE (i) showing the echo time of the ith echo, alpha showing the hydrogen proton temperature frequency coefficient,
Figure BDA0002563308500000035
phase difference before and after heating of each echo of each pixel point;
(2-2) by
Figure BDA0002563308500000041
Obtaining the temperature change before and after heating of each pixel point based on the phase diagram
Figure BDA0002563308500000042
Preferably, in step (2), a phase correction technique may be employed to improve the accuracy of the phase estimation.
In step (3), the method for calculating the signal intensity M when each pixel TE is 0 includes: according to the formula
Figure BDA0002563308500000049
Fitting the above equation using a nonlinear least squares sum method to obtain M, where S (TE (i)) isThe signal strength of the pixel point in the amplitude map of the ith echo,
Figure BDA0002563308500000043
is the significant transverse relaxation time.
In step (4), the temperature change Δ T (M) is based on
Figure BDA0002563308500000044
Wherein, T 0 Is the actual temperature of the tissue before heating, M 0 Signal intensity of the amplitude plot before heating.
In step (5), if
Figure BDA0002563308500000045
The absolute value of the difference between the sum delta T (M) and the sum delta T (M) is less than 5 degrees, and finally the sum delta T is adopted
Figure BDA0002563308500000046
The final actual tissue temperature T is T = DeltaT + T 0 Otherwise, a measurement error warning is issued.
Wherein, in the step (5),
Figure BDA0002563308500000047
and the acceptance range of the difference value of the sum delta T (M) can be flexibly adjusted according to the requirement of experiment precision.
Wherein, in step (5), if the final demand is that the temperature reaches a certain range rather than an accurate value, it may be set
Figure BDA0002563308500000048
And Δ T (M) must be in this range at the same time to meet the temperature demand.
Compared with the prior art, the invention has the main advantages that:
(1) In the prior art, only a phase diagram of multi-gradient echoes is often adopted, an amplitude diagram is abandoned, and temperature information in magnetic resonance signals is wasted. The method estimates the temperature information from the signal amplitude by reasonably calculating the amplitude diagram.
(2) The method adopts the same sequence to generate temperature images of two modes, thereby realizing the purposes of mutual correction and improving the temperature measurement accuracy.
Drawings
FIG. 1 is a flow chart of magnetic resonance temperature bimodal imaging data acquisition and real-time temperature calculation based on multiple gradient echoes;
FIG. 2 is a flow chart of temperature change before and after heating calculated based on a multi-gradient echo sequence amplitude plot;
FIG. 3 is a logic diagram of a bimodal thermographic co-fusion.
Detailed Description
The invention is described in further detail below with reference to the figures and embodiments (e.g., head imaging).
1. As shown in fig. 1, a scout image, a necessary structural image, and the like are scanned first.
2. As shown in FIG. 1, before heating, a multi-gradient-echo magnetic resonance sequence is scanned in advance for 5-10 frames, the number of echoes is N, and N is required to be more than or equal to 2.
3. As shown in fig. 1, the phase map and the amplitude map of each echo obtained by the multi-gradient echo magnetic resonance sequence scan before heating are averaged to obtain the phase map and the amplitude map of each echo before heating.
4. As shown in fig. 1, the heating apparatus is turned on and the collection of the multi-gradient-echo magnetic resonance sequence in the heating is started.
5. As shown in fig. 1, phase and amplitude maps of a multi-gradient echo magnetic resonance scan in heating are transmitted in real time.
6. As shown in FIG. 1, the phase difference before and after heating of each echo of each pixel point is calculated
Figure BDA0002563308500000051
Where i is the ith echo, starting from 1 to N,
Figure BDA0002563308500000052
to be the phase of the ith echo before heating,
Figure BDA0002563308500000053
the phase of the ith echo after heating.
7. As shown in fig. 1According to the formula
Figure BDA0002563308500000054
Calculating the temperature change Delta T (i) of each echo, wherein gamma represents the gyromagnetic ratio of hydrogen proton, B 0 Showing the static magnetic field intensity, TE (i) showing the echo time of the ith echo, and α showing the hydrogen proton temperature frequency coefficient.
8. As shown in fig. 1, by
Figure BDA0002563308500000061
Obtaining the temperature change before and after heating of each pixel point based on the phase diagram
Figure BDA0002563308500000062
9. As shown in fig. 2, according to the formula
Figure BDA0002563308500000063
Fitting the formula by adopting a nonlinear least square sum method to obtain the signal intensity M of each pixel point before heating 0 . Wherein S (TE (i)) is the signal intensity of the amplitude map of the ith echo of the pixel point before heating,
Figure BDA0002563308500000064
is a significant transverse relaxation time of the material,
Figure BDA0002563308500000065
and M 0 Are fitting parameters.
10. As shown in fig. 2, according to the formula
Figure BDA0002563308500000066
And fitting the formula by adopting a nonlinear least square sum method to obtain the signal intensity M of each pixel point in heating. Wherein S (TE (i)) is the signal intensity of the amplitude map of the ith echo in heating of the pixel point,
Figure BDA0002563308500000067
is a significant transverse relaxation time of the material,
Figure BDA0002563308500000068
and M 0 Are fitting parameters.
11. As shown in fig. 2, according to the formula
Figure BDA0002563308500000069
Calculating a temperature change before and after heating Δ T (M) based on an amplitude map, wherein T 0 Is the actual temperature of the tissue (in Kelvin) before heating, M and M 0 The signal intensities during and before heating calculated in steps 10-11, respectively.
12. As shown in fig. 3, for determining a region of interest
Figure BDA00025633085000000610
Whether or not it is satisfied.
13. As shown in FIG. 3, the temperature change of the area is determined as satisfied by step 12
Figure BDA00025633085000000611
The actual temperature T is T = DeltaT + T 0
14. As shown in fig. 3, if step 12 is not satisfied, a measurement error and termination scan command is issued.
15. As shown in fig. 3, it is determined whether the temperature change Δ T of the area reaches the termination scan request as shown in step 13.
16. As shown in fig. 3, a terminate scan command is issued, as satisfied by step 15.
17. If the instruction of terminating the scanning is not received, the multi-gradient echo magnetic resonance sequence of the next frame is scanned continuously as shown in fig. 1 and 3.
18. As shown in fig. 1, steps 5-17 are repeated.
The above-mentioned embodiments are intended to illustrate the technical solutions and advantages of the present invention, and it should be understood that the above-mentioned embodiments are only the most preferred embodiments of the present invention, and are not intended to limit the present invention, and any modifications, additions, equivalents, etc. made within the scope of the principles of the present invention should be included in the scope of the present invention.

Claims (7)

1. A bimodal magnetic resonance temperature measurement method based on a multi-gradient echo sequence is characterized by comprising the following steps:
(1) Acquiring magnetic resonance signals before and during heating through a multi-gradient echo magnetic resonance sequence, and simultaneously reserving a phase diagram and an amplitude diagram of each echo;
(2) According to the phase diagrams of different echoes before and during heating, calculating the temperature change before and after heating of each pixel point on the phase diagrams
Figure FDA0003994966560000011
(3) According to the amplitude diagram of different echoes during heating, according to the formula
Figure FDA0003994966560000012
Figure FDA0003994966560000013
Calculating the signal intensity of each pixel point TE when the signal intensity is 0, fitting the formula by adopting a nonlinear least square sum method, wherein M is the signal intensity of each pixel point TE when the signal intensity is 0, S (TE (i)) is the signal intensity of the amplitude map of the ith echo of the pixel point,
Figure FDA0003994966560000014
is the significant transverse relaxation time;
(4) Calculating temperature change delta T (M) of each pixel point based on the amplitude value according to the signal intensity of the amplitude diagram before heating and the signal intensity in the step (3);
(5) According to step (2)
Figure FDA0003994966560000015
And the delta T (M) in the step (4) judges the final actual temperature T of each pixel point or sends out error warning.
2. The bimodal magnetic resonance temperature measurement method based on the multi-gradient echo sequence as claimed in claim 1, wherein in the step (1), before heating, the multi-gradient echo magnetic resonance sequence is scanned for 5-10 frames, the number of echoes is N, and N is more than or equal to 2; and averaging the phase map and the amplitude map of each echo obtained by the multi-gradient echo magnetic resonance sequence scanning before heating to obtain the phase map and the amplitude map of each echo before heating.
3. The dual-modality magnetic resonance thermometry method based on multi-gradient echo sequence of claim 1, wherein in step (1), the sequence repetition time TR > 200ms and the flip angle < 5 degrees.
4. The dual-modality MR thermometry method based on multi-gradient echo sequence of claim 1, wherein in step (1) two echoes are used, the first echo time being within 2ms and the second echo time being close to the time constant T2 *
5. The dual-modality MR thermometry method based on multi-gradient echo sequence of claim 1, wherein in step (2), the temperature change before and after heating of each pixel point on the phase map is calculated
Figure FDA0003994966560000021
The method comprises the following steps:
(2-1) according to the formula
Figure FDA0003994966560000022
Calculating the temperature change Delta T () of each echo, wherein gamma represents the gyromagnetic ratio of hydrogen proton, B 0 Showing the static magnetic field intensity, TE (i) showing the echo time of the ith echo, alpha showing the hydrogen proton temperature frequency coefficient,
Figure FDA0003994966560000023
phase difference before and after heating of each echo of each pixel point;
(2-2) by
Figure FDA0003994966560000024
Obtaining the temperature change before and after heating of each pixel point based on the phase diagram
Figure FDA0003994966560000025
6. The dual-modality magnetic resonance thermometry method based on multiple gradient echo sequence of claim 1, wherein in step (4), the temperature change Δ T (M) is based on
Figure FDA0003994966560000026
Wherein, T 0 Is the actual temperature of the tissue before heating, M 0 Signal intensity of the amplitude plot before heating.
7. The dual-modality MR thermometry method based on multi-gradient echo sequence of claim 1, wherein in step (5), if
Figure FDA0003994966560000027
The absolute value of the difference between the sum delta T (M) and the sum delta T (M) is less than 5 degrees, and finally the sum delta T is adopted
Figure FDA0003994966560000028
The final temperature T is T = Δ T + T 0 Otherwise, a measurement error warning is issued.
CN202010621882.9A 2020-06-30 2020-06-30 Bimodal magnetic resonance temperature measurement method based on multi-gradient echo sequence Active CN111714097B (en)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN202010621882.9A CN111714097B (en) 2020-06-30 2020-06-30 Bimodal magnetic resonance temperature measurement method based on multi-gradient echo sequence
PCT/CN2021/102229 WO2022001841A1 (en) 2020-06-30 2021-06-25 Dual-mode magnetic resonance thermometry method based on multi-gradient-echo sequence

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010621882.9A CN111714097B (en) 2020-06-30 2020-06-30 Bimodal magnetic resonance temperature measurement method based on multi-gradient echo sequence

Publications (2)

Publication Number Publication Date
CN111714097A CN111714097A (en) 2020-09-29
CN111714097B true CN111714097B (en) 2023-03-21

Family

ID=72571005

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010621882.9A Active CN111714097B (en) 2020-06-30 2020-06-30 Bimodal magnetic resonance temperature measurement method based on multi-gradient echo sequence

Country Status (2)

Country Link
CN (1) CN111714097B (en)
WO (1) WO2022001841A1 (en)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111714097B (en) * 2020-06-30 2023-03-21 杭州佳量医疗科技有限公司 Bimodal magnetic resonance temperature measurement method based on multi-gradient echo sequence
CN114764133B (en) * 2021-02-08 2023-08-08 华科精准(北京)医疗科技有限公司 Ablation calculation method and ablation calculation system
CN113907871A (en) * 2021-09-13 2022-01-11 哈尔滨医科大学 Tumor ablation microwave device for multi-nuclide MRI real-time detection and using method

Family Cites Families (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE102009049520B4 (en) * 2009-10-15 2015-02-12 Siemens Aktiengesellschaft Multi-echo MR sequence with improved signal-to-noise ratio of the phase information
EP2624004A1 (en) * 2012-02-06 2013-08-07 Koninklijke Philips Electronics N.V. Temperature determination using magnetic resonance B1 field mapping
CN103284722B (en) * 2013-04-18 2015-01-21 鑫高益医疗设备股份有限公司 Temperature measurement method of permanent-magnet MRI system
CN104132736B (en) * 2014-07-31 2017-02-15 华中科技大学 Measuring method for non-intrusion type rapid temperature changes in direct-current excitation magnetic field
CN104257357B (en) * 2014-10-14 2016-05-25 深圳先进技术研究院 magnetic resonance temperature imaging method and system
CN104997512A (en) * 2015-07-07 2015-10-28 中国科学院上海技术物理研究所 MRI (magnetic resonance imaging) temperature measurement magnetic field drift phase compensation method
CN107468251A (en) * 2017-07-03 2017-12-15 中国科学技术大学 A kind of bearing calibration of Low Magnetic field MRI temperature imaging phase drift
CN111714097B (en) * 2020-06-30 2023-03-21 杭州佳量医疗科技有限公司 Bimodal magnetic resonance temperature measurement method based on multi-gradient echo sequence

Also Published As

Publication number Publication date
CN111714097A (en) 2020-09-29
WO2022001841A1 (en) 2022-01-06

Similar Documents

Publication Publication Date Title
CN111714097B (en) Bimodal magnetic resonance temperature measurement method based on multi-gradient echo sequence
US6559644B2 (en) MRI-based temperature mapping with error compensation
EP1234190B1 (en) Thermal imaging of fat and muscle
US10782373B2 (en) System and method for image processing with highly undersampled imaging data
US7906964B2 (en) Method and system for determining acquisition parameters associated with magnetic resonance imaging for a particular measurement time
US20020180436A1 (en) System and method for providing real-time motion correction by utilizing navigators
US20090030302A1 (en) Magnetic resonance imaging device
US10371779B2 (en) Apparatus and method for magnetic resonance imaging with high spatial temporal resolutions
WO2001017428A1 (en) Magnetic resonance imaging device and method therefor
CN103284722B (en) Temperature measurement method of permanent-magnet MRI system
US9207299B2 (en) Magnetic resonance imaging system and magnetic resonance imaging apparatus
US20220179023A1 (en) System and Method for Free-Breathing Quantitative Multiparametric MRI
US20170281042A1 (en) Systems and methods for magnetic resonance thermometry using balanced steady state free precession
CN104224180B (en) Temperature measurement method for in-vivo fat based on magnetic resonance imaging
CN111568390B (en) Method for reducing system temperature error in real-time magnetic resonance temperature imaging
JP4574781B2 (en) Magnetic resonance apparatus and thermotherapy apparatus
JP2000279390A (en) Magnetic resonance imaging device
JPWO2016021440A1 (en) Magnetic resonance imaging system
CN102565736B (en) Determine that layer is about the method for the position at the position of relatively its motion and magnetic resonance equipment
Das et al. Improved magnetic resonance thermal imaging by combining proton resonance frequency shift (PRFS) and apparent diffusion coefficient (ADC) data
CN114690101B (en) Quantitative detection method and device for parameters in magnetic resonance imaging and magnetic resonance scanner
JP2899649B1 (en) Magnetic resonance imaging equipment
Zhang et al. Hybrid Proton Resonance Frequency Shift and Variable Flip Angle T1 Temperature Mapping using a Golden-Angle 3D Stack-of-Radial Technique
CN112472041A (en) Magnetic resonance temperature imaging method for jointly measuring temperature of water tissue and fat tissue
EP3443371B1 (en) Parallel transmit rf power control for nmr imaging of a subject and of invasive conductive objects.

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant