一种骨水泥植入体属性定量分析方法技术

本发明专利技术实施例提供一种骨水泥植入体属性定量分析方法，涉及骨水泥分布对椎体成形术疗效影响研究技术领域。具体地，该方法从图像处理角度定量分析植入人体内的骨水泥不光滑分布状态，即统计“毛刺”数量，实现步骤为，首先将三维骨水泥模型用不同方向上的二维切片图像集合来表征，分别提取数据集在x，y，z轴方向上的切片数据，遍历每一方向上的切片数据，删除面积小于设定阈值的连通域后保存到新的数据集，然后对新生成的数据集和原输入数据集做差值处理，得到“毛刺”数据集，最后在三维坐标下统计“毛刺”数据集中的不连通区域个数，得到“毛刺”个数。所述方法可以为椎体成形术术后疗效分析提供一定理论参考。

【技术实现步骤摘要】
一种骨水泥植入体属性定量分析方法
本专利技术涉及骨水泥植入体对椎体成形术疗效影响研究
，具体而言，涉及一种骨水泥植入体属性定量分析方法。
技术介绍
骨水泥是一种用于填充骨与植入物间隙的一种具有自凝特性的生物材料。骨水泥的量被认为是影响经皮椎体成形术（PVP）和经皮椎体后凸成形术（PKP）疗效的重要因素。尽管诸多学者建议骨水泥在骨折椎体内应尽量对称分布，但目前还没有具体阐述骨水泥的分布状态与PVP或PKP疗效间相互关系的报道。描述骨水泥分布状态的常用属性有“表面积”、“体积”和“毛刺”等，其中，“毛刺”是对骨水泥表面不光滑状态的描述，即表面凸起或下凹。一般来说，医生希望植入的骨水泥表面有一定“毛刺”而非一个光滑体，因为一定数量的“毛刺”可以保证骨水泥与人体组织的紧密粘合，避免发生松动或者断裂，导致严重后果。在目前实际的临床诊断中，医生主要通过观察病人患处的多张CT断层图像并结合自己的专业知识进行主观判断，给出相应的视觉类似评分（VAS），缺乏客观定量分析。为了更好地从医学角度上去分析PVP、PKP疗效与骨水泥分布状态的关系，骨水泥属性的定量分析方法需要提出，其中，“表面积”、“体积”的计算可依据相应的计算公式，但“毛刺”数量的统计需要特定的算法。
技术实现思路
有鉴于此，本专利技术实施例从图像处理角度出发，提出一种“毛刺”数量的定量统计方法，以解决上述问题。本专利技术实施例提供的一种骨水泥植入体“毛刺”数量统计方法主要思想是将三维骨水泥模型用不同方向上的二维切片图像集合来表征。对骨水泥表面的任意一个凸起“毛刺”而言，一定在某一方向上的二维切片图像集中存在一张或连续几张“毛刺”顶端部分表现为小面积连通域的图片。所有满足上述条件的连通域集合即为表征“毛刺”的数据集。本专利技术实施例提供的一种骨水泥植入体“毛刺”数量统计方法包括：输入经过去噪、图像分割预处理后的骨水泥序列CT数据集；分别提取数据集在x，y，z轴方向上的切片数据，遍历每一方向上的切片数据，删除面积小于设定阈值的连通域后保存到新的数据集；然后对新生成的数据集和原输入数据集做差值处理，得到“毛刺”数据集；在三维坐标下统计“毛刺”数据集中的不连通区域个数，得到毛刺个数。与现有技术相比，本专利技术实施例提供的一种骨水泥“毛刺”数量统计方法是一种骨水泥植入体属性定量分析方法，可以从图像处理角度出发，定量地分析骨水泥“毛刺”数量与术后疗效问题。为使本专利技术的上述目的、特征和优点能更明显易懂，下文特举较佳实施例，并配合所附附图，作详细说明如下。附图说明为了更清楚地说明本专利技术实施例的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，应当理解，以下附图仅示出了本专利技术的某些实施例，因此不应被看作是对范围的限定，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他相关的附图。图1是本专利技术实施例中所使用的数据源骨水泥序列CT扫描图像中的任意一张。图2是本专利技术实施例中经过图像预处理后的骨水泥序列CT扫描图像中任意选取的一张切片图像。图3是本专利技术实施例提供的一种骨水泥“毛刺”数量统计方法的流程图。图4是本专利技术实施例中经过图像预处理后的骨水泥序列CT扫描图像数据集在MATLABR2016b平台下的三维视图，即植入体骨水泥的三维显示。图5是本专利技术实施例中骨水泥序列CT扫描图像数据集中z轴方向上任取一张切片图像与去除小连通域后的对比图。图6是本专利技术实施例中骨水泥删除表面毛刺顶端部分在MATLABR2016b平台下的三维显示。图7是本专利技术实施例中骨水泥的毛刺顶端部分在MATLABR2016b平台下三维显示。图8是本专利技术实施例中骨水泥三维视图上标记“毛刺”后的显示结果。具体实施方式下面将结合本专利技术实施例中附图，对本专利技术实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本专利技术一部分实施例，而不是全部的实施例。通常在此处附图中描述和展示的本专利技术实施例的显示结果可以以各种不同CT数据源来分析和设计。因此，以下对在附图中提供的本专利技术的实施例的详细描述并非旨在限制要求保护的本专利技术的范围，而是仅仅表示本专利技术的选定实施例。基于本专利技术的实施例，本领域技术人员在没有做出创造性劳动的前提下所获得的所有其他实施例，都属于本专利技术保护的范围。应注意到：相似的标号和字母在下面的附图中表示类似项，因此，一旦某一项在一个附图中被定义，则在随后的附图中不需要对其进行进一步定义和解释。如图1所示，是本专利技术实施例中所使用的数据源骨水泥序列CT扫描图像中的任意一张。本专利技术实施例所使用的数据源为骨水泥序列CT扫描图像，共29张，均来自四川省科学城医院。数据格式为DICOM,像素大小为484*484。图1所示的骨水泥CT扫描图像中白色区域部分为骨水泥。从图1可以看出本专利技术实施例所使用的骨水泥CT扫描图像中包含了一定的背景和噪声，需要经过去噪、图像分割预处理才可作为后期的骨水泥“毛刺”统计分析数据集。如图2所示，是本专利技术实施例中经过图像预处理后的骨水泥序列CT扫描图像中任意选取的一张切片图像。图像预处理，是指对初始数据源，即人体脊椎植入骨水泥部位的CT扫描图像进行预处理，经过阈值分割、中值滤波、掩盖处理去掉背景和大量噪声，生成处理后的二值图像。从图2可以看出本专利技术实施例中经过图像预处理后的骨水泥CT扫描图像是二值图像，黑色为背景，白色为骨水泥。如图3所示，是本专利技术实施例提供的一种骨水泥“毛刺”数量统计方法的流程图，请参阅图3，所述方法包括：步骤S110，输入经过去噪、图像分割预处理后的骨水泥序列CT数据集D。为了给出直观的效果，将数据集D以三维视图表示，如图4所示为数据集D的三维视图，即植入体骨水泥的三维显示。步骤S120，分别提取数据集D在x，y，z轴方向上的切片数据，在切片上删除面积小于阈值的连通域后保存到新的数据集P中。小连通域是指面积小于设定阈值的连通域。以z轴方向上提取的一张切片数据为例，如图5所示为数据集D中z轴方向上任取一张切片图像与去除小连通域后的对比图，左图部分为没有去除小连通域的切片图像，右图部分为去除小连通域的切片图像。遍历x，y，z轴上所有切片图像，将去除小连通域后的数据存放到新的数据集P中。为了给出直观的效果，将数据集P以三维视图表示，如图6所示为数据集D经过删除处理后得到的数据集P的三维视图，即骨水泥删除表面毛刺顶端部分的三维显示。步骤S130，将数据集D与数据集P做差值处理，得到“毛刺”数据集burr。为了给出直观的效果，将差值后的“毛刺”数据集burr以三维视图表示，如图7所示为数据集D与数据集P的三维视图，即骨水泥的毛刺顶端部分三维显示。步骤S140，在三维坐标下统计数据集burr中的不连通区域个数，得到毛刺个数。在三维坐标下统计数据集burr中的不连通区域个数的方法是，先取burr中z方向第一层数据，计算其连通域数目，作为毛刺初始数量；然后取出下一层数据，与上一层数据做相与操作，计算相与后图像的连通域数目，记作n；再统计该层的连通域数量记作p，那么该层单独的毛刺数量则为p–n。遍历数据集burr中z方向所有层数据，计算出“毛刺”数量。为了给出直观的效果，在骨水泥三维视图上标记“毛刺”并显示，如图8所示为骨水泥三维视图上标记“毛刺”后的显示结果。本本文档来自技高网...

