US20120249143A1

Magnetic resonance imaging apparatus

Claim Score by NHIP

Abstract

In a magnetic resonance imaging apparatus according to one embodiment, an executing unit executes a first pre-scan and a second pre-scan, each being a pre-scan in which readout gradient magnetic fields and slice direction gradient magnetic fields are applied in an identical manner to a pulse sequence for main-scanning and in which phase encode gradient magnetic fields are applied in an identical manner to the pulse sequence for main-scanning up to just before echoes used in calculating a correction amount, and each having different predetermined imaging parameters; a calculating unit calculates, as a correction amount, an amount of phase shifting by referring to phase differences present in a plurality of echo signals that are collected during the first pre-scan and the second pre-scan; and a correcting unit corrects the pulse sequence for main-scanning based on the correction amount calculated by the calculating unit.

US20120249143A1, drawing sheet 1 — Sheet 1 of 14

Term

7.4 yearsto projected expiry

Projected expiry 12 February 2034, counted from filing; an application has no term until it is granted.

20 claims: 4 independent, 16 dependent

1
A magnetic resonance imaging apparatus comprising:an executing unit configured to execute a first pre-scan and a second pre-scan, each being a pre-scan in which a readout gradient magnetic field and a slice direction gradient magnetic field are applied in an identical manner to a pulse sequence for main-scanning and in which a phase encode gradient magnetic field is applied in an identical manner to the pulse sequence for main-scanning up to just before echoes used in calculating a correction amount, and each having different predetermined imaging parameters;a calculating unit configured to calculate, as the correction amount, an amount of phase shifting by referring to phase differences present in a plurality of echo signals that are collected during the first pre-scan and the second pre-scan;and a correcting unit configured to correct the pulse sequence for main-scanning based on the correction amount calculated by the calculating unit.
18
A magnetic resonance imaging apparatus comprising:an RF coil including a plurality of coil elements;an executing unit configured to execute a pre-scan in which a readout gradient magnetic field and a slice direction gradient magnetic field are applied in an identical manner to a pulse sequence for main-scanning;a calculating unit configured to add a plurality of echo signals collected via the plurality of coil elements during the pre-scan and calculates a single phase correction amount to be used in correcting phase differences present in a plurality of echo signals collected in a main-scan;and a correcting unit configured to correct the pulse sequence for main-scanning based on the phase correction amount calculated by the calculating unit.
19
A magnetic resonance imaging apparatus comprising:an executing unit configured to execute a pre-scan in which a readout gradient magnetic field and a slice direction gradient magnetic field are applied in an identical manner to a pulse sequence for main-scanning;a calculating unit configured to, of a plurality of echo signals collected during a pre-scan, make use of data within a field of view specified in a readout direction by an operator, and calculates a phase correction amount to be used in correcting phase differences present in a plurality of echo signals collected in a main-scan;and a correcting unit configured to correct the pulse sequence for main-scanning based on the phase correction amount calculated by the calculating unit.
20
Broadest claimClaim Score 59, broad(NHIP)A magnetic resonance imaging apparatus comprising:an executing unit configured to execute a pre-scan in which a readout gradient magnetic field and a slice direction gradient magnetic field are applied in an identical manner to a pulse sequence for main-scanning;a calculating unit configured to, of a plurality of echo signals collected during a pre-scan, make use of data within a range compensated at least for magnetostatic field unevenness, and calculates a phase correction amount to be used in correcting phase differences present in a plurality of echo signals collected in a main-scan;and a correcting unit configured to correct the pulse sequence for main-scanning based on the phase correction amount calculated by the calculating unit.

Independent claims4

142 paragraphs in 5 sections, as filed

0001This application is based upon and claims the benefit of priority from Japanese Patent Application No. 2011-080111, filed on Mar. 31, 2011; and Japanese Patent Application No. 2012-045281, filed on Mar. 1, 2012, the entire contents of all of which are incorporated herein by reference.

FIELD

0002Embodiments described herein relate generally to a magnetic resonance imaging apparatus.

BACKGROUND

0003Typically, as an imaging method related to a magnetic resonance imaging apparatus, the fast spin echo (FSE) method is known. The FSE method is an imaging method in which firstly a flip pulse is applied to a person being tested and then a plurality of flop pulses is sequentially applied the person being tested so as to collect a plurality of echo signals called an echo train. Herein, a flip pulse is a radio frequency (RF) pulse used for excitation of the nucleus spin inside the person being tested; and flop pulses are RF pulses used for refocusing the phase of the nucleus.

0004In this FSE method, since a plurality of RF pulses is applied, the generation of spin echoes is accompanied by the generation of stimulated echoes. Because of the stimulated echoes, sometimes phase shifts occur in the echo signals that are collected. Such phase shifts cause degradation in the image quality in the form of unevenness in the photographic sensitivity, signal deterioration, or ghost image formation.

0005In order to prevent such degradation in the image quality; generally, a main-scan is preceded by a pre-scan so as to measure the phase differences occurring in the echo signals. Based on the phase differences measured during the pre-scan, the pulse sequence for main-scanning is corrected. In that case, for example, during a pre-scan, a pulse sequence for cancelling the stimulated echoes is implemented and only the spin echoes are collected. Then, of the spin echoes collected during the pre-scan, the first echo signal and the second echo signal are subjected to Fourier conversion in the readout direction, and a zero-dimensional phase shift and a one-dimensional phase shift between the first echo signal and the second echo signal are calculated. Subsequently, based on the zero-dimensional phase shift and the one-dimensional phase shift, a correction amount is calculated for correcting the phase shifts in the readout direction and the slice selection direction, and the pulse sequence for main-scanning is modified according to the calculated correction amount.

0006However, in the technology described above, it is not possible to correct the phase shifts that occur in the phase encode direction due to the effect of the eddy currents generated by phase encode gradient magnetic fields. Such phase shifts sometimes cause degradation in the image quality.

BRIEF DESCRIPTION OF THE DRAWINGS

0007<figref idref="DRAWINGS">FIG. 1</figref> is a diagram illustrating a configuration of an MRI apparatus according to a first embodiment;

0008<figref idref="DRAWINGS">FIG. 2</figref> is a functional block diagram illustrating a detailed configuration of a calculating system illustrated in <figref idref="DRAWINGS">FIG. 1</figref>;

0009<figref idref="DRAWINGS">FIG. 3</figref> is a diagram illustrating a pulse sequence for main-scanning according to the first embodiment;

0010<figref idref="DRAWINGS">FIG. 4</figref> is a diagram illustrating a pulse sequence for pre-scanning according to the first embodiment;

0011<figref idref="DRAWINGS">FIGS. 5 and 6</figref> are explanatory diagrams for explaining a phase shift correction operation performed by the sequence correcting unit;

0012<figref idref="DRAWINGS">FIG. 7</figref> is a flowchart for explaining a sequence of operations in the phase shift correction operation performed in the MRI apparatus according to the first embodiment;

0013<figref idref="DRAWINGS">FIG. 8</figref> is an explanatory diagram for explaining a modification example of the first embodiment;

0014<figref idref="DRAWINGS">FIG. 9A</figref> is a diagram illustrating an exemplary image captured during the conventional pre-scan;

0015<figref idref="DRAWINGS">FIG. 9B</figref> is a diagram illustrating an exemplary image captured during the pre-scan according to the first embodiment;

0016<figref idref="DRAWINGS">FIG. 10</figref> is a diagram illustrating a second pulse sequence for the second pre-scan according to a second embodiment;

0017<figref idref="DRAWINGS">FIG. 11</figref> is an explanatory diagram for explaining the calculation of phase shifts performed by the sequence correcting unit according to the second embodiment;

0018<figref idref="DRAWINGS">FIG. 12</figref> is an explanatory diagram for explaining a modification example of the first and second embodiments;

0019<figref idref="DRAWINGS">FIG. 13</figref> is a diagram illustrating the pulse sequence for pre-scanning according to a third embodiment;

0020<figref idref="DRAWINGS">FIG. 14</figref> is an explanatory diagram for explaining collection of echo signals according to the third embodiment;

0021<figref idref="DRAWINGS">FIG. 15</figref> is a diagram illustrating the second pulse sequence for the second pre-scan according to a fourth embodiment;

0022<figref idref="DRAWINGS">FIG. 16</figref> is an explanatory diagram for explaining the MRI apparatus according to a fifth embodiment;

0023<figref idref="DRAWINGS">FIG. 17</figref> is an explanatory diagram for explaining the MRI apparatus according to a sixth embodiment; and

0024<figref idref="DRAWINGS">FIG. 18</figref> is an explanatory diagram for explaining the MRI apparatus according to a seventh embodiment.

DETAILED DESCRIPTION

0025According to one embodiment, a magnetic resonance imaging apparatus includes an executing unit, a calculating unit, and a correcting unit. The executing unit executes a first pre-scan and a second pre-scan, each being a pre-scan in which readout gradient magnetic fields and slice direction gradient magnetic fields are applied in an identical manner to a pulse sequence for main-scanning and in which phase encode gradient magnetic fields are applied in an identical manner to the pulse sequence for main-scanning up to just before echoes used in calculating a correction amount, and each having different predetermined imaging parameters. The calculating unit calculates, as a correction amount, an amount of phase shifting by referring to phase differences present in a plurality of echo signals that are collected during the first pre-scan and the second pre-scan. The correcting unit corrects the pulse sequence for main-scanning based on the correction amount calculated by the calculating unit.

