图像重建方法、装置和计算机设备与流程

文档序号:17551579发布日期:2019-04-30 18:19阅读:210来源:国知局
图像重建方法、装置和计算机设备与流程

本发明涉及图像领域,特别是涉及一种图像重建方法、装置和计算机设备。



背景技术:

正电子发射计算机断层显像-电子计算机断层扫描(positronemissioncomputedtomography-computedtomography,pet-ct)检查和正电子发射计算机断层显像/核磁共振成像(positronemissioncomputedtomography/magneticresonanceimaging,pet/mri)检查是两种常见的医学检查,pet-ct和pet/mri扫描的辐射剂量有一部分来自pet药物,减少辐射的手段之一就是减少药物的注射剂量,另一方面,影响pet图像质量的一个因素是数据量的多少,采集的数据量足够多,重建出来的pet图像质量才能满足临床的需要,而采集的数据量又依赖事先注入体内的药物剂量以及采集时间,目前的pet数据采集在注射正常的剂量情况下通常需要采集3分钟每床位,如果注入一半剂量的话则至少需要采集6分钟的数据才能重建出满足临床需要的pet图像。因此,在降低pet药物注射剂量的同时,不延长pet-ct或者pet/mri的扫描时间获得满足临床需要的pet,就需要对低剂量pet图像进行重建,获取标准剂量的pet图像。

传统的基于学习类的低剂量pet图像重建方法,大多是利用传统的机器学习技术,对低剂量pet图像利用ct或者mri图像的信息,进行标准剂量pet图像的预测,进而得到标准剂量的pet图像。

但是,传统的低剂量pet图像重建方法,存在重建图像的质量低、重建时间长的问题。



技术实现要素:

基于此,有必要针对传统的基于学习类的低剂量pet图像重建方法,存在重建图像质量低、重建时间长的问题,提供一种图像重建方法、装置和计算机设备。

第一方面,本发明实施例提供一种图像重建方法,所述方法包括:

获取测试样本图像;

对所述测试样本图像进行预处理,获取预处理后的测试样本图像;

对所述预处理后的测试样本图像进行梯度计算,获取测试样本梯度图像;

对所述测试样本梯度图像进行归一化处理,获取归一化后的测试样本梯度图像;

将所述预处理后的测试样本图像和所述归一化后的测试样本梯度图像输入预设的深度学习网络模型,获得重建图像。

在其中一个实施例中,所述预处理包括重采样处理和归一化处理。

在其中一个实施例中,所述对测试样本图像进行预处理,获取预处理后的测试样本图像,包括:

判断所述测试样本图像与训练样本图像的分辨率是否相同;

若是,对所述测试样本图像进行归一化处理,得到归一化后的测试样本图像;

若否,则根据所述训练样本图像的分辨率对所述测试样本进行重采样,得到重采样后的测试样本图像,并对所述重采样后的测试样本图像进行归一化处理,得到归一化后的测试样本图像。

在其中一个实施例中,所述将所述预处理后的测试样本图像和所述归一化后的测试样本梯度图像输入预设的深度学习网络模型,获得重建图像之后,所述方法还包括:

获取用于归一化处理的测试样本图像像素值的最大值和最小值;

根据所述重建图像的像素值、所述用于归一化处理的测试样本图像像素值的最大值和最小值,获取目标重建图像。

在其中一个实施例中,所述对所述预处理后的测试样本图像进行梯度计算,获取测试样本梯度图像,包括:

采用预设的图像梯度特征算子对所述测试样本图像进行3d梯度求解,得到所述测试样本梯度图像。

在其中一个实施例中,所述将所述预处理后的测试样本图像和所述归一化后的测试样本梯度图像输入预设的深度学习网络模型,获得重建图像之前,所述方法还包括:

获取所述预处理后的测试样本图像各视角的二维切面图和所述归一化后的测试样本梯度图像各视角的二维切面图;

所述将所述预处理后的测试样本图像和所述归一化后的测试样本梯度图像输入所述深度学习网络模型,获得重建图像,包括:

将所述预处理后的测试样本图像各视角的二维切面图和所述归一化后的测试样本梯度图像各视角的二维切面图输入所述深度学习网络模型,获得所述重建图像。

在其中一个实施例中,所述深度学习网络模型包括:横断位分支模型、矢状位分支模型、冠状位分支模型和整合分支模型。

在其中一个实施例中,所述将所述预处理后的测试样本图像各视角的二维切面图和所述归一化后的测试样本梯度图像各视角的二维切面图输入所述深度学习网络模型,获得所述重建图像,包括:

将所述预处理后的测试样本图像的横断位二维切面图和所述归一化后的测试样本梯度图像的横断位的二维切面图连接,输入所述横断位分支模型,得到第一特征图像;

