CN114325528B - Magnetic resonance imaging method and related equipment - Google Patents
Magnetic resonance imaging method and related equipment Download PDFInfo
- Publication number
- CN114325528B CN114325528B CN202111605754.6A CN202111605754A CN114325528B CN 114325528 B CN114325528 B CN 114325528B CN 202111605754 A CN202111605754 A CN 202111605754A CN 114325528 B CN114325528 B CN 114325528B
- Authority
- CN
- China
- Prior art keywords
- flip angle
- image
- constructing
- flip
- data corresponding
- 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
- 238000000034 method Methods 0.000 title claims abstract description 32
- 238000002595 magnetic resonance imaging Methods 0.000 title claims abstract description 28
- 238000005070 sampling Methods 0.000 claims abstract description 30
- 238000003384 imaging method Methods 0.000 claims description 11
- 238000004590 computer program Methods 0.000 claims description 7
- 239000011159 matrix material Substances 0.000 claims description 6
- 238000005259 measurement Methods 0.000 claims description 6
- 238000005065 mining Methods 0.000 claims description 6
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 230000001133 acceleration Effects 0.000 abstract description 6
- 238000010586 diagram Methods 0.000 description 3
- 230000009977 dual effect Effects 0.000 description 3
- 238000013507 mapping Methods 0.000 description 3
- 230000003287 optical effect Effects 0.000 description 3
- 230000000644 propagated effect Effects 0.000 description 2
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000000835 fiber Substances 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 239000013307 optical fiber Substances 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 230000008707 rearrangement Effects 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000003786 synthesis reaction Methods 0.000 description 1
Landscapes
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
The embodiment of the invention discloses a magnetic resonance imaging method and related equipment, wherein the method comprises the following steps: acquiring at least two scanning data of a sample to be imaged under at least two preset flip angles; constructing a first flip angle image based on scan data corresponding to the first flip angle; constructing a low-frequency phase difference image according to center sampling data in the scanning data corresponding to the first turning angle and center sampling data in the scanning data corresponding to the second turning angle; constructing a corrected second flip angle phase image based on the low frequency phase difference image and scan data corresponding to the second flip angle; and constructing a magnetic resonance image of the sample to be imaged according to the first flip angle image and the second flip angle image. According to the technical scheme provided by the embodiment of the invention, the dependency relationship among the multiple flip angles is adopted, so that the acquisition of the multiple flip angles can obtain higher acceleration multiple under the phase constraint condition, and the magnetic resonance imaging speed is improved.
Description
Technical Field
The embodiment of the invention relates to the technical field of imaging, in particular to a magnetic resonance imaging method and related equipment.
Background
Magnetic resonance is widely used in clinic because of its advantages of non-radiation, multi-contrast and availability of multiple azimuthal images. Quantitative Magnetic Resonance Imaging (MRI) imaging has good prospects in clinical applications as it provides absolute values of physical quantities, and is also important in multicenter studies and patient disease tracking follow-up. In order to obtain quantitative images, a sequence of multiple flip angles is often used in the magnetic resonance technology to obtain data, for example, more than 2 flip angles are used to obtain a T1 mapping quantitative map, and in combination with multi-echo acquisition, a T2 mapping quantitative map or a T2 mapping can be obtained. Because different flip angle acquisitions require separate scans, this doubles the scan time. Compressed sensing (Compressed Sensing: CS) is an effective method for accelerating the magnetic resonance scan speed because it uses sparsely sampled K-space trajectories, greatly reducing the number of sampling points. The traditional acceleration method adopts a method of independently accelerating each flip angle acquisition, does not use the dependency relationship between the flip angles, is easy to generate serious aliasing artifact on the image, and has poor acceleration effect.
Disclosure of Invention
The embodiment of the invention provides a magnetic resonance imaging method and related equipment, which avoid the problem that when magnetic resonance imaging is adopted, the imaging speed is low and serious aliasing artifact is generated on an image.
In a first aspect, an embodiment of the present invention provides a magnetic resonance imaging method, including:
acquiring at least two scanning data of a sample to be imaged under at least two preset flip angles, wherein the two preset flip angles comprise a first flip angle and a second flip angle with sequence, and the scanning data comprise center data obtained by full sampling of a preset center area on the sample to be imaged and edge data obtained by sparse sampling outside the preset center area on the sample to be imaged;
constructing a first flip angle image based on scan data corresponding to the first flip angle;
constructing a low-frequency phase difference image according to the central data corresponding to the first turning angle and the central data corresponding to the second turning angle;
constructing a second flip angle image based on the low frequency phase difference image and scan data corresponding to the second flip angle;
and constructing a magnetic resonance image of the sample to be imaged according to the first flip angle image and the second flip angle image.
Optionally, the step of acquiring at least two scan data of the sample to be imaged at least two preset flip angles includes:
and acquiring at least two scanning data of the sample to be imaged through the gradient echo sequences of the at least two preset flip angles.
Optionally, the step of constructing a first flip angle image based on the scan data corresponding to the first flip angle includes:
constructing a target compressed sensing model with sparse domain constraint;
and obtaining the first flip angle image based on the target compressed sensing model.
Optionally, the step of constructing a low-frequency phase difference image according to the center data corresponding to the first flip angle and the center data corresponding to the second flip angle includes:
calculating a first low-frequency phase based on complex image data of the center data of the first flip angle;
determining a second low frequency phase based on complex image data of the center data of the second flip angle;
and obtaining the low-frequency phase difference image according to the first low-frequency phase and the second low-frequency phase.
Optionally, the step of constructing a second flip angle image based on the low frequency phase difference image and scan data corresponding to the second flip angle includes:
constructing a phase image of a second flip angle based on the scanning data corresponding to the second flip angle, and obtaining a phase difference image of two flip angles according to the low-frequency phase difference image;
correcting the target compressed sensing model based on the low-frequency phase difference image to obtain a corrected target compressed sensing model;
and generating a second flip angle image according to the corrected target compressed sensing model and the scanning data corresponding to the second flip angle.
Optionally, the step of constructing a magnetic resonance image of the sample to be imaged according to the first flip angle image and the second flip angle image includes:
generating a quantitative image of the sample to be imaged according to the first flip angle image and the second flip angle image;
a magnetic resonance image of the sample to be imaged is constructed based on the quantitative image.
Optionally, the step of constructing a target compressed sensing model with sparse domain constraint, and obtaining the first flip angle image based on the target compressed sensing model includes:
the method comprises the following steps of:
obtaining the first flip angle image based on the target compressed sensing, wherein x is a fully sampled target K space and is used for reconstructing and generating a required image; f is an undersampled measurement matrix; y is undersampled K space data, and ψ is a sparse transform domain; g is the K space convolution kernel of SPIRIT, which is obtained by the central full mining area; λ1, λ2 are regularization factors.
Optionally, the step of correcting the target compressed sensing model based on the low-frequency phase difference image to obtain a corrected target compressed sensing model includes:
correcting the target compressed sensing model based on the low-frequency phase difference image, wherein the corrected target compressed sensing model is obtained by the following steps:
wherein x is a fully sampled target K space for reconstructing a desired image; f is an undersampled measurement matrix; y is undersampled K space data, ψ is a sparse transform domain, G is a K space convolution kernel of SPIRIT, and the K space convolution kernel is obtained by a central full mining area; λ1, λ2 are regularization factors;to correct the phaseAn image.
In a second aspect, embodiments of the present application provide a magnetic resonance imaging apparatus, including:
the data acquisition module is used for acquiring at least two scanning data of a sample to be imaged under at least two preset flip angles, wherein the two preset flip angles comprise a first flip angle and a second flip angle which are in sequence, and the scanning data comprise central data obtained by fully sampling a preset central area on the sample to be imaged and edge data obtained by sparse sampling outside the preset central area on the sample to be imaged;
the conversion module is used for constructing a first flip angle image based on the scanning data corresponding to the first flip angle;
the calculating module is used for constructing a low-frequency phase difference image according to the center data corresponding to the first turning angle and the center data corresponding to the second turning angle;
a correction module for constructing a second flip angle image based on the low frequency phase difference image and scan data corresponding to the second flip angle;
and the imaging module is used for constructing a magnetic resonance image of the sample to be imaged according to the first flip angle image and the second flip angle image.
In a third aspect, embodiments of the present application provide a computer readable storage medium having a computer program stored thereon, characterized in that: the computer program, when executed by a processor, implements a magnetic resonance imaging method as described above.
The embodiments of the above invention have the following advantages or benefits:
acquiring at least two scanning data of a sample to be imaged under at least two preset flip angles, wherein the two preset flip angles comprise a first flip angle and a second flip angle with sequence, and the scanning data comprise center data obtained by full sampling of a preset center area on the sample to be imaged and edge data obtained by sparse sampling outside the preset center area on the sample to be imaged; constructing a first flip angle image based on scan data corresponding to the first flip angle; constructing a low-frequency phase difference image according to the central data corresponding to the first turning angle and the central data corresponding to the second turning angle; constructing a second flip angle image based on the low frequency phase difference image and scan data corresponding to the second flip angle; and constructing a magnetic resonance image of the sample to be imaged according to the first flip angle image and the second flip angle image. By utilizing the characteristic of signal phase consistency among different flip angles, phase constraint is introduced into a compressed sensing model to improve the scanning speed of multi-flip angle imaging, the under-sampling multiple can be increased when other flip angles acquire data, the sampling time is further reduced, and higher acceleration multiple is obtained under the condition of phase constraint, so that the scanning speed of multi-flip angle imaging is improved.
Drawings
FIG. 1 is a flow chart of a magnetic resonance imaging method provided in an embodiment of the present application;
fig. 2 is a schematic structural diagram of a magnetic resonance imaging apparatus according to an embodiment of the present application;
FIG. 3 is a schematic structural diagram of a storage medium according to an embodiment of the present application;
fig. 4 is an application scenario diagram of magnetic resonance imaging according to an embodiment of the present application.
Detailed Description
The invention is described in further detail below with reference to the drawings and examples. It is to be understood that the specific embodiments described herein are merely illustrative of the invention and are not limiting thereof. It should be further noted that, for convenience of description, only some, but not all of the structures related to the present invention are shown in the drawings.
Fig. 1 is a flowchart of a magnetic resonance imaging method provided in an embodiment of the present application, where the magnetic resonance imaging method provided in the embodiment may be applicable to other scenes scanned at different flip angles with TE, such as in the acquisition of a radio frequency power emission map. The method may be performed by a magnetic resonance imaging apparatus, which may be implemented in software and/or hardware, typically integrated in a magnetic resonance imaging device.
As shown in fig. 1, an embodiment of the present invention provides a magnetic resonance imaging method, including:
s101, acquiring at least two scanning data of a sample to be imaged under at least two preset flip angles, wherein the two preset flip angles comprise a first flip angle and a second flip angle which are sequentially arranged, the scanning data comprise center data obtained by full sampling of a preset center area on the sample to be imaged and edge data obtained by sparse sampling outside the preset center area on the sample to be imaged, and the construction process is a reconstruction process;
in one possible embodiment, the step of acquiring at least two scan data of the sample to be imaged at least two preset flip angles includes:
and acquiring at least two scanning data of the sample to be imaged through the gradient echo sequences of the at least two preset flip angles.
For example, data are acquired using gradient echo sequences of at least two flip angles, respectively, each of which may be single echo or multi-echo. Taking 2 flip angles as an example, the corresponding flip angles are FA1 and FA2, respectively. The sampling track adopts a sparse sampling mode, and the central area of the sampling track is subjected to full sampling. A typical sampling trace is the sampling pattern of the session Disk.
S102, constructing a first flip angle image based on scanning data corresponding to the first flip angle;
in one possible embodiment, the step of constructing a first flip angle image based on scan data corresponding to the first flip angle includes:
constructing a target compressed sensing model with sparse domain constraint;
and obtaining the first flip angle image based on the target compressed sensing model.
In one possible implementation manner, the step of constructing a target compressed sensing model with sparse domain constraint, and obtaining the first flip angle image based on the target compressed sensing model includes:
the method comprises the following steps of:
obtaining the first flip angle image based on the target compressed sensing, wherein x is a fully sampled target K space and is used for reconstructing and generating a required image; f is an undersampled measurement matrix; y is undersampled K-space data, ψ is a sparse transform domain, such as TV transform, wavelet transform, etc.; g is the K space convolution kernel of SPIRIT, which is obtained by the central full mining area; λ1, λ2 are regularization factors.
Exemplary, the reconstruction model adopted for reconstructing the acquired data of the first flip angle FA1 is:
or alternatively
Model
The method can be used for reconstructing single-channel data, or reconstructing data of each channel one by one for multi-channel data, and then performing phase-based channel synthesis, so that the obtained phase image is Ph1.
S103, constructing a low-frequency phase difference image according to the center data corresponding to the first turning angle and the center data corresponding to the second turning angle;
for example, when calculating the low-frequency phase image PhL1 of FA1, phL1 may perform fourier transform on the data in the central full sampling area to obtain complex data, and obtain the phase of the data. The phase image reconstructed by the FA1 can also be obtained by low-pass filtering.
S104, constructing a second flip angle image based on the low-frequency phase difference image and scanning data corresponding to the second flip angle;
in one possible embodiment, the step of constructing a low-frequency phase difference image from the center data corresponding to the first flip angle and the center data corresponding to the second flip angle includes:
calculating a first low-frequency phase based on complex image data of the center data of the first flip angle;
determining a second low frequency phase based on complex image data of the center data of the second flip angle;
and obtaining the low-frequency phase difference image according to the first low-frequency phase and the second low-frequency phase.
In one possible embodiment, the step of constructing a second flip angle image based on the low frequency phase difference image and scan data corresponding to the second flip angle includes:
correcting the target compressed sensing model based on the low-frequency phase difference image to obtain a corrected target compressed sensing model;
and generating a second flip angle image according to the corrected target compressed sensing model and the scanning data corresponding to the second flip angle.
In one possible implementation manner, the step of correcting the target compressed sensing model based on the low-frequency phase difference image to obtain a corrected target compressed sensing model includes:
correcting the target compressed sensing model based on the low-frequency phase difference image, wherein the corrected target compressed sensing model is obtained by the following steps:
wherein x is a fully sampled target K space for reconstructing a desired image; f is an undersampled measurement matrix; y is undersampled K space data, ψ is a sparse transform domain, G is a K space convolution kernel of SPIRIT, and the K space convolution kernel is obtained by a central full mining area; λ1, λ2 are regularization factors;for the corrected phase image.
Illustratively, the low frequency phase image PhL2 of FA2 is calculated by the method of calculating PhL1 as described above, and the difference image PhC of PhL2 and PhL1 is calculated, as shown in FIG. 4, by correcting Ph1 with PhC to obtain corrected phasePh1 and corrected +.>Ph2:
phase is toReconstructing the data of the second flip angle FA2 as a constraint term to obtain:
low frequency phase image PhL2 of FA2.
Illustratively, the original model is decomposed into smooth terms and non-smooth terms by a first-order original dual methodAnd (5) solving. Wherein l 1 The norm term is not a smooth term. Assume e= |.| 1 Representation l 1 Norms, and E * As a conjugate function of E, v 1 ，v 2 Respectively two of the above models 1 The dual variable of the norm. The above model is converted into an original dual problem as follows:
s105, constructing a magnetic resonance image of the sample to be imaged according to the first flip angle image and the second flip angle image.
In a possible implementation manner, the step of constructing the magnetic resonance image of the sample to be imaged according to the first flip angle image and the second flip angle image includes:
generating a quantitative image of the sample to be imaged according to the first flip angle image and the second flip angle image;
a magnetic resonance image of the sample to be imaged is constructed based on the quantitative image.
Acquiring at least two scanning data of a sample to be imaged under at least two preset flip angles, wherein the two preset flip angles comprise a first flip angle and a second flip angle with sequence, and the scanning data comprise center data obtained by full sampling of a preset center area on the sample to be imaged and edge data obtained by sparse sampling outside the preset center area on the sample to be imaged; constructing a first flip angle image based on scan data corresponding to the first flip angle; constructing a low-frequency phase difference image according to the central data corresponding to the first turning angle and the central data corresponding to the second turning angle; constructing a second flip angle image based on the low frequency phase difference image and scan data corresponding to the second flip angle; and constructing a magnetic resonance image of the sample to be imaged according to the first flip angle image and the second flip angle image. By utilizing the characteristic of signal phase consistency among different flip angles, the scanning speed of multi-flip angle imaging is improved by introducing phase constraint in a compressed sensing model, the under-sampling multiple can be increased when other flip angles acquire data, the sampling time is further reduced, and higher acceleration multiple is obtained under the phase constraint condition, so that the scanning speed of multi-flip angle imaging is improved, the prior phase constraint item is increased, and better image quality can be ensured under the high acceleration multiple.
In one possible implementation, the present embodiments provide a magnetic resonance imaging apparatus, including:
the data acquisition module 201 is configured to acquire at least two scan data of a sample to be imaged under at least two preset flip angles, where the two preset flip angles include a first flip angle and a second flip angle with a sequence, and the scan data includes center data obtained by performing full sampling on a preset center area on the sample to be imaged and edge data obtained by performing sparse sampling outside the preset center area on the sample to be imaged;
a conversion module 202, configured to construct a first flip angle image based on scan data corresponding to the first flip angle;
a calculation module 203, configured to construct a low-frequency phase difference image according to the center data corresponding to the first flip angle and the center data corresponding to the second flip angle;
a correction module 204 for constructing a second flip angle image based on the low frequency phase difference image and scan data corresponding to the second flip angle;
an imaging module 205 is configured to construct a magnetic resonance image of the sample to be imaged according to the first flip angle image and the second flip angle image.
In one possible implementation, the present embodiments provide a computer readable storage medium 400 having a computer program 411 stored thereon, characterized in that: the computer program 411 when executed by a processor implements a magnetic resonance imaging method as described above.
The computer storage media of embodiments of the invention may take the form of any combination of one or more computer-readable media. The computer readable medium may be a computer readable signal medium or a computer readable storage medium. The computer readable storage medium may be, for example, but not limited to: an electrical, magnetic, optical, electromagnetic, infrared, or semiconductor system, apparatus, or device, or a combination of any of the foregoing. More specific examples (a non-exhaustive list) of the computer-readable storage medium would include the following: an electrical connection having one or more wires, a portable computer diskette, a hard disk, a Random Access Memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EPROM or flash memory), an optical fiber, a portable compact disc read-only memory (CD-ROM), an optical storage device, a magnetic storage device, or any suitable combination of the foregoing. In this document, a computer readable storage medium may be any tangible medium that can contain, or store a program for use by or in connection with an instruction execution system, apparatus, or device.
The computer readable signal medium may include a propagated data signal with computer readable program code embodied therein, either in baseband or as part of a carrier wave. Such a propagated data signal may take any of a variety of forms, including, but not limited to, electro-magnetic, optical, or any suitable combination of the foregoing. A computer readable signal medium may also be any computer readable medium that is not a computer readable storage medium and that can communicate, propagate, or transport a program for use by or in connection with an instruction execution system, apparatus, or device.
Program code embodied on a computer readable medium may be transmitted using any appropriate medium, including but not limited to: wireless, wire, fiber optic cable, RF, etc., or any suitable combination of the foregoing.
Computer program code for carrying out operations of the present invention may be written in any combination of one or more programming languages, including an object oriented programming language such as Java, smalltalk, C ++ and conventional procedural programming languages, such as the "C" programming language or similar programming languages. The program code may execute entirely on the user's computer, partly on the user's computer, as a stand-alone software package, partly on the user's computer and partly on a remote computer or entirely on the remote computer or server. In the case of a remote computer, the remote computer may be connected to the user's computer through any kind of network, including a Local Area Network (LAN) or a Wide Area Network (WAN), or may be connected to an external computer (for example, through the Internet using an Internet service provider).
It will be appreciated by those of ordinary skill in the art that the modules or steps of the invention described above may be implemented in a general purpose computing device, they may be centralized on a single computing device, or distributed over a network of computing devices, or they may alternatively be implemented in program code executable by a computer device, such that they are stored in a memory device and executed by the computing device, or they may be separately fabricated as individual integrated circuit modules, or multiple modules or steps within them may be fabricated as a single integrated circuit module. Thus, the present invention is not limited to any specific combination of hardware and software.
Note that the above is only a preferred embodiment of the present invention and the technical principle applied. It will be understood by those skilled in the art that the present invention is not limited to the particular embodiments described herein, but is capable of numerous obvious changes, rearrangements and substitutions without departing from the scope of the invention. Therefore, while the invention has been described in connection with the above embodiments, the invention is not limited to the embodiments, but may be embodied in many other equivalent forms without departing from the spirit or scope of the invention, which is set forth in the following claims.
Claims (8)
1. A method of magnetic resonance imaging comprising:
acquiring at least two scanning data of a sample to be imaged under at least two preset flip angles, wherein the two preset flip angles comprise a first flip angle and a second flip angle with sequence, and the scanning data comprise center data obtained by full sampling of a preset center area on the sample to be imaged and edge data obtained by sparse sampling outside the preset center area on the sample to be imaged;
constructing a first flip angle image based on scan data corresponding to the first flip angle;
constructing a low-frequency phase difference image according to the central data corresponding to the first turning angle and the central data corresponding to the second turning angle;
constructing a second flip angle image based on the low frequency phase difference image and scan data corresponding to the second flip angle;
constructing a magnetic resonance image of the sample to be imaged according to the first flip angle image and the second flip angle image;
the step of constructing a first flip angle image based on scan data corresponding to the first flip angle includes:
constructing a target compressed sensing model with sparse domain constraint;
obtaining the first flip angle image based on the target compressed sensing model;
a step of constructing a second flip angle image based on the low frequency phase difference image and scan data corresponding to the second flip angle, comprising:
constructing a phase image of a second flip angle based on the scanning data corresponding to the second flip angle, and obtaining a phase difference image of two flip angles according to the low-frequency phase difference image;
correcting the target compressed sensing model based on the low-frequency phase difference image to obtain a corrected target compressed sensing model;
and generating a second flip angle image according to the corrected target compressed sensing model and the scanning data corresponding to the second flip angle.
2. The method of magnetic resonance imaging according to claim 1, characterized in that the step of acquiring at least two scan data of the sample to be imaged at least two preset flip angles comprises:
and acquiring at least two scanning data of the sample to be imaged through the gradient echo sequences of the at least two preset flip angles.
3. The method of magnetic resonance imaging according to claim 1, wherein the step of constructing a low frequency phase difference image from the center data corresponding to the first flip angle and the center data corresponding to the second flip angle comprises:
calculating a first low-frequency phase based on complex image data of the center data of the first flip angle;
determining a second low frequency phase based on complex image data of the center data of the second flip angle;
and obtaining the low-frequency phase difference image according to the first low-frequency phase and the second low-frequency phase.
4. The magnetic resonance imaging method as set forth in claim 1, wherein the step of constructing a magnetic resonance image of the sample to be imaged from the first flip angle image and the second flip angle image comprises:
generating a quantitative image of the sample to be imaged according to the first flip angle image and the second flip angle image;
a magnetic resonance image of the sample to be imaged is constructed based on the quantitative image.
5. The method of magnetic resonance imaging according to claim 1, wherein the step of constructing a target compressed sensing model with sparse domain constraints, and deriving the first flip angle image based on the target compressed sensing model, comprises:
the method comprises the following steps of:
obtaining the first flip angle image based on the target compressed sensing, wherein x is a fully sampled target K space and is used for reconstructing and generating a required image; f is an undersampled measurement matrix; y is undersampled K space data, and ψ is a sparse transform domain; g is the K space convolution kernel of SPIRIT, which is obtained by the central full mining area; λ1, λ2 are regularization factors.
6. The method of claim 1, wherein the step of correcting the target compressed sensing model based on the low frequency phase difference image to obtain a corrected target compressed sensing model comprises:
correcting the target compressed sensing model based on the low-frequency phase difference image, wherein the corrected target compressed sensing model is obtained by the following steps:
wherein x is a fully sampled target K space for reconstructing a desired image; f is an undersampled measurement matrix; y is undersampled K space data, ψ is a sparse transform domain, G is a K space convolution kernel of SPIRIT, and the K space convolution kernel is obtained by a central full mining area; λ1, λ2, λ3 are regularization factors;for the corrected phase image.
7. An imaging apparatus based on the magnetic resonance imaging method as set forth in claim 1, comprising:
the data acquisition module is used for acquiring at least two scanning data of a sample to be imaged under at least two preset flip angles, wherein the two preset flip angles comprise a first flip angle and a second flip angle which are in sequence, and the scanning data comprise central data obtained by fully sampling a preset central area on the sample to be imaged and edge data obtained by sparse sampling outside the preset central area on the sample to be imaged;
the conversion module is used for constructing a first flip angle image based on the scanning data corresponding to the first flip angle;
the calculating module is used for constructing a low-frequency phase difference image according to the center data corresponding to the first turning angle and the center data corresponding to the second turning angle;
a correction module for constructing a second flip angle image based on the low frequency phase difference image and scan data corresponding to the second flip angle;
and the imaging module is used for constructing a magnetic resonance image of the sample to be imaged according to the first flip angle image and the second flip angle image.
8. A computer-readable storage medium having stored thereon a computer program, characterized by: the computer program, when executed by a processor, implements the magnetic resonance imaging method as claimed in any one of claims 1 to 6.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111605754.6A CN114325528B (en) | 2021-12-25 | 2021-12-25 | Magnetic resonance imaging method and related equipment |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111605754.6A CN114325528B (en) | 2021-12-25 | 2021-12-25 | Magnetic resonance imaging method and related equipment |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114325528A CN114325528A (en) | 2022-04-12 |
| CN114325528B true CN114325528B (en) | 2024-01-02 |
Family
ID=81013829
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111605754.6A Active CN114325528B (en) | 2021-12-25 | 2021-12-25 | Magnetic resonance imaging method and related equipment |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114325528B (en) |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103349550A (en) * | 2013-07-04 | 2013-10-16 | 华东师范大学 | Method and device for integrating magnetic resonance imaging scanning and compressive sensing reconstruction |
| CN104569880A (en) * | 2014-12-31 | 2015-04-29 | 中国科学院深圳先进技术研究院 | Magnetic resonance fast imaging method and system |
| CN108514415A (en) * | 2018-02-09 | 2018-09-11 | 奥泰医疗系统有限责任公司 | A kind of quick magnetic susceptibility-weighted imaging scanning sequence and method |
| CN108535675A (en) * | 2018-04-08 | 2018-09-14 | 朱高杰 | A kind of magnetic resonance multichannel method for reconstructing being in harmony certainly based on deep learning and data |
| CN111263896A (en) * | 2017-08-22 | 2020-06-09 | 皇家飞利浦有限公司 | Data-driven correction of phase-dependent artifacts in magnetic resonance imaging systems |
| CN112272535A (en) * | 2019-03-29 | 2021-01-26 | 上海联影医疗科技股份有限公司 | Method and system for magnetic resonance diffusion imaging |
| CN112370040A (en) * | 2020-11-13 | 2021-02-19 | 上海东软医疗科技有限公司 | Magnetic resonance imaging method, magnetic resonance imaging apparatus, storage medium, and electronic device |
| CN113391251A (en) * | 2020-03-12 | 2021-09-14 | 上海联影医疗科技股份有限公司 | Magnetic resonance image reconstruction method, device and equipment |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US10782378B2 (en) * | 2018-05-16 | 2020-09-22 | Siemens Healthcare Gmbh | Deep learning reconstruction of free breathing perfusion |
-
2021
- 2021-12-25 CN CN202111605754.6A patent/CN114325528B/en active Active
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103349550A (en) * | 2013-07-04 | 2013-10-16 | 华东师范大学 | Method and device for integrating magnetic resonance imaging scanning and compressive sensing reconstruction |
| CN104569880A (en) * | 2014-12-31 | 2015-04-29 | 中国科学院深圳先进技术研究院 | Magnetic resonance fast imaging method and system |
| CN111263896A (en) * | 2017-08-22 | 2020-06-09 | 皇家飞利浦有限公司 | Data-driven correction of phase-dependent artifacts in magnetic resonance imaging systems |
| CN108514415A (en) * | 2018-02-09 | 2018-09-11 | 奥泰医疗系统有限责任公司 | A kind of quick magnetic susceptibility-weighted imaging scanning sequence and method |
| CN108535675A (en) * | 2018-04-08 | 2018-09-14 | 朱高杰 | A kind of magnetic resonance multichannel method for reconstructing being in harmony certainly based on deep learning and data |
| CN112272535A (en) * | 2019-03-29 | 2021-01-26 | 上海联影医疗科技股份有限公司 | Method and system for magnetic resonance diffusion imaging |
| CN113391251A (en) * | 2020-03-12 | 2021-09-14 | 上海联影医疗科技股份有限公司 | Magnetic resonance image reconstruction method, device and equipment |
| CN112370040A (en) * | 2020-11-13 | 2021-02-19 | 上海东软医疗科技有限公司 | Magnetic resonance imaging method, magnetic resonance imaging apparatus, storage medium, and electronic device |
Non-Patent Citations (2)
| Title |
|---|
| 压缩感知技术在磁共振成像技术中的应用进展分析;肖顺; 储呈晨; 王源冰; 李斌;中国医疗设备;第36卷(第11期);4-9 * |
| 探讨用磁敏感定量成像和T2*加权成像方法对地中海贫血病患者脑铁含量的检测研究;郭红宇，王雁军，黄萍，郭华;Science Discovery;第4卷(第6期);405-409 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114325528A (en) | 2022-04-12 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Bustin et al. | High‐dimensionality undersampled patch‐based reconstruction (HD‐PROST) for accelerated multi‐contrast MRI | |
| US10915990B2 (en) | Systems and methods for denoising medical images with deep learning network | |
| US10671939B2 (en) | System, method and computer-accessible medium for learning an optimized variational network for medical image reconstruction | |
| Malavé et al. | Reconstruction of undersampled 3D non‐Cartesian image‐based navigators for coronary MRA using an unrolled deep learning model | |
| CN108375746B (en) | Phase reverse winding method and equipment | |
| US8879811B2 (en) | Alternating direction of multipliers method for parallel MRI reconstruction | |
| CN110133556B (en) | Magnetic resonance image processing method, device, equipment and storage medium | |
| CN109633502B (en) | Magnetic resonance rapid parameter imaging method and device | |
| US20130279786A1 (en) | Rapid parallel reconstruction for arbitrary k-space trajectories | |
| EP2904971A1 (en) | Magnetic resonance imaging device and image processing device | |
| US20140126796A1 (en) | Mri reconstruction with motion-dependent regularization | |
| CN103027681B (en) | For the system of the parallel MRI image obtaining of reconstruct | |
| US9366741B2 (en) | Medical image imaging method, medical diagnostic apparatus using the same, and recording medium therefor | |
| Schaetz et al. | Accelerated computing in magnetic resonance imaging: real‐time imaging using nonlinear inverse reconstruction | |
| WO2022236851A1 (en) | Method for reconstructing wrap-around field-of-view magnetic resonance image, computer device, and storage medium | |
| Yu et al. | Universal generative modeling in dual domains for dynamic MRI | |
| CN116309901A (en) | Artifact removal method and device, electronic equipment and storage medium | |
| CN113298901B (en) | Method for reconstructing magnetic resonance image in convoluted field of view, computer device and storage medium | |
| CN114325528B (en) | Magnetic resonance imaging method and related equipment | |
| CN107219481B (en) | A kind of general three-dimensional undersampling trace design method | |
| CN111537930B (en) | Magnetic resonance parameter imaging method and equipment based on gradient waveform adjustment | |
| CN113050009B (en) | Three-dimensional magnetic resonance rapid parameter imaging method and device | |
| US11119173B2 (en) | Dynamic imaging based on echo planar imaging sequence | |
| JP6809870B2 (en) | Magnetic resonance imaging device and signal processing method | |
| US20230309849A1 (en) | Method for generating a magnetic resonance image |
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 |