0026Exemplary embodiments of a magnetic resonance imaging apparatus are described in detail below with reference to the accompanying drawings. In the following embodiments, the magnetic resonance image apparatus is referred to as MRI apparatus.

First Embodiment

0027<figref idref="DRAWINGS">FIG. 1</figref> is a diagram illustrating a configuration of the MRI apparatus according to a first embodiment. As illustrated in <figref idref="DRAWINGS">FIG. 1</figref>, an MRI apparatus 100 includes a magnetostatic magnet 1, a gradient coil 2, a gradient power source 3, a couch 4, a couch control unit 5, a transmitter RF coil 6, a transmitting unit 7, a receiver RF coil 8, a receiving unit 9, a sequence control unit 10, and a calculating system 20.

0028The magnetostatic magnet 1 is a hollow cylindrical magnet and forms a uniform magnetostatic field in the internal space thereof. As the magnetostatic magnet 1, for example, a permanent magnet or a superconductive magnet is used.

0029The gradient coil 2 is a hollow cylindrical coil disposed on the inside of the magnetostatic magnet 1. The gradient coil 2 is formed as a combination of three coils corresponding to the mutually-orthogonal X, Y, and Z axes. The three coils independently receive a current supply from the gradient power source 3 (described later) and generate gradient magnetic fields having varying magnetic field intensities along the X, Y, and Z axes. The gradient power source 3 supplies current to the gradient coil 2.

0030Herein, the gradient magnetic fields generated along the X, Y, and Z axes by the gradient coil 2 correspond to, for example, a slice selection gradient magnetic field Gs, a phase encode gradient magnetic field Ge, and a readout gradient magnetic field Gr. The slice selection gradient magnetic field Gs is used in arbitrarily deciding the imaging cross-section. The phase encode gradient magnetic field Ge is used in changing the phases of magnetic resonance signals according to spatial positions. The readout gradient magnetic field Gr is used in changing the frequencies of magnetic resonance signals according to spatial positions.

0031The couch 4 has a couchtop 4a on which a person being tested P is asked to lie. Under the control of a couch control unit 5 (described later), the couchtop 4a on which the person being tested P is lying is moved inside the hollow cavity (imaging opening) of the gradient coil 2. Usually, the orientation of the couch 4 is such that the longitudinal direction thereof is parallel to the magnetostatic magnet 1. The couch control unit 5 operates under the control of a control unit 26 and controls the couch 4. More particularly, the couch control unit 5 drives the couch 4 and moves the couchtop 4a in the longitudinal direction and in the vertical direction.

0032The transmitter RF coil 6 is disposed on the inside of the gradient coil 2, receives radiofrequency pulses from the transmitting unit 7, and generates high-frequency magnetic fields. The transmitting unit 7 transmits radiofrequency pulses corresponding to the Larmor frequency to the transmitter RF coil 6. The receiver RF coil 8 is disposed on the inside of the gradient coil 2 and receives magnetic resonance signals radiated from the person being tested P due to the effect of the high-frequency magnetic fields mentioned above. Upon receiving the magnetic resonance signals, the receiver RF coil 8 outputs those magnetic resonance signals to the receiving unit 9.

0033Based on the magnetic resonance signals output by the receiver RF coil 8, the receiving unit 9 generates magnetic resonance (MR) signal data. More particularly, the receiving unit 9 generates MR signal data by performing digital conversion of the magnetic resonance signals output by the receiver RF coil 8. In the MR signal data, information regarding spatial frequencies in the phase encode direction, the readout direction, and the slice encode direction is disposed in the k-space by the slice selection gradient magnetic field Gs, the phase encode gradient magnetic field Ge, and the readout gradient magnetic field Gr. Once the MR signal data is generated, the receiving unit 9 sends that MR signal data to the sequence control unit 10.

0034Based on sequence execution data sent by the calculating system 20, the sequence control unit 10 scans the person being tested P by means of driving the gradient power source 3, the transmitting unit 7, and the receiving unit 9. Herein, the sequence execution data points to the information that defines the pulse sequence indicating the sequence for scanning the person being tested P. More particularly, the sequence execution data includes the strength or the timing of the power supply from the gradient power source 3 to the gradient coil 2, the intensity or the timing of transmission of the RF signals from the transmitting unit 7 to the transmitter RF coil 6, and the timing of detection of magnetic resonance signals by the receiving unit 9. Moreover, once the gradient power source 3, the transmitting unit 7, and the receiving unit 9 are driven; when the MR signal data is sent by the receiving unit 9, the sequence control unit 10 transfers that MR signal data to the calculating system 20.

0035The calculating system 20 performs the overall control of the MRI apparatus 100. For example, the calculating system 20 drives all constituent elements of the MRI apparatus 100, performs scanning of the person being tested P, and performs image reconstruction. The calculating system 20 includes an interface unit 21, an image reconstructing unit 22, a memory unit 23, an input unit 24, a display unit 25, and the control unit 26.

0036The interface unit 21 controls the input-output of various signals communicated with the sequence control unit 10. For example, the interface unit 21 sends sequence execution data to the sequence control unit 10 and receives MR signal data from the sequence control unit 10. Upon receiving the MR signal data, the interface unit 21 stores each set of MR signal data in a corresponding manner to each person being tested P in the memory unit 23.

0037With respect to each set of MR signal data stored in the memory unit 23, the image reconstructing unit 22 performs post-processing, that is, reconstruction processing such as Fourier conversion, and generates spectrum data or image data of the desired nucleus spin inside the corresponding person being tested P.

0038The memory unit 23 stores therein a variety of data or various computer programs that are required for the operations performed by the control unit 26 (described later). For example, with respect to each person being tested P, the memory unit 23 stores therein the MR signal data received by the interface unit 21 and the spectrum data or the image data generated by the image reconstructing unit 22.

0039The input unit 24 receives various instructions or information input by the operator. As the input unit 24; it is possible to appropriately use a pointing device such as a mouse or a trackball, a selection device such as a mode switch, or an input device such as a keyboard.

0040The display unit 25 operates under the control of the control unit 26 and displays a variety of information such as the spectrum data or the image data. As the display unit 25, a display device such as a liquid crystal display device can be used.

0041The control unit 26 includes a central processing unit (CPU) and a memory, and comprehensively controls the MRI apparatus 100. For example, the control unit 26 controls the scanning operation by generating a variety of sequence execution data according to the imaging conditions, which are input by the operator via the input unit 24, and sending that sequence execution data to the sequence control unit 10. Moreover, when the sequence control unit 10 sends MR signal data as a scan result, the control unit 26 controls the image reconstructing unit 22 to perform image reconstruction based on that MR signal data.

0042Thus, the MRI apparatus 100 according to the first embodiment has the configuration described above. With such a configuration, the MRI apparatus 100 performs pre-scanning. Herein, during a pre-scan, readout gradient magnetic fields and slice direction gradient magnetic fields are applied in an identical manner to the pulse sequence for main-scanning, and phase encode gradient magnetic fields are applied in an identical manner to the pulse sequence for main-scanning up to just before the echoes used in calculating the correction amount. Then, as the correction amount, the MRI apparatus 100 calculates the amount of phase shifting by referring to the phase differences present in a plurality of echo signals obtained during the pre-scan. Moreover, the MRI apparatus 100 corrects the pulse sequence for main-scanning based on the correction amount calculated by a correction amount calculating unit 26c (described later).

0043Thus, the MRI apparatus 100 according to the first embodiment performs a pre-scan in which phase encode gradient magnetic fields are applied and, based on the phase shifts among the echo signals obtained during the pre-scan, calculates a correction amount regarding the pulse sequence for main-scanning. Hence, according to the first embodiment, it becomes possible to prevent image deterioration that may occur due to the phase shifts present in the readout direction and the slice direction as well as due to zero-dimensional phase shifts occurring due to the phase encode gradient magnetic fields. Given below is the detailed explanation of the functions of the MRI apparatus 100 according to the first embodiment.

0044<figref idref="DRAWINGS">FIG. 2</figref> is a functional block diagram illustrating a detailed configuration of the calculating system 20 illustrated in <figref idref="DRAWINGS">FIG. 1</figref>. In <figref idref="DRAWINGS">FIG. 2</figref>, from among the functions of the calculating system 20; the interface unit 21, the memory unit 23, and the control unit 26 are illustrated.

0045The memory unit 23 includes a sequence-execution-data memory unit 23a and an MR-signal-data memory unit 23b. The sequence-execution-data memory unit 23a stores therein sequence execution data that is generated by an imaging condition setting unit 26a (described later). The MR-signal-data memory unit 23b stores therein MR signal data received by the interface unit 21.

0046The control unit 26 includes the imaging condition setting unit 26a, a pre-scan executing unit 26b, the correction amount calculating unit 26c, a sequence correcting unit 26d, and a main-scan executing unit 26e.

0047The imaging condition setting unit 26a generates, based on imaging conditions input by the operator via the input unit 24, sequence execution data that defines the pulse sequence used in imaging. For example, when the operator inputs imaging conditions for the FSE method, the imaging condition setting unit 26a generates sequence execution data regarding the pulse sequence for main-scanning (described later) as well as regarding the pulse sequence for pre-scanning (described later).