将所述预处理后的测试样本图像的矢状位二维切面图和所述归一化后的测试样本梯度图像的矢状位的二维切面图连接,输入所述矢状位分支模型,得到第二特征图像;

将所述预处理后的测试样本图像的冠状位二维切面图和所述归一化后的测试样本梯度图像的冠状位的二维切面图连接,输入所述冠状位分支模型,得到第三特征图像;

将所述第一特征图像、所述第二特征图像和所述第三特征图像连接,输入所述整合分支模型,获取所述重建图像。

第二方面,本发明实施例提供一种图像重建装置,所述装置包括:

获取模块,用于获取测试样本图像;

预处理模块,用于对所述测试样本图像进行预处理,获取预处理后的测试样本图像;

梯度计算模块,用于对所述预处理后的测试样本图像进行梯度计算,获取测试样本梯度图像;

归一化处理模块,用于对所述测试样本梯度图像进行归一化处理,获取归一化后的测试样本梯度图像;

图像重建模块,用于将所述预处理后的测试样本图像和所述归一化后的测试样本梯度图像输入预设的深度学习网络模型,获得重建图像。

第三方面,本发明实施例提供一种计算机设备,包括存储器和处理器,所述存储器存储有计算机程序,所述处理器执行所述计算机程序时实现以下步骤:

获取测试样本图像;

对所述测试样本图像进行预处理,获取预处理后的测试样本图像;

对所述预处理后的测试样本图像进行梯度计算,获取测试样本梯度图像;

对所述测试样本梯度图像进行归一化处理,获取归一化后的测试样本梯度图像;

将所述预处理后的测试样本图像和所述归一化后的测试样本梯度图像输入预设的深度学习网络模型,获得重建图像。

上述实施例提供的图像重建方法、装置和计算机设备中,计算机设备可以对获取的测试样本图像进行预处理,得到预处理后的测试样本图像,然后,对预处理后的测试样本图像进行梯度计算,获取测试样本梯度图像,对测试样本梯度图像进行归一化处理,获取归一化后的测试样本梯度图像,之后,将预处理后的测试样本图像和归一化后的测试样本梯度图像输入预设的深度学习网络模型,得到重建图像。该图像重建方法中,重建图像是将预处理后的测试样本图像和归一化后的测试样本梯度图像输入预设的深度学习网络模型得到的,由于将归一化后的测试样本梯度图像作为输入可以更好的恢复图像细节,使重建出的图像更加清晰,提高了得到的重建图像的质量;另外,采用预设的深度学习模型可以直接快速的得到重建图像,减少了获取重建图像的时长。

附图说明

图1为一个实施例提供的计算机设备的内部结构示意图;

图2为一个实施例提供的图像重建方法的流程示意图;

图3(a)为一个实施例提供的测试样本图像;

图3(b)为一个实施例提供的测试样本图像的重建图像;

图4为另一个实施例提供的图像重建方法的流程示意图;

图5为另一个实施例提供的图像重建方法的流程示意图;

图6为另一个实施例提供的图像重建方法的流程示意图;

图7为一个实施例提供的获取重建图像的过程示意图;

图8为一个实施例提供的对初始深度学习网络模型进行训练的过程示意图;

图9为一个实施例提供的图像重建装置结构示意图;

图10为一个实施例提供的图像重建装置结构示意图;

图11为一个实施例提供的图像重建装置结构示意图;

图12为一个实施例提供的图像重建装置结构示意图;

图13为一个实施例提供的图像重建装置结构示意图。

具体实施方式

为了使本申请的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本申请进行进一步详细说明。应当理解,此处描述的具体实施例仅仅用以解释本申请,并不用于限定本申请。

本申请实施例提供的图像重建方法,可以适用于如图1所示的计算机设备。该计算机设备包括通过系统总线连接的处理器、存储器,该存储器中存储有计算机程序,处理器执行该计算机程序时可以执行下述方法实施例的步骤。可选的,该计算机设备还可以包括网络接口、显示屏和输入装置。其中,该计算机设备的处理器用于提供计算和控制能力。该计算机设备的存储器包括非易失性存储介质、内存储器,该非易失性存储介质存储有操作系统和计算机程序。该内存储器为非易失性存储介质中的操作系统和计算机程序的运行提供环境。该计算机设备的网络接口用于与外部的终端通过网络连接通信。可选的,该计算机设备可以是服务器,可以是pc,还可以是个人数字助理,还可以是其他的终端设备,例如pad、手机等等,还可以是云端或者远程服务器,本申请实施例对计算机设备的具体形式并不做限定。

