[发明专利]一种局部参数自适应的工业CT图像分割方法

权利要求书

1.一种局部参数自适应的工业CT图像分割方法，其特征在于包括下述步骤：

步骤(1)：基于遗传算法的局部窗口内参数自适应优化；设原图像O的长为H，宽为W，灰度级为L，i表示行，j表示列，局部窗口边长为15，当前局部窗口分割结果为R，与局部分割结果中目标对应的原图像局部窗口中的像素方差为σ_obj，与分割结果中背景对应的原图像中的像素方差为σ_bac，σ_obj和σ_bac的差的绝对值为σ_dif，步骤包括：

1)令i＝1,j＝1,σ_bac＝0,σ_obj＝0；

2)将原图像O中第i行、第j列像素设为当前局部窗口中心；

3)利用遗传算法求得使当前局部窗口内σ_dif达到最小值时的Phansalkar算法参数；

其中，计算局部窗口内σ_dif的具体步骤包括：

1)找到原图像中分别与局部窗口中背景对应的像素集合{S_bac}；

2)计算集合{S_bac}中像素的方差σ_bac；

3)找到原图像中分别与局部窗口中目标对应的像素集合{S_obj}；

4)计算集合{S_obj}中像素的方差σ_obj；

5)计算σ_dif＝|σ_bac-σ_obj|；

步骤(2)：局部窗口分割；对当前中心像素O_ij所在局部窗口内的所有像素，将步骤(1)自适应优化所得参数作为Phansalkar算法参数并计算阈值T，利用T对当前像素进行分割；

其中，窗口内用Phansalkar算法计算阈值T的具体步骤包括：

1)计算当前窗口内所有像素的灰度均值μ；

2)计算当前窗口内所有像素的灰度方差σ_local；

3)计算阈值T＝μ×(1+p×e^-q×μ)+k×((σ_lical/r)-1))；

其中，用阈值T对当前像素进行分割的步骤包括：

1)若O_ij灰度大于T，将分割结果中对应像素R_ij的灰度置为L-1，否则置为0；

步骤(3)：结束判断；若i＝H且j＝W则结束；若jW则j＝j+1，转到步骤(1)；若iH且j＝W则i＝i+1,j＝1，转到步骤(1)。