数字图像测量亚像素边沿处理方法技术

数字图像测量亚像素边沿处理方法，采用梯度模型对图像边沿进行粗处理，以梯度均值作为二值化控制阀值；基于洪水充填原理实现二值化图像最大连通区域搜索；在四个方向上搜索包含各位置点的一段梯度值超过梯度控制阀值的连续点，点数最少的方向即为法矢方向，基于样条插值细分算法计算边沿法矢方向上的梯度峰值；从收集的目标粗边沿中对拟合曲线上所包含的边沿点进行标记；经过多轮迭代计算，如果收集的目标粗边沿点全部被标记，则表明精确亚像素边沿计算完毕。本发明专利技术是基于传统梯度边沿提取、洪水充填算法及样条插值细分算法，处理过程通过粗、精两步图像处理实现。对于目标边沿实行有目的搜索，计算效率高，并且算法过程兼顾噪声抑制的要求，提高了算法的免疫力。

【技术实现步骤摘要】

本专利技术涉及一种图像测量处理方法，具体涉及一种数字图像测量亚像 素边沿处理方法。技术背景随着当前数字化信息技术的发展，数字图像处理技术正变得越来越重 要，并在很多生产生活现场都得到广泛应用。目前数字图像测量技术在光 电检观顿域也占据了非常重要的地位，它最重要的一个内容就是实现图像 目标边沿的精确定位。这点在数字图像测量应用中尤为重要，可以说在数 字图像测量应用中准确快速提取图像边缘是提取图像相关信息的基础和目前的数字图像都是二维平面图像。根据CCD感光芯片的不同，目标 图像被以不同的细分密度加以离散存储，图像上的最小尺寸单位为像素， 存储形式上存储了每个像素位置的图像色彩信息。数字图像处理就是采用 一定的计算方法对这些存储内容进行相关计算以提取图像中相关感兴趣 信息的过程。目前对图像处理采用最基本的处理算法主要是基于灰度或是梯度原理 的计算方法。它们都是通过计算图像的某个特征值(灰度或梯度)，然后 依据给定的特征值控制阀值对图形进行二值化处理。它所依据的主要思想 就是图象中目标物体的边沿与背景有明显的视觉冲击，即视觉上在物体的 边沿位置图像内容变化非常大。目前衡量这种图像内容变化的主要计算模型是灰度方差或是梯度。设图像上一以点""为中心的3X3图像区域c<table>table see original document page 4</column></row><table>每一点的灰度方差计算模型为<formula>formula see original document page 4</formula>i， j， m， n为非负整数，r(i， j)，g(i, j)，b(i, j)为点(i, j)的RGB (红、绿、蓝三原色)分量值，V (i， j)为点(i， j)的灰度方差。 梯度计算模型为<formula>formula see original document page 4</formula>其中，So6e/ =5b&," So崎O力' So崎7,0」So6e/〃力 5b崎7力 So6e/「2,0) So6e/(^，W So6e/p,"i， j， m， n为非负整数，r(m,n)，g(m，n)，b(m，n)为点(m， n)的RGB (红、绿、蓝三原色)分量值，G (i， j)为点(i, j)的梯度， Sobel(i， j)(i， j^，l，2)为Sobel算子对应(i， j)位置的元素值。目前传统的图像处理算法都是基于这两种计算模型。比如Moravec算 子就是一种利用灰度方差提取特征点的算子，此外还有Roberts、 Sobel、Laplace等梯度算子。梯度计算模型还可以通过改变计算算子中的相关元 素值达到不同的图像处理目的。另外还有一种通过图像空间向参数空间转 换实现图像目标提取的Hough变换算法。传统算法在图像处理前与处理后 一般还采用噪声滤波、平滑与锐化等辅助处理算法。这些图像处理算法对于目标图像质量非常好的对象具有较好的处理 效果，但在实际工程应用中，尤其是在实时图像处理中还存在许多的应用难点一、 由于二值化结果受特征值控制阀值影响，图像处理结果人为影响 因素太大，算法对于不同应用场合图像不具备通用性，而且不利于计算机 自动程序实现；二、 处理图像中存在多种噪声，传统算法虽然希望通过改变相关算子 或是辅助滤波算法以达到噪声抑制的目的，但对于不同图像应用效果甚 微，不利于实时图像处理；三、 对目标图像边沿定位不精确，处理结果太粗燥，难以满足精密图 像测量应用需求，同时在许多有精密图像处理要求的应用场合也存在很大 的局限性。
技术实现思路
本专利技术的目的在针对实时图像测量的应用需求，提供了一种效率高， 并且兼顾噪声抑制的要求，提高了算法的免疫力的数字图像测量亚像素边 沿处理方法。为达到上述目的，本专利技术采用的技术方案是1)首先采用梯度模型对图像边沿进行粗处理，以梯度均值作为二值 化控制阀值；2)基于洪水充填原理实现二值化图像最大连通区域搜索； 通过对起始位置向四面八方蔓延方式得到整个连通区域的大小，搜索出边沿最大连通区域，并记录这些边沿数据，而其它小的连通区域被当作噪声加以抑制；3 )在四个方向上搜索包含各位置点的一段梯度值超过梯度控制阀值 的连续点，点数最少的方向即为法矢方向，也就是精确图像边沿定位的计 算方向，图像边沿位置沿边沿的法矢方向梯度值具有单峰性，而梯度峰值 位置就是精确的边沿位置；4) 基于样条插值细分算法计算边沿法矢方向上的梯度峰值 由于物体在边缘处的梯度值最大，则对一段在边沿法矢方向且包含该边缘点的连续像素点的梯度值利用样条拟合得到一条曲线，在拟合得到的 曲线上进行10倍细分，即定位精度达到0.1个像素，比较细分后各点位 置的梯度幅值，取最大值位置即为精确边缘点；5) 从收集的目标粗边沿中对拟合曲线上所包含的边沿点进行标记；6) 经过多轮迭代计算，如果收集的目标粗边沿点全部被标记，则表 明精确亚像素边沿计算完毕。本专利技术是基于传统梯度边沿提取、洪水充填算法及样条插值细分算法， 处理过程通过粗、精两步图像处理实现。对于目标边沿实行有目的搜索， 计算效率高，并且算法过程兼顾噪声抑制的要求，提高了算法的免疫力。附图说明图1是本专利技术的流程图；图2是本专利技术基于最大连通区域搜索的目标边沿提取图，其中图2(a) 是原始图像，图2 (b)是梯度二值化处理结果的图像，图2 (c)是基于最大连通区域的目标边沿图；图3是过点0c。，y。)的图像边沿法矢方向搜索图，其中图3 (a)是边沿曲线切矢与法矢图，图3 (b)是目标边沿图像，图3 (c)是目标边沿图像法矢搜索示意图；图4是样条曲线插值细分示意图，其中横坐标与纵坐标均为像素值。具体实施方式下面结合附图对本专利技术作进一步详细说明。参见图1，本专利技术首先进行第一步采集图像数据，采集结果如图2(a)。对第一步所采集图像进行第二步梯度算子全域搜索。采用梯度算子模 型进行图像边沿粗处理，以梯度均值作为二值化控制阀值，防止二值化处 理结果中目标边沿信息丢失，同时也利于程序自动实现；利用梯度算子全域搜索的结果，对采集图像进行第三步的处理即洪水 填充区域面积搜索。由于实时图像中存在许多因素干扰目标图像，通过梯 度算子全域搜索处理后，图像存在多种随机噪声。这类噪声一般通过传统 的滤波方式难以奏效。根据图像处理的基本原理，在图像测量中一般一次 只关心某一具体的目标边沿，在只包含目标边沿的目标区域内，二值化处 理后的目标边沿应该具有最大的图像面积。通过洪水搜索算法对起始位置 向四面八方蔓延方式得到整个连通区域的大小，搜索出边沿最大连通区 域，并记录这些边沿数据集P^pi,0^i〈^nn为正整数h而其它小的连通区域 被当作噪声加以抑制。从图2 (b)与(c)的对比中可以看出噪声抑制的效果处理结果中 除目标边沿信息外，其余噪声得到完全抑制。设定不同的图像处理区域， 就可以得到所感兴趣的不同目标边沿。接着进行第四^~~选取pi，搜索pi处边缘法矢方向。前三步得到的图像边沿检测结果还只是一个粗边沿， 一般为3至5个像素宽度，这在 精密图像测量中还无法本文档来自技高网...