需要说明的是,本申请实施例提供的图像重建方法,其执行主体可以是图像重建装置,该图像重建装置可以通过软件、硬件或者软硬件结合的方式实现成为计算机设备的部分或者全部。下述方法实施例中,均以执行主体是计算机设备为例来进行说明。

下面以具体的实施例对本发明的技术方案以及本发明的技术方案如何解决上述技术问题进行详细说明。下面这几个具体的实施例可以相互结合,对于相同或相似的概念或过程可能在某些实施例中不再赘述。

图2为一个实施例提供的图像重建方法的流程示意图。图3(a)为一个实施例提供的测试样本图像,图3(b)为一个实施例提供的测试样本图像的重建图像。本实施例涉及的是计算机设备获得重建图像的具体实现过程。如图2所示,该方法可以包括:

s201,获取测试样本图像。

其中,测试样本图像为低剂量的正电子发射计算机断层显像-电子计算机断层扫描(positronemissioncomputedtomography-computedtomography,pet-ct)的pet图像,或者,低剂量的正电子发射计算机断层显像/核磁共振成像(positronemissioncomputedtomography/magneticresonanceimaging,pet/mri)的pet图像。在本实施例中,可选的,计算机设备可以从pacs(picturearchivingandcommunicationsystems,影像归档和通信系统)服务器中获取低剂量的pet图像,也可以直接从pet-ct仪或pet/mri影像设备中获取低剂量的pet图像。

s202,对所述测试样本图像进行预处理,获取预处理后的测试样本图像。

具体的,计算机设备获取到上述测试样本图像后,对上述测试样本图像进行预处理,得到预处理后的测试样本图像。其中,对测试样本图像进行预处理是指,通过已有的算法或公式,将测试样本图像处理为与训练样本图像形式相同的图像,例如,可以将测试样本图像处理为与训练样本图像分辨率相同的图像。

s203,对所述预处理后的测试样本图像进行梯度计算,获取测试样本梯度图像。

具体的,计算机设备得到预处理后的测试样本图像后,对预处理后的测试样本图像进行梯度计算,得到测试样本梯度图像。可选的,计算机设备可以采用预设的图像梯度特征算子对预处理后的测试样本图像进行3d梯度求解,得到测试样本梯度图像,例如,图像梯度特征算子可以为索贝尔算子(sobeloperator)、普利维特算子(prewitoperator)、罗伯茨交叉边缘检测算子(robertscrossoperator)等。

s204,对所述测试样本梯度图像进行归一化处理,获取归一化后的测试样本梯度图像。

具体的,计算机设备对上述获取的测试样本梯度图像,利用归一化处理公式进行归一化处理,获取归一化后的测试样本梯度图像。其中,图像归一化是指对图像进行一系列标准的处理变换,使图像变换为一固定标准形式;可选地,归一化处理可以采用公式式中i为测试样本梯度图像,max为测试样本梯度图像像素值的最大值,min为测试样本梯度图像像素值的最小值,i'为归一化后的测试样本梯度图像。

s205,将所述预处理后的测试样本图像和所述归一化后的测试样本梯度图像输入预设的深度学习网络模型,获得重建图像。

具体的,计算机设备将上述获取的预处理后的测试样本图像,和归一化后的测试样本梯度图像,输入预设的深度学习网络模型,得到重建图像。其中,重建图像为测试的低剂量样本图像对应的标准剂量的图像,预设的深度学习网络模型为预先训练的深度学习网络,用于将低剂量的pet图像重建为标准剂量的pet图像。例如,测试样本图像为pet-ct设备的低剂量pet图像,则重建图像为该pet-ct设备低剂量pet图像对应的标准剂量pet图像。如图3(a)和图3(b)所示,图3(a)为测试样本图像,图3(b)为测试样本图像的重建图像。

在本实施例中,重建图像是将预处理后的测试样本图像和归一化后的测试样本梯度图像输入预设的深度学习网络模型得到的,由于将归一化后的测试样本梯度图像作为输入可以更好的恢复图像细节,使重建出的图像更加清晰,提高了得到的重建图像的质量;另外,采用预设的深度学习模型可以直接快速的得到重建图像,减少了获取重建图像的时长。

在上述对测试样本图像进行预处理的场景中,作为一种可选的实施方式,所述预处理包括重采样处理和归一化处理。

其中,重采样处理是指将测试样本图像重采样为与训练样本图像的分辨率相同的图像;归一化处理是指指对图像进行一系列标准的处理变换,使图像变换为一固定标准形式。可选的,计算机设备可以采用最近邻内插法对测试样本图像进行重采样,也可以采用双线性内插法或三次卷积法对测试样本图像进行重采样。例如,测试样本图像的分辨率为100ppi,训练样本图像的分辨率为80ppi,则可以对测试样本图像采用双线性内插法,将测试样本图像的分辨率重采样为80ppi。可选的,计算机设备可以根据上述归一化处理公式对测试样本图像进行归一化处理。

