CN107621617B - k space motion artifact correction device - Google Patents
k space motion artifact correction device Download PDFInfo
- Publication number
- CN107621617B CN107621617B CN201710690292.XA CN201710690292A CN107621617B CN 107621617 B CN107621617 B CN 107621617B CN 201710690292 A CN201710690292 A CN 201710690292A CN 107621617 B CN107621617 B CN 107621617B
- Authority
- CN
- China
- Prior art keywords
- data
- group
- space
- frame
- corrected
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Abstract
The invention discloses a k-space motion artifact correction device, which comprises: an acquisition unit 101 adapted to continuously acquire a plurality of frames of k-space data in an imaging region; a grouping unit 102, adapted to group the data acquired by one rf excitation in each frame into a group, where there is at least one data group in each frame; the correction unit 103 is adapted to perform data correction on data of a group to be corrected by using data of a standard group to sequentially complete data correction of each frame data group, wherein one data group in a previous frame is used as the standard group, and a data group in a next frame is used as the group to be corrected; the filling unit 104 is adapted to fill the corrected data group of each frame back into the k space of the corresponding frame, so as to obtain the corrected k space data of each frame. The device can relieve the influence of unexpected movement on imaging, reduce the influence on patients caused by imaging, acquire less data and accelerate imaging speed.
Description
The application is a division of a Chinese patent application filed in the Chinese patent office on 27/4/2013, with the application number of 201310152659.4 and the name of the invention being 'k-space motion artifact correction method and device'.
Technical Field
The present invention relates to k-space artifact correction, and in particular, to a method and an apparatus for k-space motion artifact correction.
Background
The magnetic resonance imaging equipment is an imaging equipment which can detect hydrogen atom signals in human organ tissues and form medical diagnosis images. The magnetic resonance equipment provides a high-intensity polarized magnetic field, a gradient field for space encoding and a radio frequency field for enabling hydrogen atoms in human organ tissues to generate resonance phenomena to generate resonance signals in an examination area, the generated resonance signals are received by a radio frequency receiving coil, and the resonance signals are finally converted into diagnostic images capable of distinguishing different tissues through signal processing of subsequent equipment. It has better imaging effect on soft tissue, so that it is more widely applied to diagnosis of brain, muscle, heart and cancer cells than other medical imaging devices such as CT, x-ray, etc.
Most of the magnetic resonance imaging needs a long time, and in the imaging process of some organs, such as when a person breathes, internal organs also move along with the diaphragm, and such as the heart moves periodically, so that when dynamic magnetic resonance imaging of such organ tissues needs to be performed, the motion artifact is difficult to eliminate.
In the prior art, when time-dependent magnetic resonance imaging is performed, the conventional method requires the patient to hold breath to complete the image acquisition. The long imaging time brings a great burden to the patient. There are also other imaging methods that do not require the patient to hold his breath and that can alleviate the patient's pain, such as the navigator technique, PROPELLER (peripheral rolled overhead parallel Lines) technique. The principle of the navigator technique is that during magnetic resonance imaging, the navigator echo is used to detect the movement of the imaged organ tissue, the time for acquiring the magnetic resonance signal is selected to be within the time meeting a certain condition, or the navigator echo information is used to perform necessary correction on the movement of the actually acquired magnetic resonance data relative to the organ tissue. The basic idea of the PROPELLER technique is to divide the k-space sampling into several parts to be performed sequentially, as shown in fig. 1, each part is called a k-space strip, when acquiring each k-space strip, the acquisition is completed soon, so that the examinee can be considered as stationary, while the motion only occurs between the acquired k-space strips, the direction and amplitude of the motion can be estimated from the data of the overlapped sampling region between the k-space strips, and after compensating the motion, a good quality image without motion artifact interference can be reconstructed.
However, navigation techniques require attention to information other than imaging, such as diaphragm motion. The PROPELLER technology repeatedly acquires in a k space, a central area is repeatedly acquired during each acquisition, the overlapping is serious, the imaging time is influenced, and the PROPELLER technology is only suitable for single-frame image imaging and is not applied to motion artifact correction in dynamic magnetic resonance imaging.
Disclosure of Invention
The problem to be solved by the invention is to provide a k-space motion artifact correction device, which corrects rigid motion artifact of an imaging object, relieves the influence of unexpected motion on imaging (such as the influence caused by respiratory motion in cardiac scanning), relieves the influence on a patient caused by imaging (such as the pain of breath holding required in cardiac movie), reduces acquired data and shortens imaging time in a magnetic resonance and time-dependent imaging method.
In order to solve the above problem, the present invention provides a k-space motion artifact correction device, including:
the acquisition unit 101 is suitable for continuously acquiring multiple frames of k-space data in an imaging area, wherein the k-space data in each frame is acquired by at least one radio frequency excitation, and the acquisition directions of adjacent frame data are different;
a grouping unit 102, adapted to group the data acquired by one rf excitation in each frame into a group, where there is at least one data group in each frame;
the correction unit 103 is adapted to perform data correction on data of a group to be corrected by using data of a standard group to sequentially complete data correction of each frame data group, wherein one data group in a previous frame is used as the standard group, and a data group in a next frame is used as the group to be corrected;
the filling unit 104 is adapted to fill the corrected data group of each frame back into the k space of the corresponding frame, so as to obtain the corrected k space data of each frame.
Preferably, the acquisition unit acquires multiple frames of k-space data, wherein the acquisition direction of odd frame data is kx direction, and the acquisition direction of even frame data is ky direction.
Preferably, the correcting unit 103 further comprises: a k-space center correction unit 1031 adapted to perform k-space center correction; a rotation correction unit 1032 connected to the k-space center correction unit 1031 and adapted to perform rotation correction; a translation correction unit 1033 connected to the rotation correction unit 1032 and adapted to perform translation correction; one or more of the orthotic units may be omitted.
Preferably, the correcting unit 103 further comprises: a k-space center correction unit 1031 adapted to perform k-space center correction; a rotation correction unit 1032 connected to the k-space center correction unit 1031 and adapted to perform rotation correction; a translational correction unit 1033, connected to the rotational correction unit 1032, adapted to perform translational correction.
Preferably, the k-space center correction unit 1031 is adapted to perform the following steps:
converting the k-space data DataZ1 of the group to be corrected into an image domain to obtain data ImageZ 1;
windowing k-space data DataZ1 of a group to be corrected, and converting the k-space data into an image domain to obtain data ImageZ1 Window;
taking the module value of the data ImageZ1 as the module value of the corrected Image data Image;
taking the phase difference between the data ImageZ1 and the data ImageZ1Window as the phase of the corrected Image data Image;
and transforming the Image data Image to k space to obtain corrected data DataZ1 Corr.
Preferably, the rotation angle correction unit 1032 is adapted to perform the following steps:
respectively taking module values of standard group data and group data to be corrected, wherein the sizes of the standard group data and the group data to be corrected in the ky direction and the kx direction are respectively K1 and K2, the step length of the standard group data in the ky direction is delta K1, and the step length of the group data to be corrected in the kx direction is delta K2;
rotating the data of the to-be-corrected group after the modulus value is taken by an angle alpha in a K space and interpolating to a public grid space to obtain a data group MCi, wherein the size of the public grid space in the ky direction is K22, the step length is delta K22, the size of the public grid space in the kx direction is K11, the step length is delta K11, the MCi represents the data group obtained by interpolating to the public grid space after the ith rotation, i is a positive integer, and alpha is any value between 0 degree and 360 degrees;
interpolating the standard group data after modulus value extraction to the public grid space to obtain a data group MB, and calculating the correlation between the MB and the MCi;
and finding a rotation angle alpha value corresponding to the maximum correlation value, and performing rotation angle correction by taking the rotation angle value as the data of the group to be corrected.
Preferably, said translational correction unit 1033 is adapted to perform the following steps:
interpolating standard group data and group data to be corrected to the same common grid space, wherein the standard group data and the group data to be corrected are K1 and K2 in the ky direction and the kx direction respectively, the step length of the standard group data in the ky direction is delta K1, the step length of the group data to be corrected in the kx direction is delta K2, the size of the common grid space in the ky direction is K22, the step length is delta K22, the size of the common grid space in the kx direction is K11, and the step length is delta K11;
transforming the product of the standard group data and the conjugate of the group data to be corrected in the common grid space to an image domain to obtain an image, and fitting to find the position of the maximum value of the image;
and taking the offset of the maximum position of the image and the central position of the image as the offset caused by translational motion, and taking the offset as the group to be corrected to carry out translational correction.
Preferably, K11 ═ K1, K22 ═ K2, Δ K11 ═ Δ K1, and Δ K22 ═ Δ K2.
Compared with the prior art, the k-space motion artifact correction device provided by the invention can be used for correcting the motion artifact of the acquired data by utilizing the similarity of the adjacent frame data. The device can relieve the influence of unexpected movement on imaging (such as the influence caused by respiratory movement in cardiac scanning), reduce the influence on patients caused by imaging (such as the pain of breath holding required in cardiac movies), acquire less data and accelerate the imaging speed.
Drawings
FIG. 1 is a schematic diagram of a prior art PROPELLER acquisition;
FIG. 2 is a flow chart of k-space motion artifact correction according to the present invention;
FIG. 3 is a schematic diagram of a multi-frame k-space data acquisition according to an embodiment of the present invention;
FIG. 4 is a schematic diagram illustrating interpolation of a data set in a common space according to an embodiment of the present invention;
FIG. 5 is a schematic structural diagram of a k-space motion artifact correction device according to the present invention;
fig. 6 is a schematic diagram of two or more frames of k-space data acquisition according to an embodiment of the present invention.
Detailed Description
In order to make the aforementioned objects, features and advantages of the present invention comprehensible, embodiments accompanied with figures are described in detail below.
As described in the background art, in the existing k-space motion artifact correction scheme, breath holding imposes a great burden on a patient, the navigation technology needs to pay attention to information except for imaging, such as diaphragm motion and heart periodic beating, the PROPELLER technology has serious data overlap in the acquisition process, and only focuses on single-frame image correction, and the application in multi-frame magnetic resonance imaging is absent. In the implementation process of the invention, the multi-frame image motion artifact correction can be carried out, when data are acquired, the acquisition directions of adjacent frame data are different, a large amount of overlapped data does not exist, the proximity of the adjacent frame data is used, the time loss of acquiring redundant repeated data is reduced under the condition of ensuring the accuracy of the motion artifact correction data, and the imaging speed is improved.
Fig. 2 is a schematic flow chart of the k-space motion artifact correction method of the present invention, which includes the following steps:
step S101, continuously collecting multiple frames of k-space data in an imaging area, wherein the k-space data in each frame is obtained by at least one radio frequency excitation collection, and the collection directions of adjacent frame data are different;
step S102, dividing data acquired by one-time radio frequency excitation in each frame into a group, wherein at least one data group is arranged in each frame;
step S103, taking a data group in the previous frame as a standard group and a data group in the next frame as a group to be corrected, and performing data correction on the data of the group to be corrected by using the data of the standard group to sequentially finish the data correction of each frame data group;
and step S104, filling each frame data set after data correction back into the k space of the corresponding frame to obtain k space data after correction of each frame.
The above-mentioned k-space motion artifact correction method is explained in detail with specific embodiments below.
Example one
Firstly, step S101 is executed, multiple frames of k-space data are continuously acquired in an imaging area, the k-space data in each frame is acquired by at least one radio frequency excitation, and the acquisition directions of adjacent frame data are different. It should be noted here that, compared with the PROPELLER acquisition method in the prior art, the PROPELLER technology does not acquire multi-frame k-space data, and has no concept of frames, and acquires multiple sets of k-space data for one frame of image during acquisition, and each set of data is acquired by continuously rotating in multiple acquisition directions in a stepping manner. The acquisition mode of the invention is to acquire multi-frame k space data, the acquisition directions of k space data of adjacent frames are different, and the acquisition directions of each frame do not need fixed rotation angle intervals, and the acquisition has the advantage that the data are not acquired repeatedly.
In a preferred embodiment of the first embodiment of the present invention, the k-space data acquisition mode is that the data acquisition directions of the odd frames are kx directions, the data acquisition directions of the even frames are ky directions, and the k-space data acquisition mode is as shown in fig. 3. Referring to fig. 3, only the distribution of k-space data of the first four frames is illustrated, and the omitted parts are similar to the previous parts, where RO is the readout direction (also referred to as the acquisition direction), PE is the phase encoding direction, kx and ky are the spatial dimensions in k-space, and t is the time dimension, and their notations are common criteria for those skilled in the art. In addition, the thick line and the thin line in each frame are k-space data acquired through one radio frequency excitation acquisition in the frame respectively. The data acquisition interval time of one-time radio frequency excitation is short, and rigid motion does not occur between the acquired data after one-time excitation.
The odd frame data collection direction is selected to be the kx direction, and the even frame data collection direction is selected to be the ky direction, so that the odd frame data and the even frame data can be interchanged in the actual implementation process, and the implementation of the embodiment is not affected. Furthermore, the acquisition direction of each frame data can be selected at will, and because the method utilizes the adjacent frames to correct the motion artifact, the acquisition directions of k space data of the adjacent frames are different.
Then, step S102 is executed to group the data acquired by one rf excitation in each frame, where there is at least one data group in each frame. As shown in fig. 3, taking the example of correcting the k-space data of the second frame from the k-space data of the first frame, each frame of data is divided into two groups, the thick line part data is one group, and the thin line part data is one group. The grouping here is intended to be such that when acquiring k-space for high resolution images, which corresponds to a larger number of phase encodings in k-space, it is therefore necessary to excite multiple times with radio frequency pulses in order to acquire sufficient data, the data acquired after each excitation being a set of data, since the data acquisition speed within a set is sufficiently fast so that it can be assumed that no rigid body motion occurs within the data within a set. That is, the data of the thick lines or the thin lines of the first frame or the second frame in fig. 3 is a group, and no rigid motion occurs in the data of the thick lines or the solid lines. When rigid body motion correction is carried out, a group of corrected data of the previous frame and any group of uncorrected data of the next frame are taken out for correction each time.
Then, step S103 is executed, a data group in the previous frame is used as a standard group, a data group in the next frame is used as a group to be corrected, data of the group to be corrected is corrected by using data of the standard group, and data correction of each frame data group is sequentially completed.
The data correction process is divided into k-space center correction, rotation angle correction and translation correction in sequence, and the correction method is similar to the technical idea of PROPELLER. The offset of the center of its data set from the center of K-space is first corrected. Taking the first frame as an example, the k-space raw data DataZ1 is transformed into an image domain to obtain data ImageZ1, and then data DataZ1 is added with a two-dimensional triangular Window to be programmed and transformed into the image domain to obtain data ImageZ1Window, wherein the image domain data ImageZ1 and ImageZ1Window are both complex numbers. Phase subtraction of ImageZ1 with ImageZ1Window and taking the modulus value of ImageZ1, let:
abs(Image)＝abs(ImageZ1)；
angle(Image)＝angle(ImageZ1-ImageZ1Window)；
image domain data Image can be obtained, then the Image is converted back to K space, K space data DataZ1Corr is obtained, namely the corrected data, and the K space data of the second frame and the subsequent frames are subjected to data correction by the same method, namely K space center correction of each frame is completed.
Then, the rotation angle correction is performed, please refer to fig. 4. As shown in fig. 4, the upper left diagram represents the first taken data in the frame, the first frame data acquisition direction is kx direction, the solid line represents the standard group data taken out from the first frame, the upper right diagram represents the second taken data in the frame, the second frame data acquisition direction is ky direction, and the dotted line represents any group of group data to be corrected taken out from the second frame. The standard group data and the group data to be corrected are respectively K1 and K2 in the ky direction and the kx direction, the step length of the standard group data in the ky direction is delta K1, and the step length of the group data to be corrected in the kx direction is delta K2. The lower part of fig. 4 is represented as a common grid space, where each black dot represents a data point obtained by interpolation, and the common grid space has a size of K22 in the ky direction and a step size of Δ K22, and has a size of K11 in the kx direction and a step size of Δ K11.
The rotation angle correction process includes firstly, obtaining a module value from data group of two frames, then interpolating the first frame data after the module value is obtained to a common grid space represented by black dots in fig. 4, and recording the interpolated data as MB. And then, rotating the second frame data by a series of angles alpha (alpha is from 0 degree to 360 degrees), and interpolating to the common grid space in fig. 4 to obtain data MCi (here, the data to be corrected of the second frame data is interpolated to the common grid space through the ith rotation, and i is a positive integer). Then, the Correlation (Correlation) between the data sets MB and MCi is calculated, and the best Correlation value obtained in each rotation series, i.e. the rotation angle α corresponding to the maximum Correlation value, is obtained, and α is considered to be the rotation angle change of the motion artifact in the acquisition process of the two frame data sets. And correcting the data of the second frame of data to be corrected according to the rotation angle alpha, namely finishing the rotation angle correction.
And finally, performing translational correction. The method comprises the steps of interpolating standard data group data and rotation-corrected data to be corrected into a common grid space together to obtain data groups DB and DC respectively, calculating the value of F (DB × conj (DC)), wherein an operator F represents Fourier transformation, conj represents a conjugate value, namely, the data groups DB and the DC are multiplied by the conjugate, then the data groups are transformed into an image domain to obtain an image, and then the obtained image is fitted to find the position of the maximum value. The offset between the position of the maximum value and the central position of the image is the offset caused by the translational motion during the acquisition of the two frames of data. And according to the calculated offset, the translation of the data group to be corrected of the second frame relative to the data group of the standard data of the first frame is removed, and the translation correction is finished.
Here, the data correction of one data group in the second frame is completed by using one data group in the first frame as a standard group, and the other data groups in the second frame can be subjected to data correction in the same manner (using one data group in the first frame as a standard group), and then each of the other data groups can be subjected to a similar data correction method.
For convenience of calculation, K11-K1, K22-K2, Δ K11- Δ K1, and Δ K22- Δ K2 may be used for the interpolation to select the common grid space, but they may be different. And finishing the motion offset correction of the group data to be corrected selected in the second frame.
In another embodiment of this embodiment, rigid motion artifacts in the remaining data sets of the second frame may be corrected using the second frame data set that has just been corrected. The standard group and the group to be corrected selected in the data correction process should be performed by selecting adjacent frames (or similar frames) as much as possible, so that the deformation (non-rigid motion) of the moving object is not large.
And finally, executing step S104, and filling each frame data set after data correction back into the k space of the corresponding frame to obtain k space data after correction of each frame.
It should be noted that, in the above steps, if the influence of a certain motion is not great, the motion can be ignored in the steps. If the rotational motion in the cardiac cine is not significantly affected, the step of rotational correction may be omitted and the subsequent translational correction may be performed directly.
In this embodiment, a k-space motion artifact correction device corresponding to the above-mentioned k-space motion artifact correction method is further provided, as shown in fig. 5, the k-space motion artifact correction device includes:
the acquisition unit 101 is suitable for continuously acquiring multiple frames of k-space data in an imaging area, wherein the k-space data in each frame is acquired by at least one radio frequency excitation, and the acquisition directions of adjacent frame data are different;
the grouping unit 102 is connected with the acquisition unit 101 and is suitable for grouping the data acquired by the primary radio frequency excitation in each frame, and at least one data group is arranged in each frame;
the correcting unit 103 is connected with the grouping unit 102 and is suitable for taking one data group in the previous frame as a standard group and taking the data group in the next frame as a group to be corrected, and performing data correction on the data of the group to be corrected by using the data of the standard group to sequentially finish the data correction of each frame data group;
and the filling unit 104 is connected with the correcting unit 103 and is adapted to fill the corrected data group of each frame back into the k space of the corresponding frame to obtain the corrected k space data of each frame.
In this embodiment, the correcting unit 103 further includes: a k-space center correction unit 1031 adapted to perform k-space center correction; a rotation correction unit 1032 connected to the k-space center correction unit 1031 and adapted to perform rotation correction; a translational correction unit 1033, connected to the rotational correction unit 1032, adapted to perform translational correction. One or more of the correction units 1031-1033 may be omitted.
For the specific implementation of the k-space motion artifact correction apparatus, reference may be made to the implementation of the k-space motion artifact correction method, which is not described herein again.
Example two
In this embodiment, the difference from the first embodiment is that the acquisition mode of k-space data of each frame in this embodiment is as shown in fig. 6, where only the acquired data of k-space of the first frame is shown, the data is acquired through three times of radio frequency excitation acquisition, and then the acquisition mode of data of each frame is the same. Therefore, in this embodiment, each frame of data is divided into three groups, each group being data acquired by one rf excitation, and each group of data corresponds to a solid line portion, a dotted line portion, and a dotted line portion as shown in fig. 6. In the implementation of the present practical example, the specific grouping situation is determined according to the data acquisition situation, the data acquired through one excitation acquisition is divided into one group, and the number of the groups is not limited to two or three groups.
The rigid body motion correction process is the same as the method in the first embodiment, a group of data in the previous frame is taken out to be used as the standard group of data, the data group in the next frame is used as the group to be corrected, the data of the group to be corrected is corrected by using the data of the standard group, and the data correction of each frame of data group is completed in sequence. Taking the first frame to correct the second frame as an example, any group of data is taken out from the first frame as standard group data (any group of solid line, dashed line or dashed line), any group of data is taken out from the second frame as group data to be corrected (any group of solid line, dashed line or dashed line), and then k-space center correction, rotation angle correction and translation correction are sequentially performed. If the first frame takes out the solid line group as the standard group data and the second frame takes out the solid line group as the group data to be corrected, after the second frame solid line group is corrected, the second frame dotted line group and the dot-and-dash line group can be corrected according to the first frame solid line group or the second frame solid line group. After each group of data of the second frame is corrected, correcting data of a third frame, a fourth frame and the like in the same manner.
The present invention is not limited to the above embodiments, and any person skilled in the art can make modifications and improvements without departing from the spirit and scope of the present invention.
Claims (7)
1. A k-space motion artifact correction device, comprising:
the acquisition unit (101) is suitable for continuously acquiring multiple frames of k-space data in an imaging area, the k-space data in each frame is acquired by at least one radio frequency excitation, and the acquisition directions of adjacent frame data are different;
a grouping unit (102) adapted to group data acquired by one radio frequency excitation in each frame into a group, at least one data group in each frame;
the correction unit (103) is suitable for taking a data group in a previous frame as a standard group and taking a data group in a next frame as a group to be corrected, and performing data correction on the data of the group to be corrected by using the data of the standard group to sequentially finish data correction of each frame data group;
and the filling unit (104) is suitable for filling the corrected data group of each frame back to the k space of the corresponding frame to obtain the corrected k space data of each frame.
2. The k-space motion artifact correction device as claimed in claim 1, wherein said acquisition unit acquires a plurality of frames of k-space data, wherein the odd frame data acquisition direction is kx direction and the even frame data acquisition direction is ky direction.
3. The k-space motion artifact correction device as claimed in claim 1, characterized in that the correction unit (103) further comprises: a k-space center correction unit (1031) adapted to perform k-space center correction; a rotation correction unit (1032) connected to the k-space center correction unit (1031) and adapted to perform rotation correction; a translation correction unit (1033) connected to the rotation correction unit (1032) and adapted to perform a translation correction.
4. The k-space motion artifact correction apparatus as claimed in claim 3, characterized in that the k-space center correction unit (1031) is adapted to perform the steps of:
converting the k-space data DataZ1 of the group to be corrected into an image domain to obtain data ImageZ 1;
windowing k-space data DataZ1 of a group to be corrected, and converting the k-space data into an image domain to obtain data ImageZ1 Window;
taking the module value of the data ImageZ1 as the module value of the corrected Image data Image;
taking the phase difference between the data ImageZ1 and the data ImageZ1Window as the phase of the corrected Image data Image;
and transforming the Image data Image to k space to obtain corrected data DataZ1 Corr.
5. The k-space motion artifact correction device as claimed in claim 3, characterized in that the rotation angle correction unit (1032) is adapted to perform the following steps:
respectively taking module values of standard group data and group data to be corrected, wherein the sizes of the standard group data and the group data to be corrected in the ky direction and the kx direction are respectively K1 and K2, the step length of the standard group data in the ky direction is delta K1, and the step length of the group data to be corrected in the kx direction is delta K2;
rotating the data of the to-be-corrected group after the modulus value is taken by an angle alpha in a K space and interpolating to a public grid space to obtain a data group MCi, wherein the size of the public grid space in the ky direction is K22, the step length is delta K22, the size of the public grid space in the kx direction is K11, the step length is delta K11, the MCi represents the data group obtained by interpolating to the public grid space after the ith rotation, i is a positive integer, and alpha is any value between 0 degree and 360 degrees;
interpolating the standard group data after modulus value extraction to the public grid space to obtain a data group MB, and calculating the correlation between the MB and the MCi;
and finding a rotation angle alpha value corresponding to the maximum correlation value, and performing rotation angle correction by taking the rotation angle value as the data of the group to be corrected.
6. The k-space motion artifact correction device according to claim 3, characterized in that the translational correction unit (1033) is adapted to perform the steps of:
interpolating standard group data and group data to be corrected to the same common grid space, wherein the standard group data and the group data to be corrected are K1 and K2 in the ky direction and the kx direction respectively, the step length of the standard group data in the ky direction is delta K1, the step length of the group data to be corrected in the kx direction is delta K2, the size of the common grid space in the ky direction is K22, the step length is delta K22, the size of the common grid space in the kx direction is K11, and the step length is delta K11;
transforming the product of the standard group data and the conjugate of the group data to be corrected in the common grid space to an image domain to obtain an image, and fitting to find the position of the maximum value of the image;
and taking the offset of the maximum position of the image and the central position of the image as the offset caused by translational motion, and taking the offset as the group to be corrected to carry out translational correction.
7. The K-space motion artifact correction device as claimed in claim 5 or 6, characterized in that said K11-K1, K22-K2, Δ K11- Δ K1, Δ K22- Δ K2.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710690292.XA CN107621617B (en) | 2013-04-27 | 2013-04-27 | k space motion artifact correction device |
Applications Claiming Priority (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710690292.XA CN107621617B (en) | 2013-04-27 | 2013-04-27 | k space motion artifact correction device |
| CN201310152659.4A CN104122521B (en) | 2013-04-27 | 2013-04-27 | K-space motion artifacts antidote and device |
Related Parent Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201310152659.4A Division CN104122521B (en) | 2013-04-27 | 2013-04-27 | K-space motion artifacts antidote and device |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107621617A CN107621617A (en) | 2018-01-23 |
| CN107621617B true CN107621617B (en) | 2019-12-20 |
Family
ID=51768010
Family Applications (2)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201310152659.4A Active CN104122521B (en) | 2013-04-27 | 2013-04-27 | K-space motion artifacts antidote and device |
| CN201710690292.XA Active CN107621617B (en) | 2013-04-27 | 2013-04-27 | k space motion artifact correction device |
Family Applications Before (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201310152659.4A Active CN104122521B (en) | 2013-04-27 | 2013-04-27 | K-space motion artifacts antidote and device |
Country Status (1)
| Country | Link |
|---|---|
| CN (2) | CN104122521B (en) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US10672151B1 (en) * | 2019-01-07 | 2020-06-02 | Uih America, Inc. | Systems and methods for magnetic resonance image reconstruction |
| CN110673070B (en) * | 2019-09-12 | 2022-03-01 | 上海联影医疗科技股份有限公司 | Training method of magnetic resonance signal correction network and magnetic resonance signal processing method |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP0974850A2 (en) * | 1998-07-21 | 2000-01-26 | Ge Yokogawa Medical Systems, Ltd. | Magnetic resonance imaging apparatus |
| CN1515912A (en) * | 2003-08-27 | 2004-07-28 | 中国人民解放军第一军医大学 | Automatic reverse correcting method of nuclear magnetic resonance image rigid translational aliasing artifact |
| CN101232845A (en) * | 2005-07-27 | 2008-07-30 | 株式会社日立医药 | Magnetic resonance imaging device |
| CN101329389A (en) * | 2007-06-21 | 2008-12-24 | 西门子公司 | Method for correction of movement artifacts |
| JP2009160408A (en) * | 2007-12-29 | 2009-07-23 | Siemens Ag | Distortion correction method and device in magnetic resonance imaging |
| CN101669038A (en) * | 2007-04-27 | 2010-03-10 | 皇家飞利浦电子股份有限公司 | Magnetic resonance device and method |
| CN101995561A (en) * | 2010-10-22 | 2011-03-30 | 南方医科大学 | Reconstruction method of periodically rotated overlapping parallel lines with enhanced reconstruction (PROPELLER) magnetic resonance data based on image domain overlapping |
| CN102096055A (en) * | 2010-12-14 | 2011-06-15 | 南方医科大学 | Rapid and accurate reconstructing method for non-uniform sampling data of magnetic resonance imaging |
Family Cites Families (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5539312A (en) * | 1995-02-02 | 1996-07-23 | Mayo Foundation For Medical Education And Research | Detection and measurement of motion during NMR imaging using orbital navigator echo signals |
| CN1327809C (en) * | 2005-03-28 | 2007-07-25 | 南方医科大学 | Method of eliminating motion false image using PROPELLER sampling mode in T1 weighted magnetic resaonance imaging |
| US7592809B1 (en) * | 2005-11-23 | 2009-09-22 | General Electric Company | Hybrid k-t method of dynamic imaging with improved spatiotemporal resolution |
| CN102362192A (en) * | 2009-03-25 | 2012-02-22 | 皇家飞利浦电子股份有限公司 | Motion detection and correction in magnetic resonance imaging for rigid, nonrigid, translational, rotational, and through-plane motion |
| CN102540116A (en) * | 2011-12-08 | 2012-07-04 | 中国科学院深圳先进技术研究院 | Magnetic resonance imaging method and system |
-
2013
- 2013-04-27 CN CN201310152659.4A patent/CN104122521B/en active Active
- 2013-04-27 CN CN201710690292.XA patent/CN107621617B/en active Active
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP0974850A2 (en) * | 1998-07-21 | 2000-01-26 | Ge Yokogawa Medical Systems, Ltd. | Magnetic resonance imaging apparatus |
| CN1515912A (en) * | 2003-08-27 | 2004-07-28 | 中国人民解放军第一军医大学 | Automatic reverse correcting method of nuclear magnetic resonance image rigid translational aliasing artifact |
| CN101232845A (en) * | 2005-07-27 | 2008-07-30 | 株式会社日立医药 | Magnetic resonance imaging device |
| CN101669038A (en) * | 2007-04-27 | 2010-03-10 | 皇家飞利浦电子股份有限公司 | Magnetic resonance device and method |
| CN101329389A (en) * | 2007-06-21 | 2008-12-24 | 西门子公司 | Method for correction of movement artifacts |
| JP2009160408A (en) * | 2007-12-29 | 2009-07-23 | Siemens Ag | Distortion correction method and device in magnetic resonance imaging |
| CN101995561A (en) * | 2010-10-22 | 2011-03-30 | 南方医科大学 | Reconstruction method of periodically rotated overlapping parallel lines with enhanced reconstruction (PROPELLER) magnetic resonance data based on image domain overlapping |
| CN102096055A (en) * | 2010-12-14 | 2011-06-15 | 南方医科大学 | Rapid and accurate reconstructing method for non-uniform sampling data of magnetic resonance imaging |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104122521A (en) | 2014-10-29 |
| CN107621617A (en) | 2018-01-23 |
| CN104122521B (en) | 2017-12-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP3815585B2 (en) | Magnetic resonance imaging system | |
| US11143730B2 (en) | System and method for parallel magnetic resonance imaging | |
| US9797974B2 (en) | Nonrigid motion correction in 3D using autofocusing with localized linear translations | |
| US20110274331A1 (en) | Magnetic resonance imaging method for achieving water-fat separation | |
| JP4197059B2 (en) | Nuclear magnetic resonance imaging system | |
| US11039757B2 (en) | Method and system for cardiac motion corrected MR exam using deformable registration | |
| WO2013130587A1 (en) | Method and apparatus for extended phase correction in phase sensitive magnetic resonance imaging | |
| JP6523730B2 (en) | Magnetic resonance imaging system | |
| JP2001161657A (en) | Magnetic resonance imaging system | |
| EP2026084A2 (en) | Single coil parallel imaging | |
| US11002815B2 (en) | System and method for reducing artifacts in echo planar magnetic resonance imaging | |
| JP3993799B2 (en) | Magnetic resonance imaging system | |
| JP4072879B2 (en) | Nuclear magnetic resonance imaging system | |
| CN107621617B (en) | k space motion artifact correction device | |
| EP3270306A1 (en) | Method for the acquisition and processing of measurement data by a combined magnetic resonance and x-ray device | |
| KR20150045332A (en) | Magnetic resonance imaging apparatus and imaging method for magnetic resonance image thereof | |
| JP2641486B2 (en) | NMR imaging equipment | |
| US20110306869A1 (en) | Apparatus and method for producing tomographic image | |
| JP3911605B2 (en) | Magnetic resonance imaging system | |
| US10234529B2 (en) | Method for producing magnet resonance tomography recordings of cyclic movement | |
| JP4675936B2 (en) | Nuclear magnetic resonance imaging system | |
| JP4678916B2 (en) | Magnetic resonance imaging system | |
| CN111462266B (en) | Image reconstruction method, device, CT equipment and CT system | |
| US20230288509A1 (en) | Combined acquisition and reordering scheme for reconstructing images with reduced motion artifacts | |
| CN108378851B (en) | Magnetic resonance diffusion weighted imaging method and device |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CP01 | Change in the name or title of a patent holder |
Address after: 201807 No. 2258 Chengbei Road, Jiading Industrial Zone, Jiading District, Shanghai. Patentee after: Shanghai Lianying Medical Technology Co., Ltd Address before: 201807 No. 2258 Chengbei Road, Jiading Industrial Zone, Jiading District, Shanghai. Patentee before: SHANGHAI UNITED IMAGING HEALTHCARE Co.,Ltd. |
|
| CP01 | Change in the name or title of a patent holder |