0048<figref idref="DRAWINGS">FIG. 3</figref> is a diagram illustrating a pulse sequence for main-scanning according to the first embodiment. In <figref idref="DRAWINGS">FIG. 3</figref>, “RF” represents the timings at which a flip pulse for excitation and flop pulses for refocusing are applied. “Gss” represents the timings of applying the slice selection gradient magnetic fields as well as the intensities of the slice selection gradient magnetic fields. “Gro” represents the timings of applying the readout gradient magnetic fields as well as the intensities of the readout gradient magnetic fields. “Gpe” represents the timings of applying the phase encode gradient magnetic fields as well as the intensities of the phase encode gradient magnetic fields. Meanwhile, in <figref idref="DRAWINGS">FIG. 3</figref>, the pulse sequence related to only a single slice selection is illustrated. Regarding the slice encode, diagrammatic illustration is not given. Moreover, “ETS (echo train spacing)” represents echo intervals.

0049As illustrated in <figref idref="DRAWINGS">FIG. 3</figref>, the pulse sequence for main-scanning is same as the pulse sequence for the common FSE method. As illustrated in <figref idref="DRAWINGS">FIG. 3</figref>, in the pulse sequence for main-scanning, after a flip pulse fli is applied, a plurality of flop pulses flo1, flo2, . . . , flo9, flo10, flo11, . . . is sequentially applied; and a plurality of echo signals Echo1, Echo2, . . . , Echo9, Echo10, Echo11, . . . is collected. The pulse sequence illustrated in <figref idref="DRAWINGS">FIG. 3</figref> is an example of a case when the phase encode gradient magnetic field becomes zero at the echo signal Echo10 collected at the 10-th count.

0050<figref idref="DRAWINGS">FIG. 4</figref> is a diagram illustrating a pulse sequence for pre-scanning according to the first embodiment. In the pulse sequence for pre-scanning, in an identical manner to the pulse sequence for main-scanning illustrated in <figref idref="DRAWINGS">FIG. 3</figref>, the readout gradient magnetic fields and the slice direction gradient magnetic fields are applied, and the phase encode magnetic gradient fields are applied in an identical manner to the pulse sequence for main-scanning up to just before the echoes used in calculating the correction amount.

0051For example, as illustrated in <figref idref="DRAWINGS">FIG. 4</figref>, in the pulse sequence for pre-scanning, in an identical manner to the pulse sequence for main-scanning illustrated in <figref idref="DRAWINGS">FIG. 3</figref>, the readout gradient magnetic fields and the slice direction gradient magnetic fields are applied; and, of the phase encode magnetic gradient fields applied in the pulse sequence for main-scanning, at least a single phase encode magnetic gradient field is applied up to the echo in the vicinity of the k-space center. In other words, with reference to the pulse sequence for main-scanning, in the pulse sequence for pre-scanning, the phase encode gradient magnetic fields are applied up to the vicinity of the contrast TE and are not applied after the contrast TE. Meanwhile, in the example illustrated in <figref idref="DRAWINGS">FIG. 4</figref>, Echo10 corresponds to the K-space center.

0052For example, in the pulse sequence for pre-scanning, in an identical manner to the man-scan pulse sequence, the readout gradient magnetic fields and the slice selection gradient magnetic fields are applied; and, of a plurality of phase encodes used in the pulse sequence for main-scanning, the phase encode gradient magnetic field in the vicinity of the average intensity is applied.

0053Moreover, in the pre-scan sequence, the settings are such that the stimulated echoes get cancelled and only the spin echoes are collected. For example, it is possible to implement the method disclosed in U.S. Pat. No. 5,818,229. In that method, only spin echo components are extracted by adding the echo signals at the first shot, which are collected by varying the phases of the flop pulses for refocusing as π, π, π, π, . . . ; and the echo signals at the second shot, which are collected by varying the phases of the flop pulses for refocusing as π, −π, π, −π, . . . . Alternatively, only stimulated echo components are extracted by subtracting the echo signals at the second shot from the echo signals at the first shot. Then, the extracted stimulated echo components can be used in place of the spin echo components.

0054Herein, for example, in the case of a phase encode Matrix of 256 and echo collection of 19 echoes, it is necessary to have 256/19=13 shots. In that case, in the pulse sequence for pre-scanning, the setting is such that the phase encode gradient magnetic fields having the shot with the average phase encode gradient field intensity applied thereto (in the case of filling the phase encodes sequentially, the central shot (shot 7)) are applied.

0055Returning to the explanation with reference to <figref idref="DRAWINGS">FIG. 2</figref>, the pre-scan executing unit 26b executes the pulse sequence for pre-scanning.

0056More particularly, when the imaging condition setting unit 26a generates the sequence execution data for pre-scanning, the pre-scan executing unit 26b firstly reads the sequence execution data from the sequence-execution-data memory unit 23a. Then, the pre-scan executing unit 26b executes the pre-scan by sending the sequence execution data to the sequence control unit 10 via the interface unit 21. For example, the pre-scan executing unit 26b executes a pre-scan by sending the sequence execution data, which defines the pulse sequence illustrated in <figref idref="DRAWINGS">FIG. 4</figref>, to the sequence control unit 10.

0057The correction amount calculating unit 26c calculates, as the correction amount, the amount of phase shifting by referring to the phase differences present in a plurality of echo signals collected during the pre-scan. That is, from the phase differences present in a plurality of echo signals collected during the pre-scan, the correction amount calculating unit 26c calculates, as the correction amount, the amount of phase shifting (zero-dimensional and one-dimensional) occurring due to the readout gradient magnetic fields and the slice selection gradient magnetic fields, and the amount of zero-dimensional phase shifting occurring due to the phase encode gradient magnetic fields.

0058More particularly, after the pre-scan executing unit 26b performs the pre-scan, the correction amount calculating unit 26c reads, from the MR-signal-data memory unit 23b, the MR signal data related to the echo signals collected during the pre-scan. Then, the correction amount calculating unit 26c performs Fourier conversion on the MR signal data in the readout direction and calculates one-dimensional phase differences. Those one-dimensional phase differences represent the one-dimensional phase differences of the readout gradient magnetic fields.

0059The phase differences occurring on an echo-by-echo basis in the FSE method include one-dimensional phase differences serving as direct functions of location and phase differences without location dependency occurring due to a mismatch in coil arrangement. The phase differences without location dependency are called the zero-dimensional phase differences. The one-dimensional phase calculated with respect to each set of MR signal data is corrected and the zero-dimensional phase differences are calculated by obtaining the phase average. The zero-dimensional phase differences include the effect of not only the zero-dimensional phase differences of the readout gradient magnetic fields and the zero-dimensional and one-dimensional phase differences of the slice selection gradient magnetic fields, but also the zero-dimensional phase differences of the phase encode gradient magnetic fields.

0060The correction amount calculating unit 26c calculates the one-dimensional and zero-dimensional phase differences in the echo signals collected during the pre-scan. Such phase differences correspond to the phase shifts occurring due to the readout gradient magnetic fields, the slice selection gradient magnetic fields, and the phase encode gradient magnetic fields.

0061For example, assume that the pulse sequence illustrated in <figref idref="DRAWINGS">FIG. 4</figref> is executed. In that case, the correction amount calculating unit 26c calculates a phase difference p1 between the echo signal Echo10 collected at the 10-th count and the echo signal Echo11 collected at the 11-th count illustrated in <figref idref="DRAWINGS">FIG. 4</figref>. Then, as the correction amount, the correction amount calculating unit 26c calculates a phase shift p1 from the calculated phase difference.

0062Meanwhile, herein, although the 10-th echo signal and the 11-th echo signal are used in calculating the correction amount, it is also possible to further use a plurality of echo signals from the 12-th echo signal onward. For example, assuming that the phase difference between the 10-th echo signal and the 11-th echo signal is the phase difference p1, assuming that the phase difference between the 12-th echo signal and the 11-th echo signal is a phase difference p1_2, and assuming that the phase difference between the 12-th echo signal and the 13-th echo signal is a phase difference p1_3; the correction amount calculating unit 26c can set the phase shift, which occurs in the phase encode direction due to the phase encode gradient magnetic fields, as (p1+p1_2+p1_3)/3.

0063Based on the phase shift calculated by the correction amount calculating unit 26c, the sequence correcting unit 26d corrects the pulse sequence for main-scanning. More particularly, when the correction amount calculating unit 26c calculates the phase shift, the sequence correcting unit 26d refers to the calculated phase shift and corrects the sequence execution data for main-scanning stored in the sequence-execution-data memory unit 23a.

0064Herein, the sequence correcting unit 26d corrects the pulse sequence for main-scanning in such a way that the phase shifts calculated by the correction amount calculating unit 26c become zero. At that time, for example, regarding the one-dimensional phase differences, the sequence correcting unit 26d changes the pulse sequence for main-scanning in such a way that a correction gradient magnetic field is applied between the flip pulse and the flop pulses. Besides, the sequence correcting unit 26d can also add the correction gradient magnetic fields before and after the readout gradient magnetic fields in the pulse sequence for main-scanning. Moreover, for example, regarding the zero-dimensional phase differences, the sequence correcting unit 26d changes the phases of the flop pulses in such a way that the phase shifts become zero.

0065In this way, the sequence correcting unit 26d corrects the pulse sequence for main-scanning in such a way that the phase shifts observed during the pre-scan become zero. As a result, it becomes possible to obtain an image that is neither affected by the zero-dimensional eddy currents and the one-dimensional eddy currents of the readout gradient magnetic fields and the slice selection gradient magnetic fields, nor affected by the zero-dimensional eddy currents of the phase encode gradient magnetic fields.