在本实施例中,计算机设备对测试样本图像的预处理包括重采样处理和归一化处理,对测试样本进行的重采样处理,能够使测试样本图像的分辨率与训练样本图像的分辨率相同,提高了得到的重建图像的质量;另外,由于预设的深度学习网络模型是根据归一化处理后的样本图像训练得到,因此,对测试样本图像进行归一化处理,有利于测试样本图像与预设的深度学习网络模型匹配,从而提高得到的重建图像的质量。

图4为另一个实施例提供的图像重建方法的流程示意图。本实施例涉及的是对测试样本图像进行预处理,获取预处理后的测试样本图像的具体实现过程。如图4所示,在上述实施例的基础上,作为一种可选的实施方式,上述s202包括:

s401,判断所述测试样本图像与训练样本图像的分辨率是否相同,若是,则执行s402;若否,则执行s403。

具体的,计算机设备可以将测试样本图像的分辨率与训练样本图像的分辨率进行比对,判断测试样本图像与训练样本图像的分辨率是否相同。若测试样本图像与训练样本图像的分辨率相同,则对测试样本图像进行归一化处理;若测试样本图像与训练样本图像的分辨率不同,则先对测试样本进行重采样,再对重采样后的训练图像进行归一化处理。

s402,对所述测试样本图像进行归一化处理,得到归一化后的测试样本图像。

具体的,计算机设备根据上述归一化处理公式对测试样本图像进行归一化处理,得到归一化处理后的测试样本图像。

s403,根据所述训练样本图像的分辨率对所述测试样本进行重采样,得到重采样后的测试样本图像,并对所述重采样后的测试样本图像进行归一化处理,得到归一化后的测试样本图像。

具体的,计算机设备利用已有的重采样方法,根据训练样本图像的分辨率,对测试样本图像进行重采样,得到与训练样本图像的分辨率相同的重采样后的测试样本图像,并利用上述归一化处理公式对重采样后的测试样本图像进行归一化处理,得到归一化后的测试样本图像。可选的,计算机设备可以采用最近邻内插法对测试样本图像进行重采样,也可以采用双线性内插法或三次卷积法对测试样本图像进行重采样。

在本实施例中,计算机设备首先判断测试样本图像的分辨率与训练样本图像的分辨率是否相同,对与训练样本图像的分辨率相同的测试样本图像直接进行归一化处理,对与训练样本图像的分辨率不同的测试样本图像进行重采样,得到重采样后的测试样本图像,并对重采样后的测试样本图像进行归一化处理,得到归一化处理后的测试样本图像,由于计算机设备首先对测试样本图像的分辨率与训练样本图像的分辨率进行了判断,使得测试样本和训练样本的分辨率达到统一,提高了得到归一化处理后的测试样本图像的效率;另外,由于预设的深度学习网络模型是根据归一化处理后的样本图像训练得到,因此,对测试样本图像进行归一化处理,有利于测试样本图像与预设的深度学习网络模型匹配,从而提高得到的重建图像的质量。

在上述将预处理后的测试样本图像和归一化后的测试样本梯度图像输入预设的深度学习网络模型,获得重建图像的场景中,由于获得的重建图像为归一化后的图像,因此,需要对得到的重建图像再进行反归一化处理,得到目标重建图像。图5为另一个实施例提供的图像重建方法的流程示意图。本实施例涉及的是获取目标重建图像的具体实现过程。如图5所示,在上述实施例的基础上,作为一种可选的实施方式,所述方法还包括:

s501,获取用于归一化处理的测试样本图像像素值的最大值和最小值。

具体的,在上述对测试样本图像进行归一化处理时,需要获取测试样本图像像素值的最大值和最小值,在本实施例中,计算机设备首先获取对测试样本图像进行归一化处理时,得到的图像的最大值和最小值。

s502,根据所述重建图像的像素值、所述用于归一化处理的测试样本图像像素值的最大值和最小值,获取目标重建图像。

具体的,计算机设备根据上述得到的重建图像,获取重建图像的像素值,根据重建图像的像素值和上述用于归一化处理的测试样本图像像素值的最大值和最小值,获取目标重建图像。可选的,可以利用公式i=i'*(max-min)+min,获取目标重建图像,式中,i表示目标重建图像,i'表示重建图像,max为测试样本梯度图像像素值的最大值,min为测试样本梯度图像像素值的最小值。

