一种基于改进区域生长的三维断骨分割方法技术

一种基于改进区域生长的三维断骨分割方法包括：首先，读入三维

【技术实现步骤摘要】
一种基于改进区域生长的三维断骨分割方法


：
[0001]本专利技术属于医学图像处理领域，涉及一种对医学
CT
骨骼图像的三维断骨分割方法，特别运用于帮助医生进行病灶分析
、
疾病诊断与治疗等
。

技术介绍

：
[0002]随着骨科手术的发展，对于骨骼分割等医学图像处理技术的需求越来越大，其中包括颌面外科
、
脊柱外科
、
关节置换等领域
。
越来越多的医生使用医学图像处理技术来辅助诊断和治疗，提高手术的安全性和有效性
。
[0003]同时，计算机技术的迅猛发展和普及，为医学图像处理技术的发展提供了更多的可能性
。
在现有的医学图像处理技术中，骨骼分割技术是其中重要的分支之一
。
然而，现有的骨骼分割技术普遍存在一些问题
。
现有的主流的断骨分割方法分为传统方法和深度学习方法，其中传统方法主要有阈值分割
、
区域生长
、
边界跟踪等等，这些方法都存在一定的限制，例如阈值分割容易受到噪声和外部干扰的影响，区域生长方法需要人工指定生长点，生长准则难以确定，导致效果不稳定等；深度学习方法也具有数据标注难度大，训练效率低，计算成本高等问题
。
因此现有方法可能具有成本高
、
效率低
、
准确度不够
、
对于人工交互依赖大等，难以满足医生对于骨骼分割技术的需求
。
[0004]因此，提高骨骼分割的准确性和稳定性，提高分割效率以及降低人工干预的需要就成为了医学图像处理领域的重要课题之一
。
本专利技术基于改进区域生长的三维断骨分割方法，不同于传统方法和深度学习方法，直接从
CT
图像整体入手，在结合前后图像关系进行分割处理，可以提高分割的准确性和稳定性，同时可以并行生成多个断骨，大幅度提高处理效率
。
尤其在算法方面，本专利技术对区域生长做了优化，自动生成种子点，无需人工指定，降低了对人工干预的要求，提高了算法的自主智能性
。
因此本专利技术具有准确性高
、
稳定性强
、
处理效率高
、
智能化等特点，具有很高的应用价值和推广前景，可以为骨科
、
口腔科等多个医疗领域的骨骼分割工作提供更为高效和准确的处理方法
。
从而更好地提高医生的诊断准确性和治疗成功率，促进医学图像处理领域的发展，具有重要的创新性和实用性价值
。

技术实现思路

