CN111369465A - Ct动态图像增强方法及装置 - Google Patents

Abstract

本发明实施例提供一种CT动态图像增强方法及装置。本发明实施例通过通过获取待处理图像，待处理图像为CT动态图像序列对应的K空间图像序列中的一幅，对于待处理图像中的每一像素点，根据预设的时域权重值序列，确定像素点对应的在预设时长内的时域加权平均像素值，从预设的掩膜中查找出像素点对应的目标频域权重值，掩膜中第一频段对应的权重值小于第二频段对应的权重值，第一频段高于第二频段，基于目标频域权重值、时域加权平均像素值、像素点的当前像素值，获取像素点的目标像素值，将待处理图像的频域增强图像进行傅里叶逆变换，得到待处理图像对应的时域增强图像，有效减弱了拖影现象，提高了增强图像的图像质量。

Description

CT动态图像增强方法及装置

技术领域

本发明涉及医学图像处理技术领域，尤其涉及一种CT动态图像增强方法及装置。

背景技术

在应用中，经常需要对医学图像进行增强处理。例如，对CT(ComputedTomography，电子计算机断层扫描)动态图像做增强处理得到的增强图像，能够辅助医生进行更准确的诊疗。其中，CT动态图像是指通过常规CT电影扫描并重建得到的一系列时域的CT图像。

相关技术中，采用时域滤波方式对CT动态图像进行增强。然而，时域滤波会引入拖影现象，导致增强图像的质量较差。

发明内容

为克服相关技术中存在的问题，本发明提供了一种CT动态图像增强方法、装置、CT设备及CT***，提高心脏冠脉血管重建图像的图像质量。

根据本发明实施例的第一方面，提供一种CT动态图像增强方法，包括：

获取待处理图像，所述待处理图像为CT动态图像序列对应的K空间图像序列中的一幅图像；

对于所述待处理图像中的每一像素点，根据预设的时域权重值序列，确定所述像素点对应的在预设时长内的时域加权平均像素值；

从预设的掩膜中查找出所述像素点对应的目标频域权重值，所述掩膜中第一频段对应的权重值小于第二频段对应的权重值，所述第一频段高于所述第二频段；

基于所述目标频域权重值、所述时域加权平均像素值、所述像素点的当前像素值，获取所述像素点的目标像素值；

将所述待处理图像的频域增强图像进行傅里叶逆变换，得到所述待处理图像对应的时域增强图像，所述频域增强图像中的每一像素点的像素值等于所述待处理图像中相应像素点对应的目标像素值。

根据本发明实施例的第二方面，提供一种CT动态图像增强装置，包括：

图像获取模块，用于获取待处理图像，所述待处理图像为CT动态图像序列对应的K空间图像序列中的一幅图像；

确定模块，用于对于所述待处理图像中的每一像素点，根据预设的时域权重值序列，确定所述像素点对应的在预设时长内的时域加权平均像素值；

查找模块，用于从预设的掩膜中查找出所述像素点对应的目标频域权重值，所述掩膜中第一频段对应的权重值小于第二频段对应的权重值，所述第一频段高于所述第二频段；

像素值获取模块，用于基于所述目标频域权重值、所述时域加权平均像素值、所述像素点的当前像素值，获取所述像素点的目标像素值；

变换模块，用于将所述待处理图像的频域增强图像进行傅里叶逆变换，得到所述待处理图像对应的时域增强图像，所述频域增强图像中的每一像素点的像素值等于所述待处理图像中相应像素点对应的目标像素值。

本发明实施例提供的技术方案可以包括以下有益效果：

本发明实施例，通过获取待处理图像，所述待处理图像为CT动态图像序列对应的K空间图像序列中的一幅图像，对于所述待处理图像中的每一像素点，根据预设的时域权重值序列，确定所述像素点对应的在预设时长内的时域加权平均像素值，从预设的掩膜中查找出所述像素点对应的目标频域权重值，所述掩膜中第一频段对应的权重值小于第二频段对应的权重值，所述第一频段高于所述第二频段，基于所述目标频域权重值、所述时域加权平均像素值、所述像素点的当前像素值，获取所述像素点的目标像素值，将所述待处理图像的频域增强图像进行傅里叶逆变换，得到所述待处理图像对应的时域增强图像，既能保持较好的图像分辨率，又可以较好地抑制噪声，还有效减弱了拖影现象，因此大大提高了增强图像的图像质量。

应当理解的是，以上的一般描述和后文的细节描述仅是示例性和解释性的，并不能限制本说明书。

附图说明