0066<figref idref="DRAWINGS">FIGS. 5 and 6</figref> are explanatory diagrams for explaining the phase shift correction operation performed by the sequence correcting unit 26d. <figref idref="DRAWINGS">FIG. 5</figref> illustrates an example of phase shifts occurring in the phase encode direction. In <figref idref="DRAWINGS">FIG. 5</figref>, a solid line 61 represents the zero-dimensional phase shifts for the echo signals in the phase encode direction, while a dashed line 62 represents the intensity of a phase encode magnetic field. <figref idref="DRAWINGS">FIG. 6</figref> illustrates the k-space in which the echo signals illustrated in <figref idref="DRAWINGS">FIG. 5</figref> are arranged. In <figref idref="DRAWINGS">FIG. 6</figref>, the horizontal direction illustrates the phase encode direction.

0067Herein, the zero-dimensional phase shifts in the phase encode direction depend on the manner of arranging the echo signals in the k-space. In the example illustrated in <figref idref="DRAWINGS">FIG. 5</figref>, with respect to the k-space, the collected echo signals are sequentially arranged in the phase encode direction. The numbers illustrated in <figref idref="DRAWINGS">FIG. 6</figref> represent the order of the collected echo signals. When the pre-scan is executed using the pulse sequence illustrated in <figref idref="DRAWINGS">FIG. 4</figref>, the position corresponding to the echo signal Echo10 collected at the 10-th count (the position indicated by an arrow in <figref idref="DRAWINGS">FIG. 6</figref>) is corrected. That is, the zero-dimensional phase shift of the echo signal disposed in the vicinity of the k-space is corrected. Meanwhile, herein, although the explanation is given for the case of performing correction with the echo signal Echo10 collected at the 10-th count serving as the reference signal, any other signal can be used as the reference signal.

0068The main-scan executing unit 26e executes the main-scan using the pulse sequence for main-scanning corrected by the sequence correcting unit 26d. More particularly, once the sequence correcting unit 26d corrects the sequence execution data for main-scanning, the pre-scan executing unit 26e reads the post-correction sequence execution data from the sequence-execution-data memory unit 23a. Then, the main-scan executing unit 26e executes the main-scan by sending that sequence execution data to the sequence control unit 10 via the interface unit 21.

0069Explained below is a sequence of operations in the phase shift correction operation performed in the MRI apparatus 100 according to the first embodiment. <figref idref="DRAWINGS">FIG. 7</figref> is a flowchart for explaining a sequence of operations in the phase shift correction operation performed in the MRI apparatus 100 according to the first embodiment.

0070As illustrated in <figref idref="DRAWINGS">FIG. 7</figref>, in the MRI apparatus 100 according to the first embodiment, when the operator instructs the start of imaging (Yes at Step S101), the imaging condition setting unit 26a receives imaging conditions input by the operator via the input unit 24 (Step S102).

0071Subsequently, based on the imaging conditions input by the operator, the imaging condition setting unit 26a generates sequence execution data for main-scanning and sequence execution data for pre-scanning (Step S103).

0072For example, as the sequence execution data for main-scanning, the imaging condition setting unit 26a generates the sequence execution data defining the pulse sequence illustrated in <figref idref="DRAWINGS">FIG. 3</figref>. Similarly, as the sequence execution data for pre-scanning, the imaging condition setting unit 26a generates the sequence execution data defining the pulse sequence illustrated in <figref idref="DRAWINGS">FIG. 4</figref>.

0073Then, the pre-scan executing unit 26b executes a pre-scan based on the sequence execution data of pulse sequences generated by the imaging condition setting unit 26a (Step S104).

0074Subsequently, the correction amount calculating unit 26c calculates, as the correction amount, the amount of phase shifting by referring to the phase differences present in a plurality of echo signals collected during each pre-scan (Step S105). Then, the sequence correcting unit 26d corrects the sequence execution data for main-scanning based on the correction amount calculated by the correction amount calculating unit 26c (Step S106).

0075Subsequently, the main-scan executing unit 26e executes a main-scan based on the sequence execution data for main-scanning corrected by the sequence correcting unit 26d (Step S107). Then, the image reconstructing unit 22 performs image reconstruction from the MR signal data collected during the main-scan (Step S108).

0076As described above, in the first embodiment, the pre-scan executing unit 26b executes pre-scans. Herein, during a pre-scan, readout gradient magnetic fields and slice direction gradient magnetic fields are applied in an identical manner to the pulse sequence for main-scanning, and phase encode gradient magnetic fields are applied in an identical manner to the pulse sequence for main-scanning up to just before the echoes used in calculating the correction amount. For example, during a pre-scan, in an identical manner to the pulse sequence for main-scanning, the readout gradient magnetic fields and the slice direction gradient magnetic fields are applied; and, of the phase encode magnetic gradient fields applied in the pulse sequence for main-scanning, at least a single phase encode magnetic gradient field is applied up to the echo in the vicinity of the k-space center. Then, from the phase differences present in a plurality of echo signals collected during the pre-scan, the correction amount calculating unit 26c calculates, as the correction amount, the amount of phase shifting (zero-dimensional and one-dimensional) occurring due to the readout gradient magnetic fields and the slice selection gradient magnetic fields, and the amount of zero-dimensional phase shifting occurring due to the phase encode gradient magnetic field. Based on the correction amount calculated by the correction amount calculating unit 26c, the sequence correcting unit 26d corrects the pulse sequence for main-scanning. Thus, according to the first embodiment, it becomes possible to prevent image deterioration that may occur due to the zero-dimensional eddy currents and the one-dimensional eddy currents of the readout gradient magnetic fields and the slice selection gradient magnetic fields, as well as due to the zero-dimensional eddy currents of the phase encode gradient magnetic fields.

0077Moreover, in the first embodiment, in the pulse sequence for pre-scanning, phase encode gradient magnetic fields are applied in an identical manner to the pulse sequence for main-scanning either for a period identical to the period during the main-scan in which the phase encode gradient magnetic fields become zero or up to just before the echo signal that is collected in the vicinity of the period during the main-scan in which the phase encode gradient magnetic fields become zero. Thus, according to the first embodiment, the echo signal that is arranged in the vicinity of the k-space center and that greatly contributes to the image quality is subjected to phase shift correction. That enables achieving further enhancement in the image quality of the images generated by the MRI apparatus 100.

0078Meanwhile, as a modification example of the first embodiment, the pre-scan executing unit 26b can be configured to perform the pre-scan for a plurality of number of times. For example, the pre-scan executing unit 26b can be configured to perform a first pre-scan during which readout gradient magnetic fields and slice direction gradient magnetic fields are applied in an identical manner to the pulse sequence for main-scanning, and phase encode gradient magnetic fields are applied in an identical manner to the pulse sequence for main-scanning up to just before the echoes used in calculating the correction amount; and configured to perform a second pre-scan during which a second pulse sequence is used in which phase encode gradient magnetic fields are applied that are different than the phase encode gradient magnetic fields applied during the first pre-scan.

0079For example, as the first pre-scan, the pre-scan executing unit 26b applies the phase encode gradient magnetic fields of that shot, from among a plurality of shots, which applies the greatest phase encode gradient magnetic field in the leading echo; and, as the second pre-scan, the pre-scan executing unit 26b applies the phase encode gradient magnetic fields of that shot which applies the smallest phase encode gradient magnetic field in the leading echo.

0080<figref idref="DRAWINGS">FIG. 8</figref> is an explanatory diagram for explaining the modification example of the first embodiment. In <figref idref="DRAWINGS">FIG. 8</figref>, of the phase shifts in the phase encode direction illustrated in <figref idref="DRAWINGS">FIG. 5</figref>, the phase shifts from the first echo signal (Echo1) up to the fourth echo signal (Echo4) are illustrated. In an identical manner to <figref idref="DRAWINGS">FIG. 5</figref>, the solid line 61 represents the zero-dimensional phase shifts for the echo signals in the phase encode direction, while the dashed line 62 represents the intensity of a phase encode magnetic field.

0081For example, in the example illustrated in <figref idref="DRAWINGS">FIG. 8</figref>, it is illustrated that a plurality of shots a, b, c, . . . , k, l, and m are used. Moreover, for example, as illustrated in <figref idref="DRAWINGS">FIG. 8</figref>, of the shots a, b, c, . . . , k, l, and m; assume that, in the first echo signal (Echo1), the shot applying the greatest phase encode gradient magnetic field is the shot a and the shot applying the smallest phase encode gradient magnetic field is the shot m. In that case, as the first pre-scan, the pre-scan executing unit 26b applies the phase encode of the shot “a” (see white triangles illustrated in <figref idref="DRAWINGS">FIG. 8</figref>); and, as the second pre-scan, the pre-scan executing unit 26b applies the phase encode of the shot “m” (see black triangles illustrated in <figref idref="DRAWINGS">FIG. 8</figref>).

0082In this case, for example, assuming that the phase difference during the first pre-scan is the phase difference p1 and assuming that the phase difference during the second pre-scan is the phase difference p2; the correction amount calculating unit 26c can set the phase difference p1 as the correction amount with respect to the shot “a”, set the phase difference p2 as the correction amount with respect to the shot “m”, and obtain the correction amount of the shots in between by interpolation. More particularly, if a number of shots is considered to be N, the correction amount for a shot having the greatest phase encode gradient magnetic field at the i-th count is p1+(p2−p1)*(i−1)/(N−1).