在本实施例中,由于获得的重建图像为归一化后的图像,计算机设备根据重建图像的像素值、用于归一化处理的测试样本图像像素值的最大值和最小值,对重建图像进行反归一化处理,得到的目标重建图像的形式与归一化处理前的图像形式相同,提高了得到的目标重建图像的准确度;另外,计算机设备根据重建图像的像素值、用于归一化处理的测试样本图像像素值的最大值和最小值,获取目标重建图像的过程十分简单,提高了得到目标重建图像的效率。

在上述实施例的基础上,作为一种可选的实施方式,上述s205之前,所述方法还包括:获取所述预处理后的测试样本图像各视角的二维切面图和所述归一化后的测试样本梯度图像各视角的二维切面图;上述s205,包括:将所述预处理后的测试样本图像各视角的二维切面图和所述归一化后的测试样本梯度图像各视角的二维切面图输入所述深度学习网络模型,获得所述重建图像。

具体的,计算机设备获得上述预处理后的测试样本图像和测试样本梯度图像后,分别取预处理后的测试样本图像各视角的二维切面图(2dslice)和归一化后的测试样本梯度图像各视角的二维切面图(2dslice);其中,预处理后的测试样本图像各视角的二维切面图(2dslice)包括预处理后的测试样本图像横断位(axialview)的二维切面图(2dslice)、矢状位(sagittalview)的二维切面图(2dslice)和冠状位(coronalview)的二维切面图(2dslice),归一化后的测试样本梯度图像包括归一化后的测试样本梯度图像横断位(axialview)的二维切面图(2dslice)、矢状位(sagittalview)的二维切面图(2dslice)和冠状位(coronalview)的二维切面图(2dslice)。之后,计算机设备将上述获取的预处理后的测试样本图像各视角的二维切面图(2dslice)和归一化后的测试样本梯度图像各视角的二维切面图(2dslice),输入预设的深度学习网络模型,获得重建图像。

在本实施例中,计算机设备首先获取预处理后的测试样本图像各视角的二维切面图和归一化后的测试样本梯度图像各视角的二维切面图,将预处理后的测试样本图像各视角的二维切面图和归一化后的测试样本梯度图像各视角的二维切面图输入预设的深度学习网络模型,有利于深度学习网络模型充分学习预处理后的测试样本图像和归一化后的测试样本梯度图像各视角的特征,提高了得到的重建图像的质量。

在上述将预处理后的测试样本图像和归一化后的测试样本梯度图像输入预设的深度学习网络模型的场景中,所述深度学习网络模型包括:横断位分支模型、矢状位分支模型、冠状位分支模型和整合分支模型。图6为另一个实施例提供的图像重建方法的流程示意图。图7为一个实施例提供的获取重建图像的过程示意图。本实施例涉及的是将预处理后的测试样本图像各视角的二维切面图和归一化后的测试样本梯度图像各视角的二维切面图输入所述深度学习网络模型,获得重建图像的具体实现过程。如图6所示,在上述实施例的基础上,作为一种可选的实施方式,所述将所述预处理后的测试样本图像各视角的二维切面图和所述归一化后的测试样本梯度图像各视角的二维切面图输入所述深度学习网络模型,获得所述重建图像,包括:

s601,将所述预处理后的测试样本图像的横断位二维切面图和所述归一化后的测试样本梯度图像的横断位的二维切面图连接,输入所述横断位分支模型,得到第一特征图像。

具体的,计算机设备获取到预处理后的测试样本图像的横断位二维切面图(axialview2dslice)和归一化后的测试样本梯度图像的横断位二维切面图(axialview2dslice)后,如图7所示,将预处理后的测试样本图像的二维切面图(axialview2dslice)与归一化后的测试样本梯度图像的横断位二维切面图(axialview2dslice)连接,将连接的横断位二维切面图(axialview2dslice)做2.5d处理,并输入横断位分支模型(network-branch-1),得到第一特征图像(map_1)。示例性的,2.5d处理指的是预处理后的测试样本图像的5张2dslice与归一化后的测试样本梯度图像的5张2dslice,对应标准剂量pet图像的一张2dslice。例如,axialview2dslice的大小为64*64,则输入到network-branch-1的axialview2dslice的大小为64*64*10。

s602,将所述预处理后的测试样本图像的矢状位二维切面图和所述归一化后的测试样本梯度图像的矢状位的二维切面图连接,输入所述矢状位分支模型,得到第二特征图像。

