WO2012063654A1 - 磁気共鳴イメージング装置、計測空間座標補正方法、及び、画像再構成方法 - Google Patents
磁気共鳴イメージング装置、計測空間座標補正方法、及び、画像再構成方法 Download PDFInfo
- Publication number
- WO2012063654A1 WO2012063654A1 PCT/JP2011/074892 JP2011074892W WO2012063654A1 WO 2012063654 A1 WO2012063654 A1 WO 2012063654A1 JP 2011074892 W JP2011074892 W JP 2011074892W WO 2012063654 A1 WO2012063654 A1 WO 2012063654A1
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- magnetic field
- gradient magnetic
- image
- magnetic resonance
- resonance imaging
- Prior art date
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
- G01R33/48—NMR imaging systems
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
- G01R33/48—NMR imaging systems
- G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
- G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
- G01R33/565—Correction of image distortions, e.g. due to magnetic field inhomogeneities
- G01R33/56572—Correction of image distortions, e.g. due to magnetic field inhomogeneities caused by a distortion of a gradient magnetic field, e.g. non-linearity of a gradient magnetic field
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/28—Details of apparatus provided for in groups G01R33/44 - G01R33/64
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
- G01R33/48—NMR imaging systems
- G01R33/4818—MR characterised by data acquisition along a specific k-space trajectory or by the temporal order of k-space coverage, e.g. centric or segmented coverage of k-space
- G01R33/4824—MR characterised by data acquisition along a specific k-space trajectory or by the temporal order of k-space coverage, e.g. centric or segmented coverage of k-space using a non-Cartesian trajectory
Definitions
- the present invention relates to a magnetic resonance imaging (hereinafter abbreviated as “MRI”) apparatus for obtaining a tomographic image of an examination site of a subject using a nuclear magnetic resonance (hereinafter abbreviated as “NMR”) phenomenon, and in particular, a system.
- MRI magnetic resonance imaging
- NMR nuclear magnetic resonance
- the present invention relates to a technique for reducing artifacts caused by a deviation from an ideal waveform of a gradient magnetic field pulse waveform resulting from the responsiveness of.
- the gradient magnetic field is used to give positional information to the echo signal from the subject placed in the static magnetic field.
- the coordinates of the echo signal arranged in the measurement space are determined by the gradient magnetic field history in each direction applied when generating the echo signal.
- the image of the subject is reconstructed by performing inverse Fourier transform on the data of the echo signal filling the k space.
- the coordinates of this k space are determined on the assumption that the gradient magnetic field pulse has a designed waveform, and the magnitude of the gradient magnetic field pulse is determined when calculating a pulse sequence for controlling imaging.
- the echo signal is arranged at a coordinate shifted from the coordinate that should be originally arranged in the k space, and there is a problem that the image quality deteriorates.
- This problem becomes more serious in an image obtained by an imaging sequence of a non-orthogonal sampling method that samples the k space in a radial or spiral manner, compared to an imaging sequence of the orthogonal sampling method that samples the k space in parallel with the axial direction.
- the gradient magnetic field output error is the amount of gradient magnetic field pulse set at the time of sequence calculation and the amount of gradient magnetic field pulse actually output (the amount of gradient magnetic field applied to the nuclear spin in the tissue constituting the subject) This error is caused by various factors such as static magnetic field inhomogeneity, gradient magnetic field offset, eddy current, time lag of gradient magnetic field output, rise (or fall) time lag, etc. Is included.
- Non-Patent Document 1 proposes applying the latter method to the non-orthogonal sampling method. Specifically, an equivalent circuit that approximates the system response circuit of the gradient magnetic field output is obtained, and the transfer function represented by this equivalent circuit (more precisely, a function obtained by inverse Laplace transform) is convolved with the gradient magnetic field output. By doing so, the gradient magnetic field output (gradient magnetic field waveform) considering the system response of the gradient magnetic field output is calculated.
- the actual output of the gradient magnetic field pulse is often corrected by a feedback circuit or the like in order to approximate the gradient magnetic field pulse waveform to an ideal shape.
- the output gradient magnetic field pulse waveform shape cannot be approximated by the system response circuit.
- an object of the present invention is to accurately approximate the gradient magnetic field pulse waveform shape and improve the image quality.
- the MRI apparatus of the present invention includes a gradient magnetic field generating means for providing a magnetic field gradient to the static magnetic field generated by the static magnetic field generating means, and a transmission for irradiating a test object arranged in the static magnetic field with a high-frequency magnetic field
- Signal processing in an MRI apparatus comprising: means; a receiving means for receiving an echo signal emitted from the examination object by nuclear magnetic resonance; and a signal processing means for processing the echo signal to reconstruct an image to be examined
- the means generates an approximate function that approximates the gradient magnetic field pulse waveform for each section divided in time into at least two gradient magnetic field pulse waveforms to which a magnetic field gradient is given by the gradient magnetic field generation means.
- correction means for correcting the measurement space coordinates on which the echo signal is arranged using an approximation function.
- artifacts such as signal loss and distortion can be reduced in reconstructed images by various pulse sequences.
- the block diagram which shows the whole structure of an example of the MRI apparatus with which this invention is applied
- the flowchart which shows the imaging operation by 1st embodiment Time chart showing the pulse sequence of the radial method It is a figure explaining the relationship between readout gradient magnetic field and sampling, (a) shows a general sampling method, (b) shows a nonlinear sampling method. It is a figure explaining the system response which appears in a read gradient magnetic field pulse waveform, (a) is an ideal waveform shape, (b), (c) shows the waveform shape deform
- the figure explaining the approximation of the gradient magnetic field pulse waveform shape (a) shows the case where the rising / falling is divided into two sections and (b) shows the case where the rising / falling is divided into three sections. . It is a figure explaining the influence of the system response by the difference in a pulse sequence, (a) is k-space sampling of an orthogonal system pulse sequence, (b), (c) shows k-space sampling of a radial method pulse sequence, respectively. It is a figure which shows the influence of the system response given to the image of a radial method, (a) does not perform k space correction, (b) shows the case where the correction by this invention is performed.
- FIGS. 4A and 4B are diagrams for explaining evaluation judgment of an optimum parameter in a parameter search step
- FIGS. 4A and 4B are diagrams showing a phantom image reconstructed using different approximate parameters and a profile of the A-A ′ line, respectively.
- FIG. 1 is a block diagram showing an overall configuration of an example of an MRI apparatus to which the present invention is applied.
- This MRI apparatus includes a static magnetic field generation system 2, a gradient magnetic field generation system 3, a transmission system 5, a reception system 6, a signal processing system 7, a sequencer 4, and a central processing unit (CPU) 8. Yes.
- the static magnetic field generation system 2 generates a uniform static magnetic field in the space where the subject 1 is placed, and comprises a magnetic field generation means of a permanent magnet system, a normal conduction system, or a superconductivity system.
- the shape and arrangement of the magnetic field generating means differ depending on the direction of a static magnetic field such as a vertical magnetic field or a horizontal magnetic field.
- the gradient magnetic field generating system 3 includes a gradient magnetic field coil 9 wound in three orthogonal directions of X, Y, and Z, and a gradient magnetic field power source 10 for driving each gradient magnetic field coil.
- a gradient magnetic field power source 10 for driving each gradient magnetic field coil.
- the gradient magnetic fields Gx, Gy, and Gz in the X, Y, and Z directions are applied to the subject 1.
- an imaging cross section (slice) of the subject 1 can be selected, or position information can be encoded in an echo signal.
- the sequencer 4 is a control means that repeatedly applies a high-frequency magnetic field pulse (hereinafter referred to as “RF pulse”) and a gradient magnetic field pulse in a predetermined pulse sequence, and operates under the control of the CPU 8 to collect tomographic image data of the subject 1.
- RF pulse high-frequency magnetic field pulse
- Various necessary commands are sent to the transmission system 5, the gradient magnetic field generation system 3, and the reception system 6.
- the pulse sequence includes a preliminary measurement sequence for measuring the characteristics of the device.
- Various pulse sequences differ depending on the imaging method and the purpose of measurement. Are stored in advance as a program.
- the transmission system 5 irradiates an RF pulse to cause nuclear magnetic resonance to the nuclear spins of atoms constituting the biological tissue of the subject 1, and includes a high frequency oscillator 11, a modulator 12, a high frequency amplifier 13, and a transmission side And a high-frequency coil 14a.
- the high-frequency pulse output from the high-frequency oscillator 11 is amplitude-modulated by the modulator 12 at a timing according to a command from the sequencer 4, and the amplitude-modulated high-frequency pulse is amplified by the high-frequency amplifier 13 and then placed close to the subject 1.
- the subject 1 is irradiated with electromagnetic waves (RF pulses).
- the receiving system 6 detects an echo signal (NMR signal) emitted by nuclear magnetic resonance of nuclear spins constituting the biological tissue of the subject 1, and includes a high-frequency coil 14b on the receiving side, an amplifier 15, and a quadrature phase It comprises a detector 16 and an A / D converter 17.
- the response electromagnetic wave (NMR signal) of the subject 1 induced by the electromagnetic wave irradiated from the high-frequency coil 14a on the transmission side is detected by the high-frequency coil 14b arranged close to the subject 1 and amplified by the amplifier 15 Thereafter, the signals are divided into two orthogonal signals by the quadrature phase detector 16 at a timing according to a command from the sequencer 4, converted into digital quantities by the A / D converter 17, and sent to the signal processing system 7.
- the signal processing system 7 has an external storage device (storage device) such as an optical disk 19 and a magnetic disk 18 and a display 20 made up of a CRT or the like.
- an external storage device storage device
- the CPU 8 Processing such as signal processing and image reconstruction is executed, and the resulting tomographic image of the subject 1 is displayed on the display 20 and recorded on the magnetic disk 18 of the external storage device.
- the memory of the CPU 8 and the external storage device are collectively referred to as a storage device.
- the storage device records parameters (approximation parameters) of an approximate function of a gradient magnetic field pulse waveform shape to be described later. This parameter includes a parameter value obtained as a calculation result of the CPU 8 in addition to a value set in advance as an adjustment value.
- the high-frequency coils 14a and 14b on the transmission side and the reception side and the gradient magnetic field coil 9 are installed in the static magnetic field space of the static magnetic field generation system 2 arranged in the space around the subject 1. .
- a parameter (approximation parameter) of a function that approximates a gradient magnetic field pulse is stored in a storage device, and a measurement space (k space) coordinate of a measured echo signal is calculated using the approximation parameter, and gridding is performed. It is a feature.
- this embodiment will be described with reference to the flow of the imaging operation shown in FIG.
- a gradient magnetic field pulse waveform is calculated (step 201).
- An echo signal is measured using the pulse sequence calculated in step 201 (step 202).
- the pulse sequence may be either an orthogonal sampling method that samples the k space in parallel in the axial direction or a non-orthogonal sampling method that samples in a radial or spiral manner.
- Fig. 3 shows the pulse sequence of the radial method, which is a kind of non-orthogonal sampling method.
- RF, Gs, G1, G2, A / D, and echo are RF pulse, slice gradient magnetic field, readout gradient magnetic field in the first direction, readout gradient magnetic field in the second direction, AD conversion, and echo signal, respectively.
- an RF pulse 301 is applied together with a slice selection gradient magnetic field pulse 302 to excite a desired slice, then a slice rephase gradient magnetic field pulse 303 is applied, and a read gradient magnetic field pulse 304 in the biaxial direction is applied.
- the echo signal 307 is measured in the sampling window 306 which is applying the readout gradient magnetic field pulse whose polarity is reversed. This measurement is repeated at a repetition time TR, and the output of the first and second readout gradient magnetic field pulses 304 and 305 is changed little by little to acquire data necessary for one image reconstruction.
- k-space coordinates are calculated using an approximate function of the gradient magnetic field pulse waveform stored in the storage device (memory or storage device) 200 (step 203).
- gridding processing is performed to create k-space data (step 204).
- the created k-space data is two-dimensionally Fourier transformed to create an image (step 205).
- the signal processing system 7 which is a signal processing unit has a gradient magnetic field for each section divided in time into at least two gradient magnetic field pulse waveforms to which a magnetic field gradient is given by the gradient magnetic field generation unit.
- An approximate function generating means for generating an approximate function approximating the pulse waveform; and a correcting means for correcting the measurement space coordinates where the echo signal is arranged using the approximate function.
- This method is characterized in that it includes a step of correcting k-space coordinates where measured data (echo signal) is arranged using an approximate function of a gradient magnetic field pulse waveform stored as an adjustment value. Details of step 203 will be described below.
- the output G (t) of the readout gradient magnetic field pulse and the k-space coordinate k (t) have the relationship of equation (1).
- ⁇ represents a gyromagnetic ratio
- FIG. 4 shows the relationship between the readout gradient magnetic field pulse (304, 305 in FIG. 3) and the sampling 306, as shown in FIG. 4 (a), when sampling is performed when the gradient magnetic field pulse waveform becomes constant.
- the interval of k-space coordinates is constant.
- FIG. 4 (b) when sampling from the rising edge to the falling edge of the gradient magnetic field pulse 401, at the rising edge and the falling edge of the gradient magnetic field pulse 401, the interval between the k-space coordinates is as shown in FIG. It is denser than the case, but is affected by the system response.
- the sampling method shown in FIG. 4 (b) is called Ramp sampling or non-linear sampling (Nonlinear sampling) and can acquire more sampling points, which has the effect of increasing the spatial resolution or shortening the gradient magnetic field pulse application time. can get.
- the pulse sequence calculated in step 201 is calculated on the assumption that the gradient magnetic field pulse has an ideal shape as shown in FIG. 5 (a). As shown in (c) and (c), there are delays in the output response and deformation of the rising and falling shapes.
- at least two approximation functions that approximate the gradient magnetic field waveform that is actually applied are set with rising or falling inflection points as boundaries, and k-space coordinates are set using parameters of these approximation functions. to correct.
- the gradient magnetic field waveform is divided into a plurality of sections, and an approximation function for approximating the waveform is defined for each section.
- FIG. 6 shows an example of divided gradient magnetic field waveform sections. In the example shown in the figure, there are five sections: a section 601-1 at the beginning of the gradient magnetic field pulse, a section 602-1 at the end of rising, a constant section 603, a section 602-2 at the beginning of falling, and a section 601-2 at the end of falling. Dividing into sections, the sections 601-1 and 602-1 are approximated by functions expressed by the equations (2) and (3).
- D start is a time interval of the section 601-1 and can be determined to be, for example, several tens of% of the rising / falling time of the gradient magnetic field pulse.
- a, b, and c are constants.
- D end is a time interval of the section 602-1 and can be determined to be, for example, several tens of% of the rising / falling time of the gradient magnetic field pulse, similarly to D start .
- d, e, and f are constants.
- Constants a, b, c, d, e, f can also be obtained from the echo signal measured in step 202, but here the value obtained from the gradient magnetic field pulse waveform measured in advance by preliminary measurement, It is assumed that it is stored in the storage device 200 as a parameter.
- Equation (2) and Equation (3) are applied in the reverse direction.
- another approximate function in which the constants a, b, c, d, e, and f in the expressions (2) and (3) are different may be defined. Thereby, the precision of approximation can be raised.
- the gradient magnetic field waveform is divided into five sections, and the rising and falling portions are approximated using two functions of G start and G end , respectively.
- the way of dividing is not limited to the example shown in FIG. 6 (a), for example, as shown in FIG. 6 (b), the rising and falling parts are divided into three sections, and the waveform changes nonlinearly.
- step 202 When the signal measurement (imaging) executed in step 202 is a radial method pulse sequence as shown in FIG. 3, readout gradient magnetic field pulses 304 and 305 in the first direction and the second direction are used. . In that case, it is possible to use an average value for the gradient magnetic field system of each axis as the approximate parameter. However, since the actual gradient coil differs in efficiency from axis to axis due to the difference in shape, for each axis, It is preferable to define an approximate function.
- step 203 in FIG. 2 the parameters of the above formulas (2) and (3) stored in the storage device 200 are read, and the coordinates of the k space are calculated. Specifically, integration of each section is executed using a function approximating the gradient magnetic field pulse waveform for each section, as shown in Expression (4). This determines the k-space coordinates of the data sampled at time t.
- G Rise_start the same function as the G start of the formula (2)
- G rise_end is the same function as the G end
- G Fall_start is here, each function G start, G end
- G flat is a constant value.
- echo signal data is collected at regular sampling intervals.
- t ′ the sampling time of each data
- k-space coordinates where each data is arranged are obtained.
- the coordinates according to equation (4) are calculated for each axis (eg, kx, ky). A point on the k-space determined by the coordinates of these two axes becomes the coordinate of the data.
- the coordinates of k-space are represented by integers, but the k-space coordinates calculated by equation (4) are not necessarily integers.
- the data is converted into data corresponding to coordinates represented by integers using an interpolation process called gridding.
- the created k-space data is two-dimensionally Fourier transformed to create an image (step 205).
- the gradient magnetic field pulse waveform is divided into a plurality of sections, an approximate function is defined for each section, and k-space coordinates are calculated using each approximate function, so that the system response of the gradient magnetic field is uniform. Even when it cannot be approximated by two equivalent circuits, the k-space coordinates are corrected using a waveform close to the actual gradient magnetic field output, and an image in which the influence of the error due to the system response is reduced can be obtained.
- This embodiment is suitable for an imaging method that collects data at the rise and fall of the readout gradient magnetic field, as shown in FIG. 4 (b).
- the pulse sequence can be applied to both orthogonal sampling and non-orthogonal sampling imaging, but is particularly suitable for non-orthogonal sampling methods where the effect of system response is likely to appear in reconstructed images. It is.
- FIG. 7 (a) is a diagram showing the k-space data collection method using the orthogonal sampling method, (b) and (c) are each the k-space data collection method using the radial method, and (b) is centered from the outside of the k-space.
- a general radial method of sampling to the outside through (c) is a method of reading an echo signal from the rising edge of the read gradient magnetic field pulse and sampling from the center of the k space.
- the straight arrow indicates the sampling start position of the echo signal
- the end point indicates the sampling end position
- the curved arrow indicates the sampling direction.
- the influence of the deviation from the trapezoidal shape of the rising and falling of the readout gradient magnetic field is one of the deviation of the k-space interval (first influence), the other One is the k-space shift (second effect), but in the k-space data collection method in Fig. 7 (a), the first effect is the left end of the k-space as shown in the gray range in the figure. That is, it is limited to the high frequency region, and the spatial resolution of the image is only slightly reduced.
- the second effect is the peak position shift of the echo signal.
- the sampling direction is parallel to the ky axis, the peak shift occurs at the same position of each echo and does not become discontinuous between echoes. For this reason, a linear phase change occurs in the image after Fourier transform, but this is excluded from the absolute value image.
- the first effect is limited to the high frequency region of the k space, as shown as the gray range in the figure, and has little effect as in the case of orthogonal sampling.
- the second effect since the echo signal sampling direction differs for each echo, the shift of the echo peak position becomes discontinuous between echoes, and the image formability of the image is lowered.
- the sampling method of FIG. 7 (c) the first influence occurs in the k-space low-frequency region that determines the contrast of the image, resulting in a significant deterioration in image quality. Therefore, by applying this embodiment, a significant improvement in image quality can be obtained in the radial method shown in FIGS. 7B and 7C.
- FIG. 8 shows the effect of performing k-space correction using approximate parameters according to this embodiment.
- FIG. 8 is a diagram showing an image when a uniform phantom is imaged by the radial method of FIG. 7C.
- FIG. 8A shows a case where correction is not performed
- FIG. 8B shows a case where correction is performed.
- the image is doubled and the signal unevenness is generated at the edge portion.
- the correction is performed, the signal unevenness disappears and is uniform.
- a phantom image is obtained.
- at least two or more approximate functions may include different types, at least one may be an exponential function, and at least one may be a linear function. Further, at least two or more approximate functions may have the same slope between adjacent sections.
- the imaging procedure of this embodiment is the same as that of the first embodiment shown in FIG. 2, but the parameters of the approximation function approximating the gradient magnetic field waveform are changed according to the imaging parameters calculated in step 201, and the optimum It is a feature that an approximate parameter is used. That is, the parameters of at least two approximation functions are changed according to the imaging conditions.
- those related to the change of the parameter of the approximate function are specifically imaging parameters for determining the intensity of the read gradient magnetic field pulse and rise / fall times.
- the error in the gradient magnetic field pulse waveform shape reflects the characteristics of the gradient magnetic field system, and the characteristics of this gradient magnetic field system greatly depend on the gradient magnetic field pulse intensity. Further, it is conceivable that the gradient magnetic field pulse waveform shape depends not only on the gradient magnetic field pulse intensity but also on the rise / fall time of the gradient magnetic field pulse. Therefore, in the present embodiment, it is possible to perform correction with higher accuracy by changing the approximate parameter in accordance with these imaging parameters.
- the parameters of the approximate function of the gradient magnetic field pulse waveform are stored in the storage device 200, as in the first embodiment, and the approximate function of the waveform is shown in FIGS. 6 (a) and (b).
- the gradient magnetic field pulse waveform is divided into a plurality of sections, and the shape of each section is approximated by a function represented by equations (2), (3), and the like.
- equations (2), (3) and the like.
- a to f in the expressions (2) and (3) are constants.
- these a to f are values that change according to the gradient magnetic field strength and the rise time of the gradient magnetic field.
- BW is the bandwidth
- ⁇ is the magnetic rotation ratio
- FOV is the imaging field size.
- a plurality of optimum values corresponding to each gradient magnetic field pulse intensity and the rise / fall time are set in advance. Find a combination of approximate parameters. The result thus obtained is stored in the storage device 200.
- k-space coordinate calculation is performed in step 203, first, gradient magnetic field pulse intensity and rise / fall time are calculated from imaging parameters used in actual imaging (step 202). Next, an optimal approximate parameter is selected from the storage device 200 according to the calculated gradient magnetic field pulse intensity and rise time, and used for k-space coordinate calculation.
- the method for calculating k-space coordinates is the same as in the first embodiment, and after that, gridding (step 204) and image reconstruction by Fourier transform (step 205) are also the same as in the first embodiment. is there.
- parameters of each gradient magnetic field pulse intensity, rise time and waveform approximation function (first In the example of the embodiment, a relationship with each constant a to f) is obtained in advance, and an optimum approximate parameter corresponding to each gradient magnetic field pulse intensity and rise time may be calculated based on this relationship. it can.
- each gradient magnetic field pulse intensity or rise / fall time and each constant of the waveform approximation function can be obtained, for example, by the following procedure.
- This approximate function can be determined, for example, by defining a base function for each constant and optimizing the base function by the least square method using the result of the measurement point.
- a base function there are a linear function, a quadratic function, an exponential function, a logarithmic function, and the like.
- the gradient magnetic field rise / fall time obtained by empirically determining the time interval of the gradient magnetic field waveform section (D start in Equation (2) and D end in Equation (3)) Is set based on the rising / falling time of the gradient magnetic field pulse determined in step 201.
- the present embodiment is applied to an imaging method in which the readout gradient magnetic field is generated by combining two or three axes of gradient magnetic field outputs that are physically orthogonal as in oblique imaging.
- the imaging procedure of this embodiment is shown in FIG.
- the imaging procedure is the same as that of the first or second embodiment, but the calculation for coordinate transformation described below (steps 2031 and 2032) is performed in the calculation of k-space coordinates in step 203.
- the difference is that is added.
- different points will be mainly described.
- the output of the gradient magnetic field pulse at the time of oblique is given by the synthesis of each element of the pulse sequence at each time t.
- the output of the readout gradient magnetic field pulse of two axes is G r1 (t) and G r2 (t), respectively
- the physical axis of the gradient magnetic field (X / Y / The output of Z) can be calculated by equation (6).
- Equation (6) The coefficient c in Equation (6) can be generally expressed as a rotation matrix.
- the gradient magnetic field pulse waveform shape is approximated to the gradient magnetic field pulse waveforms G x (t), G y (t), and G z (t) of each axis obtained here (step 2031). That is, the gradient magnetic field pulse waveform of each axis is divided into a plurality of sections and, for example, approximate functions of equations (2) and (3) are applied. Let the waveform after approximation be G ′ x (t), G ′ y (t), and G ′ z (t). That is, the approximated waveform can be expressed by the following general formula (7).
- the gradient magnetic field pulse waveform corresponding to the two axes in the imaging plane is calculated by the following equation (8) (step 2032).
- Equation (6) the gradient magnetic field pulse waveform obtained from Equation (6) is represented as the sum of components obtained by developing G r1 (t) and G r2 (t) on each axis as represented by Equation (9).
- the approximation of the gradient magnetic field pulse waveform shape in step 2031 may be performed for each term of equation (9). That is, the results obtained by approximating G r1 (t) and G r2 (t) with respect to the physical axis (X / Y / Z) of the gradient magnetic field are respectively G r1x (t), G r1y (t), G r1z (t), G r2x (t), G r2y (t), G r2z (t), gradient magnetic field pulse waveform G ' x (t) after approximation of physical axis (X / Y / Z) of gradient magnetic field , G ′ y (t) and G ′ z (t) are as follows.
- step 2032 using the gradient magnetic field pulse waveform corresponding to the two axes in the imaging plane calculated by equation (8), the application amount for each section of the gradient magnetic field waveform is integrated, and the calculation by equation (11) is performed.
- the k-space coordinates are obtained in the same manner as in the first embodiment (step 2033).
- k-space coordinates can be corrected so as to eliminate the system error of the gradient magnetic field output by using the approximate function of the gradient magnetic field pulse waveform shape. And an image with little deterioration in image quality can be obtained.
- the optimal approximation function (approximation parameter) stored in the storage device 200 is stored, and k-space coordinates where the echo signal (sampling data) obtained by imaging is arranged.
- the approximate parameters stored in the storage device are applied as they are or after being changed according to actual imaging parameters or imaging sections.
- An MRI apparatus provided with means for calculating an approximate parameter of an optimum gradient magnetic field pulse waveform based on the shape of the gradient magnetic field pulse waveform used in the sequence will be described.
- the configuration of the apparatus and the imaging procedure are the same as those in the first to third embodiments.
- FIG. 10 shows the procedure for calculating the optimum approximate parameter according to this embodiment.
- a pulse sequence for preliminary measurement is calculated (step 101).
- the pulse sequence of the preliminary measurement may be the same as or different from the pulse sequence of the main measurement as long as the pulse sequence uses a readout gradient magnetic field similar to the pulse sequence of the main measurement for acquiring the image of the subject.
- the pulse sequence is a pulse sequence of the radial method for both the main measurements.
- the pulse sequence of the radial method as shown in FIG. 3, since the readout gradient magnetic field in the biaxial direction is used, the optimum approximation parameter is searched for each.
- step 102 Execute the pulse sequence calculated in step 101 and measure the signal (step 102). This measurement is preferably performed using a phantom with uniform contents in order to evaluate a search result to be described later.
- an optimum approximation parameter is searched for the gradient magnetic field pulse waveform in each axial direction (step 103).
- the approximate parameter retrieved in step 103 is stored in the storage device (memory or storage device) 100 (step 104).
- the subsequent procedure in the main measurement is the same as the imaging procedure in the first embodiment (steps 201 to 205 in FIG. 2), and thus the description thereof is omitted here.
- FIG. 11 is a flowchart showing details of step 103.
- step 111 set initial values of approximate parameters.
- the initial value for example, those obtained by empirically determining the constants a to f of the approximate functions of the expressions (2) and (3) described in the first embodiment can be used.
- the coordinates of the measurement space are calculated based on the approximate parameters set in step 111 and the gradient magnetic field pulse waveform created in step 101 of FIG. 10 (112). Specifically, as shown in FIG. 12, first, an approximate parameter is applied to the input gradient magnetic field pulse waveform to obtain a corrected gradient magnetic field pulse waveform (step 121). Next, in step 121, the coordinates of the measurement space are calculated from the corrected gradient magnetic field pulse waveform by equation (1) (step 122). These steps 121 and 122 are executed for the biaxial gradient magnetic field pulse waveform.
- FIG. 12 shows a calculation example for the three axes of the X axis, the Y axis, and the Z axis, but one radial method pulse sequence may be calculated for two of these three axes.
- coordinate calculation may be performed for one coincident axis. Since coordinate calculation can be executed independently for each axis, the calculation order is not limited to that shown in the figure.
- step 113 data of the measurement space is created by the gridding process.
- step 114 the data after gridding is Fourier transformed to create an image.
- FIGS. 13 (a) and 13 (b) are diagrams showing reconstructed images when k-space coordinate correction is performed using different combinations of approximate parameters, where the left is an image and the right is an AA of the image. 'Shows the signal strength profile of the line. Since this image has a uniform phantom content, ideally, the signal intensity profile has a constant signal value in the region where the phantom exists.
- the signal value of the portion where the phantom does not exist becomes zero.
- the signal lift can be confirmed at the phantom edge.
- the signal of the center part of the phantom part is high, and it becomes low as it goes outside.
- a signal rise at the edge of such an image is defined as “Overshoot”, and the uniformity of the signal inside the phantom is defined as “Uniformity”, and a value is calculated for each approximate parameter.
- the “Overshoot” value for example, the average value or maximum value of the signal within the ROI set at the edge can be used, and as the “Uniformity” value, for example, set in the area where the phantom exists The standard deviation of the signal within the ROI can be used. Then, when the value of “Overshoot” is decreased by changing the approximate parameter, and when the value of “Uniformity” is decreased, it is determined that the image quality is improved.
- the above steps 111 to 115 are repeated while changing the approximate parameter.
- the approximation parameter may be changed for each of a plurality of approximation parameters (a to f), or a plurality of approximation functions for each of the approximation functions of a plurality of approximation functions (for example, the functions of equations (2) and (3)).
- the combination of approximate parameters may be changed. It is preferable to execute the change of the approximate parameter by determining the search range of the parameter value and the interval or the number of searches of the search value. In addition, in order to shorten the search time, considering the degree of influence of the approximate parameter on the image quality, the optimum value is searched only for the parameter having a large influence, or the number of searches for the parameter having a large influence is increased. Changes may be made.
- the time interval D (for example, D start in Expression (2) or D end in Expression (3)) of each section obtained by dividing the gradient magnetic field pulse waveform shape may be added to the parameter.
- the durations D start and D end of the initial and final periods 604 and 606 of the gradient magnetic field pulse are made variable, and the relationship of Equation (14) is satisfied. You can change each value as you want.
- RiseTime is the rise / fall time of the gradient magnetic field pulse obtained from the calculation in step 101.
- step 116 it is determined whether all the approximate parameter combinations have been calculated. If the result of determination at step 116 is “No”, steps 111 to 115 are repeated again. When all combinations are calculated, the axis of the readout gradient magnetic field is made different, and the optimum parameter of the approximate function of the gradient magnetic field pulse waveform is searched for the different axis.
- the target pulse sequence in step 101 of gradient magnetic field pulse waveform calculation and step 102 of signal measurement is a pulse sequence of the radial method
- two measurements are performed to support three gradient magnetic field axes.
- Approximate parameters can be searched.
- the Z axis of the gradient magnetic field is assigned to the slice selection gradient magnetic field axis
- the remaining X and Y axes are assigned to the gradient magnetic field axes in the slice plane.
- the Y axis of the gradient magnetic field is assigned.
- the slice selection gradient magnetic field axis and the remaining X and Z axes are used as gradient magnetic field axes in the slice plane.
- the approximate parameter for the X axis and the Y axis is determined from the first measurement
- the approximate parameter for the Z axis is determined from the second measurement.
- measurement is performed for each of X, Y, and Z, and approximate parameters for each axis are determined.
- the axis search order is executed in the order of, for example, the X, Y, and Z axes.
- the search order of the axes is not limited to this, and an optimal order can be determined according to the hardware configuration of the apparatus.
- step 117 It is determined whether all the axes for searching for approximate parameters have been completed (step 117). If the result is “No”, steps 111 to 116 are repeated again. In the case of “Yes”, among the evaluation values calculated in step 115, the combination that gives the highest image quality improvement, for example, the combination that gives the smallest “Overshoot” value and “Uniformity” value is determined as the optimum approximation parameter. ,Output.
- the approximate parameter retrieved in this way is stored in the storage device 100 (FIG. 10), and is applied to k-space coordinate correction calculation in the main measurement to be executed thereafter.
- an approximate parameter for each axis of the gradient magnetic field is obtained by preliminary measurement, and reflected in the data of the main measurement, so that an image with few artifacts can be obtained even in a change in imaging section or oblique imaging. Obtainable.
- This embodiment is the same as the fourth embodiment in that the approximate parameters of each axis of the gradient magnetic field are obtained by preliminary measurement, but this embodiment is characterized in that the search is performed in a plurality of stages. .
- the approximate parameter search procedure according to this embodiment is shown in FIG. 14, steps denoted by the same reference numerals as those in FIG. 10 are processed in the same manner as the steps in FIG.
- the difference from FIG. 10 is that there are two approximate parameter search steps 141 and 142. In these two approximate parameter search steps, the parameter search ranges for calculating approximate parameters are made different.
- step 141 using the gradient magnetic field pulse waveform of the pulse sequence created in step 101 and the measurement signal measured in step 102, the optimum approximate parameter is searched for as in the fourth embodiment (i.e., FIG. 11).
- the optimum approximate parameter searched in step 1402 is set as approximate parameter 1.
- step 142 the approximate parameter 1 searched in step 141 is used as a reference, and the approximate parameter is searched at a finer pitch (amount to change the parameter for each loop in FIG. 11) than in step 141.
- the pitch used for the approximate parameter search in the second search step 142 is set to 1/10.
- the processing at this time is also the same as the processing shown in FIG.
- the optimum approximate parameter searched in the second search step 142 is set as approximate parameter 2.
- the retrieved approximate parameter 2 is recorded in the storage device 100 in step 104. Thereafter, using this approximate parameter, k-space coordinates are calculated for the actually measured measurement signal, and the image is reconstructed as in the above embodiments.
- the search step of the approximate parameter is divided into two times, and the search is performed with different pitches, so that the efficiency is reduced more efficiently than the search with a fine pitch from the beginning. It is possible to search for the optimum parameters without doing so.
- This embodiment is the same as the fourth and fifth embodiments in that the approximate parameters of each axis of the gradient magnetic field are obtained by preliminary measurement, but in this embodiment, refer to the determination when obtaining the optimal approximate parameters.
- the feature is based on images.
- the imaging procedure of this embodiment is shown in FIG. The difference from FIG. 10 is that there is a reference image measurement step.
- a pulse sequence for a reference image is calculated and signal measurement is performed (step 151).
- the acquired reference image is stored in the storage device 100 (FIG. 15).
- a pulse sequence for preliminary measurement is calculated (step 101), and signal measurement is performed (step 102).
- the optimum value of the approximate parameter is searched using the reference image acquired in step 151 and stored in the storage device 100, the measurement signal measured in step 102, and the gradient magnetic field pulse waveform calculated in step 101 (step 152).
- step 152 The internal processing of step 152 is shown in FIG.
- the difference from the processing of FIG. 11 described in the fourth embodiment is the processing content of the approximate parameter evaluation value calculation step 161.
- an image acquired by an MRI apparatus includes the effects of transmission sensitivity and reception sensitivity in addition to the effects of errors in the gradient magnetic field pulse waveform shape.
- a reference image is prepared and the deviation from the image is evaluated.
- the reference image it is possible to use a pulse sequence based on orthogonal sampling with little influence of an error of the readout gradient magnetic field, for example, an image captured by a normal gradient echo sequence.
- the degree of coincidence between such a reference image and the image obtained in step 114 is evaluated.
- the similarity can be used for evaluating the coincidence.
- the pixel value vector R in the reference image and the pixel value vector I of the image to be evaluated are defined by Expression (15).
- the similarity is calculated by the following calculation.
- the approximation parameter used in step 112 is changed, and the above steps 112 to 115 are repeated a predetermined number of times to obtain the degree of similarity and used when the degree of similarity is the highest.
- the approximate parameter is the optimum approximate parameter.
- FIG. 17 schematically shows the relationship between the similarity and the approximate parameter. From FIG. 17, the point 171 having the highest similarity can be determined as the optimum value of the approximate parameter. Further, the image reconstructed with the k-space coordinates corrected using the approximation function of the approximation parameter determined to be the optimum value is the target image of the main measurement.
- ROI region of interest
- ROI may be a line as well as a surface.
- an ROI can be provided in a region where the degree of influence is most significant. For example, when evaluating an approximate parameter that affects the uniformity of an image, an ROI is provided near the center of the image, and when evaluating an approximate parameter that affects the image formation at the edge of the image, an ROI is provided near the edge of the image. Is provided.
- FIG. 18 shows an example of a reconstructed image created according to this embodiment.
- Fig. 18 (a) Image (reference image) captured with a gradient echo sequence of normal imaging, (b) and (c) reconstruct the data obtained from the radial method gradient echo sequence using different approximation functions. It is an image.
- the graph on the right side of the figure shows a signal intensity profile on a line (indicated by a dotted line) set at the center of the image.
- the signal intensity profile of the image changes due to the difference in the approximate parameters. From this, it can be seen that an optimum approximation parameter can be obtained by evaluating the degree of coincidence between an image created using different approximation parameters and a reference image (in this example, the image of FIG. 18A).
- the optimal approximation parameter can be quantitatively evaluated by evaluating the approximation parameter based on the reference image.
- the present invention is not limited to the contents disclosed in the above embodiments, and can take various forms based on the gist of the present invention.
- the radial method does not depend on the type of the pulse sequence and can be applied to the spin echo type radial method.
- the present invention can also be applied to the radial method of asymmetric echo measurement in which sampling is performed from the center of the measurement space toward the outside. Further, the present invention can be similarly applied to the radial method of sampling from the outside of the measurement space toward the center.
- radial methods that sample in an unspecified direction of the measurement space, such as a radial method in a three-dimensional space, a radial method that samples from the center of the measurement space to the outside, and then returns to the center, etc. It is the same.
- the present invention can be similarly applied to a hybrid radial method combining a radial method and a phase encoding gradient magnetic field pulse.
- an example of an exponential function has been shown as an approximate function of a gradient magnetic field pulse waveform, the example of the function is not limited to this.
- Various forms of functions can be applied depending on the system configuration.
- the system response of the gradient magnetic field output can be considered not only for the radial method but also for all pulse sequences that can be executed by the MRI apparatus.
- image quality such as the echo planar method and fast spin echo method, which acquire multiple echo signals with a single RF irradiation.
- high-precision approximation is possible by approximating the gradient magnetic field pulse waveform shape with at least two approximation functions, and it can respond to changes in imaging conditions by using the parameters of the approximation function stored in the storage device Thus, correction calculation can be easily performed.
- 2 static magnetic field generation system 3 gradient magnetic field generation system, 4 sequencer, 5 transmission system, 6 reception system, 7 signal processing system, 8 CPU, 100, 200 storage device.
Landscapes
- Physics & Mathematics (AREA)
- Condensed Matter Physics & Semiconductors (AREA)
- General Physics & Mathematics (AREA)
- High Energy & Nuclear Physics (AREA)
- Nonlinear Science (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Engineering & Computer Science (AREA)
- Signal Processing (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
Description
本実施形態は、傾斜磁場パルスを近似する関数のパラメータ(近似パラメータ)を記憶装置に格納し、近似パラメータを用いて、計測されたエコー信号の計測空間(k空間)座標を計算し、グリッディングすることが特徴である。以下、図2に示す撮像動作のフローを参照して、本実施形態を説明する。
最後に作成したk空間データを2次元フーリエ変換して画像を作成する(ステップ205)。
本実施形態の撮像手順は、図2に示す第一の実施形態と同様であるが、傾斜磁場波形を近似する近似関数のパラメータを、ステップ201で計算される撮像パラメータに応じて変更し、最適な近似パラメータを用いることが特徴である。すなわち、少なくとも2つ以上の近似関数のパラメータは、撮像条件により変更する。
この結果から、各定数と傾斜磁場パルスについて近似関数を求める。この近似関数は、例えば、定数毎に、基盤となる関数を定義しておき、基盤となる関数を測定点の結果を用いて最小二乗法で最適化することにより、決定できる。基盤となる関数としては、1次関数や2次関数、指数関数や対数関数などがある。
本実施形態は、オブリーク撮像のように、読出し傾斜磁場が、物理的に直交する2軸或いは3軸の傾斜磁場出力の合成によって発生される撮像方法に適用される。
以下、異なる点を中心に説明する。
本実施形態による最適近似パラメータの算出手順を図10に示す。
画質向上の評価は、例えば、ステップ101でファントムを用いて計測した信号から再構成した画像のプロファイル形状から行う。画質の判定基準の例を、図13に示す。図13(a)、(b)は、それぞれ異なる近似パラメータの組み合わせを用いてk空間座標補正を行った場合の再構成画像を示す図で、図の左は画像、右は画像のA-A’ラインの信号強度プロファイルを示す。この画像は、内容物が均一なファントムなので、理想的には信号強度プロファイルは、ファントムの存在する領域では信号値が一定となる。またファントムが存在しない部分の信号値はゼロとなる。しかし、図ではファントム縁部で信号の持ち上がりが確認できる。また、ファントム部の中心部の信号が高く、外側へ向うにつれて低くなっている。このような画像の縁部での信号の持ち上がりを「Overshoot」、ファントム内部の信号の均一さを「Uniformity」と定義して、近似パラメータ毎に値を算出する。「Overshoot」の値としては、例えば、縁部に設定したROI内の信号の平均値や最大値を用いることができ、「Uniformity」の値としては、例えば、ファントムが存在する領域内に設定したROI内の信号の標準偏差を用いることができる。そして、近似パラメータを変更することによって「Overshoot」の値が減少したとき、また「Uniformity」の値が減少したときに、それぞれ、画質が向上したと判定する。
本実施形態も予備計測によって傾斜磁場の各軸の近似パラメータを求める点は、上記第四の実施形態と同じであるが、本実施形態は、検索を複数段階に分けて行うことが特徴である。
その後、この近似パラメータを用いて、実際に計測された計測信号についてk空間座標計算し、画像を再構成することは、上記各実施形態と同様である。
本実施形態も予備計測によって傾斜磁場の各軸の近似パラメータを求める点は、上記第四及び第五の実施形態と同じであるが、本実施形態では最適な近似パラメータを得る際の判定を参照画像をもとに行うことが特徴である。
Claims (15)
- 静磁場発生手段と、
前記静磁場発生手段が発生する静磁場に磁場勾配を与える傾斜磁場発生手段と、
前記静磁場の中に配置された検査対象に高周波磁場を照射する送信手段と、
前記検査対象から核磁気共鳴により放出されるエコー信号を受信する受信手段と、
前記エコー信号を処理し、前記検査対象の画像を再構成する信号処理手段と、を備えた磁気共鳴イメージング装置において、
前記信号処理手段は、前記傾斜磁場発生手段によって磁場勾配が与えられた傾斜磁場パルス波形の少なくとも2つ以上に時間的に分割されたそれぞれの区間毎に前記傾斜磁場パルス波形を近似する近似関数を生成する近似関数生成手段と、前記近似関数を用いて、前記エコー信号が配置される計測空間座標を補正する補正手段と、を備えたことを特徴とする磁気共鳴イメージング装置。 - 請求項1記載の磁気共鳴イメージング装置であって、
前記少なくとも2つ以上の近似関数は、種類の異なるものを含むことを特徴とする磁気共鳴イメージング装置。 - 請求項2記載の磁気共鳴イメージング装置であって、
前記種類の異なるものを含む近似関数は、少なくとも1つが指数関数であり、少なくとも1つが線形関数であることを特徴とする磁気共鳴イメージング装置。 - 請求項1記載の磁気共鳴イメージング装置であって、
前記近似関数生成手段は、前記エコー信号を用いて、前記近似関数を表すパラメータを算出することにより、該近似関数を近似するパラメータ算出手段を備えたことを特徴とする磁気共鳴イメージング装置。 - 請求項4記載の磁気共鳴イメージング装置であって、
前記計測したエコー信号は、傾斜磁場パルス波形形状を計測するための予備計測パルスシーケンスを複数の軸の傾斜磁場について実行することにより得たものであり、
前記近似関数生成手段は、前記複数の軸の傾斜磁場パルス波形について、それぞれ近似関数のパラメータを算出することを特徴とする磁気共鳴イメージング装置。 - 請求項4に記載の磁気共鳴イメージング装置であって、
前記少なくとも2つ以上の近似関数のパラメータは、撮像条件により変更することを特徴とする磁気共鳴イメージング装置。 - 請求項4記載の磁気共鳴イメージング装置であって、
前記パラメータ算出手段は、前記画像から前記画像の画質向上度を評価する評価値を計算する評価値算出手段と、前記評価値算出手段が算出した評価値をもとにパラメータの最適値を判定する判定手段とを備えたことを特徴とする磁気共鳴イメージング装置。 - 請求項7記載の磁気共鳴イメージング装置であって、
前記パラメータ算出手段は、前記画像に関心領域を設定し、当該関心領域の信号値をもとに前記評価値を計算し、前記評価値と前記画像を作成するときに使用したパラメータの値との関係から、パラメータの最適値を求めることを特徴とする磁気共鳴イメージング装置。 - 請求項1記載の磁気共鳴イメージング装置であって、
前記少なくとも2つ以上の近似関数は、隣り合う区間の間で互いに傾きが同じになっていることを特徴とする磁気共鳴イメージング装置。 - 磁気共鳴イメージング装置で計測したエコー信号の計測空間座標を補正する方法であって、
傾斜磁場パルス波形形状を少なくとも2つの近似関数で近似するステップと、
近似関数のパラメータを記憶手段に格納するステップと、
前記記憶手段に格納されたパラメータと計測したエコー信号とを用いて、エコー信号の計測空間座標を補正するステップと、
を含む計測空間座標補正方法。 - 請求項10に記載の計測空間座標補正方法であって、
前記近似するステップは、複数の撮像条件に対応して異なる複数の傾斜磁場パルス波形形状に対し、それぞれ近似関数を設定することを特徴とする計測空間座標補正方法。 - 請求項10に記載の計測空間座標補正方法であって、
前記近似するステップは、近似関数のパラメータを最適化するステップを含むことを特徴とする計測空間座標補正方法。 - 請求項12に記載の計測空間座標補正方法であって、
前記最適化するステップは、パラメータの初期値を用いて補正した計測空間座標と計測されたエコー信号とを用いて画像を再構成する画像再構成ステップを、前記パラメータの初期値を変えながら繰り返し、パラメータの最適値を求めることを特徴とする計測空間座標補正方法。 - 請求項13に記載の計測空間座標補正方法であって、
前記最適化するステップは、前記画像再構成ステップの後に、
前記画像から前記画像の画質向上度を評価する評価値を計算する評価値算出ステップと、
算出された評価値をもとにパラメータの最適値を判定する判定ステップと、
を含むことを特徴とする計測空間座標補正方法。 - 磁気共鳴イメージング装置で計測したエコーを計測空間座標に配置して計測空間データを作成するステップと、前記計測データを高速フーリエ変化して画像を再構成するステップとを含む画像再構成方法であって、
前記計測空間座標は、請求項11記載の計測空間座標補正方法によって補正されていることを特徴とする画像再構成方法。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US13/881,207 US9714994B2 (en) | 2010-11-09 | 2011-10-28 | Magnetic resonance imaging apparatus, and method for correcting measurement-space coordinates |
JP2012542868A JP5846450B2 (ja) | 2010-11-09 | 2011-10-28 | 磁気共鳴イメージング装置、計測空間座標補正方法、及び、画像再構成方法 |
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2010-251143 | 2010-11-09 | ||
JP2010251143 | 2010-11-09 |
Publications (1)
Publication Number | Publication Date |
---|---|
WO2012063654A1 true WO2012063654A1 (ja) | 2012-05-18 |
Family
ID=46050806
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/JP2011/074892 WO2012063654A1 (ja) | 2010-11-09 | 2011-10-28 | 磁気共鳴イメージング装置、計測空間座標補正方法、及び、画像再構成方法 |
Country Status (3)
Country | Link |
---|---|
US (1) | US9714994B2 (ja) |
JP (1) | JP5846450B2 (ja) |
WO (1) | WO2012063654A1 (ja) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2018167025A (ja) * | 2017-03-29 | 2018-11-01 | キヤノンメディカルシステムズ株式会社 | 磁気共鳴イメージング装置および医用処理装置 |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8903152B2 (en) * | 2012-06-29 | 2014-12-02 | General Electric Company | Methods and systems for enhanced tomographic imaging |
CN104297708A (zh) * | 2014-10-22 | 2015-01-21 | 浙江大学 | 磁共振成像***中的旋转倾斜梯度线圈组件 |
DE102017200446A1 (de) * | 2017-01-12 | 2018-07-12 | Siemens Healthcare Gmbh | Korrektur eines MR-Sendesignals |
JP7377231B2 (ja) * | 2021-02-02 | 2023-11-09 | 富士フイルムヘルスケア株式会社 | 磁気共鳴イメージング装置、及び、画像処理方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2010074059A1 (ja) * | 2008-12-26 | 2010-07-01 | 株式会社 日立メディコ | 磁気共鳴イメージング装置及び読み出し傾斜磁場誤差補正方法 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6969991B2 (en) * | 2002-12-11 | 2005-11-29 | The Board Of Trustees Of The Leland Stanford Junior University | Correction of the effect of spatial gradient field distortions in diffusion-weighted imaging |
DE102008061532B3 (de) * | 2008-12-10 | 2010-08-12 | Siemens Aktiengesellschaft | Verfahren und Vorrichtung zur Ermittlung von Verzeichnungskorrekturdaten |
JP2012040362A (ja) * | 2010-07-23 | 2012-03-01 | Toshiba Corp | 磁気共鳴イメージング方法、磁気共鳴イメージング装置およびその制御装置 |
WO2013002233A1 (ja) * | 2011-06-30 | 2013-01-03 | 株式会社 日立メディコ | 磁気共鳴イメージング装置および傾斜磁場波形推定方法 |
WO2015164793A1 (en) * | 2014-04-25 | 2015-10-29 | Mayo Foundation For Medical Education And Research | Integrated image reconstruction and gradient non-linearity correction for magnetic resonance imaging |
-
2011
- 2011-10-28 WO PCT/JP2011/074892 patent/WO2012063654A1/ja active Application Filing
- 2011-10-28 US US13/881,207 patent/US9714994B2/en active Active
- 2011-10-28 JP JP2012542868A patent/JP5846450B2/ja active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2010074059A1 (ja) * | 2008-12-26 | 2010-07-01 | 株式会社 日立メディコ | 磁気共鳴イメージング装置及び読み出し傾斜磁場誤差補正方法 |
Non-Patent Citations (5)
Title |
---|
ATSUSHI TAKAHASHI ET AL.: "Compensation of multi-dimensional selective excitation pulses using measured k-space trajectories", MAGNETIC RESONANCE IN MEDICINE, vol. 34, no. 3, September 1995 (1995-09-01), pages 446 - 456 * |
DANIEL M. SPIELMAN ET AL.: "Spiral imaging on a small-bore system at 4.7T", MAGNETIC RESONANCE IN MEDICINE, vol. 34, no. 4, October 1995 (1995-10-01), pages 580 - 585 * |
H. TAN ET AL.: "Estimation of k-space trajectories in spiral mri", MAGNETIC RESONANCE IN MEDICINE, vol. 61, no. 6, June 2009 (2009-06-01), pages 1396 - 1404 * |
H. TAN ET AL.: "K-space trajectory estimation in spiral imaging", PROC. INTL. SOC. MAG. RESON. MED., vol. 15, May 2007 (2007-05-01), pages 981 * |
J. Y. CHENG ET AL.: "Gradient linear system modeling using gradient characterization", PROC. INTL. SOC. MAG. RESON. MED., vol. 16, May 2008 (2008-05-01), pages 1155 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2018167025A (ja) * | 2017-03-29 | 2018-11-01 | キヤノンメディカルシステムズ株式会社 | 磁気共鳴イメージング装置および医用処理装置 |
JP7123595B2 (ja) | 2017-03-29 | 2022-08-23 | キヤノンメディカルシステムズ株式会社 | 磁気共鳴イメージング装置および医用処理装置 |
Also Published As
Publication number | Publication date |
---|---|
JPWO2012063654A1 (ja) | 2014-05-12 |
US9714994B2 (en) | 2017-07-25 |
JP5846450B2 (ja) | 2016-01-20 |
US20130207654A1 (en) | 2013-08-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP5661473B2 (ja) | 磁気共鳴イメージング装置及び読み出し傾斜磁場誤差補正方法 | |
JP5399240B2 (ja) | 磁気共鳴イメージング装置及び傾斜磁場に起因する誤差補正方法 | |
US10401456B2 (en) | Parallel MR imaging with Nyquist ghost correction for EPI | |
CN103238082B (zh) | 使用多点Dixon技术和低分辨率校准的MR成像 | |
US9766313B2 (en) | MR imaging using apt contrast enhancement and sampling at multiple echo times | |
JP6018401B2 (ja) | 拡散強調エコープラナー撮像法において高次渦電流に誘発された歪みを予測補正するためのシステムおよび方法 | |
US11067653B2 (en) | Magnetic resonance imaging method including acquiring echo signals from each k-space slice with different relaxation time weightings after echo signals from different k-space splices are acquired | |
RU2538421C2 (ru) | Картирование градиента восприимчивости | |
JP5828763B2 (ja) | 磁気共鳴イメージング装置および磁気共鳴イメージング方法 | |
US9846215B2 (en) | MRI embodiments for controlling an arrangement order of multiple echoes in a k-space | |
RU2686879C2 (ru) | Мр томография с нулевой эхо-задержкой | |
JP5846450B2 (ja) | 磁気共鳴イメージング装置、計測空間座標補正方法、及び、画像再構成方法 | |
US7689015B2 (en) | Magnetic resonance imaging apparatus and image correction estimating method | |
JP5638393B2 (ja) | 磁気共鳴イメージング装置及び方法 | |
US10126395B2 (en) | Magnetic resonance imaging apparatus and magnetic resonance imaging method | |
US10996305B2 (en) | Propeller MR imaging with relaxation weightings varying between k-space blades and k-space data lines | |
JP5808659B2 (ja) | 磁気共鳴イメージング装置及びT1ρイメージング法 | |
EP3185029A1 (en) | Mr imaging using propeller acquisition with t2 decay correction | |
US20200300948A1 (en) | Dixon type water/fat separation mr imaging with improved fat shift correction | |
JP2024506575A (ja) | Dixonタイプ水/脂肪分離MR画像 | |
JP2005028019A (ja) | 磁気共鳴イメージング装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
121 | Ep: the epo has been informed by wipo that ep was designated in this application |
Ref document number: 11840159 Country of ref document: EP Kind code of ref document: A1 |
|
ENP | Entry into the national phase |
Ref document number: 2012542868 Country of ref document: JP Kind code of ref document: A |
|
WWE | Wipo information: entry into national phase |
Ref document number: 13881207 Country of ref document: US |
|
NENP | Non-entry into the national phase |
Ref country code: DE |
|
122 | Ep: pct application non-entry in european phase |
Ref document number: 11840159 Country of ref document: EP Kind code of ref document: A1 |