0083<figref idref="DRAWINGS">FIG. 9A</figref> is a diagram illustrating an exemplary image captured during the conventional pre-scan. <figref idref="DRAWINGS">FIG. 9B</figref> is a diagram illustrating an exemplary image captured during the pre-scan according to the first embodiment. <figref idref="DRAWINGS">FIG. 9A</figref> illustrates an exemplary image collected during the conventional pre-scan by skipping the phase encode. <figref idref="DRAWINGS">FIG. 9B</figref> illustrates an exemplary image collected during the pre-scan according to the first embodiment. As illustrated in <figref idref="DRAWINGS">FIGS. 9A and 9B</figref>, it can be observed that signal deterioration, which occurs due to the zero-dimensional phase shifts caused by the phase encode gradient magnetic fields, has improved in the first embodiment.

Second Embodiment

0084In the first embodiment, the explanation is given for correcting the phase shift of the echo signal disposed in the vicinity of the k-space center. However, that is not the only possible embodiment of the MRI apparatus 100. Explained below is a second embodiment in which the phase shifts of a plurality of echo signals are corrected. According to the second embodiment, the configuration of the MRI apparatus is essentially identical to the configuration illustrated in <figref idref="DRAWINGS">FIGS. 1 and 2</figref>. However, the operations performed by the imaging condition setting unit 26a, the pre-scan executing unit 26b, the correction amount calculating unit 26c, and the sequence correcting unit 26d are different as compared to the first embodiment.

0085In the second embodiment, the imaging condition setting unit 26a generates sequence execution data regarding not only the pulse sequence for main-scanning and a first pulse sequence for the first pre-scan (i.e., the pulse sequence illustrated in <figref idref="DRAWINGS">FIG. 4</figref>), but also regarding a second pulse sequence for the second pre-scan described below.

0086<figref idref="DRAWINGS">FIG. 10</figref> is a diagram illustrating a second pulse sequence for the second pre-scan according to the second embodiment. In the second pulse sequence, as compared to the phase encode gradient magnetic fields applied during the first pre-scan, the phase encode gradient magnetic fields corresponding to a different number of echoes are applied. For example, as illustrated in <figref idref="DRAWINGS">FIG. 10</figref>, in the second pulse sequence, the phase encode gradient magnetic fields identical to the pulse sequence for main-scanning are applied for a smaller number of echoes as compared to the first pre-scan. Meanwhile, the readout gradient magnetic fields and the slice selection readout gradient magnetic fields are identical to the first pre-scan.

0087For example, as illustrated in <figref idref="DRAWINGS">FIG. 10</figref>, in the second pulse sequence, the setting is such that the phase encode gradient magnetic fields up to the third echo signal are applied.

0088In the second embodiment, the pre-scan executing unit 26b executes not only the first pre-scan, but also the second pre-scan using the second pulse sequence illustrated in <figref idref="DRAWINGS">FIG. 10</figref>.

0089Moreover, in the second embodiment, the correction amount calculating unit 26c calculates a first-type phase shift from a plurality of echo signals collected during the first pre-scan as well as calculates a second-type phase shift from the phase differences present in a plurality of echo signals collected during the second pre-scan. Then, from the first-type phase shift and the second-type phase shift, the correction amount calculating unit 26c calculates the correction amount related to a plurality of echo signals.

0090For example, assume that the first pulse sequence illustrated in <figref idref="DRAWINGS">FIG. 4</figref> and the second pulse sequence illustrated in <figref idref="DRAWINGS">FIG. 10</figref> are executed. In that case, the correction amount calculating unit 26c calculates the phase difference p1 between the 10-th echo signal Echo10 and the 11-th echo signal Echo11 illustrated in <figref idref="DRAWINGS">FIG. 4</figref> as well as calculates the phase difference p2 between the fifth echo signal Echo5 and the fourth echo signal Echo4 illustrated in <figref idref="DRAWINGS">FIG. 10</figref>.

0091Moreover, the correction amount calculating unit 26c calculates p1 as the first-type phase shift and calculates p2 as the second-type phase shift. Herein, the first-type phase shift p1 becomes a phase shift dif10 in the 10-th echo signal Echo10. Similarly, the second-type phase shift p2 becomes a phase shift dif4 in the fourth echo signal Echo4.

0092Then, from the first-type phase shift and the second-type phase shift, the correction amount calculating unit 26c calculates the correction amount for the echo signals other than the fourth echo signal and the 10-th echo signal. For example, during the main-scan, when the phase encode gradient magnetic field becomes zero in the echo signal Echo10 collected at the 10-th count, the phase shift in the echo signals gradually increases up to the echo signal Echo10 and gradually decreases subsequent to the echo signal Echo10. With that, it can be assumed that, for example, a phase shift dif16 in the 16-th echo signal is equivalent to the phase shift dif4 in the fourth echo signal Echo4.

0093There, assuming that the echo signal at which the phase encode gradient magnetic field becomes zero during the main-scan serves as the reference signal, the correction amount calculating unit 26c refers to the phase shifts in the echo signals collected prior to the reference echo signal and estimates the phase shifts in the echo signals collected subsequent to the reference echo signal. Moreover, the correction amount calculating unit 26c calculates the phase shift in a plurality of echo signals by performing linear interpolation on the phase shift in each echo signal that has been calculated.

0094<figref idref="DRAWINGS">FIG. 11</figref> is an explanatory diagram for explaining the calculation of phase shifts performed by the sequence correcting unit 26d according to the second embodiment. As illustrated in <figref idref="DRAWINGS">FIG. 11</figref>, for example, with the 10-th echo signal Echo10 as the reference signal, the sequence correcting unit 26d refers to the phase shift dif4 in the fourth echo signal Echo4 and estimates the phase shift dif16 in the 16-th echo signal Echo16. Subsequently, the correction amount calculating unit 26c calculates the phase shift of a plurality of echo signals by performing linear interpolation on the phase shift dif4 in the echo signal Echo4, the phase shift dif10 in the echo signal Echo10, and the phase shift dif16 in the echo signal Echo16.

0095Moreover, in the second embodiment, based on the phase shift calculated by the sequence correcting unit 26d, the sequence correcting unit 26d corrects the intensity of the gradient magnetic field in the phase encode direction applied in each echo signal in the pulse sequence for main-scanning. At that time, for example, the sequence correcting unit 26d changes the intensity of the phase encode gradient magnetic field applied in each echo signal in such a way that the phase shift in each echo signal becomes zero. Alternatively, the sequence correcting unit 26d can also change the intensity of a rewind gradient magnetic field so that the phase shift in each echo signal becomes zero.

0096As described above, in the second embodiment, the pre-scan executing unit 26b also performs the second pre-scan. During the second pre-scan, as compared to the phase encode gradient magnetic fields applied during the first pre-scan, the phase encode gradient magnetic fields corresponding to a different number of echoes are applied. For example, during the second pre-scan, the phase encode gradient magnetic fields identical to the pulse sequence for main-scanning are applied for a smaller number of echoes as compared to the first pre-scan. Moreover, the correction amount calculating unit 26c calculates the second-type phase shift from the phase differences present in a plurality of echo signals collected during the first pre-scan and from the phase differences present in a plurality of echo signals collected during the second pre-scan. Then, from the first-type phase shift and the second-type phase shift, the correction amount calculating unit 26c calculates the correction amount related to a plurality of echo signals. Subsequently, based on the phase shift calculated for each echo signal by the correction amount calculating unit 26c, the sequence correcting unit 26d corrects the pulse sequence for main-scanning. Thus, according to the second embodiment, the phase shift of each echo signal collected during the main-scan is corrected. As a result, image deterioration can be prevented with accuracy.

0097Moreover, in the second embodiment, assuming that the echo signal at which the phase encode gradient magnetic field becomes zero during the main-scan serves as the reference signal, the correction amount calculating unit 26c refers to the phase shifts in the echo signals collected prior to the reference echo signal and estimates the phase shifts in the echo signals collected subsequent to the reference echo signal. In contrast, it is also possible to execute the pre-scan for one more time and estimate the phase shifts in the echo signals collected subsequent to the echo signal at which the phase encode gradient magnetic field becomes zero during the main-scan. That makes it possible to more accurately prevent image deterioration that may occur due to phase shifts in the phase encode direction.

0098In the first and second embodiments, the explanation is given for a case when, in a single collection, the echo signals are arranged in the entire area of the k-space. In contrast, there are times when the k-space is divided into a plurality of areas in the phase encode direction and, in a plurality of collections, a plurality of echo signals is collected in groups. In such a case, the phase shift can be calculated on a group-by-group basis and the pulse sequence for main-scanning can be corrected on a group-by-group basis.

0099<figref idref="DRAWINGS">FIG. 12</figref> is an explanatory diagram for explaining a modification example of the first and second embodiments. <figref idref="DRAWINGS">FIG. 12</figref> illustrates an exemplary case when the k-space along the phase encode direction is divided into three areas. A vertical arrow illustrated in <figref idref="DRAWINGS">FIG. 12</figref> indicates the k-space center. Herein, for example, the echo signals are collected in two groups Gr1 and Gr2. In that case, for example, as illustrated in <figref idref="DRAWINGS">FIG. 12</figref>, with respect to the central area, the echo signals in the group Gr1 are arranged in order in the phase encode direction. Regarding the echo signals in the group Gr2, the echo signals collected in the first half are arranged in order in one of the two side areas of the k-space and the echo signals collected in the latter half are arranged in order in the other of the two side areas of the k-space.