【技术保护点】
1.一种骨水泥植入体属性定量分析方法,其特征在于,所述方法包括:从图像处理角度定量分析植入人体内的骨水泥不光滑分布状态，即统计“毛刺”数量，其中“毛刺”是对骨水泥表面不光滑状态的描述，即表面凸起或下凹。

【技术特征摘要】
1.一种骨水泥植入体属性定量分析方法,其特征在于,所述方法包括:从图像处理角度定量分析植入人体内的骨水泥不光滑分布状态，即统计“毛刺”数量，其中“毛刺”是对骨水泥表面不光滑状态的描述，即表面凸起或下凹。2.根据权利要求1中所述的方法，其特征在于，“毛刺”数量的统计，“毛刺”数量统计方法是将三维骨水泥模型用不同方向上的二维切片图像集合来表征，对骨水泥表面的任意一个“毛刺”而言，一定在某一方向上的二维切片图像集中存在一张或连续几张毛刺顶端部分表现为小面积连通域的图片，所有满足上述条件的连通域集合即为表征“毛刺”的数据集。3.根据权利要求2中所述的方法，其特征在于，将三维骨水泥模型用不同方向上的二维切片图像集合来表征，统计其“毛刺”数量的方法是首先输入经过去噪、图像分割预处理后的骨水泥CT数据集，分别提取数据集在x，y，z轴方向上的切片数据，遍历每一方向上的切片数据，删除面积小于设定阈值的连通域后保存到新的数据集；然后对新生成的数据集和输入的数据集做差值计算，得到“毛刺”数据集；在三维坐标下统计“毛刺”数据集中的不连通区域个数，得到毛...

【专利技术属性】
技术研发人员：方艳红，赵琳，祝国军，
申请(专利权)人：西南科技大学，
类型：发明
国别省市：四川,51