：
[0005]本专利技术要解决以往的断骨分割任务效率低
、
准确率低
、
计算成本高
、
需要人工介入的缺点，提出一种基于改进区域生长的三维断骨分割方法
。
[0006]具体来说本专利技术提供了一种基于改进区域生长的三维断骨分割方法，包括如下步骤
:
[0007]S1
：输入一个
CT
骨骼图像路径；
[0008]S2
：读取这一组图像，转化为三维数组，数组中的每个值代表该位置像素点的
HU
值
。
[0009]S3
：对整个三维数组按照一定阈值，进行阈值分割，并把阈值分割后的三维数组转化为非0即1的二值化的三维数组；
[0010]S31
：对数组中的每一个位置的值，如果该值大于等于该阈值，则将数组中的该位置的值赋值为1；
[0011]S32
：对数组中的每一个位置的值，如果该值小于该阈值，则将数组中的该位置的值赋值为0；
[0012]S4
：遍历新生成的三维数组，将其通过采用一定的结构化元素的形态学腐蚀操作，来去除断骨之间的粘连，操作结束后该三维数组的值也是一个非0即1的二值化的三维数组；
[0013]S5
：通过
S2
中读入的三维数组和
S4
中新生成的三维数组，来得到预处理后的有效三维数组；
[0014]S51
：遍历
S2
中的原始三维数组，如果在
S4
中新生成的三维数组里对应位置的值为1，则有效三维数组对应位置保留为
S2
位置的值；
[0015]S52
：遍历
S2
中的原始三维数组，如果在
S4
中新生成的三维数组里对应位置的值为0，则有效三维数组对应位置的值设置为一个极小值，用来代表该位置不包含骨骼；
[0016]S6
：通过
S5
中生成的有效三维数组，进行区域生长运算，最终生长得到多个连通区域，即得出了多个断骨的对应区域；
[0017]S61
：对数组染色，即声明一个颜色数组
colors
，维度和大小和
S2
中读入的三维数组大小相同，对
S4
这个二值化数组使用并查集算法得到多个连通区域，并对每个连通区域中的所有位置在颜色数组
colors
中依次编号为
1,2,3
…
；
[0018]S62
：将
S4
这个数组中对应的位置加入队列中，并按照一定的阈值规则同时对三维空间的6个方向进行区域生长，某种子点要和由他生长出来的新点的颜色编号相同，如果某生长位置已经被其他种子点生长
(
即已经被染色
)
，则不可生长该位置；
[0019]S63
：新生长出来的点加入队列，重复
S62
操作，直到队列为空；
[0020]S64
：最终该染色数组
colors
中，相同颜色编号的位置即代表是同一块断骨，不同颜色编号的位置即代表是不同断骨；
[0021]S7
：得到每个断骨包含的像素点位置，以及每个像素点位置对应的断骨编号
。
[0022]优选地，所述步骤
S3
中的一定阈值为
300HU。
[0023]优选地，所述步骤
S4
中结构化元素的大小为
[2,2]。
[0024]优选地，所述步骤
S52
中，极小值的代表参数为
‑
10000。
[0025]优选地，所述步骤
S62
中，一定的阈值规则为相邻元素的差值的小于
50。
[0026]本专利技术结合阈值分割与区域生长方法，作用于一组三维
CT
图像，用于提取三维图像的骨骼区域，在骨骼识别后，对不同断骨区域做出区分，最终识别出分离的多块分离的断骨，实现了一种效率高，错误率低的三维断骨分割方法
。
首先，读入三维
CT...

【技术保护点】

【技术特征摘要】
1.
一种三维改进区域生长的断骨分割方法，其特征在于，包括以下步骤
:S1
：输入一个三维
CT
骨骼图像路径；
S2
：读取这一组图像，转化为三维数组，数组中的每个值代表该位置像素点的
HU
值
。S3
：对整个三维数组按照一定阈值，进行阈值分割，并把阈值分割后的三维数组转化为非0即1的二值化的三维数组；
S31
：对数组中的每一个位置的值，如果该值大于等于该阈值，则将数组中的该位置的值赋值为1；
S32
：对数组中的每一个位置的值，如果该值小于该阈值，则将数组中的该位置的值赋值为0；
S4
：遍历新生成的三维数组，将其通过采用一定的结构化元素的形态学腐蚀操作，来去除断骨之间的粘连，操作结束后该三维数组的值也是一个非0即1的二值化的三维数组；
S5
：通过步骤
S2
中读入的三维数组和步骤
S4
中新生成的三维数组，来得到预处理后的有效三维数组；
S51
：遍历步骤
S2
中的原始三维数组，如果在步骤
S4
中新生成的三维数组里对应位置的值为1，则有效三维数组对应位置保留为步骤
S2
位置的值；
S52
：遍历步骤
S2
中的原始三维数组，如果在步骤
S4
中新生成的三维数组里对应位置的值为0，则有效三维数组对应位置的值设置为一个极小值，用来代表该位置不包含骨骼；
S6
：通过步骤
S5
中生成的有效三维数组，进行区域生长运算，最终生长得到多个连通区域，即得出了多个断骨的对应区域；
S61
：对数组染色，即声明一个颜色数组
colors
...

【专利技术属性】
技术研发人员：卜佳俊，赵心治，周晟，顾静军，于智，史可越，杨俊尧，
申请(专利权)人：浙江大学，
类型：发明
国别省市：