0100In such a case, for example, the correction amount calculating unit 26c can be configured to calculate the correction amount on a group-by-group basis. Moreover, in this case, the sequence correcting unit 26d corrects the pulse sequence for main-scanning on a group-by-group basis. Alternatively, for example, the correction amount calculating unit 26c can calculate the correction amount regarding one of the groups and the sequence correcting unit 26d can correct the pulse sequence for main-scanning for all groups based on that particular correction amount.

Third Embodiment

0101Explained below is an example of collection and correction of echo signals according to a third embodiment when the k-space center is positioned at the leading echo. In the third embodiment, as the pre-scan, the pre-scan executing unit 26b applies at least a single phase encode gradient magnetic field. <figref idref="DRAWINGS">FIG. 13</figref> is a diagram illustrating the pulse sequence for pre-scanning according to the third embodiment. As illustrated in <figref idref="DRAWINGS">FIG. 13</figref>, for example, the pre-scan executing unit 26b applies the phase encode magnetic fields up to the second echo signal, and makes use of the echo signals from the third echo signal onward for calculating the correction amount.

0102<figref idref="DRAWINGS">FIG. 14</figref> is an explanatory diagram for explaining collection of echo signals according to the third embodiment. As illustrated in <figref idref="DRAWINGS">FIG. 14</figref>, usually, when the k-space center is positioned at the leading echo, the echo signals are collected by dividing the k-space into two groups. In such a case, as the first pre-scan, the pre-scan executing unit 26b applies the representative phase encode gradient magnetic field of the group Gr1; and, as the second pre-scan, the pre-scan executing unit 26b applies the representative phase encode gradient magnetic field of the group Gr2. Then, in an identical manner to the first and second embodiment, the correction amount calculating unit 26c calculates the correction amount from the phase differences of a plurality of echoes in each pre-scan.

0103Meanwhile, in the pulse sequence used in the pre-scans, the gradient magnetic field such as a pre-pulse used in the main-scan is applied in an identical manner.

Fourth Embodiment

0104As a fourth embodiment, the explanation is given regarding a case in which the first to third embodiments are used in combination with a pre-scan that takes into account the phase shifts caused due to the slice encode of 3D collection sequence. <figref idref="DRAWINGS">FIG. 15</figref> is a diagram illustrating the second pulse sequence for the second pre-scan according to the fourth embodiment. In the fourth embodiment, the pre-scan executing unit 26b executes the first pre-scan using the first pulse sequence illustrated in <figref idref="DRAWINGS">FIG. 4</figref> as well as executes the second pre-scan using the second pulse sequence illustrated in <figref idref="DRAWINGS">FIG. 15</figref>.

0105Herein, for example, in the first pulse sequence, in an identical manner to the first embodiment, the representative phase encode gradient magnetic field of the main-scan is applied up to just short of the echoes in the vicinity of the k-space center. In the second pulse sequence, in addition to the phase encode gradient magnetic field applied during the first pre-scan, at least one of the slice encode phase gradient magnetic fields for the main-scan is applied. For example, of the slice encode gradient magnetic fields applied in the pulse sequence for main-scanning, the representative slice encode gradient magnetic field is applied in the second pulse sequence.

0106Then, according to the fourth embodiment, from the phase differences present in a plurality of echo signals collected during the first pre-scan and the phase differences present in a plurality of echo signals collected during the second pre-scan, the correction amount calculating unit 26c calculates, as the correction amount, the phase shift occurring in the slice encode direction due to the slice encode gradient magnetic field. Moreover, based on the correction amount calculated by the correction amount calculating unit 26c, the sequence correcting unit 26d corrects the pulse sequence for main-scanning. Hence, according to the fourth embodiment, it becomes possible to prevent image deterioration that may occur due to the phase shifts in the slice encode direction.

0107For example, when the slice encode count (slice count) is 64, the slice encode steps are −32 to 31. For example, in the second pulse sequence, assume that the slice encode corresponding to −32 is applied, assume that the phase difference (zero-dimensional or one-dimensional) between Echo10 and Echo11 during the first pre-scan is s1, and assume that the phase difference (zero-dimensional or one-dimensional) between Echo10 and Echo11 during the second pre-scan is s2. Then, the phase difference due to the i-th slice encode can be obtained by (s2−s1)*(−i)/32. By correcting the phase of the 180-degree pulse, it becomes possible to correct the phase shifts occurring due to slice encodes. Moreover, by increase the number of times of executing the pre-scan in the abovementioned case (for example, by applying the slice encode corresponding to +31), it becomes possible to enhance the accuracy.

0108As described above, according to the first to four embodiments, it becomes possible to prevent image deterioration that may be caused by zero-dimensional phase shifts of the phase encode gradient magnetic fields.

Fifth Embodiment

0109Given below is the explanation of a fifth embodiment. According to the first to fourth embodiments, it is possible to prevent image deterioration that may occur due to the zero-dimensional phase shifts of the phase encode gradient magnetic fields. However, by applying the phase encode gradient magnetic fields during the pre-scan, the phase correction amount cannot be correctly obtained due to the dependency on the RF coil or the field of view used in imaging. That may sometime cause degradation in the image quality.

0110<figref idref="DRAWINGS">FIG. 16</figref> is an explanatory diagram for explaining the MRI apparatus according to the fifth embodiment. <figref idref="DRAWINGS">FIG. 16</figref> illustrates an example when a sagittal cross-section of the lumber spine is captured with the Z-axis direction serving as the phase encode direction. In the example illustrated in <figref idref="DRAWINGS">FIG. 16</figref>, a phased array coil (PAC) is used as the RF coil and four coil elements included in the PAC coil are arranged in the Z-axis direction. In this example, as described in the abovementioned embodiments, when phase encode gradient magnetic fields are applied, one-dimensional phase shifts occur in the Z-axis direction. However, in the abovementioned embodiments, since sampling in the phase encoding direction is not performed, it is not possible to calculate and correct the one-dimensional phase shifts in the phase encode direction.

0111Herein, if the phase correction amount is obtained for each RF coil by implementing the conventional method; then, as illustrated in <figref idref="DRAWINGS">FIG. 16</figref>, at channels ch1 to ch4 connected to the four coil elements; one-dimensional phase shifts p1 to p4, respectively, are obtained. In that case, for example, the zero-dimensional phase shift amount to be corrected can be obtained as the average of the phase shifts p1 to p4 or as the weighted average of the signal intensities of the coil elements. However, for example, in the case when the phase shift p1 is −170° and the phase shift p3 is either 178°, which is in the vicinity of 180°, or −182°; it is not possible to properly perform the Unwrap operation. Moreover, depending on the type of the RF coil, the phases cannot be calculated in a correct manner. Hence, the method of obtaining the phase for each coil element may not necessarily be implemented in a proper manner.

0112Meanwhile, instead of obtaining the phase for each coil element, it is also possible to make use of data of the coil element that detects the greatest signal. However, in this method, if the channel ch1 or the channel ch4 detects the greatest signal, only the image quality in the vicinity of that coil element is enhanced. As the distance from that coil element increases, the image quality suffers from extreme degradation.

0113In order to resolve such an issue, in the fifth embodiment, the MRI apparatus 100 adds the data of all coil elements and calculates a single phase correction amount. More particularly, in the manner described below, the correction amount calculating unit 26c adds a plurality of echo signals collected via a plurality of coil elements during a pre-scan, and calculates a single phase correction amount to be used in correcting the phase differences present in a plurality of echo signals collected during the main-scan.

0114For example, regarding the calculation of phase correction amounts, the one-dimensional (inclined) phase angle (rad) is referred to as θ1 and the zero-dimensional (offset) phase angle (rad) of the i-th slice is referred to as θ0,i. Moreover, the readout point count is referred to as N, the echo count is referred to as M, the slice count is referred to as NS (0≦i<NS), the PAC channel count is referred to as L (0≦l<L), the pre-scan data is referred to as P(x, l) (however, complex data after 1D Fourier conversion: (0≦x<N)), and the score for obtaining the first-dimensional phase (e.g., four points) is referred to as dxx. Then, θ1 and θ0,i can be obtained by Equations (1) and (2) given below.

