一种低对比度图像分割方法
【技术领域】
[0001] 本发明属于图像处理领域，涉及一种低对比度图像的分割方法。
【背景技术】
[0002] 图像分割技术大致分为:基于阈值的图像分割、基于边缘的图像分割、基于区域的 图像分割、基于模型的图像分割、基于人工智能的图像分割。主要方法有阈值法、区域生长 法、活动轮廓模型法、人工神经网络法、模糊聚类法、分水岭法、边缘算子检测法等。图像分 割技术具有广泛的应用价值，可应用到模式识别、机器视觉等各行业各领域中。
[0003] 基于阈值的图像分割算法的关键在于阈值的选择，不同的阈值可能会引起完全不 同的结果。此类方差往往受噪声影响严重，且受灰度不均匀性影响严重。基于边缘的分割方 法主要基于图像灰度级的不连续性，它通过检测不同连续区域之间的边界来实现对图像的 分割，边缘检测算法有一个缺点：不能得到连续的单像素边缘。传统的边缘算子抗噪性低， 且无法得到连续的边缘，受背景灰度不均影响较大。基于区域特性的分割技术不但考虑了 像素的相似性，还考虑了空间上的邻接性，因此可以有效消除孤立噪声的干扰，具有很强的 鲁棒性。而且，无论是合并还是分裂，都能够将分割深入到像素级，因此可以保证较高的分 割精度。基于统计模式分类的分割技术将图像分割过程视为以像素为基元的模式分类过 程。这类方法，对于无法由灰度区分的复杂的纹理图像显得尤为有效。
[0004] 目前，对应用于普通图像的分割算法有大量研究，已经比较成熟，但普通图像和低 对比度图像在细节及灰度信息上具有较大差异，适用于普通图像的图像分割算法不一定适 用于低对比度图像，近年来针对低对比度图像的图像分割算法研究较少而且无法满足现实 需要。
[0005] 针对低对比度图像的特点，王思贤、刘宗义、梅建新在《W u h a η U n i v e r s i t y Journal of Natural Sciences》（2000,46(5) :641-644)的文章"低对比度图像目标分割的 研究"中提出将矩不变自动门限方法应用到低对比度的图像的自动分割中，耗时短，能克服 二值化过程中的断线、粘连等问题，具有一定实用价值，但是在有噪声情况下精度不高。针 对低对比度焊缝缺陷目标的提取，Lashkia V.在《Image&Vision Computing》（2001，19(5): 261-269)中的文章 "Defect detection in X-ray images using fuzzy reasoning" 中提 出利用区域方差和灰度对比度参数并结合模糊推理进行缺陷提取，但此方法无法广泛应用 于其它低对比度图像的分割。针对X射线图像、超声波图像以及红外图像的分析检测中，目 标与目标之间没有显著的灰度变化给区域分割带来的困难，聂守平、王鸣、刘峰在《中国激 光》(2004,31(1) :89-91)中的文章"低对比度图像分割算法研究"中提出利用空间区域方差 和灰度区域方差进行图像增强和灰度分区，能够提取出一些低对比度缺陷目标。这些方法 将区域特性和对比度结合提取缺陷目标，而针对低对比度图像，低对比度不仅存在于图像 中某个特定细节，因而以上方法的适用性不够强，需要一种对低对比度图像有较好适应性 的图像分割方法。
【发明内容】
[0006] 针对低对比度图像分割完整度差、精度低、抗噪性能较差的问题，本发明提供一种 低对比度图像分割方法，利用局部窗口内边缘像素个数信息实现局部窗口大小的自适应， 依次在局部窗口内利用局部窗口阈值对当前像素进行分割，以显著提高低对比度图像分割 的精度和抗噪性。
[0007] 本发明解决其技术问题所采用的技术方案包括以下步骤：
[0008] (1)初始化:设原图像0的灰度级为L，先计算0对应的方差图像V，再用CANNY算子计 算V的边缘，得到边缘图像E，E中边缘像素灰度值为L-1，其余像素灰度值为0。初始化窗口边 长 N=3〇
[0009] (2)当前窗口包含信息判断及当前像素分割。对原图像0进行逐像素计算，使当前 像素位于计算窗口 W的中心，步骤包括：
[0010] 1)计算边缘图像E中与原图像0对应的窗口 W中的灰度为L-1的像素的个数C，计算 边缘像素在当前窗口 W中所占的比例P，P = C/N2;
[0011] 2)当前窗口包含信息判断:若P 2(Ν-2)/Ν2则执行步骤3);若P〈(N-2)/N2则N=N+2， 返回步骤(2);
[0012] 3)当前像素分割:对0中当前窗口 W计算基于区域限制的Otsu阈值T，若当前窗口中 心像素灰度大于T则置其灰度为L-1，否则置其灰度为0。
[0013] 在上述步骤（1)中，计算V时采用逐像素计算法，以及当前像素位于3 X 3窗口中心， 在此3 X 3窗口内计算当前窗口的方差，对当前窗口方差取整，令当前窗口中心像素的灰度 值为此方差值。具体步骤为：
[0014] 1)对原图像中的所有像素，以当前像素为窗口中心，计算0中当前窗口内所有像素 的灰度均值A;
[0015] 2)计算0中当前窗口内所有像素灰度值与A的欧氏距离D1;
[0016] 3)计算当前窗口的方差5^=1￡^/9:并取整，令V的当前像素灰度为此方差值。
[0017] 在上述步骤(1)中，计算E的具体步骤为：
[0018] 1)用高斯滤波器平滑方差图像V;
[0019] 2)用一阶偏导有限差分计算方差图像V的梯度幅值和方向；
[0020] 3)对梯度幅值进行非极大值抑制；
[0021] 4)检测和连接边缘。
[0022] 在上述步骤(2)的第3)步中，计算基于区域限制的Otsu阈值T的具体步骤为：
[0023] 1)计算当前窗口 W内每个灰度级的概率和数量；
[0024] 2)计算阈值分别取各灰度级时的背景概率、目标概率、背景均值、目标均值、背景 方差、目标方差；
[0025] 3)遍历各灰度级，求得使类间方差达到最大的阈值T1;
[0026] 4)在当前窗口 W内统计[0，?\]内各灰度级的概率和数量；
[0027] 5)在当前窗口 W内，计算[0，!^]灰度级范围内阈值分别取各灰度级时的背景概率、 目标概率、背景均值、目标均值、背景方差、目标方差；
[0028] 6)遍历[0，1\]中各灰度级，求得使类间方差达到最大的阈值T。
[0029] 本发明的有益效果是:本发明提供的低对比度切片图像分割方法，可有效、完整地 对低对比度图像进行分割，并将分割误差控制在可接受范围内，抗噪性良好，且本方法无需 设置分割参数、自适应