此处的附图被并入说明书中并构成本说明书的一部分，示出了符合本说明书的实施例，并与说明书一起用于解释本说明书的原理。

图1是本发明实施例提供的CT动态图像增强方法的流程示例图。

图2是时域权重W_t的正态分布曲线图。

图3是二分类掩膜示意图。

图4是正态分布变化的mask示意图。

图5是本发明实施例提供的CT动态图像增强装置的功能方块图。

图6是本发明实施例提供的CT设备的一个硬件结构图。

具体实施方式

这里将详细地对示例性实施例进行说明，其示例表示在附图中。下面的描述涉及附图时，除非另有表示，不同附图中的相同数字表示相同或相似的要素。以下示例性实施例中所描述的实施方式并不代表与本发明相一致的所有实施方式。相反，它们仅是与如所附权利要求书中所详述的、本发明实施例的一些方面相一致的装置和方法的例子。

在本发明实施例使用的术语是仅仅出于描述特定本发明实施例的目的，而非旨在限制本发明实施例。在本发明实施例和所附权利要求书中所使用的单数形式的“一种”、“所述”和“该”也旨在包括多数形式，除非上下文清楚地表示其他含义。还应当理解，本文中使用的术语“和/或”是指并包含一个或多个相关联的列出项目的任何或所有可能组合。

应当理解，尽管在本发明实施例可能采用术语第一、第二、第三等来描述各种信息，但这些信息不应限于这些术语。这些术语仅用来将同一类型的信息彼此区分开。例如，在不脱离本发明实施例范围的情况下，第一信息也可以被称为第二信息，类似地，第二信息也可以被称为第一信息。取决于语境，如在此所使用的词语“如果”可以被解释成为“在……时”或“当……时”或“响应于确定”。

为了达到医生阅片的需求，增强图像需要满足一系列的要求，例如图像分辨率、信噪比等等。时域滤波既可以有效地抑制噪声，又能够较好地保持图像分辨率，但单一的时域滤波会引入拖影现象，这导致了增强图像质量的下降。

发明人在研究过程中发现，CT动态图像的不同帧之间在频域(即k空间)的高频区域的分布相近，而在低频区域的分布不同。基于该发现，发明人提出了一种新的基于时域和频域结合的CT动态图像增强方法。

下面通过实施例对CT动态图像增强方法进行详细说明。

图1是本发明实施例提供的CT动态图像增强方法的流程示例图。如图1所示，本实施例中，CT动态图像增强方法可以包括：

S101，获取待处理图像，所述待处理图像为CT动态图像序列对应的K空间图像序列中的一幅图像。

S102，对于所述待处理图像中的每一像素点，根据预设的时域权重值序列，确定所述像素点对应的在预设时长内的时域加权平均像素值。

S103，从预设的掩膜中查找出所述像素点对应的目标频域权重值，所述掩膜中第一频段对应的权重值小于第二频段对应的权重值，所述第一频段高于所述第二频段。

S104，基于所述目标频域权重值、所述时域加权平均像素值、所述像素点的当前像素值，获取所述像素点的目标像素值。

S105，将所述待处理图像的频域增强图像进行傅里叶逆变换，得到所述待处理图像对应的时域增强图像，所述频域增强图像中的每一像素点的像素值等于所述待处理图像中相应像素点对应的目标像素值。

本实施例中，待处理图像属于频域的CT动态图像。在应用中，对受检者的待扫描部位进行常规CT电影扫描后得到CT电影扫描数据，利用CT电影扫描数据重建得到时域的CT图像序列，对时域的CT图像序列进行傅里叶变换就可以得到CT图像序列对应的频域图像序列，也即K空间图像序列。

因此，在一个示例性的实现过程中，获取待处理图像，可以包括：

对CT动态图像序列进行傅里叶变换，得到所述CT动态图像序列对应的K空间图像序列；

按照预设提取策略，从所述K空间图像序列中提取出一幅图像，作为待处理图像。

举例说明。通过对受检者的冠脉进行常规CT电影扫描并重建，可以获得n层冠脉CT图像，记为I₀、I₁、I₂……I_n-1。I₀、I₁、I₂……I_n-1为时域图像序列。

分别对I₀、I₁、I₂……I_n-1做傅里叶变换，转化为k空间图像K₀、K₁、K₂……K_n-1。

取K₀、K₁、K₂……K_n-1中的任意一幅图像K_m(m为自然数，m的值可以为0、1、2……n-1中的任何一个)，作为待处理图像。该图像K_p中的像素点的像素值等于K_i,j,t。其中，i、j分别表示像素点在图像中的横坐标、纵坐标，t表示图像的采集时间。