0000<maths id="MATH-US-00001" num="00001"><math overflow="scroll"><mtable><mtr><mtd><mrow><msub><mi>θ</mi><mn>1</mn></msub><mo>=</mo><mrow><mi>arg</mi><mo></mo><mrow><mo>(</mo><mrow><munderover><mo>∑</mo><mrow><mi>l</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>l</mi><mo>=</mo><mrow><mi>L</mi><mo>-</mo><mn>1</mn></mrow></mrow></munderover><mo></mo><mstyle><mspace width="0.3em" height="0.3ex" /></mstyle><mo></mo><mrow><munderover><mo>∑</mo><mrow><mi>i</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>NS</mi><mo>-</mo><mn>1</mn></mrow></munderover><mo></mo><mstyle><mspace width="0.3em" height="0.3ex" /></mstyle><mo></mo><mrow><munderover><mo>∑</mo><mrow><mi>x</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>N</mi><mo>-</mo><mn>1</mn></mrow></munderover><mo></mo><mstyle><mspace width="0.3em" height="0.3ex" /></mstyle><mo></mo><mrow><mo>(</mo><mrow><mrow><mi>P</mi><mo></mo><mrow><mo>(</mo><mrow><mi>x</mi><mo>,</mo><mi>i</mi><mo>,</mo><mi>l</mi></mrow><mo>)</mo></mrow></mrow><mo></mo><mrow><mover><msub><mi>P</mi><mi>k</mi></msub><mi>_</mi></mover><mo></mo><mrow><mo>(</mo><mrow><mrow><mi>x</mi><mo>-</mo><mrow><mo></mo><mi>xx</mi></mrow></mrow><mo>,</mo><mi>i</mi><mo>,</mo><mi>l</mi></mrow><mo>)</mo></mrow></mrow></mrow><mo>)</mo></mrow></mrow></mrow></mrow><mo>)</mo></mrow></mrow></mrow></mtd><mtd><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>θ</mi><mrow><mn>0</mn><mo>,</mo><mi>i</mi></mrow></msub><mo>=</mo><mrow><mi>arg</mi><mo></mo><mrow><mo>(</mo><mrow><munderover><mo>∑</mo><mrow><mi>l</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>l</mi><mo>=</mo><mrow><mi>Lc</mi><mo>-</mo><mn>1</mn></mrow></mrow></munderover><mo></mo><mstyle><mspace width="0.3em" height="0.3ex" /></mstyle><mo></mo><mrow><munderover><mo>∑</mo><mrow><mi>x</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>N</mi><mo>-</mo><mn>1</mn></mrow></munderover><mo></mo><mstyle><mspace width="0.3em" height="0.3ex" /></mstyle><mo></mo><mrow><mo>(</mo><mrow><mrow><mi>P</mi><mo></mo><mrow><mo>(</mo><mrow><mi>x</mi><mo>,</mo><mi>i</mi><mo>,</mo><mi>l</mi></mrow><mo>)</mo></mrow></mrow><mo></mo><msup><mi></mi><mrow><mo>+</mo><mstyle><mspace width="0.3em" height="0.3ex" /></mstyle><mo></mo><mrow><mi>i</mi><mo></mo><mrow><mo>(</mo><mrow><msub><mi>θ</mi><mn>1</mn></msub><mo></mo><mrow><mo>(</mo><mrow><mi>x</mi><mo>-</mo><mfrac><mi>N</mi><mn>2</mn></mfrac></mrow><mo>)</mo></mrow></mrow><mo>)</mo></mrow></mrow></mrow></msup></mrow><mo>)</mo></mrow></mrow></mrow><mo>)</mo></mrow></mrow></mrow></mtd><mtd><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mtd></mtr></mtable></math></maths>

0115Using Equations (1) and (2) given above, the correction amount calculating unit 26c calculates θ1 and θ0,i as phase correction amounts. According to this method, it becomes possible to obtain such an average phase correction amount within the field of view which is neither affected by the difference in the signal-to-noise ratio (SNR) in each coil element nor affected by the phase differences. Besides, this method is also effective against the conventional pre-scan in which the phase encode gradient magnetic fields are not applied.

0116Thus, according to the fifth embodiment, it becomes possible to obtain the phase correction amount in a stable manner regardless of the coil elements. As a result, degradation in the image quality can be prevented.

0117Besides, although the example in which signals collected through all of the coil elements are used is explained above, it is also possible to exclude signals collected through channels on which SNR of signals are bad and/or channels on which correct phases of signals are not obtained.

Sixth Embodiment

0118Explained below is a sixth embodiment. In the fifth embodiment, during phase calculation, all of the data in the readout direction is used. However, for example, the data to be used in phase calculation can be narrowed down to the data within the field of view (FOV) specified by the user.

0119<figref idref="DRAWINGS">FIG. 17</figref> is a diagram for explaining the MRI apparatus according to the sixth embodiment. <figref idref="DRAWINGS">FIG. 17</figref> illustrates a case of capturing images of the axial cross-section of a wrist with the X-axis direction serving as the phase encode direction. In <figref idref="DRAWINGS">FIG. 17</figref>, the range indicated by a solid line represents the field of view specified by the operator and the range indicated by a dotted line represents the field of view with NoWrap (described below) taken into account.

0120Herein, NoWrap points to a function in which, a field of view wider than the field of view specified by the operator at the time of imaging planning is set as an expanded field of view and an image in the expanded field of view is generated by means of unfolding. Then, the image in the field of view specified by the operator is clipped from the expanded field of view and displayed so as to reduce the wraparound artifact.

0121As illustrated in <figref idref="DRAWINGS">FIG. 17</figref>, for example, at a position spaced apart from the center in the X-axis direction, there are times when a nonlinear phase shift occurs (see the lower diagram in <figref idref="DRAWINGS">FIG. 17</figref>). In such a case, if all of the data in the readout direction is used, then the data regarding the outside of the field of view specified by the operator also makes a contribution during phase calculation. That makes it difficult to obtain an accurate phase correction amount.

0122In order to resolve such an issue, in the sixth embodiment, the MRI apparatus 100 limits the target area for calculation in the readout direction to the area within the field of view specified by the operator. More particularly, in the manner described below, of a plurality of echo signals collected during the pre-scan, the correction amount calculation unit 26c makes use of the data within the field of view specified in the readout direction by the operator, and calculates a phase correction amount to be used in correcting the phase differences present in a plurality of echo signals collected during the main-scan.

0123For example, the field of view specified by the operator is referred to as Fu, the field of view with NoWrap taken into account is referred to as Fn, the number of calculation start points of readout is referred to as Nstart, and the number of calculation end points of readout is referred to as Nend. Then, Nstart and Nend can be obtained by Equations (3) and (4) given below.

0000 Nstart=N*(1−Nu/Nf)/2 (3)

0000 Nend=N−Nstart (4)

0124Using Equations (3) and (4), the correction amount calculating unit 26c calculates Nstart and Nend−1; and performs the addition of x, mentioned in Equations (1) and (2), from Nstart up to Nend−1. According to this method, if the phases act differently on the inside and on the outside of the field of view specified by the operator, it becomes possible to eliminate the effect of the data regarding the outside of the field of view during the calculation of the phase correction amount. Besides, this method is also effective against the conventional pre-scan.

0125Thus, according to the sixth embodiment, it becomes possible to obtain the phase correction amount in a stable manner regardless of the field of view. As a result, degradation in the image quality can be prevented.

Seventh Embodiment

0126Given below is the explanation of a seventh embodiment. In the sixth embodiment, during phase calculation, the data in the readout direction is narrowed down to the data within the field of view specified by the operator. However, for example, when the Z-axis direction serves as the readout direction, magnetostatic field unevenness may also have an adverse effect. Thus, during phase calculation, the data in the readout direction can be narrowed down to the data within a range that is less affected by magnetostatic field unevenness.

0127<figref idref="DRAWINGS">FIG. 18</figref> is an explanatory diagram for explaining the MRI apparatus according to the seventh embodiment. <figref idref="DRAWINGS">FIG. 18</figref> illustrates a case of capturing images of the abdominal region with the Z-axis direction serving as the phase encode direction. At a position spaced apart from the center in the Z-axis direction, there are times when a nonlinear phase shift occurs (see the middle diagram in <figref idref="DRAWINGS">FIG. 18</figref>) and when the signal intensity increases due the effect of magnetostatic field unevenness. In such a case, if all of the data in the readout direction is used, then the data regarding the positions spaced apart from the center in the Z-direction strongly affects the phase calculation thereby making it difficult to obtain an accurate phase correction amount.

0128In order to resolve such an issue, in the seventh embodiment, the MRI apparatus 100 limits the target area for calculation in the readout direction to such a field of view in the Z-axis direction which is compensated at least by the system. For example, of a plurality of echo signals collected during a pre-scan, the correction amount calculating unit 26c makes use of the data of such a range which is compensated at least for magnetostatic field unevenness, and calculates a phase correction amount to be used in correcting the phase differences present in a plurality of echo signals collected during the main-scan. For example, the correction amount calculating unit 26c calculates the phase correction amount by referring to the data regarding half (i.e., the area Fc indicated in <figref idref="DRAWINGS">FIG. 18</figref>) of the area of the field of view in the Z-direction that is compensated by the system (i.e., the area enclosed by dotted lines). According to this method, it becomes possible to eliminate the effect of magnetostatic field unevenness during the calculation of the phase correction amount. Besides, this method is also effective against the conventional pre-scan.

0129Thus, according to the seventh embodiment, it becomes possible to obtain the phase correction amount in a stable manner regardless of magnetostatic field unevenness. As a result, degradation in the image quality can be prevented.

0130Besides, although the example in which a field of view in the Z-axis direction is limited, it is also possible to limit a field of view in X-axis direction and/or a field of view in Y-axis direction.

0131As described above, according to the first to seventh embodiments, it becomes possible to prevent image deterioration that may occur due to phase shifts of echo signals.

0132While certain embodiments have been described, these embodiments have been presented by way of example only, and are not intended to limit the scope of the inventions. Indeed, the novel embodiments described herein may be embodied in a variety of other forms; furthermore, various omissions, substitutions and changes in the form of the embodiments described herein may be made without departing from the spirit of the inventions. The accompanying claims and their equivalents are intended to cover such forms or modifications as would fall within the scope and spirit of the inventions.

Contents5

Every citation, both ways

