Magnetic resonance imaging K spatial movement artifact is corrected parallel acquisition method for reconstructing
[technical field]
The present invention relates to magnetic resonance imaging field, particularly a kind of magnetic resonance imaging K spatial movement artifact is corrected parallel acquisition method for reconstructing.
[background technology
]
In mr imaging technique, the speed of imaging is to weigh one of formation method very major criterion.The very key factor of restriction image taking speed is data acquisition and K space-filling.General data acquisition modes will be adopted full K spatial data, then just can rebuild and obtain image.Magnetic resonance parallel gathers reconstruction technique, is the mode of utilizing coil restructuring to merge, and the data of owing sampling are filled up, and utilization is filled up complete K spatial data and rebuild.Profit in such a way, can be according to demand, only gathers a part of K spatial data, needn't adopt completely whole K space, can greatly accelerate thus the speed of imaging.
One of more conventional method for parallel reconstruction is broad sense self calibration parallel acquisition (GRAPPA, GeneRalized Autocalibrating Partially Parallel Acquisitions).As shown in Figure 1, black real point is represented as the K spatial data of actual acquisition to traditional GRAPPA algorithm; White ignore is to owe the data that sampling need to be filled up; Grey real point represents the data of the appropriate full sampling in order to calculate coil merge coefficient.GRAPPA algorithm thinks, in Fig. 1, any one white ignore can be expressed as the linear superposition of black real point around, be equivalent to the data of multiple coils to merge, and coil merge coefficient n
i j(i coil, j position) can be determined by black real point matching grey real point.At coil merge coefficient n
i jafter determining, other white ignores can be according to the coil merge coefficient n trying to achieve
i jcoil is merged to the data of plugging a gap.
But, magnetic resonance imaging (MRI), because the time of its data acquisition is long, the motion of patient or some tissues (as heart) usually causes occurring in image artifact, picture quality severe exacerbation.Therefore, there is huge demand clinically in corrective exercise artifact effectively, is also that the study hotspot of medical magnetic resonance imaging is also one of technical barrier simultaneously always.
?james G Pipe proposed PROPALLER sample mode in 1999, had good eradicating efficacy for motion artifacts.As shown in Figure 2, wherein Fig. 1 a is single blade acquisition zone to PROPALLER acquisition mode; After 1b is multiple blades acquisition zone, the complete k space scrabbling up, due to total a part of region, each blade acquisition zone, can utilize this part region to correct the variation (as rotation, phase place variation etc.) between each blade.Such mode can be used for correcting the artifact causing owing to gathering object of which movement.
But this method, because k spatial data is overlapping comparatively serious, makes the acquisition time can be long; Traditional K space acquisition mode, a way of more efficiently accelerating picking rate is to utilize parallel acquisition (GRAPPA) method.But, due to the existence of calibration data, make the multiple accelerating beat certain discount.So especially, for the acquisition mode of Fig. 1, because each number in blade acquisition zone own is just little, if add calibration data, parallel advantage is just not obvious; Certainly, in the middle of also can adopting, two total parts of blade are carried out concurrent reconstruction, and the calculating of coil merge coefficient, but because the overlapping data volume of center section is generally less can cause the coefficient that calculates inaccurate.
therefore, need to inquire into a kind of new collecting method, to overcome technical matters above.
[summary of the invention]
The problem to be solved in the present invention is to provide a kind of magnetic resonance imaging K spatial movement artifact and corrects parallel acquisition method for reconstructing, accelerates traditional PROPALLER acquisition mode.
For addressing the above problem, technical solution of the present invention provides a kind of magnetic resonance imaging K spatial movement artifact to correct parallel acquisition method for reconstructing, comprises the steps:
The object to be imaged is placed in to magnetic resonance imaging system, obtains the magnetic resonance signal for filling K space;
Divide multiple blades to gather K space, it is a certain angle of rotation that interlobate K space geometry closes;
Each blade gathering is proofreaied and correct according to PROPELLER algorithm, calculated the data variation causing due to rigid motion between blade, calculate the sports calibration coefficient of each blade;
To each blade, choose respectively other suitable blades, other blades are arrived under the cartesian coordinate system of this blade to coil merge coefficient in the direction of calculating blade missing data according to sports calibration transformation of coefficient;
Fill up the data of blade disappearance according to the coefficient calculating;
All blades were carried out after same operation, and each blade data has filled up complete, fills up complete K space according to PROPELLER algorithm, obtains image through suitable after being converted into image area.
Preferably, described each leaf packet is containing a K spatial data subset, and this data subset is, and the mode of owing sampling obtains.
Preferably, two blades form a group, and coil merge coefficient also carries out in group the data filling process lacking in blade.
Preferably, two blades in a described group are mutually vertical.
Preferably, a blade in described two blades is parallel with first direction, and the data in it are to obtain according to the mode of owing sampling collection in the second direction vertical with first direction.
Preferably, a blade in described two blades is parallel with second direction, and the data in it are to obtain according to the mode of owing sampling collection at first direction.
Preferably, gather after first blade, exchanged phase encoding and frequency coding gradient, gathered another blade.
Preferably, described K space is Fourier transform to the conversion of image area.
Preferably, described second direction is Y direction, described in owe sampling collection and refer in X-direction according to interlace mode collection, in Y direction according to gathering every row mode.
Preferably, the sports calibration coefficient of described each blade refers to rotation and the translation of this blade.
Compared with prior art, technical solution of the present invention has the following advantages: profit in this way, can, in effectively removing motion artifacts, in the situation that not gathering calibration data, be carried out parallel acquisition reconstruction, and picking rate is promoted greatly.
[accompanying drawing explanation]
Fig. 1 is the schematic diagram of the GRAPPA algorithm of prior art;
Fig. 2 is the schematic diagram of prior art PROPALLER algorithm;
Fig. 3 is individual blade acquisition mode in the embodiment of the present invention one;
Fig. 4 is the blade associated group intention that two blade acquisition zones form;
Fig. 5 is the data structure schematic diagram in blade associated group;
Fig. 6 is the schematic diagram of rebuilding K space in the magnetic resonance imaging in the embodiment of the present invention one;
Fig. 7 is the schematic diagram of rebuilding K space in the magnetic resonance imaging of the embodiment of the present invention two.
Fig. 8 is the schematic flow sheet of the method for reconstructing in K space in magnetic resonance imaging of the present invention.
[embodiment]
For above-mentioned purpose of the present invention, feature and advantage can more be become apparent, below in conjunction with accompanying drawing, the specific embodiment of the present invention is described in detail.Set forth detail in the following description so that fully understand the present invention.But the present invention can be different from alternate manner described here and implements with multiple, and those skilled in the art can do similar popularization without prejudice to intension of the present invention in the situation that.Therefore the present invention is not subject to the restriction of following public embodiment.
The present invention comprises for the MR imaging apparatus of K space reconstruction: main magnet system, gradient system, radio system, computer system, other auxiliary equipment.Main magnet system is made up of permanent magnet or superconducting magnet, for generation of a uniform and stable main field, for by imaging object magnetization, produces macroscopic magnetization vector.Gradient system comprises: gradient coil, gradient amplifier, digital to analog converter, gradient controller, gradient cooling device.Gradient coil can produce the gradient magnetic of space linearity, makes the resonant frequency difference of imaging object at space diverse location, thereby the signal of space diverse location can be made a distinction.Radio system comprises: radio-frequency sending coil, radio frequency amplifier, RF receiving coil.Radio-frequency sending coil is used for emitting electromagnetic wave, imaging object is energized, thereby launches magnetic resonance signal.RF receiving coil is used for receiving the magnetic resonance signal that imaging object emits.Computer system is processed it after receiving magnetic resonance signal, completes K space reconstruction, and by Fourier transform, finally obtains the image of tested human body (or respective organization structure).The present invention is mainly concerned with the method for reconstructing in K space, and the process flow diagram of method for reconstructing refers to Fig. 8, and in described magnetic resonance imaging, the method for reconstructing in K space comprises the steps:
S1. the object to be imaged is placed in to magnetic resonance imaging system, obtains the magnetic resonance signal for filling K space, point multiple blades gather K space, and it is a certain angle of rotation (can be inconsistent) that interlobate K space geometry closes.
S2. each blade gathering is proofreaied and correct according to PROPELLER algorithm, calculated the data variation causing due to rigid motion between blade, calculate the sports calibration coefficient of each blade.
S3. to each blade L, choose respectively other suitable blades, other blades are transformed to according to mutual relationship under the cartesian coordinate system of this blade L, calculate in the direction of blade L missing data coil merge coefficient.
S4. fill up the data of blade L disappearance according to the coefficient calculating.
S5. all blades were carried out after same operation, each blade data has filled up complete, fills up complete K space according to PROPELLER algorithm, through the suitable image area that is converted into, obtains image.
Sports calibration coefficient in above-mentioned steps S2 refers to rotation and translation.Translation refers to be offset how many pixels, and rotation also refers to rotate how many angles, and so-called calibration factor is exactly by calculating, recording kinetic side-play amount, i.e. the pixel value of translation and the angle value of rotation.
In above-mentioned steps S5, the conversion of indication is Fourier transform.
With specific embodiment, the method for reconstructing in K space in above-mentioned magnetic resonance imaging is elaborated below.
embodiment mono-
The tissue (object) that needs magnetic resonance imaging is placed in to magnetic resonance imaging system, divide multiple blades to gather K space and obtain foliated K spatial data subset (as shown in Figure 3), each K spatial data subset obtains in the mode of owing sampling, thereby quickening data acquisition, wherein, Fig. 3 is the mode that individual blade gathers, and solid line represents image data, and dotted line represents to owe sampled data (part that need fill up by calculating).Fig. 4 is collection and the reconstruction grouping schematic diagram of this programme, and in Fig. 4, two square frames of solid line are one group (heavy line adds fine line), and two square frames of dotted line are another group.Foliated K spatial data subset of Regional Representative that solid line surrounds.The scheme of introducing herein, coil merge coefficient computer data filling, occurs in group.
For convenience of description, suppose and in group, only have two blades, and the pass of blade is orthogonal (as shown in Figure 4).Organize an interior blade parallel with first direction, another one blade is parallel with second direction, and first direction is vertical with second direction.In preferential embodiment, first direction refers to X(axle) direction, second direction refers to Y(axle) direction.
In Fig. 5, the point of grey is illustrated in the data that collect on directions X, and the point of black represents the data that collect in the Y direction, and the point of white is illustrated in the data that do not gather in X and Y-direction.
As can be seen from Figure 5, Grey Point blade is continuous on directions X, is interrupted in the Y direction; Black color dots blade is contrary.
Like this, can calculate the coil merge coefficient of black color dots blade on directions X according to Grey Numbers strong point, otherwise, can calculate Grey Point blade coil merge coefficient in the Y direction according to black data point.
The window (coil merge coefficient) of Fig. 6 a for distributing according to Fig. 5 black color dots data and constructing, Fig. 6 b is grey blade in Fig. 5, Fig. 6 c is black blade in Fig. 5.Utilize the window in Fig. 6 a, in Fig. 6 b, slide, calculate according to the relation shown in Fig. 1 that is similar to, can obtain the coil merge coefficient of black data point in Fig. 5; According to this coefficient, this window that slides in grey blade, can construct the disappearance data of (representing with white point).Fill up complete by grey blade like this.Black blade can fill up complete by similar way.
Treat that all blades fill up complete, can obtain complete K spatial data, it is carried out to suitable conversion, to image area, completed whole reconstruction flow process.
embodiment bis-
If two blades do not meet vertical relation, using the blade as coefficient calculations, the lattice point transforming under the cartesian coordinate system of target blade there will be situation as shown in Figure 7.
Fig. 7 a is the mutual relationship of two blades, and Fig. 7 b is the diagram that Grey Point blade is considered separately.In this case, because grey blade point still exists the point of continuous distribution in x direction, utilize the window shown in Fig. 6 a, this part data is extracted, carry out the calculating of merge coefficient.
For the relation of the mutual data blade of calibrating is vertical as far as possible, gathering after first blade, exchange phase encoding and frequency coding gradient, (as Fig. 3 b first gathers real fine rule, then gather real thick line, one group completes to gather another blade; Gather again empty fine rule, then empty thick line ...).
In above embodiment, calculate when grouping, only chosen a blade about the calculating of merge coefficient and described, but actual can be as required, choose multiple blades and carry out merge coefficient calculating.Now, can respectively each blade gridding be arrived to target blade coordinate system grid, then extract the data of target blade missing data direction; Also can be by multiple blade data together gridding to target blade coordinate system grid, then the unified data of extracting target blade missing data direction, calculate.After gridding, can, according to the situation of gridding, add a part of weight.
Although the present invention with preferred embodiment openly as above; but it is not for limiting the present invention; any those skilled in the art without departing from the spirit and scope of the present invention; can utilize method and the technology contents of above-mentioned announcement to make possible variation and modification to technical solution of the present invention; therefore; every content that does not depart from technical solution of the present invention; any simple modification, equivalent variations and the modification above embodiment done according to technical spirit of the present invention, all belong to the protection domain of technical solution of the present invention.