本实施例中，预设提取策略可以由开发人员预先设定。例如，在一个示例中，预设提取策略可以是按照图像的采集时间顺序提取。在另一个示例中，预设提取策略也可以是根据用户的选择操作从K空间图像序列中提取。本实施例不对预设提取策略进行限制。

本实施例中，预设的时域权重值序列可以是平均分布的，也可以是正态分布(当前帧为分布的最大值，如图2所示)的。

假设时域权重值序列中的权重值为W_t，K₀、K₁、K₂……K_n-1对应的时域权重值分别为W_t0、W_t1、W_t2……W_tn-1。在平均分布的时域权重值序列中，W_t0＝W_t1＝W_t2＝……＝W_tn-1＝1/n。在正态分布的时域权重值序列中，当前帧对应的时域权重值最大。图2是时域权重W_t的正态分布曲线图。如图2所示，整个曲线为正态形式，当前图像帧的时域权重W_t＝1，从当前图像帧向两边，权重W_t的值呈现减小趋势。

在一个示例性的实现过程中，对于所述待处理图像中的每一像素点，根据预设的时域权重值序列，确定所述像素点对应的在预设时长内的时域加权平均像素值，可以包括：

从所述K空间图像序列中选择处于以目标时间点为中心的预设时长内的目标图像；

分别获取各个目标图像中与所述像素点对应的像素点的像素值，作为参考像素值；

从预设的时域权重值序列中，获取与各个参考像素值对应的匹配时域权重值，并确定各个参考像素值与对应的匹配时域权重值的乘积，作为像素值分量；

获取所述各个参考像素值对应的像素值分量之和与权重和的商，作为所述像素点对应的在预设时长内的时域加权平均像素值，所述权重和为所述各个参考像素值对应的匹配时域权重值之和。

其中，预设时长可以根据应用需要设定。本文中，称s为时间序列方向前后的步长，则预设时长＝2s+1。

假设K₀、K₁、K₂……K_n-1中当前图像帧K_m(即待处理图像)的时间为t，则目标时间点为t，目标图像为K₀、K₁、K₂……K_n-1中处于时间段[t-s，t+s]内的图像。

在应用中，可以根据实际需求来设置s值。例如，如果图片中的噪声较大，则可以设置s为较大的值，此时可以用较多的图像帧来做时间方向的平均，得到较好的去噪效果。如果图片中的噪声较小，则可以设置s为较小的值，此时可以用较少的图像帧来做时间方向的平均，减少计算量，提高处理速度。

对于K_m中的任一像素点，该像素点对应的在预设时长内的时域加权平均像素值等于

例如，假设K_m为K₃，s＝2，则K₃中的任一像素点K_i,j对应的在预设时长内的时域加权平均像素值等于W₁*K_i,j,1、W₂*K_i,j,2、W₃*K_i,j,3、W₄*K_i,j,4、W₅*K_i,j,5的和除以W₁、W₂、W₃、W₄、W₅的和所得的商。其中，K_i,j,1、K_i,j,2、K_i,j,4、K_i,j,5分别为图像K₁、K₂、K₄、K₅中与图像K₃的像素点K_i,j,对应的像素点的像素值，W₁、W₂、W₃、W₄、W₅分别为K₁、K₂、K₃、K₄、K₅对应的时域权重值。

如果待处理图像的前后步长内数据不足，则可以利用镜像方法在时域方向补充数据。在补充数据时，以K₀为界左右对称补充数据。例如待处理图像为K₁，t＝t1，s＝2，则W₁后的步长内有2幅图像K₂、K₃，数据充足，而K₁前的步长内只有一幅图像K₀，数据不足，这时可以在前补充图像K_-1，令图像K_-1＝K₁。

在CT动态图像中，图像帧之间具有很强的相关性，本实施例通过步骤S102，在频域的时域方向做加权平均的计算，能够结合待处理图像帧附近的多帧图像数据进行滤波，有效减弱了拖影现象，有助于提高图像质量。

本实施例中，第一频段为频域的高频区域，第二频段为频域的低频区域。

在步骤S103中，掩膜(mask)中第一频段对应的权重值小于第二频段对应的权重值，这表示高频区域的频域权重值小于低频区域的频域权重值。

本实施例中，掩膜是预先设置好的。掩膜可以用一个矩阵来表示，为了描述方便，本文中称这个矩阵为掩膜矩阵。掩膜矩阵的大小与待处理图像的大小相同。例如，待处理图像为256×256的矩阵，则掩膜矩阵也是256×256的矩阵。掩膜矩阵中的每个元素为待处理图像中位置相同的像素点对应的频域权重值。

在一个示例性的实现过程中，所述掩膜为如下掩膜中的任意一种：

掩膜中所述第一频段的所有像素点的权重值等于第一权重值，所述第二频段的所有像素点的权重值等于第二权重值，所述第一权重值小于所述第二权重值；

掩膜中的权重值从低频向高频呈正态分布。

例如，在一个示例中，可以设置k空间的高频区域与低频区域的mask为二分类的mask，如图3所示，图3是二分类掩膜示意图。在图3的二分类掩膜中，中间低频区域的频域权重值W_i,j为1，周围高频区域的频域权重值W_i,j为0。

例如，在另一个示例中，可以设置k空间的mask为正态分布变化的mask，如图4所示，图4是正态分布变化的mask示意图。在图4的正态分布变化的mask，中心低频区域的频域权重值为1，从低频向高频扩展，频域权重值随正态分布变化。

当然，根据应用需要，也可以设置其他类型的掩膜，本实施例对掩膜的具体形式不作限制。

在应用中，可以根据实际需要调整频域权重值W_i,j，使得序列中变化不剧烈的高频区可以用较多的帧来加权平均，而变化较大的低频区域则可以用较少的帧来加权平均。

本实施例中，采用时域(时域加权平均像素值)和频域(目标频域权重值)相结合的方式计算增强图像的像素值，既能利用时域滤波的特点保持较好的图像分辨率、较好地抑制噪声，又能利用频域滤波有效减弱时域滤波所引用的拖影现象。

在一个示例性的实现过程中，基于所述目标频域权重值、所述时域加权平均像素值、所述像素点的当前像素值，获取所述像素点的目标像素值，可以包括：

获取所述时域加权平均像素值与目标差值的乘积，作为第一值，所述目标差值为1与所述目标频域权重值的差值；

获取所述像素点的当前像素值与所述目标频域权重值的乘积，作为第二值；

将所述第一值与所述第二值的和确定为所述像素点的目标像素值。

例如，假设增强后图像中的像素值为Ks_i,j,t，则Ks_i,j,t的值可以用如下的公式(1)表示。

公式(1)中，K_i,j,t为待处理图像中像素点的像素值，W_i,j为像素点对应的频域权重值。

由于高频区域的频域权重值小于低频区域的频域权重值，因此当像素点属于高频区域时，1-W_i,j大于W_i,j，此时时域加权平均像素值在增强后像素值中的占比高于原像素值K_i,j,t在增强后像素值中的占比，从而使得高频区域的滤波程度更大，有效去除高频部分的噪声。当像素点属于低频区域时，1-W_i,j小于W_i,j，此时时域加权平均像素值在增强后像素值中的占比低于原像素值K_i,j,t在增强后像素值中的占比，从而使得低频区域的滤波程度较小。这样利用了高频区域噪声较多而低频区域噪声较少的特点，对图像的不同频率部分进行了适当的去噪，能够取得较好的增强效果。

本发明实施例提供的CT动态图像增强方法，通过获取待处理图像，所述待处理图像为CT动态图像序列对应的K空间图像序列中的一幅图像，对于所述待处理图像中的每一像素点，根据预设的时域权重值序列，确定所述像素点对应的在预设时长内的时域加权平均像素值，从预设的掩膜中查找出所述像素点对应的目标频域权重值，所述掩膜中第一频段对应的权重值小于第二频段对应的权重值，所述第一频段高于所述第二频段，基于所述目标频域权重值、所述时域加权平均像素值、所述像素点的当前像素值，获取所述像素点的目标像素值，将所述待处理图像的频域增强图像进行傅里叶逆变换，得到所述待处理图像对应的时域增强图像，既能保持较好的图像分辨率，又可以较好地抑制噪声，还有效减弱了拖影现象，因此大大提高了增强图像的图像质量。

基于上述的方法实施例，本发明实施例还提供了相应的装置、设备及存储介质实施例。

图5是本发明实施例提供的CT动态图像增强装置的功能方块图。如图5所示，本实施例中，CT动态图像增强装置可以包括：

图像获取模块510，用于获取待处理图像，所述待处理图像为CT动态图像序列对应的K空间图像序列中的一幅图像；

确定模块520，用于对于所述待处理图像中的每一像素点，根据预设的时域权重值序列，确定所述像素点对应的在预设时长内的时域加权平均像素值；

查找模块530，用于从预设的掩膜中查找出所述像素点对应的目标频域权重值，所述掩膜中第一频段对应的权重值小于第二频段对应的权重值，所述第一频段高于所述第二频段；

像素值获取模块540，用于基于所述目标频域权重值、所述时域加权平均像素值、所述像素点的当前像素值，获取所述像素点的目标像素值；

变换模块550，用于将所述待处理图像的频域增强图像进行傅里叶逆变换，得到所述待处理图像对应的时域增强图像，所述频域增强图像中的每一像素点的像素值等于所述待处理图像中相应像素点对应的目标像素值。

在一个示例性的实现过程中，所述图像获取模块510可以具体用于：

对CT动态图像序列进行傅里叶变换，得到所述CT动态图像序列对应的K空间图像序列；

按照预设提取策略，从所述K空间图像序列中提取出一幅图像，作为待处理图像。

在一个示例性的实现过程中，所述确定模块520可以具体用于：

从所述K空间图像序列中选择处于以目标时间点为中心的预设时长内的目标图像；

分别获取各个目标图像中与所述像素点对应的像素点的像素值，作为参考像素值；

从预设的时域权重值序列中，获取与各个参考像素值对应的匹配时域权重值，并确定各个参考像素值与对应的匹配时域权重值的乘积，作为像素值分量；

获取所述各个参考像素值对应的像素值分量之和与权重和的商，作为所述像素点对应的在预设时长内的时域加权平均像素值，所述权重和为所述各个参考像素值对应的匹配时域权重值之和。

在一个示例性的实现过程中，所述像素值获取模块540可以具体用于：

获取所述时域加权平均像素值与目标差值的乘积，作为第一值，所述目标差值为1与所述目标频域权重值的差值；

获取所述像素点的当前像素值与所述目标频域权重值的乘积，作为第二值；

将所述第一值与所述第二值的和确定为所述像素点的目标像素值。

在一个示例性的实现过程中，所述掩膜为如下掩膜中的任意一种：

掩膜中所述第一频段的所有像素点的权重值等于第一权重值，所述第二频段的所有像素点的权重值等于第二权重值，所述第一权重值小于所述第二权重值；

掩膜中的权重值从低频向高频呈正态分布。

本发明实施例还提供了一种CT设备。图6是本发明实施例提供的CT设备的一个硬件结构图。如图6所示，CT设备包括：内部总线601，以及通过内部总线连接的存储器602，处理器603和外部接口604，其中，所述外部接口，用于连接CT***的探测器，所述探测器包括多个探测器室及相应的处理电路；

所述存储器602，用于存储CT动态图像增强逻辑对应的机器可读指令；

所述处理器603，用于读取存储器602上的机器可读指令，并执行所述指令以实现如下操作：

获取待处理图像，所述待处理图像为CT动态图像序列对应的K空间图像序列中的一幅图像；

对于所述待处理图像中的每一像素点，根据预设的时域权重值序列，确定所述像素点对应的在预设时长内的时域加权平均像素值；

从预设的掩膜中查找出所述像素点对应的目标频域权重值，所述掩膜中第一频段对应的权重值小于第二频段对应的权重值，所述第一频段高于所述第二频段；

基于所述目标频域权重值、所述时域加权平均像素值、所述像素点的当前像素值，获取所述像素点的目标像素值；

将所述待处理图像的频域增强图像进行傅里叶逆变换，得到所述待处理图像对应的时域增强图像，所述频域增强图像中的每一像素点的像素值等于所述待处理图像中相应像素点对应的目标像素值。

在一个示例性的实现过程中，获取待处理图像，包括：

对CT动态图像序列进行傅里叶变换，得到所述CT动态图像序列对应的K空间图像序列；

按照预设提取策略，从所述K空间图像序列中提取出一幅图像，作为待处理图像。

在一个示例性的实现过程中，对于所述待处理图像中的每一像素点，根据预设的时域权重值序列，确定所述像素点对应的在预设时长内的时域加权平均像素值，包括：

从所述K空间图像序列中选择处于以目标时间点为中心的预设时长内的目标图像；

分别获取各个目标图像中与所述像素点对应的像素点的像素值，作为参考像素值；

从预设的时域权重值序列中，获取与各个参考像素值对应的匹配时域权重值，并确定各个参考像素值与对应的匹配时域权重值的乘积，作为像素值分量；

获取所述各个参考像素值对应的像素值分量之和与权重和的商，作为所述像素点对应的在预设时长内的时域加权平均像素值，所述权重和为所述各个参考像素值对应的匹配时域权重值之和。

在一个示例性的实现过程中，基于所述目标频域权重值、所述时域加权平均像素值、所述像素点的当前像素值，获取所述像素点的目标像素值，包括：

获取所述时域加权平均像素值与目标差值的乘积，作为第一值，所述目标差值为1与所述目标频域权重值的差值；

获取所述像素点的当前像素值与所述目标频域权重值的乘积，作为第二值；

将所述第一值与所述第二值的和确定为所述像素点的目标像素值。

在一个示例性的实现过程中，所述掩膜为如下掩膜中的任意一种：

掩膜中所述第一频段的所有像素点的权重值等于第一权重值，所述第二频段的所有像素点的权重值等于第二权重值，所述第一权重值小于所述第二权重值；

掩膜中的权重值从低频向高频呈正态分布。

本发明实施例还提供一种CT***，包括探测器、扫描床和CT设备，所述探测器包括多个探测器室及相应的处理电路；其中：

所述探测器室，用于在所述CT***扫描过程中，探测穿过扫描对象的X射线并转换为电信号；

所述处理电路，用于将所述电信号转换成脉冲信号，采集脉冲信号的能量信息；

所述CT设备，用于获取待处理图像，所述待处理图像为CT动态图像序列对应的K空间图像序列中的一幅图像；对于所述待处理图像中的每一像素点，根据预设的时域权重值序列，确定所述像素点对应的在预设时长内的时域加权平均像素值；从预设的掩膜中查找出所述像素点对应的目标频域权重值，所述掩膜中第一频段对应的权重值小于第二频段对应的权重值，所述第一频段高于所述第二频段；基于所述目标频域权重值、所述时域加权平均像素值、所述像素点的当前像素值，获取所述像素点的目标像素值；将所述待处理图像的频域增强图像进行傅里叶逆变换，得到所述待处理图像对应的时域增强图像，所述频域增强图像中的每一像素点的像素值等于所述待处理图像中相应像素点对应的目标像素值。

在一个示例性的实现过程中，获取待处理图像，包括：

对CT动态图像序列进行傅里叶变换，得到所述CT动态图像序列对应的K空间图像序列；

按照预设提取策略，从所述K空间图像序列中提取出一幅图像，作为待处理图像。

在一个示例性的实现过程中，对于所述待处理图像中的每一像素点，根据预设的时域权重值序列，确定所述像素点对应的在预设时长内的时域加权平均像素值，包括：

从所述K空间图像序列中选择处于以目标时间点为中心的预设时长内的目标图像；

分别获取各个目标图像中与所述像素点对应的像素点的像素值，作为参考像素值；

从预设的时域权重值序列中，获取与各个参考像素值对应的匹配时域权重值，并确定各个参考像素值与对应的匹配时域权重值的乘积，作为像素值分量；

获取所述各个参考像素值对应的像素值分量之和与权重和的商，作为所述像素点对应的在预设时长内的时域加权平均像素值，所述权重和为所述各个参考像素值对应的匹配时域权重值之和。

在一个示例性的实现过程中，基于所述目标频域权重值、所述时域加权平均像素值、所述像素点的当前像素值，获取所述像素点的目标像素值，包括：

获取所述时域加权平均像素值与目标差值的乘积，作为第一值，所述目标差值为1与所述目标频域权重值的差值；

获取所述像素点的当前像素值与所述目标频域权重值的乘积，作为第二值；

将所述第一值与所述第二值的和确定为所述像素点的目标像素值。

在一个示例性的实现过程中，所述掩膜为如下掩膜中的任意一种：

掩膜中所述第一频段的所有像素点的权重值等于第一权重值，所述第二频段的所有像素点的权重值等于第二权重值，所述第一权重值小于所述第二权重值；

掩膜中的权重值从低频向高频呈正态分布。

本发明实施例还提供一种计算机可读存储介质，其上存储有计算机程序，其中，所述程序被处理器执行时实现如下操作：

获取待处理图像，所述待处理图像为CT动态图像序列对应的K空间图像序列中的一幅图像；

对于所述待处理图像中的每一像素点，根据预设的时域权重值序列，确定所述像素点对应的在预设时长内的时域加权平均像素值；

从预设的掩膜中查找出所述像素点对应的目标频域权重值，所述掩膜中第一频段对应的权重值小于第二频段对应的权重值，所述第一频段高于所述第二频段；

基于所述目标频域权重值、所述时域加权平均像素值、所述像素点的当前像素值，获取所述像素点的目标像素值；

将所述待处理图像的频域增强图像进行傅里叶逆变换，得到所述待处理图像对应的时域增强图像，所述频域增强图像中的每一像素点的像素值等于所述待处理图像中相应像素点对应的目标像素值。

在一个示例性的实现过程中，获取待处理图像，包括：

对CT动态图像序列进行傅里叶变换，得到所述CT动态图像序列对应的K空间图像序列；

按照预设提取策略，从所述K空间图像序列中提取出一幅图像，作为待处理图像。

在一个示例性的实现过程中，对于所述待处理图像中的每一像素点，根据预设的时域权重值序列，确定所述像素点对应的在预设时长内的时域加权平均像素值，包括：

从所述K空间图像序列中选择处于以目标时间点为中心的预设时长内的目标图像；

分别获取各个目标图像中与所述像素点对应的像素点的像素值，作为参考像素值；

从预设的时域权重值序列中，获取与各个参考像素值对应的匹配时域权重值，并确定各个参考像素值与对应的匹配时域权重值的乘积，作为像素值分量；

获取所述各个参考像素值对应的像素值分量之和与权重和的商，作为所述像素点对应的在预设时长内的时域加权平均像素值，所述权重和为所述各个参考像素值对应的匹配时域权重值之和。

在一个示例性的实现过程中，基于所述目标频域权重值、所述时域加权平均像素值、所述像素点的当前像素值，获取所述像素点的目标像素值，包括：

获取所述时域加权平均像素值与目标差值的乘积，作为第一值，所述目标差值为1与所述目标频域权重值的差值；

获取所述像素点的当前像素值与所述目标频域权重值的乘积，作为第二值；

将所述第一值与所述第二值的和确定为所述像素点的目标像素值。

在一个示例性的实现过程中，所述掩膜为如下掩膜中的任意一种：

掩膜中所述第一频段的所有像素点的权重值等于第一权重值，所述第二频段的所有像素点的权重值等于第二权重值，所述第一权重值小于所述第二权重值；

掩膜中的权重值从低频向高频呈正态分布。

对于装置和设备实施例而言，由于其基本对应于方法实施例，所以相关之处参见方法实施例的部分说明即可。以上所描述的装置实施例仅仅是示意性的，其中所述作为分离部件说明的模块可以是或者也可以不是物理上分开的，作为模块显示的部件可以是或者也可以不是物理模块，即可以位于一个地方，或者也可以分布到多个网络模块上。可以根据实际的需要选择其中的部分或者全部模块来实现本说明书方案的目的。本领域普通技术人员在不付出创造性劳动的情况下，即可以理解并实施。

上述对本说明书特定实施例进行了描述。其它实施例在所附权利要求书的范围内。在一些情况下，在权利要求书中记载的动作或步骤可以按照不同于实施例中的顺序来执行并且仍然可以实现期望的结果。另外，在附图中描绘的过程不一定要求示出的特定顺序或者连续顺序才能实现期望的结果。在某些实施方式中，多任务处理和并行处理也是可以的或者可能是有利的。

本领域技术人员在考虑说明书及实践这里申请的发明后，将容易想到本说明书的其它实施方案。本说明书旨在涵盖本说明书的任何变型、用途或者适应性变化，这些变型、用途或者适应性变化遵循本说明书的一般性原理并包括本说明书未申请的本技术领域中的公知常识或惯用技术手段。说明书和实施例仅被视为示例性的，本说明书的真正范围和精神由下面的权利要求指出。

应当理解的是，本说明书并不局限于上面已经描述并在附图中示出的精确结构，并且可以在不脱离其范围进行各种修改和改变。本说明书的范围仅由所附的权利要求来限制。

以上所述仅为本说明书的较佳实施例而已，并不用以限制本说明书，凡在本说明书的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本说明书保护的范围之内。

Claims (10)

1.一种CT动态图像增强方法，其特征在于，包括：

获取待处理图像，所述待处理图像为CT动态图像序列对应的K空间图像序列中的一幅图像；

对于所述待处理图像中的每一像素点，根据预设的时域权重值序列，确定所述像素点对应的在预设时长内的时域加权平均像素值；

从预设的掩膜中查找出所述像素点对应的目标频域权重值，所述掩膜中第一频段对应的权重值小于第二频段对应的权重值，所述第一频段高于所述第二频段；

基于所述目标频域权重值、所述时域加权平均像素值、所述像素点的当前像素值，获取所述像素点的目标像素值；

将所述待处理图像的频域增强图像进行傅里叶逆变换，得到所述待处理图像对应的时域增强图像，所述频域增强图像中的每一像素点的像素值等于所述待处理图像中相应像素点对应的目标像素值。

2.根据权利要求1所述的方法，其特征在于，获取待处理图像，包括：

对CT动态图像序列进行傅里叶变换，得到所述CT动态图像序列对应的K空间图像序列；

按照预设提取策略，从所述K空间图像序列中提取出一幅图像，作为待处理图像。

3.根据权利要求1所述的方法，其特征在于，对于所述待处理图像中的每一像素点，根据预设的时域权重值序列，确定所述像素点对应的在预设时长内的时域加权平均像素值，包括：

从所述K空间图像序列中选择处于以目标时间点为中心的预设时长内的目标图像；

分别获取各个目标图像中与所述像素点对应的像素点的像素值，作为参考像素值；

从预设的时域权重值序列中，获取与各个参考像素值对应的匹配时域权重值，并确定各个参考像素值与对应的匹配时域权重值的乘积，作为像素值分量；

获取所述各个参考像素值对应的像素值分量之和与权重和的商，作为所述像素点对应的在预设时长内的时域加权平均像素值，所述权重和为所述各个参考像素值对应的匹配时域权重值之和。

4.根据权利要求1所述的方法，其特征在于，基于所述目标频域权重值、所述时域加权平均像素值、所述像素点的当前像素值，获取所述像素点的目标像素值，包括：

获取所述时域加权平均像素值与目标差值的乘积，作为第一值，所述目标差值为1与所述目标频域权重值的差值；

获取所述像素点的当前像素值与所述目标频域权重值的乘积，作为第二值；

将所述第一值与所述第二值的和确定为所述像素点的目标像素值。

5.根据权利要求1所述的方法，其特征在于，所述掩膜为如下掩膜中的任意一种：

掩膜中所述第一频段的所有像素点的权重值等于第一权重值，所述第二频段的所有像素点的权重值等于第二权重值，所述第一权重值小于所述第二权重值；

掩膜中的权重值从低频向高频呈正态分布。

6.一种CT动态图像增强装置，其特征在于，包括：

图像获取模块，用于获取待处理图像，所述待处理图像为CT动态图像序列对应的K空间图像序列中的一幅图像；

确定模块，用于对于所述待处理图像中的每一像素点，根据预设的时域权重值序列，确定所述像素点对应的在预设时长内的时域加权平均像素值；

查找模块，用于从预设的掩膜中查找出所述像素点对应的目标频域权重值，所述掩膜中第一频段对应的权重值小于第二频段对应的权重值，所述第一频段高于所述第二频段；

像素值获取模块，用于基于所述目标频域权重值、所述时域加权平均像素值、所述像素点的当前像素值，获取所述像素点的目标像素值；

变换模块，用于将所述待处理图像的频域增强图像进行傅里叶逆变换，得到所述待处理图像对应的时域增强图像，所述频域增强图像中的每一像素点的像素值等于所述待处理图像中相应像素点对应的目标像素值。

7.根据权利要求6所述的装置，其特征在于，所述图像获取模块具体用于：

对CT动态图像序列进行傅里叶变换，得到所述CT动态图像序列对应的K空间图像序列；

按照预设提取策略，从所述K空间图像序列中提取出一幅图像，作为待处理图像。

8.根据权利要求6所述的装置，其特征在于，所述确定模块具体用于：

从所述K空间图像序列中选择处于以目标时间点为中心的预设时长内的目标图像；

分别获取各个目标图像中与所述像素点对应的像素点的像素值，作为参考像素值；

从预设的时域权重值序列中，获取与各个参考像素值对应的匹配时域权重值，并确定各个参考像素值与对应的匹配时域权重值的乘积，作为像素值分量；

获取所述各个参考像素值对应的像素值分量之和与权重和的商，作为所述像素点对应的在预设时长内的时域加权平均像素值，所述权重和为所述各个参考像素值对应的匹配时域权重值之和。

9.根据权利要求6所述的装置，其特征在于，所述像素值获取模块具体用于：

获取所述时域加权平均像素值与目标差值的乘积，作为第一值，所述目标差值为1与所述目标频域权重值的差值；

获取所述像素点的当前像素值与所述目标频域权重值的乘积，作为第二值；

将所述第一值与所述第二值的和确定为所述像素点的目标像素值。

10.根据权利要求6所述的装置，其特征在于，所述掩膜为如下掩膜中的任意一种：

掩膜中所述第一频段的所有像素点的权重值等于第一权重值，所述第二频段的所有像素点的权重值等于第二权重值，所述第一权重值小于所述第二权重值；

掩膜中的权重值从低频向高频呈正态分布。

Images