具体的,计算机设备获取到预处理后的测试样本图像的矢状位二维切面图(sagittalview2dslice)和归一化后的测试样本梯度图像的矢状位的二维切面图(sagittalview2dslice)后,如图7所示,将预处理后的测试样本图像的矢状位二维切面图(sagittalview2dslice)与归一化后的测试样本梯度图像的矢状位的二维切面图(sagittalview2dslice)连接,将连接的矢状位二维切面图(sagittalview2dslice)做2.5d处理,并输入矢状位分支模型(network-branch-2),得到第二特征图像(map_2)。例如,sagittalview2dslice的大小为64*64,则输入到network-branch-2的sagittalview2dslice的大小为64*64*10。

s603,将所述预处理后的测试样本图像的冠状位二维切面图和所述归一化后的测试样本梯度图像的冠状位的二维切面图连接,输入所述冠状位分支模型,得到第三特征图像。

具体的,计算机设备获取到预处理后的测试样本图像的冠状位二维切面图(coronalview2dslice)和归一化后的测试样本梯度图像的冠状位的二维切面图(coronalview2dslice)后,如图7所示,将预处理后的测试样本图像的冠状位二维切面图(coronalview2dslice)与归一化后的测试样本梯度图像的冠状位的二维切面图(coronalview2dslice)连接,将连接的冠状位二维切面图(coronalview2dslice)做2.5d处理,并输入冠状位分支模型(network-branch-2),得到第三特征图像(map_3)。例如,coronalview2dslice的大小为64*64,则输入到network-branch-3的coronalview2dslice的大小为64*64*10。

s604,将所述第一特征图像、所述第二特征图像和所述第三特征图像连接,输入所述整合分支模型,获取所述重建图像。

具体的,计算机设备获取到上述第一特征图像(map_1)、第二特征图像(map_2)和第三特征图像(map_3)后,如图7所示,将第一特征图像(map_1)、第二特征图像(map_2)和第三特征图像(map_3)连接,输入整合分支模型(network-branch-4),得到重建图像。需要说明的是,得到的重建图像为具有横断位(axialview)、矢状位(sagittalview)和冠状位(coronalview)的三维图像。

在本实施例中,计算机设备分别将预处理后的测试样本图像和归一化后的测试样本梯度图像各视角的二维切面图连接,将各视角的二维连接切面图分别输入对应的分支模型中,有利于各分支模型充分学习各视角的特征,提高了得到的各视角的特征图像的准确度;另外,再将各视角的特征图像进行连接,将连接后的特征图像输入整合分支模型,得到重建图像,由于重建图像是根据各视角的特征图像得到的,因此,进一步的提高了得到的重建图像的准确度,提高了得到的重建图像的质量。

图8为一个实施例提供的对初始深度学习网络模型进行训练的过程示意图。在将预处理后的测试样本图像和归一化后的测试样本梯度图像输入预设的深度学习网络模型的场景中,需要对初始深度学习网络模型进行训练,在上述实施例的基础上,如图8所示,所述方法还包括:根据多幅训练样本图像和对应的多幅标准剂量pet图像,采用预设的auto-context策略,将各训练样本图像各视角的二维切面图以及对应的梯度图各视角的二维切面图连接,分别输入初始深度学习网络模型中,获得重建测试图像;将上一次训练得到的重建测试图像各视角的图像块和训练样本图像各视角的图像块以及对应的梯度图各视角的图像块连接,输入本次训练的深度学习网络模型中;将每一次训练的深度学习网络模型输出的重建测试图像的图像块和所述标准剂量pet图像的图像块输入预设的损失函数(lossfunction)中,直至所述损失函数(lossfunction)的值达到稳定值为止,将得到的深度学习网络模型确定为预设的深度学习网络模型。

具体的,如图8所示,计算机设备首先获取多幅训练样本图像,以及和多幅训练样本图像对应的多幅标准剂量pet图像;获得各训练样本图像以及对应的梯度图的各视角的二维切面图(2dslice)和各标准剂量pet图像各视角的二维切面图(2dslice);对获得的各训练样本图像以及对应的梯度图的各视角的二维切面图(2dslice)和各标准剂量pet图像各视角的二维切面图(2dslice)进行块提取,得到各训练样本图像以及对应的梯度图的各视角的图像块和各标准剂量pet图像各视角的图像块,其中,各训练样本图像以及对应的梯度图的各视角的图像块的大小为64*64*10,各标准剂量pet图像各视角的图像块的大小为64*64*1;之后,计算机设备采用预设的auto-context策略,将各训练样本图像各视角的图像块以及对应的梯度图的各视角的图像块连接,输入初始深度学习网络模型中,获得重建测试图像,初始深度学习模型也就是图8中的第一深度学习网络模型;然后,将本次得到的重建测试图像各视角的图像块和训练样本图像各视角的图像块以及对应的梯度图的各视角的图像块连接,输入下次训练的深度学习网络模型中,下次训练的深度学习网络模型也就是图8中的第二深度学习网络模型,同时,将每一次训练的深度学习网络模型输出的重建测试图像的图像块和标准剂量pet图像的图像块输入预设的损失函数(lossfunction)中,如此迭代训练,直至损失函数(lossfunction)的值达到稳定值为止,将损失函数的值达到稳定值时,得到的深度学习网络模型,确定为预设的深度学习网络模型。可选的,可以选择l1loss作为本实施例中深度学习网络模型训练的lossfunction。