Citations
Document	Relation	Office	Cited during
US9041394B2	Cited by	United States of America	Search report
JP2017000373A	Cited by	Japan	Search report
US10495716B2	Cited by	United States of America	Applicant
US10353041B2	Cited by	United States of America	Applicant
US2016363643A1	Cited by	United States of America	Pre-grant
US9046590B2	Cited by	United States of America	Search report
US2011291653A1	Cited by	United States of America	Pre-grant
US10288713B2	Cited by	United States of America	Applicant
WO2017142180A1	Cited by	World Intellectual Property Organization (WIPO)	International search
US10241183B2	Cited by	United States of America	Search report
US2016041248A1	Cited by	United States of America	Pre-grant
US10135549B2	Cited by	United States of America	Search report
WO2018097644A1	Cited by	World Intellectual Property Organization (WIPO)	International search
US2014307764A1	Cited by	United States of America	Pre-grant
US10274569B2	Cited by	United States of America	Search report
US2001017544A1	Cites	United States of America	Pre-grant
US2003102864A1	Cites	United States of America	Pre-grant
US2005218894A1	Cites	United States of America	Pre-grant
US2007238969A1	Cites	United States of America	Pre-grant
US2008265885A1	Cites	United States of America	Pre-grant
US2009134871A1	Cites	United States of America	Pre-grant
US2009177078A1	Cites	United States of America	Pre-grant
US2009212773A1	Cites	United States of America	Pre-grant
US2009300843A1	Cites	United States of America	Pre-grant
US2010164495A1	Cites	United States of America	Pre-grant
US2011105886A1	Cites	United States of America	Pre-grant
US2011112393A1	Cites	United States of America	Pre-grant
US2011245655A1	Cites	United States of America	Pre-grant
US2011291653A1	Cites	United States of America	Pre-grant
US2012313641A1	Cites	United States of America	Pre-grant
US2013057280A1	Cites	United States of America	Pre-grant
US2013169275A1	Cites	United States of America	Pre-grant
US2013293231A1	Cites	United States of America	Pre-grant
US2013314090A1	Cites	United States of America	Pre-grant
US2014084916A1	Cites	United States of America	Pre-grant
US2014111203A1	Cites	United States of America	Pre-grant
US4870361A	Cites	United States of America	Pre-grant
US4982161A	Cites	United States of America	Pre-grant
US5262725A	Cites	United States of America	Pre-grant
US5347218A	Cites	United States of America	Pre-grant
US5450010A	Cites	United States of America	Pre-grant
US5493224A	Cites	United States of America	Pre-grant
US5647362A	Cites	United States of America	Pre-grant
US5729139A	Cites	United States of America	Pre-grant
US5818229A	Cites	United States of America	Pre-grant
US5928146A	Cites	United States of America	Pre-grant
US6150815A	Cites	United States of America	Pre-grant
US6160397A	Cites	United States of America	Pre-grant
US6218834B1	Cites	United States of America	Pre-grant
US6285900B1	Cites	United States of America	Pre-grant
US6331777B1	Cites	United States of America	Pre-grant
US6369568B1	Cites	United States of America	Pre-grant
US6424152B1	Cites	United States of America	Pre-grant
US6470203B2	Cites	United States of America	Pre-grant
US6617850B2	Cites	United States of America	Pre-grant
US6781375B2	Cites	United States of America	Pre-grant
US7141970B2	Cites	United States of America	Pre-grant
US7622922B2	Cites	United States of America	Pre-grant
US7894891B2	Cites	United States of America	Pre-grant
US7898254B2	Cites	United States of America	Pre-grant
US8115485B1	Cites	United States of America	Pre-grant
US8154294B2	Cites	United States of America	Pre-grant
US8159221B2	Cites	United States of America	Pre-grant
US8473046B2	Cites	United States of America	Pre-grant
US8587305B2	Cites	United States of America	Pre-grant
US8587310B2	Cites	United States of America	Pre-grant
US8594766B2	Cites	United States of America	Pre-grant
US8624595B2	Cites	United States of America	Pre-grant

10 members in 3 offices

This recordMember, by publication datePriority claim

Priority claims10

Priority claims
Document	Office	Kind	Date
2011080111	Japan	–	2011-03-31
2011080111	Japan	A	2011-03-31
2011080111	Japan	A	2011-03-31
2012045281	Japan	–	2012-03-01
2012045281	Japan	A	2012-03-01
2012045281	Japan	A	2012-03-01
2011080111	–	–	–
2012045281	–	–	–
JP20110080111	–	–	–
JP20120045281	–	–	–

Members10

Family members, oldest first
Document	Office	Kind
US2012249143A1	United States of America	A1	2012-10-04
CN102727203A	China	A	2012-10-17
JP2012213608A	Japan	A	2012-11-08
US9046590B2	United States of America	B2	2015-06-02
CN102727203B	China	B	2015-07-29
JP5925529B2	Japan	B2	2016-05-25
JP2016104229A	Japan	A	2016-06-09
JP2016127952A	Japan	A	2016-07-14
JP6194032B2	Japan	B2	2017-09-06
JP6194033B2	Japan	B2	2017-09-06

54 transactions on the USPTO file

Allowed without a rejection on record.

Non-final rejections: 0
Final rejections: 0
RCEs: 0
Appeals: 0

Over time

Point at a mark for the transaction

Transactions

Prosecution transactions, newest first
	Event	Code
2015-06-02	Recordation of Patent Grant MailedPGM/	PGM/
2015-06-02	Patent Issue Date Used in PTA CalculationAllowedPTAC	PTAC
2015-05-15	Email NotificationEML_NTR	EML_NTR
2015-05-13	Issue Notification MailedAllowedWPIR	WPIR
2015-04-29	Dispatch to FDCD1935	D1935
2015-04-28	Email NotificationEML_NTR	EML_NTR
2015-04-28	Filing Receipt - CorrectedFLRCPT.C	FLRCPT.C
2015-04-24	Issue Fee Payment VerifiedN084	N084
2015-04-24	Application Is Considered Ready for IssuePILS	PILS
2015-04-24	Issue Fee Payment ReceivedIFEE	IFEE
2015-04-23	Email NotificationEML_NTR	EML_NTR
2015-04-23	Mail Response to 312 Amendment (PTO-271)MN271	MN271
2015-04-16	Response to Amendment under Rule 312N271	N271
2015-04-14	Printer Rush- No mailingTCPB	TCPB
2015-04-14	Printer Rush- No mailingTCPB	TCPB
2015-04-09	Pubs Case Remand to TCPUBTC	PUBTC
2015-04-06	Amendment after Notice of Allowance (Rule 312)AllowedA.NA	A.NA
2015-04-06	Workflow - Drawings FinishedDRWF	DRWF
2015-02-27	Electronic ReviewELC_RVW	ELC_RVW
2015-02-27	Email NotificationEML_NTF	EML_NTF
2015-02-27	Mail Notice of AllowanceAllowedMN/=.	MN/=.
2015-02-24	Notice of Allowance Data Verification CompletedAllowedN/=.	N/=.
2015-02-23	Reasons for AllowanceEX.R	EX.R
2015-02-23	Examiner's Amendment CommunicationEX.A	EX.A
2015-02-23	Interview Summary - Examiner Initiated - TelephonicEXET	EXET
2015-02-23	Interview Summary - Examiner InitiatedEXIE	EXIE
2014-11-17	Information Disclosure Statement consideredIDSC	IDSC
2014-11-17	Reference capture on IDSRCAP	RCAP
2014-11-17	Information Disclosure Statement (IDS) FiledM844	M844
2014-11-17	Information Disclosure Statement (IDS) FiledWIDS	WIDS
2014-09-16	Case Docketed to Examiner in GAUDOCK	DOCK
2014-07-01	Case Docketed to Examiner in GAUDOCK	DOCK
2013-10-03	Priority document has successfully retrieved via PDX/DASPD.RECVD	PD.RECVD
2012-11-04	Request for Foreign Priority (Priority Papers May Be Included)RQPR	RQPR
2012-10-05	Email NotificationEML_NTR	EML_NTR
2012-10-04	PG-Pub Issue NotificationPG-ISSUE	PG-ISSUE
2012-07-21	Case Docketed to Examiner in GAUDOCK	DOCK
2012-05-09	Case Docketed to Examiner in GAUDOCK	DOCK
2012-05-03	Application Dispatched from OIPEOIPE	OIPE
2012-04-17	Application Is Now CompleteCOMP	COMP
2012-04-17	Email NotificationEML_NTR	EML_NTR
2012-04-17	Filing ReceiptFLRCPT.O	FLRCPT.O
2012-04-13	Application Is Now CompleteCOMP	COMP
2012-04-13	Email NotificationEML_NTR	EML_NTR
2012-04-13	Filing ReceiptFLRCPT.O	FLRCPT.O
2012-04-12	Sent to Classification ContractorPGPC	PGPC
2012-04-03	Cleared by OIPE CSRL194	L194
2012-03-30	Information Disclosure Statement consideredIDSC	IDSC
2012-03-30	Reference capture on IDSRCAP	RCAP
2012-03-30	Information Disclosure Statement (IDS) FiledM844	M844
2012-03-30	Request from applicant for the USPTO to retrieve the Priority DocumentPDREQUST	PDREQUST
2012-03-30	Information Disclosure Statement (IDS) FiledWIDS	WIDS
2012-03-30	IFW Scan & PACR Auto Security ReviewSCAN	SCAN
2012-03-30	Initial Exam Team nnIEXX	IEXX

6 legal events, as the office reported them to INPADOC

Over the term

Point at a mark for the event

Events

Legal events, newest first
	Event	Code
2022-11-16	Maintenance fee paymentMAFP	MAFP
2018-11-15	Maintenance fee paymentMAFP	MAFP
2016-05-31	AssignmentAS	AS
2015-05-13	Information on status: patent grantGrantedPATENTED CASESTCF	STCF
2012-03-30	AssignmentAS	AS
2012-03-30	AssignmentAS	AS

Magnetic resonance imaging apparatus

Claim Score by NHIP

Abstract