[发明专利]一种骨水泥植入体CT图像层间插值方法

摘要

本发明提出一种骨水泥植入体CT图像层间插值方法。该方法基于形态学处理特性，首先对经二值化操作，提取出目标区域的两切片图像做或运算，得到融合图像；之后，比较融合图像与两切片图像的相似度，若与第一张切片更相似，则在融合图像上对第一张切片相较融合切片的差异部分进行一次腐蚀，再在融合图像上对第二张切片相较融合切片的差异部分进行两次腐蚀，得到一张与第一张切片更相似的插值图像；反之，则在融合图像上对第二张切片相较融合切片的差异部分进行一次腐蚀，再在融合图像上对第一张切片相较融合切片的差异部分进行两次腐蚀，得到一张与第二张切片更相似的插值图像。本发明绘制出的插值图像能描述层间变化过程，过渡自然，边界清晰。

主权项

1.一种骨水泥植入体CT图像层间插值方法，其特征在于，对相邻切片的融合图像与原切片图像的差异部分做形态学处理，使插值图像有选择地保留前后层的差异信息，包括获取融合图像、分别计算两相邻切片图像与融合图像的相似度、根据相似度关系输出插值图像三部分，第一部分包括两个步骤：步骤1，取出相邻两层CT切片图像，对其进行二值化操作，提取出目标区域，得到二值化后的相邻切片图像p1和p2；步骤2，将p1与p2两张图像进行相或运算，得到融合图像p3；第二部分包括两个步骤：步骤3，计算步骤2中p3与p1的差值，得到p3与p1的差异信息p3/ p1，统计p3/ p1中差值为0的像素个数，记为n1，得到p3与p1的相似度S1，其计算公式为：S1= n1/N其中，N 是p3的总像素点个数；步骤4，计算步骤2中p3与p2的差值，得到p3与p2的差异信息p3/ p2，统计p3/ p2中差值为0的像素个数，记为n2，得到p3与p2的相似度S2，其计算公式为：S2= n2/N其中，N 是p3的总像素点个数；第三部分包括两个步骤：步骤5，若S1> S2，先在p3上对p1相较p3的差异部分进行一次腐蚀，再在p3上对p2相较p3的差异部分进行两次腐蚀，得到与p1更相似的插值图像；步骤6，若S1< S2，先在p3上对p2相较p3的差异部分进行一次腐蚀，再在p3上对p1相较p3的差异部分进行两次腐蚀，得到与p2更相似的插值图像；步骤7，输出经步骤5或步骤6处理后的融合图像p3，得到相邻两切片的插值图像。