磁共振图像的处理方法、装置及电子设备制造方法及图纸

本发明专利技术实施例提供一种磁共振图像的处理方法、装置及电子设备。本发明专利技术实施例通过在预设方向上，将第一复数磁共振图像对应的第一K空间信号矩阵的外推空间的信号值设置为预设值，得到第二K空间信号矩阵，对第二K空间信号矩阵对应的第二复数磁共振图像的幅值图像进行图像域滤波，得到新幅值图像，利用第二复数磁共振图像的相位对新幅值图像进行相位调制，得到第三复数磁共振图像，将第一K空间信号矩阵的信号值填充到第三复数磁共振图像对应的K空间信号矩阵中的相应位置，得到目标K空间信号矩阵，根据目标K空间信号矩阵，确定目标复数磁共振图像，能够在不牺牲图像分辨率的基础上，减小磁共振图像中的吉布斯伪影，提高磁共振图像的图像质量。振图像的图像质量。振图像的图像质量。

【技术实现步骤摘要】
磁共振图像的处理方法、装置及电子设备


[0001]本专利技术涉及医学图像处理
，尤其涉及一种磁共振图像的处理方法、装置及电子设备。

技术介绍

[0002]MRI(Magnetic Resonance Imaging，磁共振成像)是目前医学领域的一项重要的成像技术，所产生的磁共振图像可以辅助医生进行诊断。核磁共振是指质子在外磁场作用下自旋方向分布满足玻尔兹曼分布，在外加特定频率的射频磁场作用下吸收能量，撤去射频磁场后发生弛豫释放能量的现象。磁共振成像主要利用该原理，结合空间编码和傅里叶变换等技术，利用检测到的核磁共振信号还原出成像物体内部结构信息。
[0003]磁共振图像通常由有限矩阵大小的像素点组成。由傅里叶变换的性质可知，有限大小的图像在频率域(K空间)的信号分布是无限大的。然而在实际MR(Magnetic Resonance，磁共振)成像过程中，只能采集有限大小的K空间信号，这就造成了K空间信号的“截断”。K空间的信号截断等价于在无限K空间信号上乘一个矩形窗；根据傅里叶变换的性质，这等价于在图像域使用Sinc函数(辛格函数)对图像进行卷积。由于Sinc函数的旁瓣效应，会在图像中产生明显的信号震荡伪影，在低分辨率图像中或图像中幅值变化剧烈的区域(如组织边界)尤为明显。这种伪影一般被称为Gibbs(吉布斯)伪影或Ringing(振铃)伪影。
[0004]吉布斯伪影的存在降低了图像质量，影响了辅助诊断的准确性。

技术实现思路

