Magnetic resonance imaging device, phase value correction method and program
View Patent ↗An object of the present invention is to suppress artifacts generated by correction of spectral distortion induced by eddy currents in MRI devices with a simple method, and thereby improve accuracy of the correction.
1. A magnetic resonance imaging device comprising:
a static magnetic field application part for applying a static magnetic field to a subject;
a gradient magnetic field application part for applying a gradient magnetic field to the subject;
a radio frequency magnetic field pulse irradiation part for irradiating a radio frequency magnetic field pulse on the subject;
a reception part for receiving magnetic resonance signals generated from the subject; and
a control part;
wherein the control part comprises:
a measurement control part for controlling operations of the gradient magnetic field application part, the radio frequency magnetic field pulse irradiation part, and the reception part to obtain a magnetic resonance signal of a desired metabolite for every measurement point,
an eddy current correction part for performing eddy current correction of the magnetic resonance signal, and
a display information generation part for generating display information from the magnetic resonance signal for every measurement point corrected by the eddy current correction part,
wherein the eddy current correction part comprises:
a phase value calculation part for calculating a phase value of an FID signal of a substance for correction showing a larger signal intensity compared with a metabolite as a measurement object for every measurement point, and
a phase value correction part for correcting a phase jump of the phase value to obtain a corrected phase value,
wherein the phase value correction part comprises:
a primary differential value calculation part for calculating a primary time differential value of the phase value for every measurement point,
a threshold value calculation part for calculating a threshold value for identifying a phase jump generation region where a phase jump is generated,
wherein the threshold value calculation part sets a redetermined region as a threshold value calculation region and calculates an absolute value of a difference of the maximum value and the minimum value of the primary time differential values in the threshold value calculation region as the threshold value,
a phase jump generation region identification part for identifying the phase jump generation region of the phase value using the threshold value and the primary time differential value, and
a phase jump correction part for correcting the phase jump of the phase value by correcting the primary time differential value of the phase jump generation region, and
wherein the eddy current correction part performs the eddy current correction by using the phase jump-corrected phase value.
2. The magnetic resonance imaging device according to claim 1 , wherein: the threshold value calculation region is a region from a measurement start time to a predetermined time.
3. The magnetic resonance imaging device according to claim 2 , wherein the predetermined time is a time when absolute value of signal intensity of the FID signal of the substance for correction first becomes the predetermined value.
4. The magnetic resonance imaging device according to claim 1 , wherein the threshold value calculation part divides a sequence of the primary time differential values into a plurality of small regions in the direction of time, calculates standard deviation of the primary time differential values included in each small region for every small region, identifies the smallest standard deviation out of all the calculated standard deviations, and calculates a value corresponding to the identified smallest standard deviation multiplied with a predetermined coefficient as the threshold value.
5. The magnetic resonance imaging device according to claim 1 , wherein the phase jump generation region identification part calculates an absolute value of a difference of the maximum value and the minimum value of the primary time differential values included in the evaluation value calculation region, which is a predetermined time width of which center is a predetermined evaluation point, as an evaluation value, compares the evaluation value calculated for every measurement point as the evaluation point with the threshold value, and identifies measurement points at which the evaluation value is larger than the threshold value as the phase jump generation region, and the other measurement points as the non-phase jump generation region.
6. The magnetic resonance imaging device according to claim 5 , wherein the phase jump correction part connects the primary time differential values of the non-phase jump generation regions by interpolation to correct the phase jump.
7. The magnetic resonance imaging device according to claim 6 , wherein the interpolation is linear interpolation for connecting the primary time differential values at the ends of the adjacent non-phase jump generation regions with a straight line.
8. The magnetic resonance imaging device according to claim 6 , wherein the phase jump correction part gives a power sum of the primary time differential values obtained after the interpolation as the phase jump-corrected phase value.
9. The magnetic resonance imaging device according to claim 6 , wherein the phase jump correction part gives a value of a function obtained by integrating the primary time differential values obtained by fitting of the primary time differential values obtained after the interpolation corresponding to each measurement point as a phase jump-corrected phase value.
10. A phase value correction method for correcting a phase jump of a phase value in a magnetic resonance imaging device, comprising:
applying, by a static magnetic field application part, a static magnetic field to a subject,
applying, by a gradient magnetic field application part, a gradient magnetic field to the subject,
irradiating, by a radio frequency magnetic field pulse irradiation part, a radio frequency magnetic field pulse on the subject,
receiving, by a reception part, magnetic resonance signals generated from the subject,
controlling, by a measurement control part, operations of the gradient magnetic field application part, the radio frequency magnetic field pulse irradiation part, and the reception part to obtain a magnetic resonance signal of a desired metabolite for every measurement point,
performing, by an eddy current correction part, eddy current correction of the magnetic resonance signal by using a phase value of an FID signal of a substance for correction showing a larger signal intensity compared with a metabolite as a measurement object for every measurement point, and
generating, by a display information generation part, display information from the magnetic resonance signal for every measurement point corrected by the eddy current correction part, which comprises:
a primary differential value calculation step of calculating a primary time differential value of the phase value for every measurement point,
a threshold value calculation step of calculating a threshold value for identifying a phase jump generation region where a phase jump is generated,
wherein the threshold value calculation step sets a predetermined region as a threshold value calculation region and calculates an absolute value of a difference of the maximum value and the minimum value of the primary time differential values in the threshold value calculation region as the threshold value,
a phase jump generation region identification step of identifying the phase jump generation region of the phase value using the threshold value and the primary time differential value,
a phase jump correction step of correcting the phase jump of the phase value by correcting the primary time differential value of the phase jump generation region, and
a corrected phase value calculation step of obtaining a corrected phase value from the corrected primary time differential value.
11. A program for operating a computer of a magnetic resonance imaging device comprising:
a static magnetic field application part for applying a static magnetic field to a subject;
a gradient magnetic field application part for applying a gradient magnetic field to the subject;
a radio frequency magnetic field pulse irradiation part for irradiating a radio frequency magnetic field pulse on the subject;
a reception part for receiving magnetic resonance signals generated from the subject;
a measurement control part for controlling operations of the gradient magnetic field application part, the radio frequency magnetic field pulse irradiation part, and the reception part to obtain a magnetic resonance signal of a desired metabolite for every measurement point;
an eddy current correction part for performing eddy current correction of the magnetic resonance signal by using a phase value of an FID signal of a substance for correction showing a larger signal intensity compared with a metabolite as a measurement object for every measurement point; and
a display information generation part for generating display information from the magnetic resonance signal for every measurement point corrected by the eddy current correction part, as:
a primary differential value calculation part for calculating a primary time differential value of the phase value for every measurement point,
a threshold value calculation part for calculating a threshold value for identifying a phase jump generation region where a phase jump is generated,
wherein the threshold value calculation part sets a and calculates an absolute value of a difference of the maximum value and the minimum value of the primary time differential values in the threshold value calculation region as the threshold value,
a phase jump generation region identification part for identifying the phase jump generation region of the phase value using the threshold value and the primary time differential value, and
a phase jump correction part for correcting the phase jump of the phase value by correcting the primary time differential value of the phase jump generation region and obtaining a corrected phase value from the corrected primary time differential value.