【技术保护点】
数字图像测量亚像素边沿处理方法，其特征在于：　１）首先采用梯度模型对图像边沿进行粗处理，以梯度均值作为二值化控制阀值；　２）基于洪水充填原理实现二值化图像最大连通区域搜索；　通过对起始位置向四面八方蔓延方式得到整个连通区域的大小，搜索出边沿最大连通区域，并记录这些边沿数据，而其它小的连通区域被当作噪声加以抑制；　３）在四个方向上搜索包含各位置点的一段梯度值超过梯度控制阀值的连续点，点数最少的方向即为法矢方向，也就是精确图像边沿定位的计算方向，图像边沿位置沿边沿的法矢方向梯度值具有单峰性，而梯度峰值位置就是精确的边沿位置；　４）基于样条插值细分算法计算边沿法矢方向上的梯度峰值　由于物体在边缘处的梯度值最大，则对一段在边沿法矢方向且包含该边缘点的连续像素点的梯度值利用样条拟合得到一条曲线，在拟合得到的曲线上进行１０倍细分，即定位精度达到０．１个像素，比较细分后各点位置的梯度幅值，取最大值位置即为精确边缘点；　５）从收集的目标粗边沿中对拟合曲线上所包含的边沿点进行标记；　６）经过多轮迭代计算，如果收集的目标粗边沿点全部被标记，则表明精确亚像素边沿计算完毕。

【技术特征摘要】

【专利技术属性】
技术研发人员：丁建军，蒋庄德，李兵，郭俊杰，张飞，
申请(专利权)人：西安交通大学，
类型：发明
国别省市：87[中国|西安]