The invention provides a magnetic resonance image reconstruction method and a device thereof. Including the magnetic resonance image reconstruction method: magnetic resonance image is decomposed into smooth and structured texture; the smooth part and the texture part regularization, establish the smooth part and the estimation of the texture component; according to the acquired magnetic resonance imaging scanner in the sampled K spatial data. For the estimation of the iterative smoothing part and texture; according to the iterative solution obtained smooth and texture reconstruction of magnetic resonance images. This method can effectively improve the imaging quality of magnetic resonance image reconstruction.
【技术实现步骤摘要】
本专利技术涉及磁共振成像领域,具体而言,涉及一种磁共振图像重建方法及装置。
技术介绍
磁共振成像(MagneticResonanceImaging,简称MRI)是医学成像的热点问题,被广泛应用于医学诊断中,这是由于磁共振成像具有非侵入性以及没有电离辐射等问题,而且磁共振成像还能够提供人体内部器官和结构的重要解剖和功能信息,同时具有非常高的软组织对比度,对于疾病的检测和诊断有着巨大的作用。因此,磁共振成像的质量也尤为重要。在磁共振成像中,成像的质量是衡量成像方法的一个重要标准。磁共振成像过程中,K-空间数据和磁共振图像之间是傅里叶变换和傅里叶逆变换的关系,采集全部的K-空间数据比较慢,只部分采集磁共振成像设备K-空间数据,会带来较多伪影,图像质量较差。现有的磁共振成像方法大多是引入一些先验信息(如稀疏性)来通过一个优化函数,求解重建出清晰图像。这些方法将图像视为一个整体来进行重建,无法克服在大多数磁共振图像重建方法中使用全变分对平滑区域带来的阶梯效应,从而影响了成像的质量。
技术实现思路
有鉴于此,本专利技术实施例提供了一种磁共振图像重建方法及装置,将磁共振图像结构化分解为平滑部分和纹理部分,再针对两个部分进行不同的正则化并建立估算式,然后进行求解得到需要重建的磁共振图像,克服了现有方法中将磁共振图像视为一个整体所带来的对图像中平滑区域的阶梯效应,并且能较好的恢复纹理部分的信息,提高了成像的质量。为了实现上述目的,本专利技术采用的技术方案如下:一种磁共振图像重建方法,所述方法包括:将磁共振图像结构化分解为平滑部分和纹理部分;对所述平滑部分以及所述纹理部分正则化 ...
【技术保护点】
一种磁共振图像重建方法,其特征在于,所述方法包括:将磁共振图像结构化分解为平滑部分和纹理部分;对所述平滑部分以及所述纹理部分正则化,建立所述平滑部分和所述纹理部分的估算式;根据磁共振图像扫描仪中获取的下采样K‑空间数据,分别对所述估算式中的平滑部分和纹理部分迭代求解;根据迭代求解获得的平滑部分以及纹理部分重建磁共振图像。
【技术特征摘要】
1.一种磁共振图像重建方法,其特征在于,所述方法包括:将磁共振图像结构化分解为平滑部分和纹理部分;对所述平滑部分以及所述纹理部分正则化,建立所述平滑部分和所述纹理部分的估算式;根据磁共振图像扫描仪中获取的下采样K-空间数据,分别对所述估算式中的平滑部分和纹理部分迭代求解;根据迭代求解获得的平滑部分以及纹理部分重建磁共振图像。2.根据权利要求1所述的方法,其特征在于,所述对所述平滑部分以及所述纹理部分正则化,建立所述平滑部分和所述纹理部分的估算式包括:运用各项同性的二阶全变分对所述平滑部分正则化,得到所述平滑部分对应的正则化ISOTV(us),运用非局部的全变分对所述纹理部分正则化,得到所述纹理部分对应的正则化NLTV(ut),us和ut分别表示平滑部分和纹理部分;根据下采样K-空间数据与磁共振图像之间的关系模型y=Au+ε、ISOTV(us)以及NLTV(ut)建立us和ut的估算式:其中,u=us+ut,y表示下采样K-空间数据,A表示下采样傅里叶算子,ε表示附加噪声,||·||2表示L2范数,λ和μ表示正则化参数。3.根据权利要求2所述的方法,其特征在于,所述运用非局部的全变分对所述纹理部分正则化,得到所述纹理部分对应的正则化NLTV(ut)包括:获取所述纹理部分中一个像素点对应的非局部离散梯度表达式;根据所述非局部离散梯度表达式得到所述纹理部分对应的正则化NLTV(ut)。4.根据权利要求1所述的方法,其特征在于,所述根据磁共振图像扫描仪中获取的下采样K-空间数据,分别对所述估算式中的平滑部分和纹理部分迭代求解包括:利用交替最小化方法,根据所述估算式建立所述平滑部分对应的迭代公式以及所述纹理部分对应的迭代公式;根据磁共振图像扫描仪中获取的下采样K-空间数据对所述平滑部分对应的迭代公式求解得到所述平滑部分对应的解,以及对所述纹理部分对应的迭代公式求解得到所述纹理部分对应的解;对所述平滑部分对应的解与所述纹理部分对应的解进行求和,当求和结果收敛时,迭代求解过程结束。5.根据权利要求3所述的方法,其特征在于,所述根据磁共振图像扫描仪中获取的下采样K-空间数据对所述平滑部分对应的迭代公式求解得到所述平滑部分对应的解包括:运用半二次方分裂最小化方法,增加辅助变量,根据所述平滑部分对应的迭代公式得到关于所述辅助变量以及所述平滑部分的估算式;利用交替最小化方法得到所述辅助变量对应的迭代公式以及所述平滑部分对应的求解迭代式;根据磁共振图像扫描仪中获取的下采样K-空间数据对所述辅助变量对应的迭代公式以及所述平滑部分对应的求解迭代式进行迭代求解,得到所述平滑部分。6.一种磁共振图像重建装置,其特征在于,所述装置包括:图像分解模块、估算式建立模块、迭代求解模块以及图...
【专利技术属性】
技术研发人员:杨晓梅,郭朝云,刘凯,曾晓东,
申请(专利权)人:四川大学,
类型:发明
国别省市:四川;51
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。