Summary of the invention
Recycle in order to solve existing magnetic resonance imaging local coil data in the process of rebuilding that to cause taking computational resource too much, the problem that computing velocity declines, the invention provides a kind of MR image reconstruction method.
A kind of MR image reconstruction method, comprising:
Obtain the hyperchannel K spatial data K of full sampling
1full;
In K space, to described data K
1fullproofread and correct the data K that is optimized
full, described optimization data K
full=N*K
1full, N is local coil-VTC coil correction parameter;
By described optimization data K
fullgo to image area and obtain image.
Preferably, obtain by the following method described local coil-VTC coil correction parameter N:
Obtain the data K that prescan imaging region local coil and VTC coil gather respectively
1, K
vtc;
Utilize described data K
1, K
vtctry to achieve respectively the prescan K spatial data K in a certain aspect
1d, K
vtcd;
Merge the data K of all passages
vtcd, obtain the uniform K spatial data of signal K
vctd0;
Choose convolution kernel and pass through mapping relations K
1dn=K
vtcd0calculate described local coil-VTC coil correction parameter N.
Preferably, the size of described convolution kernel is 2 × 2 × CH.
Preferably, directly obtain the hyperchannel K spatial data K of full sampling by full sampling
1full; Or obtain the multi-channel data of owing sampling, by algorithm, the described multi-channel data of owing sampling is filled up to the hyperchannel K spatial data K that obtains described full sampling
1full.
Preferably, described algorithm is GRAPPA algorithm or SPIRIT algorithm.
Preferably, by Fourier transform by described optimization data K
fullgo to image area and obtain image.
The present invention also provides a kind of MR image reconstruction device, comprising:
Data acquisition unit, for obtaining the hyperchannel K spatial data K of full sampling
1full;
K spatial data correcting unit, in K space to described data K
1fullproofread and correct the data K that is optimized
full, described optimization data K
full=N*K
1full, N is local coil-VTC coil correction parameter;
Image reconstruction unit, for by described optimization data K
fullgo to image area and obtain image.
Preferably, also comprise correction parameter acquiring unit, for obtaining in the following way local coil-VTC coil correction parameter N:
Obtain the data K that prescan imaging region local coil and VTC coil gather respectively
1, K
vtc;
Utilize described data K
1, K
vtctry to achieve respectively the prescan K spatial data K in a certain aspect
1d, K
vtcd;
Merge the data K of all passages
vtcd, obtain the uniform K spatial data of signal K
vctd0;
Choose convolution kernel and pass through mapping relations K
1dn=K
vtcd0calculate described local coil-VTC coil correction parameter N.
Preferably, described data acquisition unit is for directly obtaining the hyperchannel K spatial data K of full sampling by full sampling
1full; Or obtain the multi-channel data of owing sampling, by algorithm, the described multi-channel data of owing sampling is filled up to the hyperchannel K spatial data K that obtains described full sampling
1full.
Preferably, described image reconstruction unit by Fourier transform by described optimization data K
fullgo to image area and obtain image.
MR image reconstruction method provided by the invention, by K space to described data K
1fullproofread and correct the data K that is optimized
fullmethod, the process that prescan homogenization is proofreaied and correct is mentioned K space from image area, and passage is merged rapidly, in follow-up image area process of reconstruction, has reduced calculative data volume, and the speed of rebuilding is provided.The present invention also provides a kind of MR image reconstruction device.
Embodiment
A lot of details are set forth in the following description so that fully understand the present invention.But the present invention can implement to be much different from alternate manner described here, and those skilled in the art can do similar popularization without prejudice to intension of the present invention in the situation that, and therefore the present invention is not subject to the restriction of following public concrete enforcement.Secondly, the present invention utilizes schematic diagram to be described in detail, and in the time that the embodiment of the present invention is described in detail in detail, for ease of explanation, described schematic diagram is example, and it should not limit the scope of protection of the invention at this.
A kind of MR image reconstruction method, comprises the steps:
Obtain the hyperchannel K spatial data K of full sampling
1full; In K space, to described data K
1fullproofread and correct the data K that is optimized
full, described optimization data K
full=N*K
1full, N is local coil-VTC coil correction parameter; By described optimization data K
fullgo to image area and obtain image.
Explain the process of image reconstruction below in conjunction with accompanying drawing.As shown in Figure 1, comprising:
Step S10, obtains the hyperchannel K spatial data K of full sampling
1full.
Here can directly obtain the hyperchannel K spatial data K of full sampling by full sampling
1full; Or obtain the multi-channel data of owing sampling, by algorithm, the described multi-channel data of owing sampling is filled up to the hyperchannel K spatial data K that obtains described full sampling
1full.Those skilled in that art know or should know concrete acquisition method, and therefore not to repeat here.
Step S20, in K space, to described data K
1fullproofread and correct the data K that is optimized
full, described optimization data K
full=N*K
1full, N is local coil-VTC coil correction parameter.
Obtain by the following method described local coil-VTC coil correction parameter N:
Obtain the data K that prescan imaging region local coil and VTC coil gather respectively
1, K
vtc; Utilize described data K
1, K
vtctry to achieve respectively the prescan K spatial data K in a certain aspect
1d, K
vtcd.Particularly, by K
1, K
vtctransform to respectively image area I
m1, I
mvtc; Then according to aspect position, interpolation obtains I
m1d, I
mvtcd, be transformed back to respectively K space and obtain data K
1d, K
vtcd.As shown in Figure 2, wherein black lattice point represents pre-scan images numeric field data I to Interpolation Process
m1, I
mvtc, solid black lines represents current imaging aspect, the numerical value I in its solid line
m1d, I
mvtcddrawn by the black lattice point interpolation of being close to.
Merge the data K of all passages
vtcd, obtain the uniform K spatial data of signal K
vctd0.
Here merge K
vtcdmeaning be to obtain a single passage as object space, in this space, signal evenly and only have a passage, is set up mapping relations afterwards.
Choose convolution kernel and pass through mapping relations K
1d* N=K
vtcd0calculate described local coil-VTC coil correction parameter N.The algorithm (as GRAPPA, SPIRIT etc.) of rebuilding with reference to parallel acquisition, chooses convolution kernel, calculates the coefficient N in mapping relations A*N=B, and wherein A is from K
1d, B is from K
vtcd0; It determines that A is consistent with the method for B and parallel acquisition reconstruction, and the convolution kernel that slides, sets it as window, selects data set A and data set B.
First the size of determining convolution kernel, as shown in Figure 3, the size of convolution kernel is 2*2*CH, the quantity that CH is passage; Real point data representation K in figure
1dspace, ignore data representation K
vctd0space, wherein square expression data set A, arrow represents linear mapping coefficient N (merge coefficient), trigpoint represents data set B.It should be noted that and can have multiple convolution kernel here, but ultimate principle is identical, there is variation in position and the quantity of square data just, mapping relations are all also similar, the meaning of this mapping coefficient and parallel acquisition is very similar on mathematics, parallel acquisition is also broadly to shine upon doing in fact, is only interchannel and shines upon mutually, and the technical program is mapped to outside passage.
Step S30, by described optimization data K
fullgo to image area and obtain image.In one embodiment, by Fourier transform by described optimization data K
fullgo to image area and obtain image.
As shown in Figure 4, the figure on the left side is the image that available technology adopting conventional method is rebuild, the figure on the right is the image that utilizes technical scheme provided by the invention to rebuild, and can find out and utilize the quality of the image that the technical program obtains to reach the picture quality that conventional method is obtained; Utilize the technical program simultaneously, by K space to described data K
1fullproofread and correct the data K that is optimized
fullmethod, the process that prescan homogenization is proofreaied and correct is mentioned K space from image area, and passage is merged rapidly, in follow-up image area process of reconstruction, has reduced calculative data volume, has improved the speed of rebuilding.
As shown in Figure 5, the invention provides a kind of MR image reconstruction device, comprising:
Data acquisition unit 10, for obtaining the hyperchannel K spatial data K of full sampling
1full.
Described data acquisition unit 10 is for directly obtaining the hyperchannel K spatial data K of full sampling by full sampling
1full; Or obtain the multi-channel data of owing sampling, by algorithm, the described multi-channel data of owing sampling is filled up to the hyperchannel K spatial data K that obtains described full sampling
1full.
K spatial data correcting unit 20, in K space to described data K
1fullproofread and correct the data K that is optimized
full, described optimization data K
full=N*K
1full, N is local coil-VTC coil correction parameter.Described local coil-VTC coil correction parameter obtains by correction parameter acquiring unit 40, particularly, obtains in the following way:
Obtain the data K that prescan imaging region local coil and VTC coil gather respectively
1, K
vtc; Utilize described data K
1, K
vtctry to achieve respectively the prescan K spatial data K in a certain aspect
1d, K
vtcd; Merge the data K of all passages
vtcd, obtain the uniform K spatial data of signal K
vctd0; Choose convolution kernel and pass through mapping relations K
1d* N=K
vtcd0calculate described local coil-VTC coil correction parameter N.
Image reconstruction unit 30, for by described optimization data K
fullgo to image area and obtain image.One preferred embodiment in, described image reconstruction unit 30 by Fourier transform by described optimization data K
fullgo to image area and obtain image.
MR image reconstruction method provided by the invention, by K space to described data K
1fullproofread and correct the data K that is optimized
fullmethod, the process that prescan homogenization is proofreaied and correct is mentioned K space from image area, and passage is merged rapidly, in follow-up image area process of reconstruction, has reduced calculative data volume, has improved the speed of rebuilding.The present invention also provides a kind of MR image reconstruction device.
Although the present invention discloses as above, the present invention is not defined in this.Any those skilled in the art, without departing from the spirit and scope of the present invention, all can make various changes or modifications, and therefore protection scope of the present invention should be as the criterion with claim limited range.