CN102488519B - 弥散张量成像方法及系统 - Google Patents
弥散张量成像方法及系统 Download PDFInfo
- Publication number
- CN102488519B CN102488519B CN201110387526.6A CN201110387526A CN102488519B CN 102488519 B CN102488519 B CN 102488519B CN 201110387526 A CN201110387526 A CN 201110387526A CN 102488519 B CN102488519 B CN 102488519B
- Authority
- CN
- China
- Prior art keywords
- space
- sampling
- phase
- readout direction
- module
- 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
- Image Processing (AREA)
Abstract
一种弥散张量成像方法,包括以下步骤:在K空间在读出方向上进行划分,将K空间划分成多个相互独立的子空间;在K空间读出方向上对所述多个相互独立的子空间进行激发,并对K空间在相位编码方向上进行采样,得到K空间的数据;对所述K空间数据进行重建得到重建图像。此外,还提供了弥散张量成像系统。
Description
【技术领域】
本发明涉及磁共振技术,特别是涉及一种弥散张量成像方法及系统。
【背景技术】
弥散张量成像(DTI),是在弥散加权成像(DWI)基础上发展起来的新方法,是磁共振成像(MRI)的特殊形式,是近几年发展迅速的一项新的磁共振成像技术。弥散张量成像技术是利用水分子的弥散各向异性进行成像,可从微观的领域评价组织结构的完整性,可以在细胞及分子水平给出疾病状况,是功能磁共振成像的一个重要组成部分。弥散张量成像是目前唯一无创性活体研究脑白质结构及白质束形态的检查技术,也是目前唯一反映活体组织空间组成信息及病理状态下各组织成分之间水分子交换功能状况的检查方法。基于弥散张量成像的高空间分辨率、非入侵性的优点,弥散张量成像主要应用于脑部各方向白质纤维和白质纤维束的评价,且逐渐扩展到人体的其他部位(比如心脏、肾脏、骨骼肌等),可以为疾病的诊断、治疗提供更多信息。
在磁共振领域里,k空间是寻常空间在傅利叶转换下的对偶空间,主要应用在磁共振造影的成像分析。
弥散张量成像常使用单次激发-平面回波成像(single shot-echo planarimaging,SE-EPI)成像。单次激发采集由于回波链过长,链尾的信号衰减较快,在磁化率差异大的组织边界会产生严重的几何变形,这样会造成弥散张量成像易于出现伪影,使得图像质量较差,限制了其临床应用。
【发明内容】
基于此,有必要提供一种图像质量高的弥散张量成像方法。
一种弥散张量成像方法,包括以下步骤:在K空间在读出方向上进行划分,将K空间划分成多个相互独立的子空间;在K空间读出方向上对所述多个相互独立的子空间进行激发,并对K空间在相位编码方向上进行采样,得到K空间的数据;对所述K空间数据进行重建得到重建图像。
进一步地,所述将K空间读出方向上划分出多个子空间的步骤具体为:将K空间读出方向上中间区域划分出多个相互独立的子空间;所述在K空间读出方向上对所述多个子空间进行激发,并对K空间在相位编码方向上进行采样,得到K空间的数据的步骤具体为:在K空间读出方向上对中部区域的所述多个相互独立的子空间进行激发,并对K空间在相位编码方向上进行采样。
进一步地,对K空间在相位编码方向上进行采样的方式为欠采样。
进一步地,所述欠采样具体为随机欠采样、等间距欠采样及变密度欠采样中的一种。
进一步地,所述对所述K空间数据进行重建得到重建图像的步骤具体为:
对所述K空间数据进行压缩感知重建得到重建图像。
此外,还有必要提供一种快速的弥散张量成像系统。
一种弥散张量成像系统,包括处理模块、采样模块及重建模块;处理模块,用于在K空间在读出方向上进行划分,将K空间划分成多个相互独立的子空间;采样模块,与所述处理模块连接,所述采样模块用于在K空间读出方向上对所述多个相互独立的子空间进行激发,并对K空间在相位编码方向上进行采样,得到K空间的数据;重建模块,与所述采样模块连接,用于对所述K空间数据进行重建得到重建图像。
进一步地,所述处理模块具体用于将K空间读出方向上中间区域划分出多个相互独立的子空间;所述采样模块具体用于在K空间读出方向上对中部区域的所述多个相互独立的子空间进行激发,并对K空间在相位编码方向上进行采样。
进一步地,对K空间在相位编码方向上进行采样的方式为欠采样。
进一步地,所述欠采样具体为随机欠采样、等间距欠采样及变密度欠采样中的一种。
进一步地,所述重建模块对所述K空间数据进行压缩感知重建得到重建图像。
上述弥散张量成像方法及系统中,通过在K空间在读出方向上进行划分,将K空间划分成多个相互独立的子空间,并分别激发,多次激发缩短了单次激发后的回波链长度,增加激发次数,解决了仅对K空间读出方向单次激发所带来的弥散张量成像易于出现伪影的问题,大大提高了图像的质量。
进一步地,仅在K空间读出方向上中间区域划分出多个相互独立的子空间,并激发,做到了在K空间读出方向上的欠采样,在保证成像质量的同时节省了大量的时间,大大加快了弥散张量成像的速度。
【附图说明】
图1为弥散张量成像领域K空间示意图;
图2为一个实施例中弥散张量成像的流程图;
图3为一个实施例中在读出方向将K空间划分的示意图;
图4为一个实施例中弥散张量成像的示意图;
图5为另一实施例在读出方向将K空间划分的示意图;
图6为另一实施例中弥散张量成像的示意图;
图7为一个实施例中弥散张量成像系统的结构示意图;
【具体实施方式】
为了解决在传统的弥散张量成像常使用单次激发-平面回波成像过程中出现伪影,导致图像质量较差的问题,提出了一种图像质量高的弥散张量成像方法。
请参阅图1,在弥散张量成像领域中,K空间包括Kx和Ky两个方向,Kx为相位编码方向(phase-encoding direction),Ky为读出方向(readout direction)。
请参阅图2,一种弥散张量成像方法,包括以下步骤:
步骤S10,在K空间的读出方向上进行划分,将K空间划分成多个相互独立的子空间。请参阅图3,在K空间读出方向Ky上将K空间划分成多个相互独立的子空间。
步骤S20,在K空间读出方向上对多个相互独立的子空间进行激发,并对K空间在相位编码方向上进行采样,得到K空间的数据。
请参阅图4,对K空间读出方向上划分出来的多个相互独立的子空间,分别使用平面回波信号进行激发,在K空间读出方向上形成多次激发。同时,在相位编码方向上进行采样,得到每一个子空间所对应的K空间数据,将这些子空间所对应的K空间数据相加,即得到一个完整的K空间数据。
需要指出的是,在K空间相位编码方向上可以采用欠采样方式对数据进行采样。欠采样的方式可以为随机欠采样、等间距欠采样和变密度欠采样等方式中的一种。随机欠采样即在K空间相位编码方向上进行随机采样以替代全采样。等间距采样即在K空间相位编码方向上进行隔行采样。变密度欠采样即根据自适应或非自适应函数,对K空间相位编码方向上进行欠采样。在K空间相位编码方向欠采样,减少了所需采集的数据量,大大缩短了采样所需的时间,提高了整个弥散张量成像的速度。
步骤S30,对K空间数据进行重建得到重建图像。对所采集到的K空间数据进行反傅里叶变换等过程,最终完成重建得到重建图像。
需要注意的是,如果K空间数据时欠采样得到的数据,则可以采用压缩感知重建,以快速得到重建图像。压缩感知是近几年新兴起的一种能实现快速成像的理论。压缩感知基于信号或图像的稀疏性,突破奈奎斯特采样定理的限制,通过欠采样得到的极少量的采样点或观测点恢复出原始信号或图像。但必须满足以下条件:(1)数据信号是稀疏的,或者可以被稀疏表示;(2)由于欠采样所得的K空间数据所导致的混叠伪影是不一致的。使用压缩感知能大大加快图像重建的速度,缩短了弥散张量成像的时间。
上述弥散张量成像方法中,通过在K空间在读出方向上进行划分,将K空间划分成多个相互独立的子空间,并分别激发,多次激发缩短了单次激发后的回波链长度,增加激发次数,解决了仅对K空间读出方向单次激发所带来的弥散张量成像易于出现伪影的问题,大大提高了图像的质量。
在另一实施例中,步骤S10具体为:将K空间读出方向上中间区域划分出多个相互独立的子空间。
步骤S20具体为在K空间读出方向上对中部区域的多个相互独立的子空间进行激发,并对K空间在相位编码方向上进行采样。
如图5所示,仅在K空间读出方向中间区域划分出多个相互独立的子空间。由于图像中的大部分数据都集中在读出方向的中部区域,如图6所示,发射平面回波,对K空间读出方向中间区域划分出的多个相互独立的子空间进行激发,同时对K空间相位编码方向上进行采样,即能得到K空间读出方向上多个子空间所对应的数据,最后通过压缩感知或其它重建方式,得到完整的重建图像。在K空间读出方向中间区域进行数据欠采样,减少了有用信息的丢失,保证了图像的质量。
仅在K空间读出方向上中间区域划分出多个相互独立的子空间,并激发,做到了在K空间读出方向上的欠采样,在保证了成像质量的同时节省了大量的时间,大大加快了弥散张量成像的速度。
请参阅图7,还提供了一种弥散张量成像系统,包括处理模块100、采样模块200及重建模块300。
处理模块100用于在K空间在读出方向上进行划分,将K空间划分成多个相互独立的子空间。在K空间读出方向Ky上将K空间划分成多个相互独立的子空间。
采样模块200与处理模块100连接,采样模块200用于在K空间读出方向上对多个相互独立的子空间进行激发,并对K空间在相位编码方向上进行采样,得到K空间的数据。
整个系统对K空间读出方向上划分出来的多个相互独立的子空间,分别使用平面回波信号进行激发,在K空间读出方向上形成多次激发。同时,在相位编码方向上进行采样,得到每一个子空间所对应的K空间数据,将这些子空间所对应的K空间数据相加,即得到一个完整的K空间数据。
具体的,在本实施例中,采样模块200用于在K空间相位编码方向上采用欠采样方式对数据进行采样。欠采样的方式可以为随机欠采样、等间距欠采样及变密度欠采样等方式中的一种。随机欠采样即在K空间相位编码方向上进行随机采样以替代全采样。等间距采样即在在K空间相位编码方向上进行隔行采样。变密度欠采样即根据自适应或非自适应函数,对K空间相位编码方向上进行欠采样。在K空间相位编码方向欠采样,减少了所需采集的数据量,大大缩短了采样所需的时间,提高了整个弥散张量成像的速度。当然,为了保证数据的高度完整性,可以在K空间相位编码方向上进行全采样。
重建模块300与采样模块200连接,用于对所述K空间数据进行重建得到重建图像。重建模块300对所采集到的K空间数据进行反傅里叶变换等过程,最终完成重建得到重建图像。
需要注意的是,如果K空间数据时欠采样得到的数据,重建模块300可以采用压缩感知重建,以快速得到重建图像。
上述弥散张量成像系统中,通过在K空间在读出方向上进行划分,将K空间划分成多个相互独立的子空间,并分别激发,多次激发缩短了单次激发后的回波链长度,增加激发次数,解决了仅对K空间读出方向单次激发所带来的弥散张量成像易于出现伪影的问题,大大提高了图像的质量。
在另一实施例中,处理模块100具体用于将K空间读出方向上中间区域划分出多个相互独立的子空间。
采样模块200具体用于在K空间读出方向上对中部区域的多个相互独立的子空间进行激发,并对K空间在相位编码方向上进行采样。
仅在K空间读出方向上中间区域划分出多个相互独立的子空间,并激发,做到了在K空间读出方向上的欠采样,在保证了成像质量的同时节省了大量的时间,大大加快了弥散张量成像的速度。
以上所述实施例仅表达了本发明的几种实施方式,其描述较为具体和详细,但并不能因此而理解为对本发明专利范围的限制。应当指出的是,对于本领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干变形和改进,这些都属于本发明的保护范围。因此,本发明专利的保护范围应以所附权利要求为准。
Claims (4)
1.一种弥散张量成像方法,其特征在于,包括以下步骤:
在K空间在读出方向上中间区域进行划分,将K空间读出方向上中间区域划分成多个相互独立的子空间;
在K空间读出方向上对中部区域的所述多个相互独立的子空间进行激发,并对K空间在相位编码方向上进行采样,得到K空间的数据;对K空间在相位编码方向上进行采样的方式为欠采样;所述欠采样具体为变密度欠采样;所述变密度欠采样根据自适应函数,对K空间相位编码方向上进行欠采样;
对所述K空间数据进行重建得到重建图像。
2.根据权利要求1所述的弥散张量成像方法,其特征在于,所述对所述K空间数据进行重建得到重建图像的步骤具体为:
对所述K空间数据进行压缩感知重建得到重建图像。
3.一种弥散张量成像系统,其特征在于,包括:
处理模块,用于在K空间在读出方向上中间区域进行划分,将K空间读出方向上中间区域划分成多个相互独立的子空间;
采样模块,与所述处理模块连接,所述采样模块用于在K空间读出方向上对中部区域的所述多个相互独立的子空间进行激发,并对K空间在相位编码方向上进行采样,得到K空间的数据;对K空间在相位编码方向上进行采样的方式为欠采样;所述欠采样具体为变密度欠采样;所述变密度欠采样根据自适应函数,对K空间相位编码方向上进行欠采样;及
重建模块,与所述采样模块连接,用于对所述K空间数据进行重建得到重建图像。
4.根据权利要求3所述的弥散张量成像系统,其特征在于,所述重建模块对所述K空间数据进行压缩感知重建得到重建图像。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110387526.6A CN102488519B (zh) | 2011-11-29 | 2011-11-29 | 弥散张量成像方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110387526.6A CN102488519B (zh) | 2011-11-29 | 2011-11-29 | 弥散张量成像方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102488519A CN102488519A (zh) | 2012-06-13 |
CN102488519B true CN102488519B (zh) | 2014-07-16 |
Family
ID=46180430
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201110387526.6A Active CN102488519B (zh) | 2011-11-29 | 2011-11-29 | 弥散张量成像方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102488519B (zh) |
Families Citing this family (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103163496A (zh) * | 2011-12-12 | 2013-06-19 | 中国科学院深圳先进技术研究院 | 平面回波成像方法及系统 |
CN103529414B (zh) * | 2012-07-04 | 2016-06-29 | 上海联影医疗科技有限公司 | 磁共振成像方法与装置、k空间的重建方法与装置 |
CN103529413B (zh) * | 2012-07-04 | 2016-03-02 | 上海联影医疗科技有限公司 | 磁共振成像方法与装置、k空间的重建方法与装置 |
KR101663229B1 (ko) * | 2013-10-23 | 2016-10-06 | 삼성전자주식회사 | 자기 공명 영상 장치 및 그에 따른 자기 공명 영상의 이미징 방법 |
EP2866045A1 (en) | 2013-10-23 | 2015-04-29 | Samsung Electronics Co., Ltd | Magnetic resonance imaging apparatus and method |
CN105259525B (zh) * | 2015-10-28 | 2019-08-13 | 北京大学 | 基于邻域共享压缩感知的动态对比增强磁共振快速成像方法 |
CN108459289B (zh) * | 2018-01-30 | 2021-03-12 | 奥泰医疗系统有限责任公司 | 一种基于数据一致性的多次激发弥散加权磁共振成像方法 |
CN109358303B (zh) | 2018-12-07 | 2021-04-02 | 上海联影医疗科技股份有限公司 | 磁共振成像的扫描方法和弥散张量成像系统 |
CN109738840B (zh) * | 2018-12-29 | 2022-06-14 | 佛山瑞加图医疗科技有限公司 | 一种磁共振成像系统和方法 |
US10672151B1 (en) * | 2019-01-07 | 2020-06-02 | Uih America, Inc. | Systems and methods for magnetic resonance image reconstruction |
CN112782629B (zh) * | 2020-12-30 | 2022-05-03 | 深圳市联影高端医疗装备创新研究院 | 一种磁共振扫描控制方法、存储介质及系统 |
US11965947B2 (en) * | 2021-11-23 | 2024-04-23 | Shanghai United Imaging Intelligence Co., Ltd. | Multiplex MRI image reconstruction |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1102320A (zh) * | 1993-07-09 | 1995-05-10 | 菲利浦电子有限公司 | 梯度和自旋回波磁共振成像中的重影抑制 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP5569714B2 (ja) * | 2009-06-16 | 2014-08-13 | 株式会社国際電気通信基礎技術研究所 | 磁気共鳴画像化装置および磁気共鳴画像化方法 |
JP5023127B2 (ja) * | 2009-09-30 | 2012-09-12 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | 磁気共鳴イメージング装置 |
JP5613065B2 (ja) * | 2010-01-21 | 2014-10-22 | 株式会社東芝 | 磁気共鳴イメージング装置 |
-
2011
- 2011-11-29 CN CN201110387526.6A patent/CN102488519B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1102320A (zh) * | 1993-07-09 | 1995-05-10 | 菲利浦电子有限公司 | 梯度和自旋回波磁共振成像中的重影抑制 |
Non-Patent Citations (2)
Title |
---|
JP特开2011-167509A 2011.09.01 |
JP特开2011-184A 2011.01.06 |
Also Published As
Publication number | Publication date |
---|---|
CN102488519A (zh) | 2012-06-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102488519B (zh) | 弥散张量成像方法及系统 | |
US8878535B2 (en) | Method and apparatus for capturing magnetic resonance image | |
Summers et al. | A preliminary study of the effects of trigger timing on diffusion tensor imaging of the human spinal cord | |
US10852376B2 (en) | Magnetic resonance imaging method and device | |
CN1683939A (zh) | 在磁共振断层造影中融合加权成像的运动校正的多点方法 | |
CN103477238A (zh) | 采用来自先验采集的约束的压缩感测mr图像重建 | |
US10684345B2 (en) | Reconstructing magnetic resonance images for contrasts | |
US10219787B2 (en) | Respiratory mode (“R-Mode”)—acquisition and display of cardiovascular images to show respiratory effects | |
CN106443533B (zh) | 基于多次激发的导航磁共振扩散成像方法及装置 | |
CN105572613B (zh) | 磁共振化学交换饱和转移成像方法和系统 | |
CN102772207B (zh) | 磁共振成像方法和装置 | |
CN110346743A (zh) | 一种磁共振弥散加权成像方法和装置 | |
CN103033782B (zh) | 并行磁共振成像装置及其成像的方法 | |
CN103163496A (zh) | 平面回波成像方法及系统 | |
CN102890255B (zh) | 一种用于缩短回波时间的由内而外平面回波成像方法 | |
CN105005012A (zh) | 基于压缩感知的腹部器官动态对比增强磁共振成像方法 | |
Matsuda et al. | Imaging of a large collection of human embryo using a super-parallel MR microscope | |
US10267881B2 (en) | Multi-echo parametric variation saturation (MePaVARS CEST imaging and other MR imaging | |
CN205251561U (zh) | 一种嵌入型可分离式pet/mr扫描仪 | |
CN106772168B (zh) | 磁共振成像方法及装置 | |
CN101509964B (zh) | 校正磁共振成像系统中主磁场不稳定的方法 | |
CN103048632B (zh) | 一种基于压缩传感的快速高角分辨率扩散成像方法 | |
CN102727206A (zh) | 颅脑、颈动脉和主动脉三位一体的扫描方法及扫描系统 | |
CN106137198A (zh) | 一种磁共振成像方法及装置 | |
CN106646299B (zh) | 梯度回波三回波水脂分离方法及应用它的磁共振成像系统 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |