CT图片的配准方法与流程

本发明涉及ct图片的处理方法，尤其涉及一种ct图片的配准方法。

背景技术：

每次做肺部ct都会得到几百张ct图片，这些图片代表身体不同切片层的图像，由于每次做ct扫描时扫描的起始位置不同，扫描的层厚也不同，所以要从不同时期的ct图片中根据身体的部位进行一一匹配的难度很大。

技术实现要素：

本发明要解决的技术问题是：为了解决现有技术中对于不同时期的ct图片配准难度大，准确性低的问题，本发明提供了一种ct图片的配准方法来解决上述问题。

本发明解决其技术问题所采用的技术方案是：一种ct图片的配准方法，包括以下步骤：

s1、从第一ct图片集中选取一个切片层的1张ct图片和前后2张ct图片，从第二ct图片集中选取同一切片层的1张ct图片和前后2张ct图片，以上10张ct图片作为一个训练样本；按照上述方式从不同的切片层获取若干个训练样本组成训练集；

s2、构造卷积神经网络，所述卷积神经网络包括共享权值的第一通道和第二通道；所述卷积神经网络的损失函数为：

其中m为预设值，d为第一通道输出的第一特征向量与第二通道输出的第二特征向量的欧氏距离；

s3、将所述训练集输入到卷积神经网络中进行正样本训练和负样本训练，得到训练好的卷积神经网络：

正样本训练：从同一个训练样本中选取两张ct图并且分别输入到第一通道和第二通道中，将第一特征向量和第二特征向量输入损失函数中，其中l为1，根据损失函数的值反向传播，更新卷积神经网络的权值；

负样本训练：从不同的两个训练样本中分别选取一张ct图并且分别输入到第一通道和第二通道中，将第一特征向量和第二特征向量输入损失函数中，其中l为0，根据损失函数的值反向传播，更新卷积神经网络的权值。

作为优选，还包括步骤：将dicom格式的ct图片转换成维度为227*227*1并且格式为jpg的ct图片后再作为训练样本。

作为优选，所述第一通道包括输入层、第一卷积层、第一relu层、第二卷积层、第二relu层、第三卷积层、第三relu层、第四卷积层、第四relu层、第五卷积层、第五relu层、第一全连接层、第二全连接层和输出层；所述第二通道与第一通道的结构相同；

所述步骤s3中，将ct图片输入第一通道和第二通道后进行以下步骤：

s301、在第一卷积层中，维度为227*227*1的ct图片通过96个173*173*1的卷积核进行卷积，步数为1，得到55*55*96的三维特征图；

s302、第一relu层中，通过函数relu对三维特征图进行非线性映射；

s303、在第二卷积层中，维度为55*55*96的三维特征图通过256个29*29*96的卷积核进行卷积，步数为1，得到27*27*256的三维特征图；

s304、第二relu层中，通过函数relu对三维特征图进行非线性映射；

s305、在第三卷积层中，维度为27*27*256的三维特征图通过384个15*15*256的卷积核进行卷积，步数为1，得到13*13*384的三维特征图；

s306、第三relu层中，通过函数relu对三维特征图进行非线性映射；

s307、在第四卷积层中，维度为13*13*384的三维特征图通过384个3*3*384的卷积核进行卷积，步数为1，补零后得到13*13*384的三维特征图；

s308、第四relu层中，通过函数relu对三维特征图进行非线性映射；

s309、在第五卷积层中，维度为13*13*384的三维特征图通过256个3*3*384的卷积核进行卷积，步数为1，补零后得到13*13*256的三维特征图；

s310、第五relu层中，通过函数relu对三维特征图进行非线性映射；

s311、在第一全连接层中，将13*13*256的三维特征图转换为4096维的特征向量v1；

s312、在第二全连接层中，将4096维的特征向量v1转换为4096维的特征向量v2；

其中，所述第一通道的特征向量v2为所述第一特征向量，所述第二通道的特征向量v2为所述第二特征向量。

作为优选，还包括步骤:

s4、从第三ct图片集中选取一张ct图片a输入训练好的卷积神经网络的第一通道，得到第一特征向量v’，将第四ct图片集的n张ct图片一张张输入训练好的卷积神经网络的第二通道，得到第二特征向量vi’，其中i＝1,2,3…n；

分别计算第二特征向量vi’与第一特征向量v’之间的欧氏距离d并代入损失函数中计算，令l＝0，最小的损失函数值所对应的ct图片为第四图片集中与ct图片a匹配度最高的ct图片。

本发明的有益效果是，这种ct图片的配准方法利用双通道的卷积神经网络对不同时期的ct图片进行配准，匹配效率和准确性较高。

附图说明

下面结合附图和实施例对本发明进一步说明。

图1是本发明的一种ct图片的配准方法的最优实施例的流程图。

具体实施方式

下面详细描述本发明的实施例，所述实施例的示例在附图中示出，其中自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。下面通过参考附图描述的实施例是示例性的，仅用于解释本发明，而不能理解为对本发明的限制。

在本发明的描述中，需要理解的是，术语“中心”、“纵向”、“横向”、“长度”、“宽度”、“厚度”、“上”、“下”、“前”、“后”、“左”、“右”、“竖直”、“水平”、“顶”、“底”“内”、“外”、“轴向”、“径向”、“周向”等指示的方位或位置关系为基于附图所示的方位或位置关系，仅是为了便于描述本发明和简化描述，而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作，因此不能理解为对本发明的限制。

此外，术语“第一”、“第二”等仅用于描述目的，而不能理解为指示或暗示相对重要性。在本发明的描述中，需要说明的是，除非另有明确的规定和限定，术语“相连”、“连接”应做广义理解，例如，可以是固定连接，也可以是可拆卸连接，或一体地连接；可以是机械连接，也可以是电连接；可以是直接相连，也可以通过中间媒介间接相连。对于本领域的普通技术人员而言，可以具体情况理解上述术语在本发明中的具体含义。此外，在本发明的描述中，除非另有说明，“多个”的含义是两个或两个以上。

流程图中或在此以其他方式描述的任何过程或方法描述可以被理解为，表示包括一个或更多个用于实现特定逻辑功能或过程的步骤的可执行指令的代码的模块、片段或部分，并且本发明的优选实施方式的范围包括另外的实现，其中可以不按所示出或讨论的顺序，包括根据所涉及的功能按基本同时的方式或按相反的顺序，来执行功能，这应被本发明的实施例所属技术领域的技术人员所理解。

如图1所示，是一种ct图片的配准方法的实施例，包括以下步骤：

s1、选取不同时期ct扫描的第一ct图片集和第二ct图片集，从第一ct图片集中选取一个切片层的1张ct图片和前后2张ct图片，从第二ct图片集中选取同一切片层的1张ct图片和前后2张ct图片，将上述10张ct图片从dicom格式转换成维度为227*227*1并且格式为jpg的ct图片，然后放入一个文件夹中，作为一个训练样本；按照上述方式从不同的切片层获取若干个训练样本组成训练集。

s2、构造卷积神经网络，卷积神经网络包括共享权值的第一通道和第二通道，第一通道和第二通道的训练结构相同，包括输入层、第一卷积层、第一relu层、第二卷积层、第二relu层、第三卷积层、第三relu层、第四卷积层、第四relu层、第五卷积层、第五relu层、第一全连接层、第二全连接层和输出层；卷积神经网络的损失函数为：

其中，d为第一通道输出的第一特征向量与第二通道输出的第二特征向量的欧氏距离，m为预设值；从第一ct图片集挑选一张ct图片p1，从第二ct图片集挑选中挑选5张与ct图片p1差别最大的ct图片q1～q5，将ct图片p1输入第一通道中，得到第一特征向量vp1，将5张ct图片q1～q5依次输入到第二通道中，得到第二特征向量vq1～vq5，分别计算第二特征向量vq1～vq5与第一特征向量vp1之间的欧氏距离，选取最大的值，重复几次上述步骤并且每次都重新挑选ct图片，m的值大于几个最大的欧氏距离；

s3、将训练集输入到卷积神经网络中进行正样本训练和负样本训练，得到训练好的卷积神经网络：

正样本训练：从同一个训练样本中选取两张ct图并且分别输入到第一通道和第二通道中，将第一特征向量和第二特征向量输入损失函数中，令l为1，根据损失函数的值反向传播，更新卷积神经网络的权值；

负样本训练：从不同的两个训练样本中分别选取一张ct图并且分别输入到第一通道和第二通道中，将第一特征向量和第二特征向量输入损失函数中，其中l为0，根据损失函数的值反向传播，更新卷积神经网络的权值。

将ct图片输入第一通道和第二通道后具体的训练过程如下：

s301、在第一卷积层中，维度为227*227*1的ct图片通过96个173*173*1的卷积核进行卷积，步数为1，得到55*55*96的三维特征图；

s302、第一relu层中，通过函数relu对三维特征图进行非线性映射；

s303、在第二卷积层中，维度为55*55*96的三维特征图通过256个29*29*96的卷积核进行卷积，步数为1，得到27*27*256的三维特征图；

s304、第二relu层中，通过函数relu对三维特征图进行非线性映射；

s305、在第三卷积层中，维度为27*27*256的三维特征图通过384个15*15*256的卷积核进行卷积，步数为1，得到13*13*384的三维特征图；

s306、第三relu层中，通过函数relu对三维特征图进行非线性映射；

s307、在第四卷积层中，维度为13*13*384的三维特征图通过384个3*3*384的卷积核进行卷积，步数为1，补零后得到13*13*384的三维特征图；

s308、第四relu层中，通过函数relu对三维特征图进行非线性映射；

s309、在第五卷积层中，维度为13*13*384的三维特征图通过256个3*3*384的卷积核进行卷积，步数为1，补零后得到13*13*256的三维特征图；

s310、第五relu层中，通过函数relu对三维特征图进行非线性映射；

s311、在第一全连接层中，将13*13*256的三维特征图转换为4096维的特征向量v1；

s312、在第二全连接层中，将4096维的特征向量v1转换为4096维的特征向量v2；

其中，第一通道的特征向量v2为第一特征向量，第二通道的特征向量v2为第二特征向量。

s4、训练好的卷积神经网络可以用来进行图片的配比：收集不同时期ct扫描得到的第三ct图片集和第四ct图片集；从第三ct图片集中选取一张ct图片a输入训练好的卷积神经网络的第一通道，得到第一特征向量v’，将第四ct图片集的n张ct图片一张张输入训练好的卷积神经网络的第二通道，得到第二特征向量vi’，其中i＝1,2,3…n；

分别计算第二特征向量vi’与第一特征向量v’之间的欧氏距离d并代入损失函数中计算，令l＝0，最小的损失函数值所对应的ct图片为第四图片集中与ct图片a匹配度最高的ct图片。

在本说明书的描述中，参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中，对所述术语的示意性表述不一定指的是相同的实施例或示例。而且，描述的具体特征、结构、材料或者特点可以在任何的一个或多个实施例或示例中以合适的方式结合。

以上述依据本发明的理想实施例为启示，通过上述的说明内容，相关工作人员完全可以在不偏离本项发明技术思想的范围内，进行多样的变更以及修改。本项发明的技术性范围并不局限于说明书上的内容，必须要根据权利要求范围来确定其技术性范围。