CN104931904B - 一种ppi的多对比度磁共振图像的联合重构方法 - Google Patents
一种ppi的多对比度磁共振图像的联合重构方法 Download PDFInfo
- Publication number
- CN104931904B CN104931904B CN201510299052.8A CN201510299052A CN104931904B CN 104931904 B CN104931904 B CN 104931904B CN 201510299052 A CN201510299052 A CN 201510299052A CN 104931904 B CN104931904 B CN 104931904B
- Authority
- CN
- China
- Prior art keywords
- image
- magnetic resonance
- formula
- contrast magnetic
- ppi
- 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
Landscapes
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
本发明涉及磁共振成像领域,旨在提供一种PPI的多对比度磁共振图像的联合重构方法。该PPI的多对比度磁共振图像的联合重构方法包括下述步骤:采集需要进行图像重构的多对比度磁共振图像,采用空间灵敏度编码技术对采集的图像进行重建,利用模型进行图像重构,最后通过计算得到重构图像。本发明通过建立合理的数学模型,设计出快速、高效的重构算法应用与磁共振图像的联合重构问题上,从而达到减少扫描时间、提高成像质量、减少患者痛苦和治疗费用的目的。
Description
技术领域
本发明是关于磁共振成像领域,特别涉及一种PPI的多对比度磁共振图像的联合重构方法。
背景技术
磁共振成像是一种利用核磁共振原理进行人体成像的技术,可精确获得人体组织和器官的生理功能、解剖结构和病变信息,成为当今医学诊断的重要手段之一。
然而,受Nyquist采样定理限制,要想得到好的成像效果,需要获取大量的K空间(即傅立叶空间)数据,这将花费比较长的扫描时间。扫描时间过长会造成下述几个方面的问题:1、给病人造成额外的痛苦;2、某些需短时扫描的情况,由于器官的运动(譬如呼吸、眨眼、吞咽等)会造成图像模糊和失真;3、无法满足心脏实时成像等高端临床应用的要求等。
为了减少扫描时间,人们常常采用减少K空间数据的扫描量,然后通过图像重构技术,重构出具有高质量的MR图像。其中并行成像(PI)法,特别是最近受人们愈来愈多关注的部分并行成像(PPI)法,是一种十分有效的方法。然而,由于PPI是运用欠采样方式,但是欠采样方式的采样频率不满足Nyquist采样频率,从而导致了混叠伪影。无法满足重构图像的要求。
发明内容
本发明的主要目的在于克服现有技术中的不足,提供一种能够在不降低重构图像质量的前提下增加欠采样比例,从而减少扫描时间的联合重构方法。为解决上述技术问题,本发明的解决方案是:
提供一种PPI的多对比度磁共振图像的联合重构方法,用于对欠采样的多对比度磁共振图像进行图像重构,所述PPI的多对比度磁共振图像的联合重构方法具体包括下述步骤:
步骤A:采集需要进行图像重构的多对比度磁共振图像,所述多对比度磁共振图像包括T1、T2和PD加权图像,其中T1表示纵向弛豫图像,T2表示横向弛豫图像,PD表示质子密度图像,且多对比度磁共振图像都是独立重构图像;
步骤B:采用空间灵敏度编码技术对采集的图像进行重建,用于消除因PPI采用欠采样方式而产生的图像中的伪影,空间灵敏度编码技术(SENSE)基于下面的方程实现:
MF(Sj⊙u)=fj;
其中,所述⊙是两个向量之间的Hadamard(或分量)积;所述fj是第j个线圈所带的采集器测量到的K空间的信号,将信号按规则排列成一个向量;所述M是描述采集轨迹的mask,它是一个在采集数据的地方取1,其它地方取0的二值矩阵;所述F是Fourier变换;所述u指要重建的图像;所述Sj是第j个线圈的灵敏度映射(sensitivity map),能反映线圈与物体之间的距离对图像灰度的影响;
步骤C:将T1图像记为u1,T2图像记为u2,PD图像记为u3,利用式(1)的模型进行图像重构:
其中,所述||(MF(Sju))i-(fj)i||表示在范数意义下重构信号与观察信号的差;所述λ表示根据经验取的正数;所述i是1至3的自然数,表示图像序号;所述j是1至l的自然数,表示扫的线圈序号;所述TV(u)表示u的TV范数;
且T1,T2和PD图像的边界信息相同,灵敏度映射Sj在T1,T2和PD图像中不变;
步骤D:利用对于适当的线性算子K和适当的矩阵A,令F(Ku)=TV(u),A和b分别取下述形式:
将式(1)写成下述式(2)的形式:
minu∈XF(Ku)+λ||Au-b||2 式(2);
其中,u=(u1,u2,u3)T,即向量(u1,u2,u3)的转置,表示重构图像,且式(2)的u和步骤B的方程中的u表示同一个需要重建的目标图片;
步骤E:先用算子分裂算法对式(2)进行变量分离,得到
其中,所述z用于代替Kx,然后加个约束使得z=Kx,这里的x表示集合X中的一个元素;所述subject to是使得的意思,把上述无约束问题拆分成两个变量的问题再加个约束,以方便计算;
然后写出式(3)的增广拉格朗日公式:
其中,ρ指惩罚项的系数;第k步迭代并更新x、z过程是:
λk+1=λk+ρ(Kuk+1-zk+1);
其中,所述k是指迭代到第k步;
因为A不容易求逆,对A作线性化展开:
再用ADMM算法交替迭代更新z和u,直至达到预先设定的最大迭代步数或者满足事先设定的迭代终止条件,得到的u即为最后重构的图像。
与现有技术相比,本发明的有益效果是:
通过建立合理的模型,设计出快速、高效的重构算法应用与磁共振图像的联合重构问题上,从而达到减少扫描时间、提高成像质量、减少患者痛苦和治疗费用的目的。
具体实施方式
下面结合具体实施方式对本发明作进一步详细描述:
一种PPI的多对比度的磁共振图像的联合重构模型,包括以下步骤:
步骤A,图像采集;
步骤B,采用空间灵敏度编码技术;
步骤C,提出图像重构的数学模型;
步骤D,将模型看成非光滑凸优化问题;
步骤E,交替迭代求解方程,迭代出最终的重构结果。
上述步骤的具体实施过程如下:
步骤A:采集需要进行图像重构的多对比度磁共振图像,所述多对比度磁共振图像包括T1、T2和PD加权图像,其中T1表示纵向弛豫图像,T2表示横向弛豫图像,PD表示质子密度图像,且多对比度磁共振图像都是独立重构图像;
步骤B:PPI采用欠采样方式,但是欠采样频率不满足Nyquist采样频率,这导致了混叠伪影。采用空间灵敏度编码技术进行重建,可以消除图像中的伪影。灵敏度技术(SENSE)是PPI中常用的图像重构方法之一。它基于下面的方程
MF(Sj⊙u)=fj
其中⊙是两个向量之间的Hadamard(或分量)积,fj是第j个线圈所带的采集器测量到的K空间的信号,我们可以把它按一定的规则排列成一个向量,M是描述采集轨迹的mask,它是一个在采集数据的地方取1,其它地方取0的二值矩阵,F是Fourier变换,u是要重建的图像,Sj是第j个线圈的灵敏度映射(sensitivity map),它反映了线圈与物体之间的距离对图像灰度的影响。
步骤C:将T1图像记为u1,T2图像记为u2,PD图像记为u3,利用式(1)的模型进行图像重构:
其中,所述||(MF(Sju))i-(fj)i||表示在范数意义下重构信号与观察信号的差;λ表示根据经验取的正数;i是1至3的自然数,表示图像序号;j是1至l的自然数,表示扫的线圈序号。TV(u)表示u的TV范数。且灵敏度映射Sj在T1,T2和PD图像中不变;边界信息相同,尽管T1,T2和PD加权图像的对比度(灰度值分布)不同,但由于这三种不同对比度的图像的扫描是同一个物体,所以它们理应共享相同的轮廓和边界。从数学角度上描述为这些多对比度图像的梯度具有相似稀疏模式(similar sparsity pattern)
步骤D:利用对于适当的线性算子K和适当的矩阵A,令F(Ku)=TV(u),A和b分别取下述形式
式(1)很容易就写成式(2)的形式:
minu∈XF(Ku)+λ||Au-b||2 式(2)
其中,u=(u1,u2,u3)T,即向量(u1,u2,u3)的转置,表示重构图像。这是一类非光滑凸优化问题;
步骤E:先用算子分裂算法对式(2)进行变量分离,得到
然后写出式(3)的增广拉格朗日公式,
迭代计算并更新x、z过程是:
λk+1=λk+ρ(Kuk+1-zk+1)
因为A不容易求逆,我们先对A作线性化展开
再用ADMM算法就能交替迭代更新z和u,在达到预先设定的最大迭代步数或者满足事先设定的迭代终止条件的时候得到的u就是我们最后重构的图像。
最后,需要注意的是,以上列举的仅是本发明的具体实施例。显然,本发明不限于以上实施例,还可以有很多变形。本领域的普通技术人员能从本发明公开的内容中直接导出或联想到的所有变形,均应认为是本发明的保护范围。
Claims (1)
1.一种PPI的多对比度磁共振图像的联合重构方法,用于对欠采样的多对比度磁共振图像进行图像重构,其特征在于,所述PPI的多对比度磁共振图像的联合重构方法具体包括下述步骤:
步骤A:采集需要进行图像重构的多对比度磁共振图像,所述多对比度磁共振图像包括T1、T2和PD加权图像,其中T1表示纵向弛豫图像,T2表示横向弛豫图像,PD表示质子密度图像,且多对比度磁共振图像都是独立重构图像;
步骤B:采用空间灵敏度编码技术对采集的图像进行重建,用于消除因PPI采用欠采样方式而产生的图像中的伪影,空间灵敏度编码技术基于下面的方程实现:
MF(Sj⊙u)=fj;
其中,所述⊙是两个向量之间的Hadamard积;所述fj是第j个线圈所带的采集器测量到的K空间的信号,将信号按规则排列成一个向量;所述M是描述采集轨迹的mask,它是一个在采集数据的地方取1,其它地方取0的二值矩阵;所述F是Fourier变换;所述u指要重建的图像;所述Sj是第j个线圈的灵敏度映射,能反映线圈与物体之间的距离对图像灰度的影响;
步骤C:将T1图像记为u1,T2图像记为u2,PD图像记为u3,利用式(1)的模型进行图像重构:
其中,所述||(MF(Sju))i-(fj)i||表示在范数意义下重构信号与观察信号的差;所述λ表示根据经验取的正数;所述i是1至3的自然数,表示图像序号;所述j是1至l的自然数,表示扫的线圈序号;所述TV(u)表示u的TV范数;
且T1,T2和PD图像的边界信息相同,灵敏度映射Sj在T1,T2和PD图像中不变;
步骤D:利用对于适当的线性算子K和适当的矩阵A,令F(Ku)=TV(u),A和b分别取下述形式:
将式(1)写成下述式(2)的形式:
minu∈XF(Ku)+λ||Au-b||2 式(2);
其中,u=(u1,u2,u3)T,即向量(u1,u2,u3)的转置,表示重构图像,且式(2)的u和步骤B的方程中的u表示同一个需要重建的目标图片;
步骤E:先用算子分裂算法对式(2)进行变量分离,得到
其中,所述z用于代替Kx,然后加个约束使得z=Kx,这里的x表示集合X中的一个元素;所述subject to是使得的意思,把上述无约束问题拆分成两个变量的问题再加个约束,以方便计算;
然后写出式(3)的增广拉格朗日公式:
其中,ρ指惩罚项的系数;第k步迭代并更新x、z过程是:
λk+1=λk+ρ(Kuk+1-zk+1);
其中,所述k是指迭代到第k步;
因为A不容易求逆,对A作线性化展开:
再用ADMM算法交替迭代更新z和u,直至达到预先设定的最大迭代步数或者满足事先设定的迭代终止条件,得到的u即为最后重构的图像。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510299052.8A CN104931904B (zh) | 2015-01-27 | 2015-06-03 | 一种ppi的多对比度磁共振图像的联合重构方法 |
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2015100405604 | 2015-01-27 | ||
CN201510040560 | 2015-01-27 | ||
CN201510299052.8A CN104931904B (zh) | 2015-01-27 | 2015-06-03 | 一种ppi的多对比度磁共振图像的联合重构方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104931904A CN104931904A (zh) | 2015-09-23 |
CN104931904B true CN104931904B (zh) | 2018-10-30 |
Family
ID=54119157
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510299052.8A Active CN104931904B (zh) | 2015-01-27 | 2015-06-03 | 一种ppi的多对比度磁共振图像的联合重构方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104931904B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106772167B (zh) * | 2016-12-01 | 2019-05-07 | 中国科学院深圳先进技术研究院 | 核磁共振成像方法及装置 |
CN107242873B (zh) * | 2017-07-05 | 2020-06-09 | 成都信息工程大学 | 一种基于功能磁共振心理生理交互的脑网络构建方法 |
CN107561467A (zh) * | 2017-07-31 | 2018-01-09 | 上海东软医疗科技有限公司 | 磁共振多对比度图像重建方法和装置 |
CN107576925B (zh) * | 2017-08-07 | 2020-01-03 | 上海东软医疗科技有限公司 | 磁共振多对比度图像重建方法和装置 |
CN111784788A (zh) * | 2020-06-04 | 2020-10-16 | 深圳深透医疗科技有限公司 | 一种基于深度学习的pet快速成像方法和系统 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102362191A (zh) * | 2009-03-25 | 2012-02-22 | 皇家飞利浦电子股份有限公司 | 具有经运动校正的线圈灵敏度的磁共振部分并行成像(ppi) |
CN103403567A (zh) * | 2011-03-01 | 2013-11-20 | 皇家飞利浦有限公司 | 涉及图像比约束重建的加速mr温度测定绘图 |
CN103430038A (zh) * | 2011-03-17 | 2013-12-04 | 皇家飞利浦有限公司 | 没有图像劣化的较快的逐个信道重建的mri方法 |
CN103767705A (zh) * | 2012-10-23 | 2014-05-07 | 三星电子株式会社 | 磁共振成像系统和磁共振成像方法 |
CN103810731A (zh) * | 2014-01-20 | 2014-05-21 | 浙江大学 | 一种基于tv范数的pet图像重建方法 |
CN104156994A (zh) * | 2014-08-14 | 2014-11-19 | 厦门大学 | 一种压缩感知磁共振成像的重建方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8989465B2 (en) * | 2012-01-17 | 2015-03-24 | Mayo Foundation For Medical Education And Research | System and method for medical image reconstruction and image series denoising using local low rank promotion |
-
2015
- 2015-06-03 CN CN201510299052.8A patent/CN104931904B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102362191A (zh) * | 2009-03-25 | 2012-02-22 | 皇家飞利浦电子股份有限公司 | 具有经运动校正的线圈灵敏度的磁共振部分并行成像(ppi) |
CN103403567A (zh) * | 2011-03-01 | 2013-11-20 | 皇家飞利浦有限公司 | 涉及图像比约束重建的加速mr温度测定绘图 |
CN103430038A (zh) * | 2011-03-17 | 2013-12-04 | 皇家飞利浦有限公司 | 没有图像劣化的较快的逐个信道重建的mri方法 |
CN103767705A (zh) * | 2012-10-23 | 2014-05-07 | 三星电子株式会社 | 磁共振成像系统和磁共振成像方法 |
CN103810731A (zh) * | 2014-01-20 | 2014-05-21 | 浙江大学 | 一种基于tv范数的pet图像重建方法 |
CN104156994A (zh) * | 2014-08-14 | 2014-11-19 | 厦门大学 | 一种压缩感知磁共振成像的重建方法 |
Non-Patent Citations (1)
Title |
---|
Computational Acceleration for MR Image Reconstruction in Partially Parallel Imaging;Xiaojing Ye et al;《IEEE Transactions on Medical Imaging》;20110531;第30卷(第5期);1055-1063 * |
Also Published As
Publication number | Publication date |
---|---|
CN104931904A (zh) | 2015-09-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104931904B (zh) | 一种ppi的多对比度磁共振图像的联合重构方法 | |
CN111081354B (zh) | 用于通过深度学习网络对医疗图像进行去噪的系统和方法 | |
Cordero-Grande et al. | Sensitivity encoding for aligned multishot magnetic resonance reconstruction | |
Knoll et al. | Joint MR-PET reconstruction using a multi-channel image regularizer | |
Lingala et al. | Deformation corrected compressed sensing (DC-CS): a novel framework for accelerated dynamic MRI | |
US9140770B2 (en) | Method for reconstructing a signal from distorted experimental measurements and device for its implementation | |
CN110148215B (zh) | 一种基于平滑约束和局部低秩约束模型的四维磁共振图像重建方法 | |
CN113436290A (zh) | 用于使用深度神经网络从图像中选择性地去除条纹伪影和噪声的方法和系统 | |
Gramfort et al. | Denoising and fast diffusion imaging with physically constrained sparse dictionary learning | |
Ma et al. | Accelerated cardiac diffusion tensor imaging using joint low-rank and sparsity constraints | |
CN103278784B (zh) | 一种多约束滑动窗的磁共振并行成像方法 | |
Dar et al. | Synergistic reconstruction and synthesis via generative adversarial networks for accelerated multi-contrast MRI | |
Rasch et al. | Dynamic MRI reconstruction from undersampled data with an anatomical prescan | |
Hampe et al. | Investigating the challenges and generalizability of deep learning brain conductivity mapping | |
Sarma et al. | Accelerating dynamic magnetic resonance imaging (MRI) for lung tumor tracking based on low-rank decomposition in the spatial–temporal domain: A feasibility study based on simulation and preliminary prospective undersampled MRI | |
CN109118554A (zh) | 基于低秩稀疏分解的电阻抗图像重建方法 | |
US9274198B2 (en) | Method of reconstructing a signal in medical imaging on the basis of perturbed experimental measurements, and medical imaging device implementing this method | |
Jacob et al. | Improved model-based magnetic resonance spectroscopic imaging | |
Wang et al. | Accelerating MR imaging via deep Chambolle-Pock network | |
CN109920017A (zh) | 基于特征向量的自一致性的联合全变分Lp伪范数的并行磁共振成像重构方法 | |
US11748921B2 (en) | Learnable extrapolation for deep learning-based medical image reconstruction | |
CN106093814A (zh) | 一种基于多尺度低秩模型的心脏磁共振成像方法 | |
US20130027037A1 (en) | Self consistent parallel imaging with temporal sensitivity estimation in real-time magnetic resonance imaging | |
Yi et al. | Fast and Calibrationless low-rank parallel imaging reconstruction through unrolled deep learning estimation of multi-channel spatial support maps | |
Liu et al. | Accelerating the 3D T1ρ mapping of cartilage using a signal-compensated robust tensor principal component analysis model |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
CB03 | Change of inventor or designer information |
Inventor after: Hong Yuan Inventor before: Kong Dexing |
|
CB03 | Change of inventor or designer information | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CP03 | Change of name, title or address |
Address after: 310012 Room 709, 710, 7-storey East Building, No. 90 Wensan Road, Xihu District, Hangzhou City, Zhejiang Province Patentee after: Zhejiang Deshang Yunxing Medical Technology Co., Ltd. Address before: 310027 18F-F (1806) No. 173 Yugu Road, Xihu District, Hangzhou City, Zhejiang Province Patentee before: ZHEJIANG DESHANG YUNXING IMAGE SCIENCE & TECHNOLOGY CO., LTD. |
|
CP03 | Change of name, title or address |