[0005]为克服相关技术中存在的问题，本专利技术提供了一种磁共振图像的处理方法、装置及电子设备，提高磁共振图像的图像质量。
[0006]根据本专利技术实施例的第一方面，提供一种磁共振图像的处理方法，包括：
[0007]在预设方向上，将第一复数磁共振图像对应的第一K空间信号矩阵的外推空间的信号值设置为预设值，得到第二K空间信号矩阵；
[0008]对第二K空间信号矩阵对应的第二复数磁共振图像的幅值图像进行图像域滤波，得到新幅值图像；
[0009]利用所述第二复数磁共振图像的相位对所述新幅值图像进行相位调制，得到第三复数磁共振图像；
[0010]将所述第一K空间信号矩阵的信号值填充到所述第三复数磁共振图像对应的K空间信号矩阵中的相应位置，得到目标K空间信号矩阵；
[0011]根据所述目标K空间信号矩阵，确定目标复数磁共振图像。
[0012]根据本专利技术实施例的第二方面，提供一种磁共振图像的处理装置，包括：
[0013]设置模块，用于在预设方向上，将第一复数磁共振图像对应的第一K空间信号矩阵的外推空间的信号值设置为预设值，得到第二K空间信号矩阵；
[0014]滤波模块，用于对第二K空间信号矩阵对应的第二复数磁共振图像的幅值图像进行图像域，得到新幅值图像；
[0015]相位调制模块，用于利用所述第二复数磁共振图像的相位对所述新幅值图像进行相位调制，得到第三复数磁共振图像；
[0016]填充模块，用于将所述第一K空间信号矩阵的信号值填充到所述第三复数磁共振图像对应的K空间信号矩阵中的相应位置，得到目标K空间信号矩阵；
[0017]确定模块，用于根据所述目标K空间信号矩阵，确定目标复数磁共振图像。
[0018]根据本专利技术实施例的第三方面，提供一种电子设备，包括：内部总线，以及通过内部总线连接的存储器、处理器和外部接口；其中：
[0019]所述存储器，用于存储磁共振图像的处理逻辑对应的机器可读指令；
[0020]所述处理器，用于读取所述存储器上的所述机器可读指令，并执行如下操作：
[0021]在预设方向上，将第一复数磁共振图像对应的第一K空间信号矩阵的外推空间的信号值设置为预设值，得到第二K空间信号矩阵；
[0022]对第二K空间信号矩阵对应的第二复数磁共振图像的幅值图像进行图像域滤波，得到新幅值图像；
[0023]利用所述第二复数磁共振图像的相位对所述新幅值图像进行相位调制，得到第三复数磁共振图像；
[0024]将所述第一K空间信号矩阵的信号值填充到所述第三复数磁共振图像对应的K空间信号矩阵中的相应位置，得到目标K空间信号矩阵；
[0025]根据所述目标K空间信号矩阵，确定目标复数磁共振图像。
[0026]本专利技术实施例提供的技术方案可以包括以下有益效果：
[0027]本专利技术实施例，通过在预设方向上，将第一复数磁共振图像对应的第一K空间信号矩阵的外推空间的信号值设置为预设值，得到第二K空间信号矩阵，对第二K空间信号矩阵对应的第二复数磁共振图像的幅值图像进行图像域滤波，得到新幅值图像，利用所述第二复数磁共振图像的相位对所述新幅值图像进行相位调制，得到第三复数磁共振图像，将所述第一K空间信号矩阵的信号值填充到所述第三复数磁共振图像对应的K空间信号矩阵中的相应位置，得到目标K空间信号矩阵，根据所述目标K空间信号矩阵，确定目标复数磁共振图像，能够在不牺牲图像的原有分辨率的基础上，快速、鲁棒地减小磁共振图像中的吉布斯伪影，提高磁共振图像的图像质量，进而提高基于磁共振图像的辅助诊断的准确性。
[0028]应当理解的是，以上的一般描述和后文的细节描述仅是示例性和解释性的，并不能限制本说明书。
附图说明
[0029]此处的附图被并入说明书中并构成本说明书的一部分，示出了符合本说明书的实施例，并与说明书一起用于解释本说明书的原理。
[0030]图1是本专利技术实施例提供的磁共振图像的处理方法的流程示例图。
[0031]图2是吉布斯伪影产生和减小的原理示意图。
[0032]图3是原始磁共振图像与经过本实施例方法处理后的磁共振图像的对比示例图。
[0033]图4是本专利技术实施例提供的磁共振图像的处理装置的功能方块图。
[0034]图5是本专利技术实施例提供的电子设备的一个硬件结构图。
具体实施方式
[0035]这里将详细地对示例性实施例进行说明，其示例表示在附图中。下面的描述涉及附图时，除非另有表示，不同附图中的相同数字表示相同或相似的要素。以下示例性实施例中所描述的实施方式并不代表与本专利技术相一致的所有实施方式。相反，它们仅是与如所附权利要求书中所详述的、本专利技术实施例的一些方面相一致的装置和方法的例子。
[0036]在本专利技术实施例使用的术语是仅仅出于描述特定本专利技术实施例的目的，而非旨在限制本专利技术实施例。在本专利技术实施例和所附权利要求书中所使用的单数形式的“一种”、“所述”和“该”也旨在包括多数形式，除非上下文清楚地表示其他含义。还应当理解，本文中使用的术语“和/或”是指并包含一个或多个相关联的列出项目的任何或所有可能组合。
[0037]应当理解，尽管在本专利技术实施例可能采用术语第一、第二、第三等来描述各种信息，但这些信息不应限于这些术语。这些术语仅用来将同一类型的信息彼此区分开。例如，在不脱离本专利技术实施例范围的情况下，第一信息也可以被称为第二信息，类似地，第二信息也可以被称为第一信息。取本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种磁共振图像的处理方法，其特征在于，包括：在预设方向上，将第一复数磁共振图像对应的第一K空间信号矩阵的外推空间的信号值设置为预设值，得到第二K空间信号矩阵；对第二K空间信号矩阵对应的第二复数磁共振图像的幅值图像进行图像域滤波，得到新幅值图像；利用所述第二复数磁共振图像的相位对所述新幅值图像进行相位调制，得到第三复数磁共振图像；将所述第一K空间信号矩阵的信号值填充到所述第三复数磁共振图像对应的K空间信号矩阵中的相应位置，得到目标K空间信号矩阵；根据所述目标K空间信号矩阵，确定目标复数磁共振图像。2.根据权利要求1所述的方法，其特征在于，在预设方向上，将第一复数磁共振图像对应的第一K空间信号矩阵的外推空间的信号值设置为预设值，得到第二K空间信号矩阵之前，还包括：确定所述第一K空间信号矩阵的外推空间。3.根据权利要求2所述的方法，其特征在于，确定所述第一K空间信号矩阵的外推空间，包括：根据所述第一K空间信号矩阵的尺寸，确定所述第一K空间信号矩阵的外推空间对应的目标尺寸；若所述第一K空间信号矩阵对应的采样方式为满采，确定从所述第一K空间信号矩阵的边界向外扩展所述目标尺寸的区域为所述第一K空间信号矩阵的外推空间。4.根据权利要求2所述的方法，其特征在于，确定所述第一K空间信号矩阵的外推空间，包括：根据所述第一K空间信号矩阵的尺寸，确定所述第一K空间信号矩阵的外推空间对应的目标尺寸；若所述第一K空间信号矩阵对应的采样方式为欠采，若欠采区域的尺寸大于或等于所述目标尺寸，确定所述欠采区域为所述第一K空间信号矩阵的外推空间；若欠采区域的尺寸小于所述目标尺寸，确定所述欠采区域和从所述第一K空间信号矩阵的边界向外扩展第一尺寸的区域为所述第一K空间信号矩阵的外推空间；所述第一尺寸等于所述目标尺寸与所述欠采区域的尺寸之差。5.根据权利要求1所述的方法，其特征在于，根据所述目标K空间信号矩阵，确定目标复数磁共振图像，包括：对所述目标K空间信号矩阵进行傅里叶逆变换，得到中间复数磁共振图像；若所述目标K空间信号矩阵与所述第一K空间信号矩阵尺寸相同，将所述中间复数磁共振图像确定为目标复数磁共振图像。6.根据权利要求1所述的方法，...

【专利技术属性】
技术研发人员：熊昱辉，郭红宇，
申请(专利权)人：上海东软医疗科技有限公司，
类型：发明
国别省市：