CN103630860B - 对核磁共振图像的灰度失真进行修正的方法 - Google Patents
对核磁共振图像的灰度失真进行修正的方法 Download PDFInfo
- Publication number
- CN103630860B CN103630860B CN201210309108.XA CN201210309108A CN103630860B CN 103630860 B CN103630860 B CN 103630860B CN 201210309108 A CN201210309108 A CN 201210309108A CN 103630860 B CN103630860 B CN 103630860B
- Authority
- CN
- China
- Prior art keywords
- image
- magnetic resonance
- nuclear magnetic
- distortion
- prime
- 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
- 238000005481 NMR spectroscopy Methods 0.000 title claims abstract description 21
- 238000000034 method Methods 0.000 title claims abstract description 19
- 238000012937 correction Methods 0.000 claims abstract description 9
- 238000013507 mapping Methods 0.000 claims abstract description 4
- 238000003702 image correction Methods 0.000 claims description 3
- 238000003384 imaging method Methods 0.000 abstract description 10
- 238000003759 clinical diagnosis Methods 0.000 abstract description 4
- 238000013461 design Methods 0.000 description 4
- 238000013421 nuclear magnetic resonance imaging Methods 0.000 description 4
- 210000004556 brain Anatomy 0.000 description 3
- 238000013459 approach Methods 0.000 description 2
- 238000002595 magnetic resonance imaging Methods 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 208000024806 Brain atrophy Diseases 0.000 description 1
- 230000001154 acute effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000000747 cardiac effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 210000001640 nerve ending Anatomy 0.000 description 1
- 230000010355 oscillation Effects 0.000 description 1
- 238000002672 stereotactic surgery Methods 0.000 description 1
- 230000000638 stimulation Effects 0.000 description 1
- 238000004804 winding Methods 0.000 description 1
Landscapes
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
本发明揭示了一种对核磁共振图像的灰度失真进行修正的方法,该方法通过对非线性的梯度场进行球谐函数展开,获得线性项及高阶的非线性项,并通过常规的快速傅里叶变换FFT求出梯度场的非均匀性,以及图像几何位置和亮度的非均匀映射关系,进而求出梯度场非线性与图像灰度失真之间的加权系数,从而完成图像灰度的校正。相较于现有技术,本发明所述对核磁共振图像的灰度失真进行修正的方法可以快速实现对核磁共振图像的灰度失真进行修正,从而提高成像质量,提高临床诊断的可靠性,在核磁共振成像领域具有重要意义。
Description
技术领域
本发明涉及核磁共振成像,尤其涉及一种对核磁共振图像的灰度失真进行修正的方法。
背景技术
在核磁共振成像***中,图像质量依赖于成像对象的空间编码的准确性。脑立体定位手术、脑萎缩和脑功能成像纵向测量方面等,都会把图像变形和相应地磁场均匀性作为决定性要素。能造成图像灰度不均匀性的磁场非均匀性包括电磁涡流,梯度非均匀性,以及B0和B1的非均匀性。
近年来梯度的发展趋势具有短腔设计和高速翻转的特点,但是这些设计遇到了由大翻转电流引起的神经末梢热刺激问题。尤其是在某些应用如心脏的EPI(echo planar imaging)和大脑的DTI(diffusiontensor imaging)方面,这个瓶颈问题更加突出。作为妥协的措施,设计人员只能想办法减少梯度线圈的圈数,因而增加了梯度的非线性。
通常我们假设梯度场是线性的,因此k空间信号可以线性采样,然后通过FFT(fast Fourier transformation)进行图像重建。但是,实际上梯度场是非线性的。由于对非线性梯度场进行线性采样,从而造成图像空间的错误编码,进而使得成像质量下降。一个可行的办法是,用非线性傅里叶变换来做图像重建。但是,重建的计算时间比通常的FFT办法将增加N/log2N倍,使得实时图象重建变成不可能。
鉴于上述问题,有必要提供一种对核磁共振图像的灰度失真进行修正的方法,以便解决上述问题。
发明内容
针对现有技术的不足,本发明解决的技术问题是提供一种对核磁共振图像的灰度失真进行修正的方法,所述对核磁共振图像的灰度失真进行修正的方法可以快速实现对核磁共振图像的灰度失真进行修正,从而提高成像质量,提高临床诊断的可靠性,在核磁共振成像领域具有重要意义。
为解决上述技术问题,本发明的技术方案是这样实现的:
一种对核磁共振图像的灰度失真进行修正的方法,包括如下步骤:a:定义非线性梯度场其中L代表梯度线性分量,N代表非线性分量,v是任意坐标轴,并对非线性的梯度场进行球谐函数展开,从而获得线性项及高阶的非线性项,其中
b:通过常规的快速傅里叶变换FFT求出梯度场的非均匀性;
c:求出图像几何位置和亮度的非均匀映射关系,定义:
其中
d:求出梯度场非线性与图像灰度失真之间的加权系数;
e:完成图像灰度的校正,其中,设磁场中心点的标准亮度为V0,任一点的亮度为Vi,则:
图像校正为:ρ(x′,y′,z′)=ρ(x,y,z),
灰度校正为:
本发明的有益效果是:本发明所述对核磁共振图像的灰度失真进行修正的方法可以快速实现对核磁共振图像的灰度失真进行修正,从而提高成像质量,提高临床诊断的可靠性,在核磁共振成像领域具有重要意义。
具体实施方式
为了使本发明的目的、技术方案和优点更加清楚,下面对本发明对核磁共振图像的灰度失真进行修正的方法进行详细描述。
在核磁共振成像理论中,图像空间(r=(x,y,z))和k空间(k=(kx,ky,kz))通过傅里叶变换发生联系。扫描对象经过rf脉冲激发以后,加上三维梯度空间编码,产生震荡信号。用三维傅里叶变换可以写成:
s(k)=∫ρ(r)e-i2πk·rd3r
其中,s(k)是k空间中的信号,ρ(r)是图像空间中的扫描对象自旋密度,是最终得到的图像空间中的图像灰度。k的3个时间函数分量可由各自的梯度分量积分得到:
其中,是常量,(Gx,Gy,Gz)是梯度分量。
在标准的核磁共振成像过程中,时间域信号以均匀速度采集,而且k空间是基于梯度场均匀的假设,因而形成了以频率编码的图像空间。由于梯度场是非线性的,从而造成频率编码的非线性,进而造成图像的几何形变。因此,要解决图像变形的问题,首先要知道梯度场的真实轮廓。
描述梯度场最常用的方法是用球谐函数作基本函数:
BV(n,m)(r,θ,φ)=rn[aV(n,m)cos(mφ)+bV(n,m)sin(mφ)]P(n,m)(cosθ)
其中V代表X,Y或Z梯度场,BV(n,m)是n阶m度的球谐函数扩展。a,b是基于特定梯度场的常参量,r是距离磁场中心的半径,P是关联的拉格朗日函数。球谐函数的有限项之和,可以作为真实梯度场B2的近似。
用反卷积的方法以得到球谐函数的参量。一般来说,对DSV(diameter spherical volume)的多平面采样,可以保证磁场的过采样以及关于非均匀性的准确估计。通常的经验是,对主磁体的磁场使用24平面,得到23阶系数。而依据梯度***设计的理论数据,梯度场非均匀性使用至少8平面7阶系数。即使针对特定的梯度设计,梯度线圈的绕线误差会引入更多磁场误差,使得单纯使用简单的理论梯度场扩展也会造成不准确。更好的办法是测量实际的梯度产生的磁场,从中求出能够表达这个磁场的函数。
重新定义梯度场:
其中L代表梯度线性分量,N代表非线性分量,v是任意坐标轴,Bz是成像中起作用的磁场总和。
根据球谐函数,我们有:
和
再定义直角坐标系中非线性部分:
其中
图像校正可以写成:
ρ(x′,y′,z′)=ρ(x,y,z)
其中
求解过程如下:
每一步迭代都计算误差项:
当迭代之间的误差项小于0.01时,迭代停止。
几何坐标校正之后,还需要进行像素点亮度的校正。
假设磁场中心点的标准亮度V0,欲求出任一点的亮度Vi:
或者
综上所述,本发明对核磁共振图像的灰度失真进行修正的方法包括如下步骤:
a:对非线性的梯度场进行球谐函数展开,获得线性项及高阶的非线性项;
b:通过常规的快速傅里叶变换FFT求出梯度场的非均匀性;
c:求出图像几何位置和亮度的非均匀映射关系;
d:求出梯度场非线性与图像灰度失真之间的加权系数;
e:完成图像灰度的校正。
相较于现有技术,本发明所述对核磁共振图像的灰度失真进行修正的方法可以快速实现对核磁共振图像的灰度失真进行修正,从而提高成像质量,提高临床诊断的可靠性,在核磁共振成像领域具有重要意义。
特别需要指出,对于本领域的普通技艺人员来说,在本发明的教导下所作的针对本发明的等效变化,仍应包含在本发明申请专利范围所主张的范围中。
Claims (1)
1.一种对核磁共振图像的灰度失真进行修正的方法,其特征在于:
所述对核磁共振图像的灰度失真进行修正的方法包括如下步骤:
a:定义非线性梯度场其中L代表梯度线性分量,N代表非线性分量,v是任意坐标轴,并对非线性的梯度场进行球谐函数展开,从而获得线性项及高阶的非线性项,其中
b:通过常规的快速傅里叶变换FFT求出梯度场的非均匀性;
c:求出图像几何位置和亮度的非均匀映射关系,定义:
其中
d:求出梯度场非线性与图像灰度失真之间的加权系数;
e:完成图像灰度的校正,其中,设磁场中心点的标准亮度为V0,任一点的亮度为Vi,则:
图像校正为:ρ(x′,y′,z′)=ρ(x,y,z),
灰度校正为:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210309108.XA CN103630860B (zh) | 2012-08-28 | 2012-08-28 | 对核磁共振图像的灰度失真进行修正的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210309108.XA CN103630860B (zh) | 2012-08-28 | 2012-08-28 | 对核磁共振图像的灰度失真进行修正的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103630860A CN103630860A (zh) | 2014-03-12 |
CN103630860B true CN103630860B (zh) | 2016-12-21 |
Family
ID=50212111
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210309108.XA Active CN103630860B (zh) | 2012-08-28 | 2012-08-28 | 对核磁共振图像的灰度失真进行修正的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103630860B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105989371B (zh) * | 2015-03-03 | 2019-08-23 | 香港中文大学深圳研究院 | 一种核磁共振影像的灰度标准化方法和装置 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1678080A (zh) * | 2004-04-02 | 2005-10-05 | 乐金电子(沈阳)有限公司 | 混合灰度校正方法 |
CN101915901A (zh) * | 2010-08-17 | 2010-12-15 | 中国科学院深圳先进技术研究院 | 磁共振成像方法及装置 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1978080A (zh) * | 2005-11-29 | 2007-06-13 | 中冶赛迪工程技术股份有限公司 | 一种带钢生产工艺-esp |
EP2402780A1 (en) * | 2010-06-23 | 2012-01-04 | Koninklijke Philips Electronics N.V. | Method of reconstructing a magnetic resonance image of an object considering higher-order dynamic fields |
-
2012
- 2012-08-28 CN CN201210309108.XA patent/CN103630860B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1678080A (zh) * | 2004-04-02 | 2005-10-05 | 乐金电子(沈阳)有限公司 | 混合灰度校正方法 |
CN101915901A (zh) * | 2010-08-17 | 2010-12-15 | 中国科学院深圳先进技术研究院 | 磁共振成像方法及装置 |
Non-Patent Citations (2)
Title |
---|
A novel phantom and method for comprehensive 3-dimensional measurement and correction of geometric distortion in magnetic resonance imaging;Deming Wang et al.;《Magnetic Resonance Imaging》;20040531;第22卷(第4期);全文 * |
基于图像灰度的SSDA匹配算法;刘晓龙等;《航空计算技术》;20100131;第40卷(第1期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN103630860A (zh) | 2014-03-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Buonincontri et al. | MR fingerprinting with simultaneous B1 estimation | |
Janke et al. | Use of spherical harmonic deconvolution methods to compensate for nonlinear gradient effects on MRI images | |
US7088099B2 (en) | Correction of magnetic resonance images | |
US10874326B2 (en) | System and method for fast imaging in magnetic resonance imaging | |
US10185012B2 (en) | System and method for magnetic resonance imaging | |
CN107607895B (zh) | 磁共振成像方法和*** | |
Koch et al. | Rapid calculations of susceptibility-induced magnetostatic field perturbations for in vivo magnetic resonance | |
US9448289B2 (en) | Background field removal method for MRI using projection onto dipole fields | |
US9329252B2 (en) | Apparatus for real-time phase correction for diffusion-weighted magnetic resonance imaging using adaptive RF pulses | |
JP5868160B2 (ja) | 磁気共鳴イメージング装置 | |
US9529068B2 (en) | System and method for image warp correction for magnetic resonance imaging | |
WO2010035569A1 (ja) | 磁気共鳴イメージング装置 | |
Xu et al. | Prospective and retrospective high order eddy current mitigation for diffusion weighted echo planar imaging | |
CN104337516A (zh) | 磁共振控制序列的确定 | |
CN109143135B (zh) | 在梯度运动传感成像用途中校正梯度不均匀性的方法 | |
US11204408B2 (en) | System and method for magnetic resonance imaging | |
Haskell et al. | Off‐resonance artifact correction for MRI: A review | |
Meineke et al. | Data consistency‐driven determination of‐fluctuations in gradient‐echo MRI | |
US20160003929A1 (en) | Magnetic resonance imaging with different frequency encoding patterns | |
WO2018184056A1 (en) | Magnetic resonance imaging method and apparatus | |
Wang et al. | Stochastic optimization of 3D non-cartesian sampling trajectory (SNOPY) | |
US11022667B2 (en) | System and method for image reconstruction | |
CN103630860B (zh) | 对核磁共振图像的灰度失真进行修正的方法 | |
Li et al. | Comparison of uniform‐density, variable‐density, and dual‐density spiral samplings for multi‐shot DWI | |
Luengviriya et al. | Combination of prospective and retrospective motion correction for multi-channel MRI |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
CB02 | Change of applicant information |
Address after: Suzhou City, Jiangsu province Kunshan 215347 Yushan Zhen Yu Yang Road No. 299, room 9 Applicant after: Jiangsu Maigesi Frequency Instrument Co., Ltd. Address before: 215347, 13 -1304, South Gate Road, Yushan Town, Kunshan City, Jiangsu, Suzhou province 1699 Applicant before: JIANGSU JINGYU JINGYU DETECTOR TECHNOLOGY CO., LTD. |
|
COR | Change of bibliographic data | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |