CN104714198B - The Susceptibility effect minimizing technology of layer directional compensation gradient is selected in adaptive change - Google Patents
The Susceptibility effect minimizing technology of layer directional compensation gradient is selected in adaptive change Download PDFInfo
- Publication number
- CN104714198B CN104714198B CN201310693631.1A CN201310693631A CN104714198B CN 104714198 B CN104714198 B CN 104714198B CN 201310693631 A CN201310693631 A CN 201310693631A CN 104714198 B CN104714198 B CN 104714198B
- Authority
- CN
- China
- Prior art keywords
- gradient
- magnetic field
- compensation gradient
- signal
- compensation
- 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.)
- Expired - Fee Related
Links
- 230000000694 effects Effects 0.000 title abstract description 25
- 230000003044 adaptive effect Effects 0.000 title abstract description 4
- 238000000034 method Methods 0.000 claims abstract description 26
- 230000008557 oxygen metabolism Effects 0.000 claims abstract description 7
- 238000012545 processing Methods 0.000 claims abstract description 5
- 238000007476 Maximum Likelihood Methods 0.000 claims description 2
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 abstract description 6
- 229910052760 oxygen Inorganic materials 0.000 abstract description 6
- 239000001301 oxygen Substances 0.000 abstract description 6
- 238000003384 imaging method Methods 0.000 abstract description 4
- 238000002592 echocardiography Methods 0.000 abstract description 3
- 230000010354 integration Effects 0.000 abstract description 2
- 238000005457 optimization Methods 0.000 abstract description 2
- 238000010586 diagram Methods 0.000 description 5
- 210000004556 brain Anatomy 0.000 description 3
- INGWEZCOABYORO-UHFFFAOYSA-N 2-(furan-2-yl)-7-methyl-1h-1,8-naphthyridin-4-one Chemical compound N=1C2=NC(C)=CC=C2C(O)=CC=1C1=CC=CO1 INGWEZCOABYORO-UHFFFAOYSA-N 0.000 description 2
- 108010002255 deoxyhemoglobin Proteins 0.000 description 2
- 235000013399 edible fruits Nutrition 0.000 description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 2
- 241001269238 Data Species 0.000 description 1
- 241000282414 Homo sapiens Species 0.000 description 1
- 239000002616 MRI contrast agent Substances 0.000 description 1
- 239000008280 blood Substances 0.000 description 1
- 210000004369 blood Anatomy 0.000 description 1
- 210000000476 body water Anatomy 0.000 description 1
- 230000000052 comparative effect Effects 0.000 description 1
- 238000005094 computer simulation Methods 0.000 description 1
- 238000002059 diagnostic imaging Methods 0.000 description 1
- 229940124642 endogenous agent Drugs 0.000 description 1
- 239000011259 mixed solution Substances 0.000 description 1
- 230000001575 pathological effect Effects 0.000 description 1
- 230000035479 physiological effects, processes and functions Effects 0.000 description 1
- 238000012805 post-processing Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Landscapes
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
The invention discloses one kind to be based on the asymmetric Rapid Imaging with Spin Echo technology (MASE) of more echoes, by selecting layer direction to apply the new method that the Susceptibility effect of adaptive change compensation gradient (Z-shimming Gradient) removes.Including data acquisition module, Δ B signal attenuation model processing module, removal three parts of Susceptibility effect module.The system well gets up the Z-Shimming gradient method of optimization and the decaying sinc model integration of MR signal caused by Susceptibility effect, adaptively apply Z-Shimming gradient according to the different Δ TE times, and efficiency finally is improved with pattern fitting method, the data of estimation Δ B are obtained in 24 seconds, and the data of the oxygen metabolism indexes such as data and estimation oxygen intake score OEF that can only obtain estimation Δ B in new method simultaneously with single pass using Z-MASE sequence, there is no geometry distributional difference between both data, so that Susceptibility effect removal processing is more convenient, accurately, robust.
Description
Technical field
It is specifically a kind of that layer directional compensation is selected based on adaptive change the invention belongs to magnetic resonance medical imaging technology field
The Susceptibility effect minimizing technology of gradient (Z-shimming Gradient).
Background technique
Biological tissue's magnetic susceptibility is not both function with its own, and physiology or pathological state are closely related, this machine
It makes and forms useful comparative information in MRI image.By BOLD contrast mechanism, it is known that deoxyhemoglobin is compared as endogenous
Agent causes signal decaying (Ogawa et al, Biophysical Journal1993 on T2* weighted image；Thulborn
Et al, Biochim Biophys Acta1982).But in tissue other than deoxyhemoglobin, cause there are also other
Near the reason of signal is decayed, such as interface between air tissue or the uneven distribution of main field itself, magnetic is formed
Sensitive artifact, leads to Magnetic field inhomogeneity, will affect the quantitative survey to oxygen metabolism indexes such as oxygen metabolism indexes such as oxygen intake score OEF
Amount.
Method proposed in this paper can remove influence caused by this Susceptibility effect, to obtain to tissue oxygen metabolism
Index, such as the more accurately estimation of oxygen intake score OEF.This method is in the asymmetric Rapid Imaging with Spin Echo skill of more echoes
On the basis of art (MASE), newly applies the choosing layer directional compensation gradient (ZShimming) changed with 180 degree shift time, selecting
MRI signal is encoded on layer gradient direction, using these information, Magnetic field inhomogeneity distribution map Δ B is estimated, finally, Δ B quilt
For removing the influence of the Susceptibility effect in original signal.
Summary of the invention
Present invention combination Z-Shimming method and the respective advantage of post-processing approach, with the asymmetric spin echo of more echoes
Based on fast imaging techniques (MASE), a kind of new method for removing Susceptibility effect is proposed.By being applied on selecting layer direction
Addend compensates gradient (Z-Shimming Gradient), followed by collected signal, in conjunction with its signal according to sinc letter
Several changing rules obtains accurate Magnetic field inhomogeneity distribution map Δ B, and finally removes macroscopic magnetic field and unevenly refer to oxygen metabolism
Mark the influence of estimation.
The present invention can be designed by magnetic resonance sequences, apply several compensation gradients on selecting layer direction, and after utilization
Reason method combines its signal according to the changing rule of sinc function, completes following task:
Layer is selected to return due to incomplete 1. overcoming using multiple and different choosing layer rephasing gradients is applied on selecting layer direction
Gradient echo image caused by poly- mechanism is selecting signal decaying namely Z-Shimming method on layer direction；
2. the relationship between MR signal decaying caused by couple macroscopical non-uniform magnetic field Δ B is analyzed, it is known that difference compensation
MR signal under gradient is compensation gradient dGzjWith Δ TEiSinc function；
3. utilizing limited quantity under the premise of known compensation gradient and its application time by above-mentioned sinc equation
The MR signal that gradient obtains is compensated, the profile Δ B of macroscopical non-uniform magnetic field is estimated.
To achieve the above object, the present invention takes following technical scheme: based on selecting layer directional compensation gradient (Z-shimming
Gradient Susceptibility effect minimizing technology), comprising the following three steps:
1. compensating ladder for different Z-Shimming is adaptively applied according to the different Δ TE time on selecting layer direction
It spends square (gradient moment)；
2. establishing corresponding macroscopical non-uniform magnetic field variable signal sinc model, and accurate using greatest hope estimation method
Estimate macroscopical non-uniform magnetic field profile Δ B；
3. further removing magnetic susceptibility in original MR signal using macroscopic view non-uniform magnetic field profile Δ B obtained in step 2
The influence of artifact, the oxygen intake score OEF result for estimating MASE reach best.
The invention is characterized in that above technical scheme, pseudo- by the Z-Shimming gradient method of optimization and magnetic susceptibility well
The decaying sinc model integration of MR signal caused by shadow gets up, and adaptively applies Z-Shimming gradient according to the different Δ TE times,
And efficiency finally is improved with pattern fitting method, the data of estimation Δ B are obtained in 24 seconds, and Z- is utilized in new method
MASE sequence can only obtain the oxygen metabolism indexes such as data and the estimation oxygen intake score OEF of estimation Δ B simultaneously with single pass
Data, there is no geometry distributional difference between both data, so that Susceptibility effect removal processing is more convenient, accurately.
The present invention is in computer simulation data, Susceptibility effect hydraulic model trial and 10 healthy volunteer's experimental datas
Did test；The bottom of cylindrical body water mould is tightly attached to using the mixed solution of MRI contrast agent and water as Susceptibility effect signal source
Portion is selected far from magnetic susceptibility source suitable distance, and at axis location image, the asymmetric spin echo image center the MR of Δ TE=62 occurs one
Signal attenuation region caused by a Susceptibility effect.The spin echo image of identical imaging position restores ladder since 180 degree is overturn
The presence of degree, is not influenced by Susceptibility effect, is existed without signal attenuation region.Using it is proposed that new method estimation
Water mould Δ B Distribution of Magnetic Field shows that reach maximum value be about 0.08ppm to Δ B at image center location, after correcting Susceptibility effect
Signal decaying is resumed in the asymmetric spin echo image of MR.
By the analysis of 10 volunteer's regions of interest data the result shows that, before removing Susceptibility effect, 4.38 ±
0.61% (full brain)；But after removing Susceptibility effect, DBV is reduced to 3.59 ± 0.58% (full brains), and P value is <
0.00001 (full brain) shows before and after removing Susceptibility effect that there are significant differences by venous blood volume DBV.The knot of this and forefathers
Fruit is consistent (An et al, MRM2002).
Detailed description of the invention
Fig. 1 is three module diagrams of this system.
Fig. 2 is three echo Z-MASE sequence Z-Shimming gradient allocation plan schematic diagrames.
Fig. 3 is analogue data Δ B estimation result schematic diagram.
Fig. 4 is Susceptibility effect hydraulic model trial result schematic diagram.
Fig. 5 is that normal health volunteer removes oxygen intake score OEF, DBV, R2 and R2 ' estimation knot before and after Susceptibility effect
Fruit compares schematic diagram.
Specific embodiment
The present invention is described in further detail below in conjunction with the accompanying drawings, and the present invention includes three main modulars altogether, such as Fig. 1 institute
Show, specific as follows:
1. data acquisition module: during actual scanning, at the MASE sequence maximum Δ TE moment, adaptively obtaining Z-
Shimming gradient Gzmax, while keeping Δ TE constant, by the way that maximum Z-Shimming gradient square is put down in positive negative direction
It distributes, according to [0, -1/4, -2/4, -3/4, -4/4,1/4,2/4,3/4] × GzmaxRule apply, and then it is average to apply this
Compensation gradient after distribution, to obtain multiple magnetic resonance signals compensated by Z-Shimming gradient.
2. Δ B signal attenuation model processing module: compensation gradient and original macroscopic view on layer direction are selected in any one application
Caused MR signal decaying meets sinc function model under non-uniform magnetic field Δ B collective effect, by the model, in known compensation
Gradient and under the premise of the time, the MR signal obtained using the compensation gradient of limited quantity, can pass through maximum likelihood parameter
Estimation method estimates the profile Δ B of macroscopical non-uniform magnetic field.
3. removing Susceptibility effect module: after estimating macroscopical non-uniform magnetic field Δ B, can according to following equation,
The influence of Susceptibility effect Δ B is removed in original signal.
Claims (4)
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201310693631.1A CN104714198B (en) | 2013-12-17 | 2013-12-17 | The Susceptibility effect minimizing technology of layer directional compensation gradient is selected in adaptive change |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201310693631.1A CN104714198B (en) | 2013-12-17 | 2013-12-17 | The Susceptibility effect minimizing technology of layer directional compensation gradient is selected in adaptive change |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104714198A CN104714198A (en) | 2015-06-17 |
| CN104714198B true CN104714198B (en) | 2019-05-17 |
Family
ID=53413689
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201310693631.1A Expired - Fee Related CN104714198B (en) | 2013-12-17 | 2013-12-17 | The Susceptibility effect minimizing technology of layer directional compensation gradient is selected in adaptive change |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104714198B (en) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106772168B (en) * | 2017-02-24 | 2019-07-02 | 深圳先进技术研究院 | Magnetic resonance imaging method and device |
| CN107861080B (en) * | 2017-10-25 | 2019-11-29 | 北京大学 | A kind of method of dynamic measurement oxygen uptake rate |
| CN118229821B (en) * | 2024-05-22 | 2024-08-09 | 山东奥新医疗科技有限公司 | Magnetic sensitive artifact correction method, device, equipment and medium |
Family Cites Families (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP3525007B2 (en) * | 1996-03-28 | 2004-05-10 | 株式会社日立メディコ | Magnetic resonance imaging system |
| US7561909B1 (en) * | 2002-09-16 | 2009-07-14 | The United States Of America As Represented By The Department Of Health And Human Services | MRI navigator methods and systems |
| CN101470180B (en) * | 2007-12-29 | 2016-01-20 | 西门子（中国）有限公司 | The method and apparatus of distortion calibration in magnetic resonance imaging |
| EP2506026A1 (en) * | 2011-03-29 | 2012-10-03 | Universitätsklinikum Freiburg | Method of dynamically compensating for magnetic field heterogeneity in magnetic resonance imaging |
| CN102749601B (en) * | 2011-04-22 | 2016-02-10 | 株式会社东芝 | Image processing apparatus, image processing method and MR imaging apparatus |
| JP5808659B2 (en) * | 2011-12-12 | 2015-11-10 | 株式会社日立メディコ | Magnetic resonance imaging apparatus and T1ρ imaging method |
| CN103278786B (en) * | 2013-03-29 | 2015-08-19 | 深圳先进技术研究院 | A kind of rapid magnetic resonance imaging method and system |
-
2013
- 2013-12-17 CN CN201310693631.1A patent/CN104714198B/en not_active Expired - Fee Related
Also Published As
| Publication number | Publication date |
|---|---|
| CN104714198A (en) | 2015-06-17 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Lutti et al. | Using high-resolution quantitative mapping of R1 as an index of cortical myelination | |
| Kundu et al. | Multi-echo fMRI: A review of applications in fMRI denoising and analysis of BOLD signals | |
| Caan et al. | MP2RAGEME: T1, T2*, and QSM mapping in one sequence at 7 tesla | |
| Jeurissen et al. | Multi-tissue constrained spherical deconvolution for improved analysis of multi-shell diffusion MRI data | |
| Li et al. | Quantitative susceptibility mapping of human brain reflects spatial variation in tissue composition | |
| Kuehn et al. | Body topography parcellates human sensory and motor cortex | |
| Shams et al. | A comparison of in vivo MRI based cortical myelin mapping using T1w/T2w and R1 mapping at 3T | |
| Li et al. | A subspace approach to spectral quantification for MR spectroscopic imaging | |
| Forstmann et al. | Towards a mechanistic understanding of the human subcortex | |
| Bazin et al. | A computational framework for ultra-high resolution cortical segmentation at 7 Tesla | |
| Hutton et al. | A comparison between voxel-based cortical thickness and voxel-based morphometry in normal aging | |
| Metzler-Baddeley et al. | How and how not to correct for CSF-contamination in diffusion MRI | |
| Eippert et al. | Denoising spinal cord fMRI data: Approaches to acquisition and analysis | |
| Hammelrath et al. | Morphological maturation of the mouse brain: An in vivo MRI and histology investigation | |
| Shmueli et al. | Magnetic susceptibility mapping of brain tissue in vivo using MRI phase data | |
| Chuang et al. | An MRI-based atlas and database of the developing mouse brain | |
| US9940712B2 (en) | Quantitating disease progression from the MRI images of multiple sclerosis patients | |
| Bao et al. | Quantitative Susceptibility Mapping Using Structural Feature Based Collaborative Reconstruction<? Pub _newline=""?>(SFCR) in the Human Brain | |
| Corbett-Detig et al. | 3D global and regional patterns of human fetal subplate growth determined in utero | |
| Xiao et al. | Multicontrast multiecho FLASH MRI for targeting the subthalamic nucleus | |
| Raoult et al. | Arterial spin labeling for motor activation mapping at 3T with a 32-channel coil: reproducibility and spatial accuracy in comparison with BOLD fMRI | |
| CN104714198B (en) | The Susceptibility effect minimizing technology of layer directional compensation gradient is selected in adaptive change | |
| Rocca et al. | Pyramidal tract lesions and movement-associated cortical recruitment in patients with MS | |
| Boillat et al. | Surface-based characteristics of the cerebellar cortex visualized with ultra-high field MRI | |
| Van Veluw et al. | FLAIR images at 7 Tesla MRI highlight the ependyma and the outer layers of the cerebral cortex |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190517 Termination date: 20191217 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |