CN114325528B - 一种磁共振成像方法及相关设备 - Google Patents
一种磁共振成像方法及相关设备 Download PDFInfo
- Publication number
- CN114325528B CN114325528B CN202111605754.6A CN202111605754A CN114325528B CN 114325528 B CN114325528 B CN 114325528B CN 202111605754 A CN202111605754 A CN 202111605754A CN 114325528 B CN114325528 B CN 114325528B
- Authority
- CN
- China
- Prior art keywords
- flip angle
- image
- constructing
- flip
- data corresponding
- 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
- 238000000034 method Methods 0.000 title claims abstract description 32
- 238000002595 magnetic resonance imaging Methods 0.000 title claims abstract description 28
- 238000005070 sampling Methods 0.000 claims abstract description 30
- 238000003384 imaging method Methods 0.000 claims description 11
- 238000004590 computer program Methods 0.000 claims description 7
- 239000011159 matrix material Substances 0.000 claims description 6
- 238000005259 measurement Methods 0.000 claims description 6
- 238000005065 mining Methods 0.000 claims description 6
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 230000001133 acceleration Effects 0.000 abstract description 6
- 238000010586 diagram Methods 0.000 description 3
- 230000009977 dual effect Effects 0.000 description 3
- 238000013507 mapping Methods 0.000 description 3
- 230000003287 optical effect Effects 0.000 description 3
- 230000000644 propagated effect Effects 0.000 description 2
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000000835 fiber Substances 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 239000013307 optical fiber Substances 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 230000008707 rearrangement Effects 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000003786 synthesis reaction Methods 0.000 description 1
Landscapes
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
本发明实施例公开了一种磁共振成像方法及相关设备，该方法包括：在至少两个预设翻转角下获取待成像样本的至少两个扫描数据；基于对应所述第一翻转角的扫描数据构建第一翻转角图像；根据对应所述第一翻转角的扫描数据内的中心采样数据和对应所述第二翻转角的扫描数据内的中心采样数据构建低频相位差图像；基于所述低频相位差图像和对应所述第二翻转角的扫描数据构建校正第二翻转角相位图像；根据所述第一翻转角图像和所述第二翻转角图像构建所述待成像样本的磁共振图像。通过本发明实施例的技术方案，采用多翻转角之间的依赖关系，使得多翻转角的采集可以在相位约束条件下获得更高的加速倍数，提高了磁共振成像的速度。
Description
技术领域
本发明实施例涉及成像技术领域，尤其涉及一种磁共振成像方法及相关设备。
背景技术
磁共振由于它的无辐射、多对比度以及可获得多个方位图像的优势，在临床中被广泛使用。定量磁共振成像(MRI)成像由于提供了物理量的绝对值，可以作为生物标志物，在临床应用中具有很好的前景，在多中心研究和患者疾病跟踪随访中也非常重要。为了获得定量图像，磁共振技术中经常采用多个翻转角的序列来获取数据，比如采用2个以上的翻转角的数据来获得T1 mapping定量图，结合多回波采集，又可以获得T2 mapping定量图或者T2*mapping。因为不同的翻转角采集都需要单独扫描，这样会翻倍增加扫描时间。压缩感知(Compressed Sensing:CS)是一种用来加速磁共振扫描速度的有效方法，因为它采用稀疏采样的K空间轨迹，大大减少采样点数。而传统加速方法都是采用每个翻转角采集进行独立加速的方法，并未使用他们之间的依赖关系，易在图像上产生严重卷叠伪影的问题，加速效果不佳。
发明内容
本发明实施例提供了一种磁共振成像方法及相关设备，避免了在采用磁共振成像时，成像速度慢，图像上产生严重卷叠伪影的问题。
第一方面，本发明实施例提供了一种磁共振成像方法，包括：
在至少两个预设翻转角下获取待成像样本的至少两个扫描数据，其中，所述两个预设翻转角包括具有先后顺序的第一翻转角和第二翻转角，所述扫描数据包括对待成像样本上预设中心区域进行全采样获取的中心数据和对待成像样本上预设中心区域外进行稀疏采样获取的边缘数据；
基于对应所述第一翻转角的扫描数据构建第一翻转角图像；
根据对应所述第一翻转角的中心数据和对应所述第二翻转角的中心数据构建低频相位差图像；
基于所述低频相位差图像和对应所述第二翻转角的扫描数据构建第二翻转角图像；
根据所述第一翻转角图像和所述第二翻转角图像构建所述待成像样本的磁共振图像。
可选的，所述在至少两个预设翻转角下获取待成像样本的至少两个扫描数据的步骤，包括：
通过所述至少两个预设翻转角的梯度回波序列获取待成像样本的至少两个扫描数据。
可选的，所述基于对应所述第一翻转角的扫描数据构建第一翻转角图像的步骤，包括：
构建带有稀疏域约束的目标压缩感知模型；
基于所述目标压缩感知模型求得所述第一翻转角图像。
可选的，所述根据对应所述第一翻转角的中心数据和对应所述第二翻转角的中心数据构建低频相位差图像的步骤，包括：
基于所述第一翻转角的中心数据的复数图像数据求出第一低频相位；
基于所述第二翻转角的中心数据的复数图像数据求出第二低频相位；
根据所述第一低频相位和所述第二低频相位得到所述低频相位差图像。
可选的，基于所述低频相位差图像和对应所述第二翻转角的扫描数据构建第二翻转角图像的步骤，包括：
基于对应所述第二翻转角的扫描数据，构建第二翻转角的相位图像，根据所述低频相位差图像得到两翻转角的相位差图像；
基于所述低频相位差图像校正所述目标压缩感知模型，得到校正后的目标压缩感知模型；
根据所述校正后的目标压缩感知模型和对应所述第二翻转角的扫描数据生成第二翻转角图像。
可选的，所述根据所述第一翻转角图像和所述第二翻转角图像构建所述待成像样本的磁共振图像的步骤，包括：
根据所述第一翻转角图像和所述第二翻转角图像生成所述待成像样本的定量图像；
基于所述定量图像构建所述待成像样本的磁共振图像。
可选的，所述构建带有稀疏域约束的目标压缩感知模型,基于所述目标压缩感知模型求得所述第一翻转角图像的步骤，包括：
构建带有稀疏域约束的目标压缩感知模型为：
基于所述目标压缩感知求得所述第一翻转角图像，其中，x为全采样的目标K空间，用来重建生成需要的图像；F为欠采样测量矩阵；y为欠采样的K空间数据，Ψ为稀疏变换域；G为SPIRiT的K空间卷积核,由中心满采区域求得；λ1，λ2为正则化因子。
可选的，所述基于所述低频相位差图像校正所述目标压缩感知模型，得到校正后的目标压缩感知模型的步骤，包括：
基于所述低频相位差图像校正所述目标压缩感知模型，得到校正后的目标压缩感知模型为：
其中，x为全采样的目标K空间，用来重建生成需要的图像；F为欠采样测量矩阵；y为欠采样的K空间数据，Ψ为稀疏变换域，G为SPIRiT的K空间卷积核，由中心满采区域求得；λ1，λ2为正则化因子；为进行校正后的相位图像。
第二方面，本申请实施例提供了一种磁共振成像装置，包括：
数据获取模块，用于在至少两个预设翻转角下获取待成像样本的至少两个扫描数据，其中，所述两个预设翻转角包括具有先后顺序的第一翻转角和第二翻转角，所述扫描数据包括对待成像样本上预设中心区域进行全采样获取的中心数据和对待成像样本上预设中心区域外进行稀疏采样获取的边缘数据；
转换模块，用于基于对应所述第一翻转角的扫描数据构建第一翻转角图像；
计算模块，用于根据对应所述第一翻转角的中心数据和对应所述第二翻转角的中心数据构建低频相位差图像；
校正模块，用于基于所述低频相位差图像和对应所述第二翻转角的扫描数据构建第二翻转角图像；
成像模块，用于根据所述第一翻转角图像和所述第二翻转角图像构建所述待成像样本的磁共振图像。
第三方面，本申请实施例提供了一种计算机可读存储介质，其上存储有计算机程序，其特征在于：所述计算机程序被处理器执行时实现如上述的磁共振成像方法。
上述发明中的实施例具有如下优点或有益效果：
通过在至少两个预设翻转角下获取待成像样本的至少两个扫描数据，其中，所述两个预设翻转角包括具有先后顺序的第一翻转角和第二翻转角，所述扫描数据包括对待成像样本上预设中心区域进行全采样获取的中心数据和对待成像样本上预设中心区域外进行稀疏采样获取的边缘数据；基于对应所述第一翻转角的扫描数据构建第一翻转角图像；根据对应所述第一翻转角的中心数据和对应所述第二翻转角的中心数据构建低频相位差图像；基于所述低频相位差图像和对应所述第二翻转角的扫描数据构建第二翻转角图像；根据所述第一翻转角图像和所述第二翻转角图像构建所述待成像样本的磁共振图像。利用不同翻转角间信号相位一致性的特性，在压缩感知模型中引入相位约束来提高多翻转角成像的扫描速度，其它翻转角采集数据时可以加大欠采样的倍数，进而减少采样时间，在相位约束条件下获得更高的加速倍数，从而提高多翻转角成像的扫描速度。
附图说明
图1是本申请实施例提供的一种磁共振成像方法的流程图；
图2是本申请实施例提供的一种磁共振成像装置的结构示意图；
图3是本申请实施例提供的一种存储介质的结构示意图；
图4是本申请实施例提供的一种磁共振成像的运用场景图。
具体实施方式
下面结合附图和实施例对本发明作进一步的详细说明。可以理解的是，此处所描述的具体实施例仅仅用于解释本发明，而非对本发明的限定。另外还需要说明的是，为了便于描述，附图中仅示出了与本发明相关的部分而非全部结构。
图1为本申请实施例提供的一种磁共振成像方法的流程图，本实施例提供的磁共振成像方法可适用于同TE不同翻转角扫描的其它场景中，比如射频功率发射图的采集中。该方法可以由磁共振成像装置来执行，该装置可以由软件和/或硬件的方式来实现，通常集成于磁共振成像设备。
如图1所示，本发明实施例提供了一种磁共振成像方法，包括：
S101、在至少两个预设翻转角下获取待成像样本的至少两个扫描数据，其中，所述两个预设翻转角包括具有先后顺序的第一翻转角和第二翻转角，所述扫描数据包括对待成像样本上预设中心区域进行全采样获取的中心数据和对待成像样本上预设中心区域外进行稀疏采样获取的边缘数据，所述构建过程为重建过程；
在一种可能的实施方式中，所述在至少两个预设翻转角下获取待成像样本的至少两个扫描数据的步骤，包括：
通过所述至少两个预设翻转角的梯度回波序列获取待成像样本的至少两个扫描数据。
示例性的，分别采用至少两个翻转角的梯度回波序列采集数据，每个翻转角的序列可以是单回波或多回波。以2个翻转角为例进行说明，对应翻转角分别是FA1和FA2。采样轨迹采用稀疏采样的方式，其中心区域要进行全采样。典型的采样轨迹是Possion Disk的采样方式。
S102、基于对应所述第一翻转角的扫描数据构建第一翻转角图像；
在一种可能的实施方式中，所述基于对应所述第一翻转角的扫描数据构建第一翻转角图像的步骤，包括：
构建带有稀疏域约束的目标压缩感知模型；
基于所述目标压缩感知模型求得所述第一翻转角图像。
在一种可能的实施方式中，所述构建带有稀疏域约束的目标压缩感知模型,基于所述目标压缩感知模型求得所述第一翻转角图像的步骤，包括：
构建带有稀疏域约束的目标压缩感知模型为：
基于所述目标压缩感知求得所述第一翻转角图像，其中，x为全采样的目标K空间，用来重建生成需要的图像；F为欠采样测量矩阵；y为欠采样的K空间数据，Ψ为稀疏变换域，比如TV变换，小波变换等；G为SPIRiT的K空间卷积核,由中心满采区域求得；λ1，λ2为正则化因子。
示例性的，第一翻转角FA1采集数据的重建采用的重建模型为：
或者
模型
可以用来重建单通道数据，或者对于多通道数据可以逐个重建每个通道的数据然后进行基于相位的通道合成，进而求得的相位图像为Ph1。
S103、根据对应所述第一翻转角的中心数据和对应所述第二翻转角的中心数据构建低频相位差图像；
示例性的，计算FA1的低频相位图像PhL1时，PhL1可以对中心全采样区域的数据进行傅里叶变换得到复数数据，取数据的相位得到。也可以用对FA1重建的相位图像进行低通滤波得到。
S104、基于所述低频相位差图像和对应所述第二翻转角的扫描数据构建第二翻转角图像；
在一种可能的实施方式中，所述根据对应所述第一翻转角的中心数据和对应所述第二翻转角的中心数据构建低频相位差图像的步骤，包括：
基于所述第一翻转角的中心数据的复数图像数据求出第一低频相位；
基于所述第二翻转角的中心数据的复数图像数据求出第二低频相位；
根据所述第一低频相位和所述第二低频相位得到所述低频相位差图像。
在一种可能的实施方式中，所述基于所述低频相位差图像和对应所述第二翻转角的扫描数据构建第二翻转角图像的步骤，包括：
基于所述低频相位差图像校正所述目标压缩感知模型，得到校正后的目标压缩感知模型；
根据所述校正后的目标压缩感知模型和对应所述第二翻转角的扫描数据生成第二翻转角图像。
在一种可能的实施方式中，所述基于所述低频相位差图像校正所述目标压缩感知模型，得到校正后的目标压缩感知模型的步骤，包括：
基于所述低频相位差图像校正所述目标压缩感知模型，得到校正后的目标压缩感知模型为：
其中，x为全采样的目标K空间，用来重建生成需要的图像；F为欠采样测量矩阵；y为欠采样的K空间数据，Ψ为稀疏变换域，G为SPIRiT的K空间卷积核，由中心满采区域求得；λ1，λ2为正则化因子；为进行校正后的相位图像。
示例性的，采用同上述计算PhL1的方法计算FA2的低频相位图像PhL2，计算PhL2与PhL1的差图像PhC，如图4所示，用PhC校正Ph1得到校正后的相位Ph1和校正后的/>以及Ph2：
将相位作为约束项对第二翻转角FA2的数据进行重建，得到：
FA2的低频相位图像PhL2。
示例性的，采用一阶原始对偶方法将原模型分解为光滑项和非光滑项进行求解。其中l1范数项时非平滑项。假定E＝|.|1表示l1范数，并且E*为E的共轭函数，v1，v2分别是上述模型中两个l1范数的对偶变量。将上述模型转换成一个原始对偶问题，如下所示：
S105、根据所述第一翻转角图像和所述第二翻转角图像构建所述待成像样本的磁共振图像。
在一种可能的实施方式中，所述根据所述第一翻转角图像和所述第二翻转角图像构建所述待成像样本的磁共振图像的步骤，包括：
根据所述第一翻转角图像和所述第二翻转角图像生成所述待成像样本的定量图像；
基于所述定量图像构建所述待成像样本的磁共振图像。
通过在至少两个预设翻转角下获取待成像样本的至少两个扫描数据，其中，所述两个预设翻转角包括具有先后顺序的第一翻转角和第二翻转角，所述扫描数据包括对待成像样本上预设中心区域进行全采样获取的中心数据和对待成像样本上预设中心区域外进行稀疏采样获取的边缘数据；基于对应所述第一翻转角的扫描数据构建第一翻转角图像；根据对应所述第一翻转角的中心数据和对应所述第二翻转角的中心数据构建低频相位差图像；基于所述低频相位差图像和对应所述第二翻转角的扫描数据构建第二翻转角图像；根据所述第一翻转角图像和所述第二翻转角图像构建所述待成像样本的磁共振图像。利用不同翻转角间信号相位一致性的特性，在压缩感知模型中引入相位约束来提高多翻转角成像的扫描速度，其它翻转角采集数据时可以加大欠采样的倍数，进而减少采样时间，在相位约束条件下获得更高的加速倍数，从而提高多翻转角成像的扫描速度，增加了先验的相位约束项，在高加速倍下仍然能够保证较好的图像质量。
在一种可能的实施方式中，本申请实施例提供了一种磁共振成像装置，包括：
数据获取模块201，用于在至少两个预设翻转角下获取待成像样本的至少两个扫描数据，其中，所述两个预设翻转角包括具有先后顺序的第一翻转角和第二翻转角，所述扫描数据包括对待成像样本上预设中心区域进行全采样获取的中心数据和对待成像样本上预设中心区域外进行稀疏采样获取的边缘数据；
转换模块202，用于基于对应所述第一翻转角的扫描数据构建第一翻转角图像；
计算模块203，用于根据对应所述第一翻转角的中心数据和对应所述第二翻转角的中心数据构建低频相位差图像；
校正模块204，用于基于所述低频相位差图像和对应所述第二翻转角的扫描数据构建第二翻转角图像；
成像模块205，用于根据所述第一翻转角图像和所述第二翻转角图像构建所述待成像样本的磁共振图像。
在一种可能的实施方式中，本申请实施例提供了一种计算机可读存储介质400，其上存储有计算机程序411，其特征在于：所述计算机程序411被处理器执行时实现如上述的磁共振成像方法。
本发明实施例的计算机存储介质，可以采用一个或多个计算机可读的介质的任意组合。计算机可读介质可以是计算机可读信号介质或者计算机可读存储介质。计算机可读存储介质例如可以是但不限于：电、磁、光、电磁、红外线、或半导体的***、装置或器件，或者任意以上的组合。计算机可读存储介质的更具体的例子(非穷举的列表)包括：具有一个或多个导线的电连接、便携式计算机磁盘、硬盘、随机存取存储器(RAM)、只读存储器(ROM)、可擦式可编程只读存储器(EPROM或闪存)、光纤、便携式紧凑磁盘只读存储器(CD-ROM)、光存储器件、磁存储器件、或者上述的任意合适的组合。在本文件中，计算机可读存储介质可以是任何包含或存储程序的有形介质，该程序可以被指令执行***、装置或者器件使用或者与其结合使用。
计算机可读的信号介质可以包括在基带中或者作为载波一部分传播的数据信号，其中承载了计算机可读的程序代码。这种传播的数据信号可以采用多种形式，包括但不限于电磁信号、光信号或上述的任意合适的组合。计算机可读的信号介质还可以是计算机可读存储介质以外的任何计算机可读介质，该计算机可读介质可以发送、传播或者传输用于由指令执行***、装置或者器件使用或者与其结合使用的程序。
计算机可读介质上包含的程序代码可以用任何适当的介质传输，包括但不限于：无线、电线、光缆、RF等等，或者上述的任意合适的组合。
可以以一种或多种程序设计语言或其组合来编写用于执行本发明操作的计算机程序代码，上述程序设计语言包括面向对象的程序设计语言，诸如Java、Smalltalk、C++，还包括常规的过程式程序设计语言—诸如“C”语言或类似的程序设计语言。程序代码可以完全地在用户计算机上执行、部分地在用户计算机上执行、作为一个独立的软件包执行、部分在用户计算机上部分在远程计算机上执行、或者完全在远程计算机或服务器上执行。在涉及远程计算机的情形中，远程计算机可以通过任意种类的网络，包括局域网(LAN)或广域网(WAN)，连接到用户计算机，或者，可以连接到外部计算机(例如利用因特网服务提供商来通过因特网连接)。
本领域普通技术人员应该明白，上述的本发明的各模块或各步骤可以用通用的计算装置来实现，它们可以集中在单个计算装置上，或者分布在多个计算装置所组成的网络上，可选地，他们可以用计算机装置可执行的程序代码来实现，从而可以将它们存储在存储装置中由计算装置来执行，或者将它们分别制作成各个集成电路模块，或者将它们中的多个模块或步骤制作成单个集成电路模块来实现。这样，本发明不限制于任何特定的硬件和软件的结合。
注意，上述仅为本发明的较佳实施例及所运用技术原理。本领域技术人员会理解，本发明不限于这里上述的特定实施例，对本领域技术人员来说能够进行各种明显的变化、重新调整和替代而不会脱离本发明的保护范围。因此，虽然通过以上实施例对本发明进行了较为详细的说明，但是本发明不仅仅限于以上实施例，在不脱离本发明构思的情况下，还可以包括更多其他等效实施例，而本发明的范围由所附的权利要求范围决定。
Claims (8)
1.一种磁共振成像方法，其特征在于，包括：
在至少两个预设翻转角下获取待成像样本的至少两个扫描数据，其中，所述两个预设翻转角包括具有先后顺序的第一翻转角和第二翻转角，所述扫描数据包括对待成像样本上预设中心区域进行全采样获取的中心数据和对待成像样本上预设中心区域外进行稀疏采样获取的边缘数据；
基于对应所述第一翻转角的扫描数据构建第一翻转角图像；
根据对应所述第一翻转角的中心数据和对应所述第二翻转角的中心数据构建低频相位差图像；
基于所述低频相位差图像和对应所述第二翻转角的扫描数据构建第二翻转角图像；
根据所述第一翻转角图像和所述第二翻转角图像构建所述待成像样本的磁共振图像；
所述基于对应所述第一翻转角的扫描数据构建第一翻转角图像的步骤，包括：
构建带有稀疏域约束的目标压缩感知模型；
基于所述目标压缩感知模型求得所述第一翻转角图像；
基于所述低频相位差图像和对应所述第二翻转角的扫描数据构建第二翻转角图像的步骤，包括：
基于对应所述第二翻转角的扫描数据，构建第二翻转角的相位图像，根据所述低频相位差图像得到两翻转角的相位差图像；
基于所述低频相位差图像校正所述目标压缩感知模型，得到校正后的目标压缩感知模型；
根据所述校正后的目标压缩感知模型和对应所述第二翻转角的扫描数据生成第二翻转角图像。
2.根据权利要求1所述的磁共振成像方法，其特征在于，所述在至少两个预设翻转角下获取待成像样本的至少两个扫描数据的步骤，包括：
通过所述至少两个预设翻转角的梯度回波序列获取待成像样本的至少两个扫描数据。
3.根据权利要求1所述的磁共振成像方法，其特征在于，所述根据对应所述第一翻转角的中心数据和对应所述第二翻转角的中心数据构建低频相位差图像的步骤，包括：
基于所述第一翻转角的中心数据的复数图像数据求出第一低频相位；
基于所述第二翻转角的中心数据的复数图像数据求出第二低频相位；
根据所述第一低频相位和所述第二低频相位得到所述低频相位差图像。
4.根据权利要求1所述的磁共振成像方法，其特征在于，所述根据所述第一翻转角图像和所述第二翻转角图像构建所述待成像样本的磁共振图像的步骤，包括：
根据所述第一翻转角图像和所述第二翻转角图像生成所述待成像样本的定量图像；
基于所述定量图像构建所述待成像样本的磁共振图像。
5.根据权利要求1所述的磁共振成像方法，其特征在于，所述构建带有稀疏域约束的目标压缩感知模型,基于所述目标压缩感知模型求得所述第一翻转角图像的步骤，包括：
构建带有稀疏域约束的目标压缩感知模型为：
基于所述目标压缩感知求得所述第一翻转角图像，其中，x为全采样的目标K空间，用来重建生成需要的图像；F为欠采样测量矩阵；y为欠采样的K空间数据，Ψ为稀疏变换域；G为SPIRiT的K空间卷积核,由中心满采区域求得；λ1，λ2为正则化因子。
6.根据权利要求1所述的磁共振成像方法，其特征在于，所述基于所述低频相位差图像校正所述目标压缩感知模型，得到校正后的目标压缩感知模型的步骤，包括：
基于所述低频相位差图像校正所述目标压缩感知模型，得到校正后的目标压缩感知模型为：
其中，x为全采样的目标K空间，用来重建生成需要的图像；F为欠采样测量矩阵；y为欠采样的K空间数据，Ψ为稀疏变换域，G为SPIRiT的K空间卷积核，由中心满采区域求得；λ1，λ2,λ3都为正则化因子；为进行校正后的相位图像。
7.一种基于权利要求1所述的磁共振成像方法的成像装置，其特征在于，包括：
数据获取模块，用于在至少两个预设翻转角下获取待成像样本的至少两个扫描数据，其中，所述两个预设翻转角包括具有先后顺序的第一翻转角和第二翻转角，所述扫描数据包括对待成像样本上预设中心区域进行全采样获取的中心数据和对待成像样本上预设中心区域外进行稀疏采样获取的边缘数据；
转换模块，用于基于对应所述第一翻转角的扫描数据构建第一翻转角图像；
计算模块，用于根据对应所述第一翻转角的中心数据和对应所述第二翻转角的中心数据构建低频相位差图像；
校正模块，用于基于所述低频相位差图像和对应所述第二翻转角的扫描数据构建第二翻转角图像；
成像模块，用于根据所述第一翻转角图像和所述第二翻转角图像构建所述待成像样本的磁共振图像。
8.一种计算机可读存储介质，其上存储有计算机程序，其特征在于：所述计算机程序被处理器执行时实现如权利要求1至6中任一项所述的磁共振成像方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111605754.6A CN114325528B (zh) | 2021-12-25 | 2021-12-25 | 一种磁共振成像方法及相关设备 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111605754.6A CN114325528B (zh) | 2021-12-25 | 2021-12-25 | 一种磁共振成像方法及相关设备 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114325528A CN114325528A (zh) | 2022-04-12 |
| CN114325528B true CN114325528B (zh) | 2024-01-02 |
Family
ID=81013829
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111605754.6A Active CN114325528B (zh) | 2021-12-25 | 2021-12-25 | 一种磁共振成像方法及相关设备 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114325528B (zh) |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103349550A (zh) * | 2013-07-04 | 2013-10-16 | 华东师范大学 | 一种磁共振成像扫描与压缩传感重建一体化方法及其装置 |
| CN104569880A (zh) * | 2014-12-31 | 2015-04-29 | 中国科学院深圳先进技术研究院 | 一种磁共振快速成像方法及*** |
| CN108514415A (zh) * | 2018-02-09 | 2018-09-11 | 奥泰医疗***有限责任公司 | 一种快速磁敏感加权成像扫描序列及方法 |
| CN108535675A (zh) * | 2018-04-08 | 2018-09-14 | 朱高杰 | 一种基于深度学习和数据自洽的磁共振多通道重建方法 |
| CN111263896A (zh) * | 2017-08-22 | 2020-06-09 | 皇家飞利浦有限公司 | 在磁共振成像***中对相位相关的伪影的数据驱动的校正 |
| CN112272535A (zh) * | 2019-03-29 | 2021-01-26 | 上海联影医疗科技股份有限公司 | 磁共振弥散成像的方法和*** |
| CN112370040A (zh) * | 2020-11-13 | 2021-02-19 | 上海东软医疗科技有限公司 | 磁共振成像方法、装置、存储介质及电子设备 |
| CN113391251A (zh) * | 2020-03-12 | 2021-09-14 | 上海联影医疗科技股份有限公司 | 磁共振图像重建方法、装置和设备 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US10782378B2 (en) * | 2018-05-16 | 2020-09-22 | Siemens Healthcare Gmbh | Deep learning reconstruction of free breathing perfusion |
-
2021
- 2021-12-25 CN CN202111605754.6A patent/CN114325528B/zh active Active
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103349550A (zh) * | 2013-07-04 | 2013-10-16 | 华东师范大学 | 一种磁共振成像扫描与压缩传感重建一体化方法及其装置 |
| CN104569880A (zh) * | 2014-12-31 | 2015-04-29 | 中国科学院深圳先进技术研究院 | 一种磁共振快速成像方法及*** |
| CN111263896A (zh) * | 2017-08-22 | 2020-06-09 | 皇家飞利浦有限公司 | 在磁共振成像***中对相位相关的伪影的数据驱动的校正 |
| CN108514415A (zh) * | 2018-02-09 | 2018-09-11 | 奥泰医疗***有限责任公司 | 一种快速磁敏感加权成像扫描序列及方法 |
| CN108535675A (zh) * | 2018-04-08 | 2018-09-14 | 朱高杰 | 一种基于深度学习和数据自洽的磁共振多通道重建方法 |
| CN112272535A (zh) * | 2019-03-29 | 2021-01-26 | 上海联影医疗科技股份有限公司 | 磁共振弥散成像的方法和*** |
| CN113391251A (zh) * | 2020-03-12 | 2021-09-14 | 上海联影医疗科技股份有限公司 | 磁共振图像重建方法、装置和设备 |
| CN112370040A (zh) * | 2020-11-13 | 2021-02-19 | 上海东软医疗科技有限公司 | 磁共振成像方法、装置、存储介质及电子设备 |
Non-Patent Citations (2)
| Title |
|---|
| 压缩感知技术在磁共振成像技术中的应用进展分析;肖顺; 储呈晨; 王源冰; 李斌;中国医疗设备;第36卷(第11期);4-9 * |
| 探讨用磁敏感定量成像和T2*加权成像方法对地中海贫血病患者脑铁含量的检测研究;郭红宇，王雁军，黄萍，郭华;Science Discovery;第4卷(第6期);405-409 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114325528A (zh) | 2022-04-12 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Bustin et al. | High‐dimensionality undersampled patch‐based reconstruction (HD‐PROST) for accelerated multi‐contrast MRI | |
| US10671939B2 (en) | System, method and computer-accessible medium for learning an optimized variational network for medical image reconstruction | |
| US10915990B2 (en) | Systems and methods for denoising medical images with deep learning network | |
| CN108375746B (zh) | 一种相位反卷折方法与设备 | |
| US9482732B2 (en) | MRI reconstruction with motion-dependent regularization | |
| US8879811B2 (en) | Alternating direction of multipliers method for parallel MRI reconstruction | |
| JP5926285B2 (ja) | 任意ｋ空間トラジェクトリの高速並行再構成 | |
| US10274567B2 (en) | Magnetic resonance imaging apparatus and image processing apparatus | |
| CN110133556B (zh) | 一种磁共振图像处理方法、装置、设备及存储介质 | |
| Malavé et al. | Reconstruction of undersampled 3D non‐Cartesian image‐based navigators for coronary MRA using an unrolled deep learning model | |
| CN105143904A (zh) | 用于多信道发射机的多频带射频/磁共振成像脉冲设计 | |
| Maier et al. | CG‐SENSE revisited: Results from the first ISMRM reproducibility challenge | |
| US9366741B2 (en) | Medical image imaging method, medical diagnostic apparatus using the same, and recording medium therefor | |
| US20190033409A1 (en) | Reconstructing magnetic resonance images with different contrasts | |
| US20200209332A1 (en) | System and method for magnetic resonance imaging | |
| CN113298902B (zh) | 卷褶视野磁共振图像的重建方法、计算机设备及存储介质 | |
| CN113298901B (zh) | 卷褶视野磁共振图像的重建方法、计算机设备及存储介质 | |
| CN114325528B (zh) | 一种磁共振成像方法及相关设备 | |
| CN111537930B (zh) | 基于梯度波形调整的磁共振参数成像方法与设备 | |
| US11119173B2 (en) | Dynamic imaging based on echo planar imaging sequence | |
| Yu et al. | Universal generative modeling in dual domains for dynamic MRI | |
| JP6809870B2 (ja) | 磁気共鳴イメージング装置、および、信号処理方法 | |
| US20230309849A1 (en) | Method for generating a magnetic resonance image | |
| US11893662B2 (en) | Density compensation function in filtered backprojection | |
| US20240219501A1 (en) | System and Method for Distortion and Motion Artifact-Free Diffusion Imaging Using Single-Shot Diffusion-Prepared Turbo-Spin-Echo Sequence With Spiral-Ring Readouts and Magnitude Stabilizers |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |