CN105467339A - Quick multilayer magnetic resonance imaging method and device - Google Patents
Quick multilayer magnetic resonance imaging method and device Download PDFInfo
- Publication number
- CN105467339A CN105467339A CN201511030684.0A CN201511030684A CN105467339A CN 105467339 A CN105467339 A CN 105467339A CN 201511030684 A CN201511030684 A CN 201511030684A CN 105467339 A CN105467339 A CN 105467339A
- Authority
- CN
- China
- Prior art keywords
- imaging
- rebuild
- reconstruction
- sampling
- data
- 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.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
- G01R33/48—NMR imaging systems
Landscapes
- Physics & Mathematics (AREA)
- High Energy & Nuclear Physics (AREA)
- Condensed Matter Physics & Semiconductors (AREA)
- General Physics & Mathematics (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
The invention discloses a quick multilayer magnetic resonance imaging method. The method comprises steps: calculation of a sampling template is carried out, and the sampling template which gives consideration to requirements of both parallel imaging and compressive sensing imaging is generated; data are acquired according to the sampling template; the sampling data are rebuilt, wherein the rebuilding comprises compressive sensing rebuilding and parallel imaging rebuilding. The invention also discloses a quick multilayer magnetic resonance imaging device. According to the multilayer imaging sampling mode and the corresponding rebuilding method, an original image of a high quality can be rebuilt when only few samples need to be acquired, the acquisition line number in k space is reduced, the scanning time is reduced, and quick multilayer imaging can be realized.
Description
Technical field
The application relates to magnetic resonance imaging arts, is specifically related to a kind of Quick multi-layer MR imaging method and device.
Background technology
In order to improve magnetic resonance image (MRI) picking rate, parallel imaging technique is widely used in magnetic resonance imaging.This technology mainly utilizes the spatial sensitivities difference of single receiving coil in phased-array coil to carry out encodes spatial information, is reduced to as necessary phase encoding step number, obtains sweep velocity faster.Traditional reconstruction algorithm mainly contains SENSE (sensitivityencoding), GRAPPA (Generalizedautocalibratingpartiallyparallelacquisitions) etc.
It similarly is the another kind of fast imaging method utilizing coil sensitivities information that multilayer is excited into.Multilayer is excited into picture and excites multi-layer image at every turn, obtains the image of aliasing, and the difference of recycling interlayer coil sensitivities, by the image of aliasing separately, obtains multi-layer image.
Summary of the invention
The application provides a kind of Quick multi-layer MR imaging method, comprising: sample template calculates, and generates the sample template taking into account parallel imaging and compressed sensing imaging requirements; According to sample template image data; Described sampled data is rebuild, comprises and sequentially carry out compressed sensing reconstruction and parallel imaging reconstruction.
Above-mentionedly to comprise according to sampling masterplate image data: the multilayer excitation pulse adopting many plane phase to be shifted carrys out imaging, and the excitation pulse of adjacent layer has the phase differential of PI, the adjacent layer image shift FOV/2 after imaging.
Above-mentioned sample template calculates, and generates the sample template taking into account parallel imaging and compressed sensing imaging requirements and comprises: carry out variable density sampling evenly falling on the basis of adopting, form final sample template.
Above-mentioned reconstruction described compressed sensing comprises: according to reconstruction formula
rebuild the K space data gathered, wherein x is the image rebuild, and y is the k-space data collected, F
pfor the fourier descriptor of lack sampling, Ψ is sparse transformation, and λ is bound term coefficient.
Above-mentioned parallel imaging of stating is rebuild and is comprised: the view data obtained after rebuilding described compressed sensing to K space through inverse fourier transform, is rebuild the data concurrent reconstruction algorithm after described conversion, and is separated by multi-layer image.
According to the second aspect of the application, a kind of Quick multi-layer MR imaging apparatus is provided, comprises: template generation module, calculate for sample template, generate the sample template taking into account parallel imaging and compressed sensing imaging requirements; Sampling module, for according to described sample template image data; Rebuild module, for rebuilding described sampled data, comprising and sequentially carrying out compressed sensing reconstruction and parallel imaging reconstruction.
The multilayer excitation pulse of above-mentioned sampling module also for adopting many plane phase to be shifted carrys out imaging, and the excitation pulse of adjacent layer has the phase differential of PI, the adjacent layer image shift FOV/2 after imaging.
Above-mentioned template generation module also for carrying out variable density sampling evenly falling on the basis of adopting, forms final sample template.
Above-mentioned reconstruction module is also for according to reconstruction formula
rebuild the K space data gathered, wherein x is the image rebuild, and y is the k-space data collected, F
pfor the fourier descriptor of lack sampling, Ψ is sparse transformation, and λ is bound term coefficient.
Above-mentioned reconstruction module also for the view data that obtains after rebuilding described compressed sensing through inverse fourier transform to K space, the data concurrent reconstruction algorithm after described conversion is rebuild, and multi-layer image is separated.
Owing to have employed above technical scheme, the beneficial effect that the application is possessed is:
In the embodiment of the application, due to compressed sensing lack sampling and corresponding method for reconstructing, only need gather a small amount of sample and can reconstruct original image by high-quality, reduce the gathering line number of k-space, reduce sweep time, can Quick multi-layer imaging.
Accompanying drawing explanation
Fig. 1 is the process flow diagram according to the application's method embodiment;
Fig. 2 is the sampling masterplate according to the application's method embodiment;
Fig. 3 is the structural representation according to the application's device embodiment.
Embodiment
By reference to the accompanying drawings the present invention is described in further detail below by embodiment.
Embodiment one:
Fig. 1 is the process flow diagram according to the application's method embodiment, comprising:
Step 102: sample template calculates, generates the sample template taking into account parallel imaging and compressed sensing imaging requirements.A kind of embodiment, sampling masterplate is calculated as follows: entirely adopt along frequency coding direction (kx) direction, first evenly owe to adopt (accelerating multiple=R1) according to parallel imaging mode along phase-encoding direction, namely often R1 is capable gathers a phase encoding line; Then the basis of adopting carries out variable density sampling (accelerating multiple=R2) evenly falling, i.e. the line number=R2 of the line number/should gather of actual acquisition, variable density sampling pattern meets the requirement of compressed sensing imaging theory.Through calculating, obtain the phase encoding line position that finally need gather, i.e. sample template, as shown in Figure 2, wherein solid dot is the data point gathered, and hollow dots is the data do not gathered in parallel imaging, and plus sige is the data do not gathered in compressed sensing imaging.
Step 104: according to sample template image data.Adopt the multilayer excitation pulse of POMP (the many plane phase displacements of PhaseoffsetMultiplanar) technology to carry out imaging, the excitation pulse of adjacent layer has the phase differential of PI, the adjacent layer image shift FOV/2 after imaging.
Step 106: rebuild sampled data, comprises and sequentially carries out compressed sensing reconstruction and parallel imaging reconstruction.A kind of implementation, this step comprises further:
A, directly to rebuild owing to adopt data by compressed sensing, reconstructing overlapping multilayer and exciting image.
A kind of embodiment, adopts compressed sensing method for reconstructing, and directly rebuild the K space data gathered, reconstruction formula is as follows:
wherein x is the image rebuild, and y is the k-space data collected, F
pfor the fourier descriptor of lack sampling, Ψ is sparse transformation, and λ is bound term coefficient, adopts Nonlinear conjugate gradient descent algorithm to solve.
B, utilize SENSE method by the image of overlap separately, form two width images.
A kind of embodiment, is designated as the view data that obtains in A to K space through inverse fourier transform
will
rebuild with concurrent reconstruction algorithm (SENSE), obtain final reconstructed results.POMP pulse is adopted during owing to exciting, more abundant to the coil sensitivities Information Pull of interlayer, the noise in picture centre region can be reduced.
Embodiment two:
Fig. 3 is the structural representation according to the application's device embodiment, comprising: template generation module, sampling module and reconstruction module.
Template generation module, calculates for sample template, generates the sample template taking into account parallel imaging and compressed sensing imaging requirements.A kind of embodiment, carries out variable density sampling evenly falling on the basis of adopting, forms final sample template.Entirely sample along frequency coding direction, along phase-encoding direction according to parallel imaging mode uniform subsampling, carry out stochastic sampling according to variable density sampling pattern to the data collected, variable density sampling pattern meets the requirement of compressed sensing imaging theory.
Sampling module, for according to described sample template image data.A kind of embodiment, the multilayer excitation pulse adopting many plane phase to be shifted carrys out imaging, and the excitation pulse of adjacent layer has the phase differential of PI, the adjacent layer image shift FOV/2 after imaging.
Rebuild module, for rebuilding sampled data, comprising and sequentially carrying out compressed sensing reconstruction and parallel imaging reconstruction.A kind of embodiment, according to reconstruction formula
rebuild the K space data gathered, wherein x is the image rebuild, and y is the k-space data collected, F
pfor the fourier descriptor of lack sampling, Ψ is sparse transformation, and λ is bound term coefficient.The view data obtained after rebuilding compressed sensing to K space through inverse fourier transform, is rebuild the data concurrent reconstruction algorithm after described conversion, and is separated by multi-layer image.
Compressive sensing theory utilizes the openness of signal, only need gather a small amount of sample and can reconstruct raw data by high-quality.The application utilizes this theory, can reconstruct original image from the k-space of owing to adopt, thus reduces the gathering line number of k-space, reduces sweep time, reaches the object of fast imaging.
Above content is in conjunction with concrete embodiment further description made for the present invention, can not assert that specific embodiment of the invention is confined to these explanations.For general technical staff of the technical field of the invention, without departing from the inventive concept of the premise, some simple deduction or replace can also be made.
Claims (10)
1. a Quick multi-layer MR imaging method, is characterized in that, comprising:
Sample template calculates, and generates the sample template taking into account parallel imaging and compressed sensing imaging requirements;
According to described sample template image data;
Described sampled data is rebuild, comprises and sequentially carry out compressed sensing reconstruction and parallel imaging reconstruction.
2. the method for claim 1, is characterized in that, describedly comprises according to described sampling masterplate image data:
The multilayer excitation pulse adopting many plane phase to be shifted carrys out imaging, and the excitation pulse of adjacent layer has the phase differential of PI, the adjacent layer image shift FOV/2 after imaging.
3. the method for claim 1, is characterized in that, described sample template calculates, and generates the sample template taking into account parallel imaging and compressed sensing imaging requirements, comprising:
Carry out variable density sampling evenly falling on the basis of adopting, form final sample template.
4. the method for claim 1, is characterized in that, described compressed sensing is rebuild and comprised:
According to reconstruction formula
rebuild the K space data gathered, wherein x is the image rebuild, and y is the k-space data collected, F
pfor the fourier descriptor of lack sampling, Ψ is sparse transformation, and λ is bound term coefficient.
5. the method for claim 1, is characterized in that, described parallel imaging is rebuild and comprised:
The view data obtained after rebuilding described compressed sensing to K space through inverse fourier transform, is rebuild the data concurrent reconstruction algorithm after described conversion, and is separated by multi-layer image.
6. a Quick multi-layer MR imaging apparatus, is characterized in that, comprising:
Template generation module, calculates for sample template, generates the sample template taking into account parallel imaging and compressed sensing imaging requirements;
Sampling module, for according to described sample template image data;
Rebuild module, for rebuilding described sampled data, comprising and sequentially carrying out compressed sensing reconstruction and parallel imaging reconstruction.
7. device as claimed in claim 6, it is characterized in that, the multilayer excitation pulse of described sampling module also for adopting many plane phase to be shifted carrys out imaging, and the excitation pulse of adjacent layer has the phase differential of PI, the adjacent layer image shift FOV/2 after imaging.
8. device as claimed in claim 6, is characterized in that, described template generation module also for carrying out variable density sampling evenly falling on the basis of adopting, forms final sample template.
9. device as claimed in claim 6, is characterized in that, described reconstruction module is also for according to reconstruction formula
rebuild the K space data gathered, wherein x is the image rebuild, and y is the k-space data collected, F
pfor the fourier descriptor of lack sampling, Ψ is sparse transformation, and λ is bound term coefficient.
10. device as claimed in claim 6, it is characterized in that, described reconstruction module also for the view data that obtains after rebuilding described compressed sensing through inverse fourier transform to K space, the data concurrent reconstruction algorithm after described conversion is rebuild, and multi-layer image is separated.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201511030684.0A CN105467339A (en) | 2015-12-31 | 2015-12-31 | Quick multilayer magnetic resonance imaging method and device |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201511030684.0A CN105467339A (en) | 2015-12-31 | 2015-12-31 | Quick multilayer magnetic resonance imaging method and device |
Publications (1)
Publication Number | Publication Date |
---|---|
CN105467339A true CN105467339A (en) | 2016-04-06 |
Family
ID=55605252
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201511030684.0A Pending CN105467339A (en) | 2015-12-31 | 2015-12-31 | Quick multilayer magnetic resonance imaging method and device |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105467339A (en) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106093814A (en) * | 2016-06-02 | 2016-11-09 | 浙江理工大学 | A kind of cardiac magnetic resonance imaging method based on multiple dimensioned low-rank model |
CN106339982A (en) * | 2016-08-24 | 2017-01-18 | 深圳先进技术研究院 | Fast magnetic resonance heart real-time cine imaging method and fast magnetic resonance heart real-time cine imaging system |
CN106597336A (en) * | 2016-11-29 | 2017-04-26 | 广东工业大学 | Scanning trajectory design method for MRI (magnetic resonance imaging) and device thereof |
CN106997034A (en) * | 2017-04-25 | 2017-08-01 | 清华大学 | Based on the magnetic resonance diffusion imaging method that reconstruction is integrated by example of Gauss model |
CN108564542A (en) * | 2018-04-04 | 2018-09-21 | 中国科学院长春光学精密机械与物理研究所 | A kind of control method of parallelly compressed perception imaging system |
CN108825205A (en) * | 2018-04-09 | 2018-11-16 | 中国石油大学(北京) | Downhole NMR spectroscopic signal compressed sensing acquisition method and device |
CN113534032A (en) * | 2020-04-14 | 2021-10-22 | 上海联影医疗科技股份有限公司 | Magnetic resonance imaging method and system |
US11327132B2 (en) | 2017-06-29 | 2022-05-10 | Shanghai United Imaging Healthcare Co., Ltd. | System and method for magnetic resonance imaging acceleration |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20120299590A1 (en) * | 2011-05-27 | 2012-11-29 | Riederer Stephen J | Method for Self-Calibrated Parallel Magnetic Resonance Image Reconstruction |
CN103064046A (en) * | 2012-12-25 | 2013-04-24 | 深圳先进技术研究院 | Image processing method based on sparse sampling magnetic resonance imaging |
CN103349550A (en) * | 2013-07-04 | 2013-10-16 | 华东师范大学 | Method and device for integrating magnetic resonance imaging scanning and compressive sensing reconstruction |
CN103505206A (en) * | 2012-06-18 | 2014-01-15 | 山东大学威海分校 | Fast and parallel dynamic MRI method based on compressive sensing technology |
CN104077791A (en) * | 2014-05-22 | 2014-10-01 | 南京信息工程大学 | Joint reconstruction method for multiple dynamic contrast enhancement nuclear magnetic resonance images |
-
2015
- 2015-12-31 CN CN201511030684.0A patent/CN105467339A/en active Pending
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20120299590A1 (en) * | 2011-05-27 | 2012-11-29 | Riederer Stephen J | Method for Self-Calibrated Parallel Magnetic Resonance Image Reconstruction |
CN103505206A (en) * | 2012-06-18 | 2014-01-15 | 山东大学威海分校 | Fast and parallel dynamic MRI method based on compressive sensing technology |
CN103064046A (en) * | 2012-12-25 | 2013-04-24 | 深圳先进技术研究院 | Image processing method based on sparse sampling magnetic resonance imaging |
CN103349550A (en) * | 2013-07-04 | 2013-10-16 | 华东师范大学 | Method and device for integrating magnetic resonance imaging scanning and compressive sensing reconstruction |
CN104077791A (en) * | 2014-05-22 | 2014-10-01 | 南京信息工程大学 | Joint reconstruction method for multiple dynamic contrast enhancement nuclear magnetic resonance images |
Non-Patent Citations (3)
Title |
---|
BO LIU ET AL.: "《SPARSESENSE:APPLICATION OF COMPRESSED SENSING IN PARALLEL MRI》", 《INTERNATIONAL CONFERENCE ON IEEE》 * |
FELIX A. BREUER ET AL.: "《Controlled aliasing in parallel imaging results in higher acceleration》", 《MAGNETIC RESONANCE IN MEDICINE》 * |
何珊: "《基于部分K空间数据的并行磁共振成像》", 《中国优秀硕士学位论文全文数据库》 * |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106093814A (en) * | 2016-06-02 | 2016-11-09 | 浙江理工大学 | A kind of cardiac magnetic resonance imaging method based on multiple dimensioned low-rank model |
CN106339982A (en) * | 2016-08-24 | 2017-01-18 | 深圳先进技术研究院 | Fast magnetic resonance heart real-time cine imaging method and fast magnetic resonance heart real-time cine imaging system |
CN106339982B (en) * | 2016-08-24 | 2019-12-24 | 深圳先进技术研究院 | Rapid magnetic resonance heart real-time film imaging method and system |
CN106597336A (en) * | 2016-11-29 | 2017-04-26 | 广东工业大学 | Scanning trajectory design method for MRI (magnetic resonance imaging) and device thereof |
CN106997034A (en) * | 2017-04-25 | 2017-08-01 | 清华大学 | Based on the magnetic resonance diffusion imaging method that reconstruction is integrated by example of Gauss model |
US11327132B2 (en) | 2017-06-29 | 2022-05-10 | Shanghai United Imaging Healthcare Co., Ltd. | System and method for magnetic resonance imaging acceleration |
CN108564542A (en) * | 2018-04-04 | 2018-09-21 | 中国科学院长春光学精密机械与物理研究所 | A kind of control method of parallelly compressed perception imaging system |
CN108564542B (en) * | 2018-04-04 | 2022-02-22 | 中国科学院长春光学精密机械与物理研究所 | Control method of parallel compressed sensing imaging system |
CN108825205A (en) * | 2018-04-09 | 2018-11-16 | 中国石油大学(北京) | Downhole NMR spectroscopic signal compressed sensing acquisition method and device |
CN108825205B (en) * | 2018-04-09 | 2020-09-22 | 中国石油大学(北京) | Method and device for compressed sensing acquisition of underground nuclear magnetic resonance spectrum signals |
CN113534032A (en) * | 2020-04-14 | 2021-10-22 | 上海联影医疗科技股份有限公司 | Magnetic resonance imaging method and system |
CN113534032B (en) * | 2020-04-14 | 2023-01-31 | 上海联影医疗科技股份有限公司 | Magnetic resonance imaging method and system |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105467339A (en) | Quick multilayer magnetic resonance imaging method and device | |
US9229082B2 (en) | Method for diffusion magnetic resonance imaging | |
CN107182216B (en) | A kind of rapid magnetic resonance imaging method and device based on depth convolutional neural networks | |
CN104068859B (en) | For determining method and the magnetic resonance equipment of multiple magnetic resonance image (MRI) | |
CN103185878B (en) | Magnetic resonance parallel image acquisition and image reconstruction method | |
CN1252488C (en) | Magnetic resonance imaging method with sub-sampling | |
CN107064845A (en) | One-dimensional division Fourier's parallel MR imaging method based on depth convolution net | |
CN104714199B (en) | A kind of MR imaging method and device | |
RU2013103063A (en) | DYNAMIC CONTRAST IMPROVED MR VISUALIZATION WITH RECONSTRUCTION OF COMPRESSED MEASUREMENT | |
CN105232045A (en) | Single-scanning quantitative magnetic resonance diffusion imaging method based on dual echoes | |
CN102540116A (en) | Magnetic resonance imaging method and system | |
CN103142228A (en) | Compressed sensing magnetic resonance fast imaging method | |
CN102937706B (en) | A kind of MRI parallel imaging method utilizing and magnetic resonance imager | |
CN102914754A (en) | Method to generate magnetic resonance exposures | |
CN101975935A (en) | Partial echo compressed sensing-based quick magnetic resonance imaging method | |
CN104749538A (en) | Phase processing method for parallel magnetic resonance imaging | |
CN105548927A (en) | Multi-excitation magnetic resonance diffusion imaging method based on multilayer simultaneous excitation | |
CN104597419A (en) | Method for correcting motion artifacts in combination of navigation echoes and compressed sensing | |
CN105738846A (en) | K-space data acquisition method and magnetic resonance imaging method | |
CN103027681A (en) | System used for reconstructing and parallelly obtaining mri image | |
CN104635188A (en) | K-space reconstruction method and magnetic resonance imaging method | |
CN102590773A (en) | Magnetic resonance imaging method and system | |
CN105931242B (en) | Dynamic nuclear magnetic resonance (DNMR) image rebuilding method based on dictionary learning and time gradient | |
CN105678822B (en) | A kind of three canonical magnetic resonance image reconstructing methods based on Split Bregman iteration | |
CN113971706A (en) | Rapid magnetic resonance intelligent 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 | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20160406 |