[go: up one dir, main page]

CN104166113A - Magnetic resonance image reconstruction method and device - Google Patents

Magnetic resonance image reconstruction method and device Download PDF

Info

Publication number
CN104166113A
CN104166113A CN201310753856.1A CN201310753856A CN104166113A CN 104166113 A CN104166113 A CN 104166113A CN 201310753856 A CN201310753856 A CN 201310753856A CN 104166113 A CN104166113 A CN 104166113A
Authority
CN
China
Prior art keywords
data
full
vtc
coil
image
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.)
Granted
Application number
CN201310753856.1A
Other languages
Chinese (zh)
Other versions
CN104166113B (en
Inventor
翟人宽
张卫国
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shanghai United Imaging Healthcare Co Ltd
Original Assignee
Shanghai United Imaging Healthcare Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Shanghai United Imaging Healthcare Co Ltd filed Critical Shanghai United Imaging Healthcare Co Ltd
Priority to CN201310753856.1A priority Critical patent/CN104166113B/en
Publication of CN104166113A publication Critical patent/CN104166113A/en
Application granted granted Critical
Publication of CN104166113B publication Critical patent/CN104166113B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Magnetic Resonance Imaging Apparatus (AREA)

Abstract

A magnetic resonance image reconstruction method provided by the invention comprises the following steps: fully-sampled multi-channel K-space data K1full is acquired; the data K1full is corrected in the K space to obtain optimized data Kfull, wherein the optimized data Kfull=N*K1full, and N is local coil-VTC coil correction parameter; and the optimized data Kfull is transformed to an image domain to acquire an image. By adopting the method in which the data K1full is corrected in the K space to obtain the optimized data Kfull, the process of pre-scanning homogenization correction is moved from the image domain to the K space, and channels are enabled to be merged quickly. Therefore, the amount of data needing calculation in the subsequent process of image domain reconstruction is reduced, and the speed of image reconstruction is improved. The invention further provides a magnetic resonance image reconstruction device.

Description

MR image reconstruction method and device
Technical field
The present invention relates to magnetic resonance imaging field, relate in particular to a kind of MR image reconstruction method and device.
Background technology
In mr imaging technique, because magnetic resonance adopts hyperchannel, data are gathered, therefore the inhomogeneous meeting of coil sensitivities causes the bright dark mistake of image, and clinical diagnosis meeting is affected.So the bright dark mistake of image is proofreaied and correct, seem particularly important in magnetic resonance.In conventional technology, relatively more outstanding is prescan homogenization alignment technique, utilize the region that imaging is wanted in scanning in advance, when scanning, utilize two kinds of coils of local coil and VTC coil to scan, because VTC coil is that signal is uniform, so utilize VTC and the bright dark difference of local coil to carry out image rectification.But because this correction occurs in image area, so in process of reconstruction, a large amount of local coil data will be recycled, finally just can obtain the data after merging, affect so on the one hand computing velocity.
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.
Brief description of the drawings
Fig. 1 is the process flow diagram of MR image reconstruction method provided by the invention;
Fig. 2 is that the present invention obtains the schematic diagram in local coil-VTC coil correction parameter N process;
Fig. 3 is the schematic diagram that the present invention obtains convolution algorithm in local coil-VTC coil correction parameter N process;
Fig. 4 is the comparison schematic diagram that adopts the image that the image that obtains of technical solution of the present invention and prior art obtain;
Fig. 5 is the structural representation of MR image reconstruction device provided by the invention.
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.

Claims (10)

1. a MR image reconstruction method, is characterized in that, 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.
2. MR image reconstruction method as claimed in claim 1, is characterized in that, obtains 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 1d* N=K vtcd0calculate described local coil-VTC coil correction parameter N.
3. MR image reconstruction method as claimed in claim 2, is characterized in that, the size of described convolution kernel is 2*2*CH.
4. MR image reconstruction method as claimed in claim 1, is characterized in that, directly obtains 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.
5. MR image reconstruction method as claimed in claim 4, is characterized in that, described algorithm is GRAPPA algorithm or SPIRIT algorithm.
6. MR image reconstruction method as claimed in claim 1, is characterized in that, by Fourier transform by described optimization data K fullgo to image area and obtain image.
7. a MR image reconstruction device, is characterized in that, 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 fulldescribed 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.
8. MR image reconstruction device as claimed in claim 7, is characterized in that, also comprises 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 1d* N=K vtcd0calculate described local coil-VTC coil correction parameter N.
9. MR image reconstruction device as claimed in claim 7, is characterized in that, 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.
10. MR image reconstruction device as claimed in claim 7, is characterized in that, described image reconstruction unit by Fourier transform by described optimization data K fullgo to image area and obtain image.
CN201310753856.1A 2013-12-31 2013-12-31 magnetic resonance image reconstruction method and device Active CN104166113B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310753856.1A CN104166113B (en) 2013-12-31 2013-12-31 magnetic resonance image reconstruction method and device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310753856.1A CN104166113B (en) 2013-12-31 2013-12-31 magnetic resonance image reconstruction method and device

Publications (2)

Publication Number Publication Date
CN104166113A true CN104166113A (en) 2014-11-26
CN104166113B CN104166113B (en) 2017-02-15

Family

ID=51910014

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310753856.1A Active CN104166113B (en) 2013-12-31 2013-12-31 magnetic resonance image reconstruction method and device

Country Status (1)

Country Link
CN (1) CN104166113B (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109345605A (en) * 2017-09-15 2019-02-15 上海联影医疗科技有限公司 System and method for reducing Nyquist Ghost artifacts
CN110226100A (en) * 2017-01-25 2019-09-10 上海联影医疗科技有限公司 System and method for magnetic resonance imaging

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6560353B1 (en) * 2000-03-21 2003-05-06 Magnetic Resonance Innovations, Inc. Method of MRI image reconstruction from partially acquired data in two or more dimensions using a multidimensional inverse transform technique
CN101216541A (en) * 2008-01-15 2008-07-09 新奥博为技术有限公司 Magnetic resonance image-forming system gradient field spherical harmonic coefficient obtaining method
CN103185878A (en) * 2011-12-27 2013-07-03 上海联影医疗科技有限公司 Magnetic resonance parallel image acquisition and image reconstruction method
CN103229069A (en) * 2010-12-02 2013-07-31 皇家飞利浦电子股份有限公司 MR imaging using a multi-point dixon technique

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6560353B1 (en) * 2000-03-21 2003-05-06 Magnetic Resonance Innovations, Inc. Method of MRI image reconstruction from partially acquired data in two or more dimensions using a multidimensional inverse transform technique
CN101216541A (en) * 2008-01-15 2008-07-09 新奥博为技术有限公司 Magnetic resonance image-forming system gradient field spherical harmonic coefficient obtaining method
CN103229069A (en) * 2010-12-02 2013-07-31 皇家飞利浦电子股份有限公司 MR imaging using a multi-point dixon technique
CN103185878A (en) * 2011-12-27 2013-07-03 上海联影医疗科技有限公司 Magnetic resonance parallel image acquisition and image reconstruction method

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
付正博 等: "《感应加热与节能—感应加热器(炉)的设计与应用》", 31 October 2008 *

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110226100A (en) * 2017-01-25 2019-09-10 上海联影医疗科技有限公司 System and method for magnetic resonance imaging
US11204408B2 (en) 2017-01-25 2021-12-21 Shanghai United Imaging Healthcare Co., Ltd. System and method for magnetic resonance imaging
CN110226100B (en) * 2017-01-25 2022-03-15 上海联影医疗科技股份有限公司 System and method for magnetic resonance imaging
CN109345605A (en) * 2017-09-15 2019-02-15 上海联影医疗科技有限公司 System and method for reducing Nyquist Ghost artifacts
CN109345605B (en) * 2017-09-15 2022-11-25 上海联影医疗科技股份有限公司 Systems and methods for reducing Nyquist ghost artifacts

Also Published As

Publication number Publication date
CN104166113B (en) 2017-02-15

Similar Documents

Publication Publication Date Title
US20220036516A1 (en) Training methods of a denoising model and image denoising methods and apparatuses
Sun et al. Histogram-based normalization technique on human brain magnetic resonance images from different acquisitions
CN102144923B (en) Method and magnetic resonance system to reduce distortions in diffusion imaging
CN104569882B (en) Phase correction method of fast spin echo of magnetic resonance imaging system and application thereof
US11842427B2 (en) Method and system of motion correction for magnetic resonance imaging
St-Jean et al. Non Local Spatial and Angular Matching: Enabling higher spatial resolution diffusion MRI datasets through adaptive denoising
CN102279375B (en) For reducing method and the magnetic resonance equipment of the distortion in diffusion imaging
JP7663212B2 (en) Optical fiber image calibration method for optical fiber image calibration device, optical fiber image calibration device, and computer-readable storage medium
RU2016142190A (en) CORRECTION OF EXTERNAL EPI ECHO SIGNALS
DE102014219782A1 (en) Operation of an imaging medical examination device with a plurality of subsystems
DE102014219784B4 (en) Operation of an imaging medical examination device with a plurality of subsystems
CN105869192A (en) Technology for reconstructing MRI fingerprint identification based on sliding window
Bastin et al. On the use of water phantom images to calibrate and correct eddy current induced artefacts in MR diffusion tensor imaging
RU2014119867A (en) MAGNETIC RESONANCE (MR) TOMOGRAPHY USING GENERAL INFORMATION FOR IMAGES WITH DIFFERENT CONTRAST
CN108305221B (en) Magnetic resonance parallel imaging method and device
CN104635188A (en) K-space reconstruction method and magnetic resonance imaging method
DE102014219778A1 (en) Operation of an imaging medical examination device with a plurality of subsystems
KR20130101476A (en) Method for determining a set of b1 field maps
CN104166112B (en) MR image reconstruction method and device
CN109345605B (en) Systems and methods for reducing Nyquist ghost artifacts
CN111190131A (en) Method and magnetic resonance system for magnetic resonance image reconstruction
CN104166113A (en) Magnetic resonance image reconstruction method and device
CN102521809A (en) Regularization correction method of magnetic resonance phase array coil image uniformity
JP2023530845A (en) Methods for Post-Processing Samples of WASAB1 Acquisition Sequences
JP2020175194A (en) Magnetic resonance imaging apparatus and magnetic resonance imaging method

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CP01 Change in the name or title of a patent holder

Address after: 201815 No. 1180 Xingxian Road, Jiading Industrial Zone, Jiading District, Shanghai

Patentee after: Shanghai Lianying Medical Technology Co., Ltd

Address before: 201815 No. 1180 Xingxian 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
CP02 Change in the address of a patent holder

Address after: 201807 2258 Chengbei Road, Jiading District, Shanghai

Patentee after: Shanghai Lianying Medical Technology Co.,Ltd.

Address before: 201815 No. 1180 Xingxian Road, Jiading Industrial Zone, Jiading District, Shanghai

Patentee before: Shanghai Lianying Medical Technology Co.,Ltd.

CP02 Change in the address of a patent holder