在本实施例中,计算机设备采用auto-context策略对初始深度学习网络模型进行迭代训练,充分学习了各训练样本图像的特征,提高了确定的预设的深度学习网络模型的质量,进而提高了根据预设的深度学习网络模型得到的重建图像的质量。

应该理解的是,虽然图2-8的流程图中的各个步骤按照箭头的指示依次显示,但是这些步骤并不是必然按照箭头指示的顺序依次执行。除非本文中有明确的说明,这些步骤的执行并没有严格的顺序限制,这些步骤可以以其它的顺序执行。而且,图2-8中的至少一部分步骤可以包括多个子步骤或者多个阶段,这些子步骤或者阶段并不必然是在同一时刻执行完成,而是可以在不同的时刻执行,这些子步骤或者阶段的执行顺序也不必然是依次进行,而是可以与其它步骤或者其它步骤的子步骤或者阶段的至少一部分轮流或者交替地执行。

图9为一个实施例提供的图像重建装置结构示意图。如图9所示,该装置可以包括:第一获取模块10、预处理模块11、梯度计算模块12、归一化处理模块13和图像重建模块14。

具体的,第一获取模块10,用于获取测试样本图像;

预处理模块11,用于对所述测试样本图像进行预处理,获取预处理后的测试样本图像;

梯度计算模块12,用于对所述预处理后的测试样本图像进行梯度计算,获取测试样本梯度图像;

归一化处理模块13,用于对所述测试样本梯度图像进行归一化处理,获取归一化后的测试样本梯度图像;

图像重建模块14,用于将所述预处理后的测试样本图像和所述归一化后的测试样本梯度图像输入预设的深度学习网络模型,获得重建图像。

可选的,所述预处理包括重采样处理和归一化处理。

可选的,所述深度学习网络模型包括:横断位分支模型、矢状位分支模型、冠状位分支模型和整合分支模型。

本实施例提供的图像重建装置,可以执行上述方法实施例,其实现原理和技术效果类似,在此不再赘述。

图10为一个实施例提供的图像重建装置结构示意图。在上述实施例的基础上,可选的,如图10所示,预处理模块11包括判断单元111、第一处理单元112和第二处理单元113。

具体的,判断单元111,用于判断所述测试样本图像与训练样本图像的分辨率是否相同;

第一处理单元112,用于对所述测试样本图像进行归一化处理,得到归一化后的测试样本图像;

第二处理单元113,用于根据所述训练样本图像的分辨率对所述测试样本进行重采样,得到重采样后的测试样本图像,并对所述重采样后的测试样本图像进行归一化处理,得到归一化后的测试样本图像。

本实施例提供的图像重建装置,可以执行上述方法实施例,其实现原理和技术效果类似,在此不再赘述。

图11为一个实施例提供的图像重建装置结构示意图。在上述实施例的基础上,可选的,如图11所示,所述装置还包括:第二获取模块15和目标重建图像模块16。

具体的,第二获取模块15,用于获取用于归一化处理的测试样本图像像素值的最大值和最小值;

目标重建图像模块16,用于根据所述重建图像的像素值、所述用于归一化处理的测试样本图像像素值的最大值和最小值,获取目标重建图像。

本实施例提供的图像重建装置,可以执行上述方法实施例,其实现原理和技术效果类似,在此不再赘述。

在上述实施例的基础上,可选的,梯度计算模块12,具体用于采用预设的图像梯度特征算子对所述测试样本图像进行3d梯度求解,得到所述测试样本梯度图像。

本实施例提供的图像重建装置,可以执行上述方法实施例,其实现原理和技术效果类似,在此不再赘述。

图12为一个实施例提供的图像重建装置结构示意图。在上述实施例的基础上,可选的,如图12所示,所述装置还包括:第三获取模块17。

具体的,第三获取模块17,用于获取所述预处理后的测试样本图像各视角的二维切面图和所述归一化后的测试样本梯度图像各视角的二维切面图;

图像重建模块14,具体用于将所述预处理后的测试样本图像各视角的二维切面图和所述归一化后的测试样本梯度图像各视角的二维切面图输入所述深度学习网络模型,获得所述重建图像。

