The present application claims priority under 35 U.S.C. § 119 to Japanese Patent Application No. 2023-150407, filed Sep. 15, 2023, the content of which is hereby expressly incorporated by reference, in its entirety, into the present application.
The present invention relates to a magnetic resonance imaging apparatus (hereinafter, referred to as an MRI apparatus), and particularly, to a technology of controlling a data sampling pattern of the MRI apparatus.
An MRI apparatus encodes nuclear magnetic resonance signals by collecting the nuclear magnetic resonance signals while applying gradient magnetic fields having different magnitudes to a static magnetic field and assigns positional information to the nuclear magnetic resonance signals, which are a sum of nuclear spins having spatial distributions. The encoded nuclear magnetic resonance signals are collected as k-space data in a frequency domain, and the k-space data is converted into an image through calculations such as FFT.
The k-space data is 2D or 3D data having axes corresponding to gradient magnetic fields in three axial directions in a real space, and in general imaging, Cartesian sampling is performed to collect the k-space data along the axial directions. In addition, there is also radial sampling in which sampling is radially performed from a center of a k-space.
Since the radial sampling is a technique of frequently sampling a central region of the k-space, that is, a low-frequency region contributing to image contrast, the radial sampling is an imaging method that is robust against motion and can suppress blood flow artifacts. However, there are techniques such as sampling where an angle changes for each radial line or using a plurality of sampling lines for each angle for measurement (blade type radial measurement). A plurality of lines at the same angle are called blades.
A general problem of MRI is to reduce artifacts and to shorten an imaging time, and various techniques have been proposed to solve the problem. For example, Magnetic Resonance in Medicine 57:1086-1098 (2007), “Undersampled Radial MRI with Multiple Coils. Iterative Image Reconstruction Using a Total Variation Constraint” discloses a technology of reducing artifacts by performing iterative reconstruction to minimize a difference between each channel data and channel data after sensitivity correction in radial sampling MRI in which a single line is rotated. In addition, JP2021-029777A discloses that k-space sampling is controlled to reduce body movement artifacts, full sampling is performed in a center region, and thinning-out sampling is performed in a peripheral region.
For the blade type radial measurement, reducing artifacts while shortening the imaging time is required.
In the blade type radial measurement, measurements are repeatedly performed for each blade by varying the angle of the blade, and a circular region centered on an origin of the k-space is measured, and a sampling pattern is determined by the number of blades, that is, the number of shots, and the number of lines constituting each blade, that is, echo train length (ETL). Here, in order to accelerate imaging, it is necessary to set the ELT or the number of shots to be smaller and to reduce the radial sampling lines, but reducing the sampling lines may lead to a gap (a region where data cannot be collected) occurring in a high-frequency region, causing streak artifacts.
To address such a problem of the blade type radial measurement, the technology described in Magnetic Resonance in Medicine 57:1086-1098 (2007), “Undersampled Radial MRI with Multiple Coils. Iterative Image Reconstruction Using a Total Variation Constraint” employs iterative reconstruction to reduce artifacts, but image quality may change depending on the sampling pattern, and there is a possibility that artifact reduction may not be achieved. In addition, the technology described in JP2021-029777A is based on Cartesian sampling as its fundamental approach, so that optimization of the sampling pattern in the blade type radial measurement in which sampling patterns differ cannot be applied.
An object of the present invention is to provide an MRI apparatus capable of performing optimal sampling from the viewpoint of improving image quality and reducing artifacts in a case of performing a blade type radial measurement.
In order to achieve the above-described object, the present invention makes it possible to suppress artifacts while reducing the number of blades by controlling a sampling pattern in a blade in the blade type radial measurement.
That is, according to an aspect of the present invention, there is provided an MRI apparatus comprising: an imaging unit that collects k-space data by measuring a nuclear magnetic resonance signal generated by a subject in accordance with a predetermined pulse sequence; and a controller that controls sampling of the nuclear magnetic resonance signal. The controller performs control of radially sampling a k-space by varying an angle of a blade configured with a plurality of parallel sampling lines, and in the control, sets interline intervals between the plurality of parallel sampling lines to uneven intervals for at least one of a plurality of the blades having different angles.
In addition, according to another aspect of the present invention, there is provided a control method for an MRI apparatus including an imaging unit that collects k-space data by measuring a nuclear magnetic resonance signal generated by a subject in accordance with a predetermined pulse sequence, the control method comprising: performing control of radially sampling a k-space by varying an angle of a blade configured with a plurality of parallel sampling lines, and in the control, setting interline intervals between the plurality of parallel sampling lines to uneven intervals for at least one of a plurality of the blades having different angles.
According to the aspects of the present invention, by setting the interline intervals in the blade to uneven intervals, it is possible to mitigate a significant difference in sampling density between a high-frequency region and a low-frequency region of the k-space and to optimize the overall sampling density. As a result, it is possible to reduce streak artifacts caused by the sparse high-frequency region even in a case in which the number of blades is reduced.
Hereinafter, embodiments of an MRI apparatus according to the present invention will be described with reference to the accompanying drawings.
The MRI apparatus comprises, in a broad sense, an imaging unit that collects data for reconstructing an image of a subject by measuring a nuclear magnetic resonance signal generated from the subject, and a computer that controls an operation of the imaging unit and that performs various calculations using the data collected by the imaging unit. In the MRI apparatus of the embodiment of the present invention, a configuration of the imaging unit is similar to that of a general MRI apparatus, but the MRI apparatus of the embodiment of the present invention is characterized in that a function of controlling and optimizing sampling of a k-space in a case in which the imaging unit performs a blade type radial measurement is provided as a main function of the computer.
First, an outline of a configuration of the MRI apparatus will be described. As shown in
A configuration of the imaging unit 10 is similar to that of a general MRI apparatus and comprises a static magnetic field magnet 101 that forms a uniform static magnetic field in the imaging space, gradient magnetic field coils 102 that provide magnetic field gradients in three orthogonal axial directions to the static magnetic field, an RF transmission coil 103 for irradiating the subject 50 with a high-frequency magnetic field, an RF reception coil 104 that receives a nuclear magnetic resonance signal generated from the subject 50, and a sequencer 108.
The gradient magnetic field coils 102 in the three axial directions are each connected to a gradient magnetic field power supply 105 and generate gradient magnetic fields corresponding to power supplied from the gradient magnetic field power supply 105. In addition, by combining the gradient magnetic fields in the three axial directions, it is possible to generate a gradient magnetic field of any magnitude in any direction.
The RF transmission coil 103 is connected to a transmitter 106 provided with a high-frequency magnetic field generator, a high-frequency amplifier, and the like and irradiates the subject with a high-frequency magnetic field having a predetermined frequency as an RF pulse in response to an output from the transmitter 106.
The RF reception coil 104 is connected to a receiver 107 provided with an amplifier, a quadrature detector, an A/D converter, and the like. The receiver 107 samples the nuclear magnetic resonance signal by using the A/D converter and passes the nuclear magnetic resonance signal to a signal processing system (computer 30) as digital data. The RF reception coil 104 may be a multi-channel coil consisting of a plurality of small coils, and in that case, data is collected for each channel of each small coil.
The sequencer 108 operates the gradient magnetic field power supply 105, the transmitter 106, and the receiver 107 in accordance with a predetermined pulse sequence. Various pulse sequences are stored in the apparatus, and in imaging, the sequencer 108 calculates a pulse sequence to be used for the imaging by using an imaging condition and an imaging parameter set via the UI unit 40 or the like and performs the imaging based on the calculated pulse sequence. In the present embodiment, imaging based on the blade type radial measurement is performed.
The computer 30 can be configured as a computer provided with a CPU and a memory and implements a function as a control unit, such as an imaging control unit 310 and a display control unit 320, and a function as a calculation unit, such as an image reconstruction unit 330. The imaging control unit 310 controls an operation of the imaging unit 10 via the sequencer 108. In addition, in the MRI apparatus 1 of the present embodiment, the imaging control unit 310 controls a sampling pattern of each blade in a case of the blade type radial measurement.
As shown in the drawing, the imaging control unit 310 is provided with a sampling pattern setting unit 340. The sampling pattern setting unit 340 determines the number of blades (the number of shots), a rotation angle, a blade width, and the number of measurement lines in the blade from the imaging condition input by the user and sets a sampling pattern using uneven intervals for interline intervals, in the blade type radial measurement.
The image reconstruction unit 330 performs image reconstruction on k-space data obtained in the blade type radial measurement by using an appropriate image reconstruction method. For example, image reconstruction using fast Fourier transformation or image reconstruction using iterative calculations such as compressed sensing is performed. In addition, in a case in which the RF reception coil is a multi-channel coil, combination of the data collected for each channel, parallel imaging reconstruction using a sensitivity distribution of each small coil, or the like is performed.
Further, the computer 30 may be provided with a function of evaluating the sampling pattern of the blade type radial measurement (evaluation unit 350). For example, the evaluation unit 350 calculates an evaluation value based on a predetermined evaluation criterion for the sampling pattern and generates a pre-image such as a simulation image in a case of using the sampling pattern.
The display control unit 320 displays a GUI for accepting a user setting for the imaging condition of the blade type radial measurement, particularly, for the sampling pattern, on the UI unit 40 and accepts the user setting through the GUI to pass the user setting to the imaging control unit 310. In a case in which the evaluation unit 350 is provided, an evaluation result by the evaluation unit 350 or the pre-image generated by a pre-image generation unit 353 is displayed on the UI unit 40.
The functions of the control unit and the calculation unit of the computer 30 described above are stored in advance in the memory as programs and are implemented by the CPU uploading and executing the program of each function. Note that a part of the functions of the computer 30 can also be executed by a programmable IC.
Next, the blade type radial measurement and the sampling pattern thereof will be described.
First, a pulse sequence of the blade type radial measurement and k-space data collected by the blade type radial measurement will be described.
The pulse sequence used for imaging of the blade type radial measurement is not limited, but a pulse sequence in which a plurality of nuclear magnetic resonance signals are continuously acquired after one excitation RF pulse (hereinafter, referred to as an excitation pulse) is applied is executed. Examples of such a pulse sequence include a sequence of a fast spin echo (FSE) method and a gradient echo-based pulse sequence.
By executing a similar pulse sequence by varying axes (combination of two axes) of a readout gradient magnetic field (Gr) and the phase-encoding gradient magnetic field (Gp), data with different blade angles on the same slice plane can be collected. Although
A k-space arrangement of data obtained by such a blade type radial measurement, that is, the sampling pattern, is shown in
The line intervals in the blade are determined by how an encoding step of the phase-encoding gradient magnetic field pulse 203 is changed, and in the general sampling pattern, the change in the encoding step indicated by a broken line in
The aspect of the uneven intervals is not limited, and an example of the uneven intervals in which data in a peripheral region is thinned out by making a center region of the blade dense is shown on a lower side of
In a case of iterative calculations involving image reconstruction with transformation into a sparse space such as compressed sensing, it is preferable that the data is random. Therefore, random sampling is suitable. In that case, as the patterns for the uneven intervals, instead of a pattern in which the line intervals gradually increase/decrease as shown in
In addition, in the sampling pattern, all the blades may have uneven intervals, or some blades may have even intervals and some other blades may have uneven intervals. Further, in a case in which a plurality of blades have uneven intervals, the patterns for the plurality of blades may be the same or different.
Among the sampling patterns exemplified above, which sampling pattern to use may be set in advance or may be set by accepting a user setting or user adjustment via the UI unit 40.
Hereinafter, a flow of the imaging control in the MRI apparatus having the above-described configuration will be described. As shown in
In a case in which there is no user designation (S3), the imaging control unit 310 controls the imaging unit 10 in the determined encoding step to perform the blade type radial measurement (S4). That is, echoes are measured, for example, in accordance with the FSE pulse sequence as shown in
The image reconstruction unit 330 reconstructs the image through calculations using the collected k-space data (S5). The technique of image reconstruction is the same as the reconstruction of a general radial scan, and the image is reconstructed through fast Fourier transformation (FFT) after gridding (data rearrangement onto grid points in the k-space) or iterative calculations. In that case, sensitivity correction of the reception coil or channel combination is performed as necessary. The optimal reconstruction method may be automatically selected according to the set sampling pattern.
In a case in which the adjustment or the designation by the user is accepted for the sampling pattern set by the sampling pattern setting unit 340 (S3), the sampling pattern setting unit 340 determines the sampling pattern again based on the condition designated by the user, and the process transitions to step S4.
According to the present embodiment, in the blade type radial measurement in which the k-space data is collected by varying the angle of the blade, by controlling the interline intervals in the blade to uneven intervals, it is possible to suppress the occurrence of artifacts that are likely to occur in a case of the even intervals. In addition, even in a case in which the k-space is thinned out and measured, the sampling pattern can be optimized, and the imaging can be accelerated while suppressing the decrease in image quality.
In the present embodiment, an adjustment function based on the evaluation of the sampling pattern and the evaluation result is added to the control unit of Embodiment 1.
As shown in
In addition, in the present embodiment, as the functions of the imaging control unit 310, a line interval adjustment unit 341 that adjusts line intervals in the blade may be provided. A function of adjusting the line intervals will be described below.
The evaluation of the sampling pattern is to evaluate how the sampling pattern affects the point spread function (PSF). Examples of the indicator for evaluating the influence on the PSF include a maximum value of sidelobes of the PSF in the sampling pattern and a sampling density.
The PSF in the sampling pattern can be created by gridding simulated data in which a sampling position is set to 1 and the other positions are set to 0, and then performing Fourier inverse transformation. From the sidelobe shape of the PSF, the maximum value of the sidelobes is calculated, and the value is used as the evaluation value.
The sidelobes and the maximum value thereof will be described in detail by using
As shown in (A) of
As shown in (D) of
In a case in which the sidelobes are clearly divided, the maximum value is large, which causes the occurrence of the artifacts. Therefore, by knowing the magnitude of the maximum value, it is possible to evaluate to what extent the artifacts are suppressed.
Regarding the sampling density, the sampling pattern is divided into regions, the density of the sampling point is calculated for each region, and it is evaluated whether the sampling density of each region is within an appropriate range. For example, predetermined threshold values (Th1, Th2) are set for the peripheral region and the center region, and then the evaluation value for the density is calculated based on a criterion such as (1) the peripheral region is not completely sparse (for example, the sampling density of the peripheral region is equal to or greater than the threshold value Th1) and (2) the center region is relatively dense (for example, the sampling density of the center region is the threshold value Th2).
For the pre-image, for example, an image in a case in which the blade type radial measurement is performed for each sampling pattern by using a digital phantom can be generated in advance. In a case in which the sampling pattern is set, the pre-image generated for the pattern is displayed. Since the images differ depending on the number of lines in the blade and the thinning-out rate, the pre-image is generated by varying these conditions for each sampling pattern, and then the pre-image is selected and displayed in conformity with the conditions for the set sampling pattern.
The display control unit 320 may display the images generated by the pre-image generation unit 353 as the sequentially changed pre-images in a thumbnail manner each time the sampling pattern, which is the processing target of the evaluation unit 350, and the conditions thereof are changed.
The user can determine whether the measurement may proceed with the set sampling pattern by referring to the above-described evaluation value and pre-image and can also change the pattern or the condition to perform the selection of the pattern, the adjustment of the number of lines or the acceleration factor, and the like again.
Next, a flow of processing of the blade type radial measurement of the present embodiment will be described with reference to
In the imaging control unit 310, in a case in which the blade type radial measurement is selected as the imaging method and the imaging condition is set (S1), the sampling pattern setting unit 340 sets the sampling pattern based on the imaging condition (S2). The evaluation unit 350 evaluates the sampling pattern set by the sampling pattern setting unit 340 (S21). In
As described above, for the evaluation, at least one of the generation of the pre-image or the calculation of the evaluation value based on the maximum value of the sidelobes or the sampling density is performed. The display control unit 320 displays the evaluation result of the evaluation unit 350 (S22). The user determines whether or not to make a new designation based on the displayed result. In a case in which there is no user designation (S3), the process transitions to step S4, the blade type radial measurement is executed with the set sampling pattern (S4), and the reconstructed image is displayed (S5).
Based on the result of the evaluation, the line intervals of the sampling pattern may be automatically adjusted on an apparatus side.
In this case, as shown in
The display control unit 320 may also display the evaluation result during adjustment or after adjustment, the sampling pattern, the condition thereof or the pre-image at that time, and the like on the UI unit 40 even in a case in which the line interval adjustment is automatically performed in the apparatus.
By automatically adjusting the line intervals in this way, the user can avoid complicated work such as resetting the imaging condition. In addition, by displaying the result, it is possible to confirm that the sampling pattern is a desired sampling pattern.
According to the present embodiment, in addition to the same effects as in Embodiment 1, by performing the evaluation of the sampling pattern and the sampling pattern setting reflecting the result thereof, it is possible to further optimize the sampling with higher accuracy, and it is possible to suppress artifacts while achieving the acceleration.
In addition, by presenting the set sampling pattern to the user and displaying the GUI for accepting the condition related to the setting of the sampling pattern, the degree of freedom of the user in the setting of the sampling pattern by the user can be enhanced, enabling an appropriate setting. In particular, by providing a parameter for the user to adjust the set uneven intervals as the imaging condition, the user can clearly understand the adjustment target.
In the block 1201 for the imaging condition, as conditions specialized in the blade type radial measurement, the number of blades (the number of shots), ETL (the number of lines in the blade), selection of even intervals (even) or uneven intervals (uneven) for sampling, and the like are displayed. A value set in advance as a default is displayed, and the user can change the value. In a case in which the uneven intervals for sampling are selected, for example, a new window 1204 may be opened, and the sampling pattern (Type B, Type C, and the like) with the uneven intervals may be selected, the sampling density may be set, or an uneven interval degree and the like may be set.
In a case in which the function of the evaluation unit 350 is provided, the user operates each button of the operation button block 1203 to perform the evaluation, such as the calculation of the evaluation value and the generation of the pre-image, and the display of the evaluation result. The pre-image is displayed on the image display block 1202. The evaluation value may be displayed together with the image.
In addition, the set sampling pattern and condition are decided on through an operation of a determination button by the user.
The GUI of
By providing such a user interface, the convenience of user-involved sampling pattern design is improved.
Number | Date | Country | Kind |
---|---|---|---|
2023-150407 | Sep 2023 | JP | national |