CN104434108B - 一种基于相邻层信息和欠采样的多层磁共振快速成像方法 - Google Patents
一种基于相邻层信息和欠采样的多层磁共振快速成像方法 Download PDFInfo
- Publication number
- CN104434108B CN104434108B CN201410745312.5A CN201410745312A CN104434108B CN 104434108 B CN104434108 B CN 104434108B CN 201410745312 A CN201410745312 A CN 201410745312A CN 104434108 B CN104434108 B CN 104434108B
- Authority
- CN
- China
- Prior art keywords
- big
- signaling point
- scanning
- lack sampling
- magnetic resonance
- 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
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/05—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves
- A61B5/055—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves involving electronic [EMR] or nuclear [NMR] magnetic resonance, e.g. magnetic resonance imaging
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Engineering & Computer Science (AREA)
- Medical Informatics (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- High Energy & Nuclear Physics (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Radiology & Medical Imaging (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
本发明公开了一种基于相邻层信息和欠采样的多层磁共振快速成像方法。该方法首先从与待扫描层相邻并且已经完成扫描的层面中,获取其K空间的局部大信号点和全局大信号点的分布位置;其次,据此规划待扫描层的欠采样扫描轨迹,并完成扫描;最后,通过非线性共轭梯度迭代法,对欠采样数据在稀疏变换域内以最小化一范数为目标进行优化的方式重建图像。本方法通过采样点的减少缩短了扫描时间,并且通过相邻层信息的指导提高了成像的质量和稳定性。
Description
技术领域
本发明涉及磁共振成像(Magnetic Resonance Imaging,MRI)技术领域,具体涉及一种基于相邻层信息和欠采样的多层磁共振快速成像方法。
背景技术
磁共振成像是利用物质中的原子核在磁场中受到射频脉冲激励而发生核磁共振现象,产生磁共振信号,经过电子计算机处理,重建出物质内部结构图像的成像技术。磁共振成像作为一项新的医学影像诊断技术,近年来发展十分迅速。磁共振成像无放射性,一般来说所提供的信息量较其他影像技术更为丰富,能够更好的为疾病诊断提供支持。
但是和其他医学影像技术相比,磁共振成像技术成像速度相对较慢,这限制了磁共振成像技术更广泛的应用。图像压缩技术的研究表明,磁共振图像的可压缩比很高。这表明磁共振信号采集过程中有很多冗余成分。通过欠采样技术可以在采样过中节省采样时间,同时获得与满采样相比图像质量相近的图像。但是通过随机的方式决定欠采样点位置的采样方法,难以保证重建图像质量的稳定性。
对于磁共振连续多层成像来说,其相邻层面图像的K空间具有较高的相似度。本发明公开的方法,利用了连续多层扫描过程中,已完成扫描层面的K空间提供的信息来辅助规划新扫描层的欠采样扫描轨迹,在加快扫描速度的同时提高重建图像的质量和稳定性。
发明内容
本发明的目的是针对现有技术存在的上述磁共振成像速度较慢的问题,提出一种基于相邻层信息和欠采样的多层磁共振快速成像方法。
为了实现上述目的,本发明采用以下技术措施:
一种基于相邻层信息和欠采样的多层磁共振快速成像方法,包括以下步骤:
步骤1、对扫描层面序列的第一层,进行满采样扫描,并通过傅里叶变换重建图像;
步骤2、进入下一层的扫描,对与当前扫描层面相邻的前一层的K空间数据进行分析,获取其中局部大信号点和全局大信号点的位置信息。挑选局部大信号点的标准为在K空间中,数据幅值比其四邻域或八邻域中其他点都大的点;挑选全局大信号点的标准为在K空间中,数据幅值从大到小进行排列,选取预定比例数据幅值较大的点。作为优选方案,在采样过程中挑选出的全局大信号点的数量应和局部大信号点的数量相近或相同。通过这种方式可以实现图像轮廓和细节之间的平衡,获得视觉效果更佳的图像;
步骤3、根据步骤2获得的全局大信号点和局部大信号点位置信息,规划当前层面的欠采样扫描轨迹。统计笛卡尔扫描轨迹中每行轨迹所包含的全局大信号点和局部大信号点数量,在笛卡尔扫描轨迹中按照预定行数选取包含全局大信号点和局部大信号点数量较大的行,从而确定欠采样扫描轨迹,预定行数由欠采样比率所确定;
步骤4、根据步骤3得到的欠采样扫描轨迹,对当前层面进行扫描,得到当前层面的欠采样K空间数据;
步骤5、将步骤4中得到的欠采样K空间数据通过零填充和傅里叶逆变换得到当前层面重建图像I0 ;
步骤6、以步骤5中的重建图像I0为重建初始值,通过非线性共轭梯度法,以重建图像If的稀疏变换的1范数最小化为目标,以重建图像If在傅里叶逆变换后与步骤4中获得的欠采样K空间数据之差的2范数小于设定值为约束条件,进行优化迭代运算,从而得到最终的重建图像If。可以使用稀疏变换包括离散余弦变换、小波变换、有限差分变换等。磁共振图像在上述几种稀疏变换域内都有较好的稀疏化效果,但是不同身体部位的磁共振图像,在不同的稀疏变换域内的稀疏程度不同,需要根据具体情况合理选取;
步骤7、通过对重建图像If的傅里叶逆变换,得到当前层面的完整K空间数据;返回步骤2直至扫描层面序列中所有扫描及图像重建完成。
本发明具有以下优点和积极效果:
本发明可以加快多层磁共振扫描的速度,提高MRI设备的使用效率,节省病人的使用成本,并为MRI新方法的研究提供技术基础。
附图说明
图1发明方法扫描流程图;
图2实施例2中其中一层K空间大点分布情况;
图3实施例2中其中一层K空间欠采样扫描轨迹;
图4实施例2中其中一层迭代初始图像;
图5欠采样图像重建算法流程;
图6(a)为利用相邻层信息指导的欠采样重建图像;
图6(b)为利用随机欠采样后进行非线性共轭梯度重建算法所得的图像。
具体实施方式
下面结合附图和实施示例对本发明进一步说明:
实施例1:
一种基于相邻层信息和欠采样的多层磁共振快速成像方法,包括以下步骤:
步骤1、对需要扫描层面序列的第一层,进行满采样扫描,并通过傅里叶变换重建图像;
步骤2、进入下一层的扫描,对与当前扫描层面相邻的前一层的K空间数据进行分析,获取其中局部大信号点和全局大信号点的位置信息;
步骤3、根据步骤2获得的全局大信号点和局部大信号点位置信息,规划当前层面的欠采样扫描轨迹;
步骤4、根据步骤3得到的欠采样扫描轨迹,对当前层面进行扫描,得到当前层面的欠采样K空间数据;
步骤5、将步骤4中得到的欠采样K空间数据通过零填充和傅里叶逆变换得到当前层面重建图像I0 ;
步骤6、以步骤5中的重建图像I0为重建初始值,通过非线性共轭梯度法,以重建图像If的稀疏变换的1范数最小化为目标,以重建图像If在傅里叶逆变换后与步骤4中获得的欠采样K空间数据之差的2范数小于设定值为约束条件,进行优化迭代运算,从而得到最终的重建图像If;
步骤7、通过对重建图像If的傅里叶逆变换,得到当前层面的完整K空间数据;返回步骤2直至扫描层面序列中所有扫描及图像重建完成。
步骤2中局部大信号点为在K空间中,数据幅值比其四邻域或八邻域中其他点都大的点;全局大信号点为在K空间中,数据幅值从大到小进行排列,选取预定比例数据幅值较大的点。
步骤3中规划当前层面的欠采样扫描轨迹包括以下步骤:
统计笛卡尔扫描轨迹中每行轨迹所包含的全局大信号点和局部大信号点数量,在笛卡尔扫描轨迹中按照预定行数选取包含全局大信号点和局部大信号点数量较大的行,从而确定欠采样扫描轨迹,预定行数由欠采样比率所确定;
步骤6中稀疏变换为离散余弦变换或小波变换或有限差分变换等。
实施例2:
步骤1、小鼠头部扫描层数为5层,每层厚度为0.4mm,分辨率为128*128。 首先对需要扫描层面序列的第一层,采用自旋回波(SE)序列,设定TR=2s,TE=8.667ms,采样轨迹为笛卡尔采样方式,进行满采样扫描,得到K空间数据,通过傅里叶变换重建图像;
步骤2、进入下一层的扫描,对与当前扫描层面相邻的前一层的K空间数据求取幅值,统计得出前一层的K空间中的八邻域局部大点位置以及全局大点位置。经过统计,K空间中八邻域局部大信号点数量占总点数的15%,选取全局大信号点的个数与局部大信号点一致;前15%的全局大信号点和局部大信号点之间重复的点数占总点数的6%,因此两者之和共占总点数的24%。将全局大点和局部大点位置信息融合得到K空间大点分布,如图2所示;
步骤3、按照笛卡尔采样轨迹,统计K空间中每一行所包含的大信号点(包括全局大信号点和局部大信号点)个数,按从大到小的顺序进行排序。设定欠采样比率为50%,选择前64行为所要扫描的行,确定欠采样扫描轨迹,如图3所示;
步骤4、根据步骤3得到的欠采样扫描轨迹,对当前层进行扫描,得到当前层欠采样的K空间数据;
步骤5、将步骤4中得到的欠采样K空间数据通过零填充,得到大小为128*128的欠采样K空间数据。零填充后对K空间数据做傅里叶变换得到当前层面重建图像I0,如图4所示;
步骤6、以步骤5中的重建图像I0为重建初始值并通过非线性共轭梯度法进行图像重建,该算法的流程图如图5所示。具体方法为:首先,对重建图像I0进行小波变换得到其在小波变换域内128*128的矩阵J0,计算矩阵J0的1范数,以矩阵J0的1范数最小化为目标,以重建图像I0在傅里叶逆变换后与步骤4中获得的欠采样K空间数据之差的2范数小于设定值为约束条件进行优化运算,得到重建图像I0优化后的矩阵I1以及矩阵I1在小波变换空间128*128的矩阵J1,重复进行上述的迭代运算,直到重建图像收敛于If或者达到最大迭代步数,从而得到最终的重建图像If,如图6(a)所示。图6(b)则是对该层进行随机欠采样后进行非线性共轭梯度重建算法所得的图像;
步骤7、通过对重建图像If的傅里叶逆变换,得到当前层面的完整K空间数据;返回步骤2直至扫描层面序列中所有扫描及图像重建完成。
本文中所描述的具体实施示例仅仅是对本发明精神作举例说明。本发明所属技术领域的技术人员可以对所描述的具体实施示例做各种各样的修改或补充或采用类似的方式替代,但并不会偏离本发明的精神或者超越所附权利要求书所定义的范围。
Claims (3)
1.一种基于相邻层信息和欠采样的多层磁共振快速成像方法,其特征在于,包括以下步骤:
步骤1、对需要扫描层面序列的第一层,进行满采样扫描,并通过傅里叶变换重建图像;
步骤2、进入下一层的扫描,对与当前扫描层面相邻的前一层的K空间数据进行分析,获取其中局部大信号点和全局大信号点的位置信息;
步骤3、根据步骤2获得的全局大信号点和局部大信号点位置信息,规划当前层面的欠采样扫描轨迹;
步骤4、根据步骤3得到的欠采样扫描轨迹,对当前层面进行扫描,得到当前层面的欠采样K空间数据;
步骤5、将步骤4中得到的欠采样K空间数据通过零填充和傅里叶逆变换得到当前层面重建图像I0;
步骤6、以步骤5中的重建图像I0为重建初始值,通过非线性共轭梯度法,以重建图像If的稀疏变换的1范数最小化为目标,以重建图像If在傅里叶逆变换后与步骤4中获得的欠采样K空间数据之差的2范数小于设定值为约束条件,进行优化迭代运算,从而得到最终的重建图像If;
步骤7、通过对重建图像If的傅里叶逆变换,得到当前层面的完整K空间数据;返回步骤2直至扫描层面序列中所有扫描及图像重建完成,
所述的步骤3中规划当前层面的欠采样扫描轨迹包括以下步骤:
统计笛卡尔扫描轨迹中每行轨迹所包含的全局大信号点和局部大信号点数量,在笛卡尔扫描轨迹中按照预定行数选取包含全局大信号点和局部大信号点数量较大的行,从而确定欠采样扫描轨迹,预定行数由欠采样比率所确定。
2.根据权利要求1所述的一种基于相邻层信息和欠采样的多层磁共振快速成像方法,其特征在于,所述的步骤2中局部大信号点为在K空间中,数据幅值比其四邻域或八邻域中其他点都大的点;全局大信号点为在K空间中,数据幅值从大到小进行排列,选取预定比例数据幅值较大的点。
3.根据权利要求1所述的一种基于相邻层信息和欠采样的多层磁共振快速成像方法,其特征在于,所述的步骤6中稀疏变换为离散余弦变换或小波变换或有限差分变换。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410745312.5A CN104434108B (zh) | 2014-12-09 | 2014-12-09 | 一种基于相邻层信息和欠采样的多层磁共振快速成像方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410745312.5A CN104434108B (zh) | 2014-12-09 | 2014-12-09 | 一种基于相邻层信息和欠采样的多层磁共振快速成像方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104434108A CN104434108A (zh) | 2015-03-25 |
CN104434108B true CN104434108B (zh) | 2017-02-01 |
Family
ID=52881413
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410745312.5A Active CN104434108B (zh) | 2014-12-09 | 2014-12-09 | 一种基于相邻层信息和欠采样的多层磁共振快速成像方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104434108B (zh) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2019090447A1 (en) * | 2017-11-10 | 2019-05-16 | Shenzhen United Imaging Healthcare Co., Ltd. | System and method for image reconstruction |
CN110133557B (zh) * | 2018-02-09 | 2021-08-31 | 深圳先进技术研究院 | 一种新型的非线性并行重建的磁共振成像方法、装置及介质 |
CN109239631B (zh) | 2018-08-03 | 2021-02-23 | 上海东软医疗科技有限公司 | 一种磁共振成像方法和装置 |
CN109613460B (zh) * | 2018-12-28 | 2021-03-23 | 上海联影医疗科技股份有限公司 | 磁共振成像方法、装置、计算机设备和存储介质 |
CN110082696B (zh) * | 2019-04-30 | 2021-06-08 | 上海东软医疗科技有限公司 | 一种多回波序列的成像方法和装置 |
CN112927313B (zh) * | 2019-12-05 | 2022-12-06 | 上海联影医疗科技股份有限公司 | 磁共振图像重建方法、装置、计算机设备及可读存储介质 |
CN112526423B (zh) * | 2020-12-10 | 2023-06-06 | 川北医学院 | 基于共轭和层间信息的并行磁共振成像算法 |
CN115932687A (zh) * | 2021-08-12 | 2023-04-07 | 安徽华米健康科技有限公司 | 磁共振成像mri的扫描方法、装置和电子设备 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5073752A (en) * | 1990-04-19 | 1991-12-17 | Picker International, Inc. | Discrete fourier transform imaging |
AU2002366452A1 (en) * | 2001-12-17 | 2003-06-30 | Koninklijke Philips Electronics N.V. | Magnetic resonance method |
ATE528733T1 (de) * | 2007-12-20 | 2011-10-15 | Wisconsin Alumni Res Found | Verfahren für dynamische eingeschränkte bildrekonstruktion mit vorhergehendem bild |
US8310233B2 (en) * | 2009-02-18 | 2012-11-13 | Mayo Foundation For Medical Education And Research | Method for image reconstruction from undersampled medical imaging data |
CN103027681B (zh) * | 2011-10-06 | 2016-05-04 | 西门子公司 | 用于重构并行获取的mri图像的系统 |
-
2014
- 2014-12-09 CN CN201410745312.5A patent/CN104434108B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN104434108A (zh) | 2015-03-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104434108B (zh) | 一种基于相邻层信息和欠采样的多层磁共振快速成像方法 | |
CN103477238B (zh) | 采用来自先验采集的约束的压缩感测mr图像重建 | |
US9396562B2 (en) | MRI reconstruction with incoherent sampling and redundant haar wavelets | |
JP4646015B2 (ja) | 磁気共鳴イメージング(mri)システムの作動方法 | |
CN101484823B (zh) | 磁共振成像设备和方法 | |
CN105074491B (zh) | 具有使用压缩感知的图像重建的动态mri | |
CN108957375A (zh) | 磁共振成像方法、装置、设备及存储介质 | |
CN103380384B (zh) | 使用交错的水参考扫描的具有自动相位和b0校正的磁共振波谱成像方法和装置 | |
DE102007045054A1 (de) | Vorrichtung und Verfahren zur beschleunigten Propeller-Bildgebung | |
CN106574954A (zh) | 针对epi的具有奈奎斯特伪影校正的并行mr成像 | |
CN105143906A (zh) | 金属抗mr成像 | |
CN104067137A (zh) | 具有对流动伪影的抑制的mr成像 | |
US11002815B2 (en) | System and method for reducing artifacts in echo planar magnetic resonance imaging | |
CN110133556A (zh) | 一种磁共振图像处理方法、装置、设备及存储介质 | |
WO2023093842A1 (zh) | 一种肝脏多参数定量成像方法 | |
Chen et al. | Self‐calibrating wave‐encoded variable‐density single‐shot fast spin echo imaging | |
CN106981051A (zh) | 用于运动解析的mri的系统和方法 | |
CN107110945A (zh) | 用于使用在两种不同磁场强度下采集的图像对图像细节进行分割的mri协议 | |
JP2006223869A (ja) | 連続的テーブル移動により取得された磁気共鳴画像におけるアーチファクトの除去方法 | |
Chao et al. | Fast diffusion imaging with high angular resolution | |
Wang et al. | High‐fidelity direct contrast synthesis from magnetic resonance fingerprinting | |
HashemizadehKolowri et al. | Jointly estimating parametric maps of multiple diffusion models from undersampled q‐space data: a comparison of three deep learning approaches | |
US6377831B1 (en) | Real-time MR image subtraction and reconstruction | |
CN106137199B (zh) | 扩散磁共振成像中的广义球面去卷积 | |
US10928475B2 (en) | Dynamic contrast enhanced magnetic resonance imaging with flow encoding |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |