一种CT图像重建方法及其系统与流程

本发明涉及医学成像技术领域，尤其涉及的是一种ct图像重建方法及其系统。

背景技术：

ct医学成像系统自20世纪70年代发明了之后经过了长足的进步，扫描速度从开始要几分钟到现在的0.2秒。探测器排数也从开始的单排双排，到现在的64排，128排，甚至256排。这其中的变化不单单是系统硬件的升级换代，系统的图像重建技术也带来了革命性的变化。初期的ct系统由于只有一排探测器，所以x射线束是扇形束，所用到的重建技术也都是二维的扇形束重建技术。因为每次只能扫描一层，整个扫描需要很长的时间，后来多排ct引入就是为了加快扫描的速度，比如16排，32排的系统。这时候x射线也变成了三维的锥形束，和之前的几何结构有区别，所以必须要用锥形束重建技术来重建图像。

虽然锥形束精确重建技术的数学理论已经被提出来很久，但是由于算法的复杂性，还没有能够普遍的应用到系统中。目前主流ct产品中所应用的锥形束重建技术都是基于fdk算法的近似算法来解析图像。

但是由于重建误差和锥角(正比于探测器的排数)的平方成正比，所以在锥角(探测器排数)比较小的时候还能够得到比较好的效果(小于等于64排)。但是当探测器排数增大到128甚至256排的时候就会带来很大的误差。如图4和图5所示，图4是采用常规处理方法对锥形束伪影图像重建后图像，图5是目标图像。现有迭代算法虽然能够改善图像质量减少锥形束伪影对图像质量带来的影响，但是算法的计算量很大，很难满足临床的实际应用。

因此，现有技术还有待于改进和发展。

技术实现要素：

本发明要解决的技术问题在于，针对现有技术的上述缺陷，提供一种ct图像重建方法及其系统，旨在解决现有技术中的迭代算法在减少锥形束伪影对图像质量带来的影响的过程中计算量很大，很难满足临床的实际应用的问题。

本发明解决技术问题所采用的技术方案如下：

一种ct图像重建方法，其中，包括步骤：

通过图像重建技术对圆周扫描得到的初始锥形束投影数据重建得到初始图像；

对初始图像进行非线性处理得到矫正图像；

对矫正图像进行数值模拟产生矫正锥形束投影数据；

对矫正锥形束投影数据和初始锥形束投影数据进行计算得到差异投影数据；

通过图像重建技术对差异投影数据处理得到差异图像；

差异图像和矫正图像融合得到最终图像。

所述ct图像重建方法，其中，所述通过图像重建技术对圆周扫描得到的初始锥形束投影数据重建得到初始图像的步骤具体包括：

通过圆周扫描采集锥形束投影数据；

采用fdk算法对圆周扫描得到的锥形束投影数据重建得到初始图像。

所述ct图像重建方法，其中，所述对初始图像进行非线性处理得到矫正图像的步骤中，所述非线性处理为减弱伪影的非线性函数处理，对于减弱伪影的非线性函数处理过程可以是基于不同组织的图像处理函数，或者为训练好的神经网络处理，易于理解的是用于减弱初始图像中伪影的算法过程均可应用在本过程中，应均为本方案的相似方案；本实施例中的优选方案为：减弱伪影的非线性函数处理为训练好的神经网络处理，所述训练好的神经网络处理是通过优化目标图像和输出图像的差异来使得输出图像尽可能接近理想图像。

所述ct图像重建方法，其中，所述对初始图像进行非线性处理得到矫正图像的步骤中的具体步骤为：

构建图像数据库，图像数据库的样本包括图像中无伪影的目标图像，以及存在锥形束伪影的输出图像；

构建多尺度卷积神经网络模型；

构建模型的目标函数；

利用图像数据库和模型按照目标函数训练神经网络；

利用训练好的神经网络处理初始图像得到矫正图像。

所述ct图像重建方法，其中，在所述构建多尺度卷积神经网络的步骤中，所述多尺度卷积神经网络模型为改进的resnet或改进的unet。

所述ct图像重建方法，其中，所述目标函数为：

其中imgerr为图像的误差，img为网络的输出图像，imgtarg是数据库里面的目标图像；imgk为输出图像中坐标为k的像素点，imgtarg,k为数据库里面的目标图像中坐标为k的像素点。

所述ct图像重建方法，其中，在所述对矫正图像进行数值模拟产生矫正锥形束投影数据的过程中：

数值模拟过程为模拟ct系统的成像过程，从矫正图像得到矫正锥形束投影数据。

所述ct图像重建方法，其中，在所述通过图像重建技术对差异投影数据处理得到差异图像的步骤中，其图像重建技术与步骤通过图像重建技术对初始锥形束投影数据重建得到初始图像中的图像重建技术为同一图像重建技术。

所述ct图像重建方法，其中，在所述对矫正锥形束投影数据和初始锥形束投影数据进行计算得到差异投影数据的步骤中，通过对矫正锥形束投影数据和初始锥形束投影数据直接相减得到的误差投影数据。

一种ct图像重建系统，其中，包括：处理器，以及与所述处理器连接的存储器；

所述存储器存储有圆周扫描ct锥形束图像重建程序，所述圆周扫描ct锥形束图像重建程序被所述处理器执行时实现如上所述的ct图像重建方法。

本发明所提供的一种ct图像重建方法及其系统，所述方法包括步骤：通过图像重建技术对圆周扫描得到的初始锥形束投影数据重建得到初始图像；对初始图像进行非线性处理得到矫正图像；通过对矫正图像进行数值模拟产生矫正锥形束投影数据；对矫正锥形束投影数据和初始锥形束投影数据进行计算得到差异投影数据；通过图像重建技术对差异投影数据处理得到差异图像；差异图像和矫正图像融合得到最终图像。从而使带伪影的差异图像和矫正伪影的矫正图像结合实现锥形束伪影图像重建；通过在现有的重建技术的基础上只加入了投影模拟和非线性变换，在伪影图像重建过程中实现简单而较少的计算量，相比迭代算法要更加简单，同时还能取得好的图像重建效果，满足临床的实际应用，具有通用性。

附图说明

图1是本发明中ct图像重建方法的实施例的步骤流程图。

图2是本发明中ct图像重建方法较佳实施例的步骤流程图。

图3是本发明中ct图像重建方法的实施例流程示意图。

图4是采用常规处理方法对锥形束伪影图像重建后图像。

图5是目标图像。

图6是本发明中ct图像重建系统较佳实施例的功能原理框图。

具体实施方式

为使本发明的目的、技术方案及优点更加清楚、明确，以下参照附图并举实施例对本发明进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。

请同时参阅图1-图3，本发明提供了一种ct图像重建方法的一些实施例。

如图1所示，一种ct图像重建方法，包括以下步骤：

步骤s100、通过图像重建技术对圆周扫描得到的初始锥形束投影数据重建得到初始图像。

具体地，通过ct医学成像系统采用圆周扫描方式采集初始锥形束投影数据，这一过程通过现有ct成像系统完成，采集的初始锥形束投影数据经过数据重建后得到初始图像，这个初始图像由于现有探测器排数增大到128甚至256排的时候就会带来很大的误差，导致图像不能够满足临床诊断的要求，需要对图像进行处理。

如图2所示，所述步骤s100具体包括：

步骤s110、通过圆周扫描采集锥形束投影数据。

具体地，可通过模拟ct系统或ct医学成像系统直接采用圆周扫描的方式采集锥形束投影数据。

步骤s120、采用fdk算法对圆周扫描得到的锥形束投影数据重建得到初始图像。

本发明中，对于初始图像的重建算法可采用很多种，除fdk算法外，还可采用art，sart等算法，这些算法均可应用在图像识别上，实现把锥形束投影数据重建得到初始图像。本发明中优选fdk算法或其改进算法，其对锥束几何圆形扫描轨迹实现三维图像重建，实现图像快速重建。

步骤s200、对初始图像进行非线性处理得到矫正图像。

具体地，得到初始图像后，对初始图像进行矫正，矫正过程所采用的是非线性变化的方式，使处理后的矫正图像更接近于需要的目标图像，即不带伪影的图像。非线性变化使图像接近ct所拍摄的物体本身的图像，是量化处理图像的重要方法之一。本实施例中的所述非线性处理为采用多层神经网络处理，所述多层神经网络处理通过优化目标图像和输入图像的差异来使得输出图像尽可能接近理想图像。

多层神经网络可以是基于近年来人工智能的神经网络训练方法来处理图像。输出图像就是带有锥束伪影的重建图像，目标图像就是没有锥形束伪影的原始图像。可以通过优化目标图像和输出图像的差异来使得输出图像尽可能的接近理想图像，输出图像在接近理想图像的过程中对神经网络训练。

采用训练后的神经网络，能使初始图像快速得到矫正图像，且能保证矫正图像的成像质量。

如图2所示，所述步骤s200具体包括：

步骤s210、构建图像数据库，图像数据库的样本包括图像中无伪影的目标图像，以及存在锥形束伪影的输出图像。

具体地，建立图像数据库，图像数据库中包括很多样本，每一个样本包括目标图像和输出图像，其中目标图像通过实际系统或仿真的投影数据生成，目标图像是理想的重建图像，图像中无伪影。输出图像通过采集得到，输出图像存在锥束伪影。

步骤s220、构建多尺度卷积神经网络模型。

具体地，所述多尺度卷积神经网络模型包括一层或多层卷积层且还包括池化层，多尺度卷积神经网络模型为改进的resnet、改进的unet或者其他应用在图像优化处理上的多尺度卷积神经网络模型。resnet或unet作为解析算法，需要根据图像的锥束伪影应用进行与之相关联的改进，因此改进的resnet、改进的unet是根据本发明中的图像处理而进行改进后适用于本发明中的resnet、unet。本实施中优先采用resnet解析算法，随着网络的加深，出现了训练集准确率下降的现象，resnet可以解决随着网络加深，准确率不下降的问题，应用在解决锥束伪影的多尺度卷积神经网络模型上对图像处理有更好的效果。

步骤s230、构建模型的目标函数。

本发明中，所述目标函数为：

其中imgerr为图像的误差，img为网络的输出图像，imgtarg是数据库里面的目标图像；imgk为输出图像中坐标为k的像素点，imgtarg,k为数据库里面的目标图像中坐标为k的像素点。

通过目标函数，使输出图像与目标图像进行对比，训练优化多尺度卷积神经网络模型。

步骤s240、利用图像数据库和模型按照目标函数训练神经网络。

具体地，通过图像数据库中的输出图像和目标图像，利用梯度下降法针对目标函数，对神经网络的参数进行优化。训练完成后，整个神经网络的参数被保留下来，供后续过程使用。

步骤s250、利用训练好的神经网络处理初始图像得到矫正图像。

具体地，初始图像通过训练好的神经网络处理后得到矫正图像，此时的矫正图像在成像效果上尽可能的接近理想图像。

步骤s300、对矫正图像进行数值模拟产生矫正锥形束投影数据。

具体地，对矫正图像进行数值化处理，数值模拟产生矫正锥形束投影数据，数值模拟过程为模拟ct系统的成像过程，这一过程与实际的数值模拟过程相符合，从而矫正图像得到矫正锥形束投影数据更接近实际值。

步骤s400、对矫正锥形束投影数据和初始锥形束投影数据进行计算得到差异投影数据。

具体地，通过对矫正锥形束投影数据和初始锥形束投影数据融合得到的误差投影数据，融合计算包括直接相减、相加、部分相减、部分相加、线性组合等，本实施例中的融合过程只要可以对比得到误差投影数据的方式均可适应，本实施例中优先采用直接相减的方式实现融合。误差投影数据得出矫正的部分数据。

步骤s500、通过图像重建技术对差异投影数据处理得到差异图像。

具体地，其图像重建技术与步骤通过图像重建技术对初始锥形束投影数据重建得到初始图像中的图像重建技术为同一图像重建技术，可均采用fdk算法重建。这样采用同一图像重建技术对图像的处理结果具有一致性，便于后续的图像处理。

步骤s600、差异图像和矫正图像融合得到最终图像。

具体地，差异图像和矫正图像融合得到最终图像，实现图像的重建，融合过程包括直接相减、相加、部分相减、部分相加、线性组合等，使伪影在矫正过程中消除，本实施例中通过差异图像和矫正图像进行相加得到更接近实际的最终图像，即直接把两幅图像加到一起，相当于带有伪影的差异图像加上矫正伪影的矫正图像，得到伪影重建的图像。

如图2、图3所示，本发明方法在现有的重建技术的基础上只是加入了投影模拟和非线性变换，计算量相对于原来的方法增加不多，相比迭代算法要简单的多。同时还能取得好的图像重建效果，具有通用性。

如图6所示，本发明还提出一种ct图像重建系统，其中，包括有处理器10，以及与所述处理器连接的存储器20。

所述存储器20存储有圆周扫描ct锥形束伪影图像重建程序，所述圆周扫描ct锥形束伪影图像重建程序被所述处理器10执行时实现如上所述的ct图像重建方法。

综上所述，本发明所提供的一种ct图像重建方法及其系统，所述方法包括步骤：通过图像重建技术对圆周扫描得到的初始锥形束投影数据重建得到初始图像；对初始图像进行非线性处理得到矫正图像；通过对矫正图像进行数值模拟产生矫正锥形束投影数据；对矫正锥形束投影数据和初始锥形束投影数据进行计算得到差异投影数据；通过图像重建技术对差异投影数据处理得到差异图像；差异图像和矫正图像融合得到最终图像。从而使带伪影的差异图像和矫正伪影的矫正图像结合实现锥形束伪影图像重建；通过在现有的重建技术的基础上只加入了投影模拟和非线性变换，在伪影图像重建过程中实现简单而较少的计算量，相比迭代算法要更加简单，同时还能取得好的图像重建效果，满足临床的实际应用，具有通用性。

应当理解的是，本发明的应用不限于上述的举例，对本领域普通技术人员来说，可以根据上述说明加以改进或变换，所有这些改进和变换都应属于本发明所附权利要求的保护范围。