本实施例提供的图像重建装置,可以执行上述方法实施例,其实现原理和技术效果类似,在此不再赘述。

图13为一个实施例提供的图像重建装置结构示意图。在上述实施例的基础上,可选的,如图13所示,图像重建模块14包括横断位重建单元141、矢状位重建单元142、冠状位重建单元143和整合单元144。

具体的,横断位重建单元141,用于将所述预处理后的测试样本图像的横断位二维切面图和所述归一化后的测试样本梯度图像的横断位的二维切面图连接,输入所述横断位分支模型,得到第一特征图像;

矢状位重建单元142,用于将所述预处理后的测试样本图像的矢状位二维切面图和所述归一化后的测试样本梯度图像的矢状位的二维切面图连接,输入所述矢状位分支模型,得到第二特征图像;

冠状位重建单元143,用于将所述预处理后的测试样本图像的冠状位二维切面图和所述归一化后的测试样本梯度图像的冠状位的二维切面图连接,输入所述冠状位分支模型,得到第三特征图像;

整合单元144,用于将所述第一特征图像、所述第二特征图像和所述第三特征图像连接,输入所述整合分支模型,获取所述重建图像。

本实施例提供的图像重建装置,可以执行上述方法实施例,其实现原理和技术效果类似,在此不再赘述。

关于图像重建装置的具体限定可以参见上文中对于图像重建方法的限定,在此不再赘述。上述图像重建装置中的各个模块可全部或部分通过软件、硬件及其组合来实现。上述各模块可以硬件形式内嵌于或独立于计算机设备中的处理器中,也可以以软件形式存储于计算机设备中的存储器中,以便于处理器调用执行以上各个模块对应的操作。

在一个实施例中,提供了一种计算机设备,包括存储器和处理器,存储器中存储有计算机程序,该处理器执行计算机程序时实现以下步骤:

获取测试样本图像;

对所述测试样本图像进行预处理,获取预处理后的测试样本图像;

对所述预处理后的测试样本图像进行梯度计算,获取测试样本梯度图像;

对所述测试样本梯度图像进行归一化处理,获取归一化后的测试样本梯度图像;

将所述预处理后的测试样本图像和所述归一化后的测试样本梯度图像输入预设的深度学习网络模型,获得重建图像。

上述实施例提供的计算机设备,其实现原理和技术效果与上述方法实施例类似,在此不再赘述。

在一个实施例中,提供了一种计算机可读存储介质,其上存储有计算机程序,计算机程序被处理器执行时实现以下步骤:

获取测试样本图像;

对所述测试样本图像进行预处理,获取预处理后的测试样本图像;

对所述预处理后的测试样本图像进行梯度计算,获取测试样本梯度图像;

对所述测试样本梯度图像进行归一化处理,获取归一化后的测试样本梯度图像;

将所述预处理后的测试样本图像和所述归一化后的测试样本梯度图像输入预设的深度学习网络模型,获得重建图像。

上述实施例提供的计算机可读存储介质,其实现原理和技术效果与上述方法实施例类似,在此不再赘述。

本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程,是可以通过计算机程序来指令相关的硬件来完成,所述的计算机程序可存储于一非易失性计算机可读取存储介质中,该计算机程序在执行时,可包括如上述各方法的实施例的流程。其中,本申请所提供的各实施例中所使用的对存储器、存储、数据库或其它介质的任何引用,均可包括非易失性和/或易失性存储器。非易失性存储器可包括只读存储器(rom)、可编程rom(prom)、电可编程rom(eprom)、电可擦除可编程rom(eeprom)或闪存。易失性存储器可包括随机存取存储器(ram)或者外部高速缓冲存储器。作为说明而非局限,ram以多种形式可得,诸如静态ram(sram)、动态ram(dram)、同步dram(sdram)、双数据率sdram(ddrsdram)、增强型sdram(esdram)、同步链路(synchlink)dram(sldram)、存储器总线(rambus)直接ram(rdram)、直接存储器总线动态ram(drdram)、以及存储器总线动态ram(rdram)等。

以上所述实施例的各技术特征可以进行任意的组合,为使描述简洁,未对上述实施例中的各个技术特征所有可能的组合都进行描述,然而,只要这些技术特征的组合不存在矛盾,都应当认为是本说明书记载的范围。

以上所述实施例仅表达了本发明的几种实施方式,其描述较为具体和详细,但并不能因此而理解为对发明专利范围的限制。应当指出的是,对于本领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干变形和改进,这些都属于本发明的保护范围。因此,本发明专利的保护范围应以所附权利要求为准。

当前第1页1 2 
网友询问留言 已有0条留言
  • 还没有人留言评论。精彩留言会获得点赞!
1