主权项 |
1、一种子像素级阶跃边缘快速提取方法,其特征在于,阶跃边缘轮廓模型定义为:<math> <mrow> <mi>y</mi> <mo>=</mo> <mfrac> <mi>a</mi> <mrow> <mn>1</mn> <mo>+</mo> <msup> <mi>e</mi> <mfrac> <mrow> <mi>x</mi> <mo>-</mo> <mi>b</mi> </mrow> <mi>σ</mi> </mfrac> </msup> </mrow> </mfrac> <mo>+</mo> <mi>c</mi> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mo>[</mo> <mn>1</mn> <mo>]</mo> </mrow> </math> 提取子像素级阶跃边缘的步骤如下:(1)使用Gauss卷积核对原图像进行平滑滤波处理;(2)从上到下、从左到右按行/列扫描图像,若当前扫描点的灰度值比此前扫描点的灰度值的下降幅度>30,则当前点即为阶跃边缘的稳态灰度最大值边界点Ps;然后在Ps点之后继续扫描,若当前扫描点的灰度值比此前扫描点的灰度值的变化幅度<5,则当前点即为阶跃边缘的稳态灰度最小值的边界点Pe,灰度值的取值范围为0~255,于是有Ps点的灰度值ys和横坐标xs,Pe点的灰度值ye和横坐标xe满足下式:ys=ymax,xs=xmin ye=ymin,xe=xmax [5](3)计算边缘模型参数的初始值a0,b0,c0;按照下式计算得到a、b和c的初始值分别为:a0=ymax-ymin=ys-ye [6]b0=(xmax +xmin)/2=(xs+xe)/2 [7]c0=ymin=ye [8](4)计算a,b,c的最优值,确定阶跃边缘,其步骤为:①在点Ps和Pe所在的行或列之间按行或列搜索灰度值ye≤yi≤ys的点Pi(xi,yi),i=1,2,...,m,其中m为搜索到的点的个数,m≥3;②由式[1],令<math> <mrow> <mi>f</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>)</mo> </mrow> <mo>=</mo> <mi>y</mi> <mo>-</mo> <mfrac> <mi>a</mi> <mrow> <mn>1</mn> <mo>+</mo> <msup> <mi>e</mi> <mfrac> <mrow> <mi>x</mi> <mo>-</mo> <mi>b</mi> </mrow> <mi>σ</mi> </mfrac> </msup> </mrow> </mfrac> <mo>+</mo> <mi>c</mi> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mo>[</mo> <mn>9</mn> <mo>]</mo> </mrow> </math> 取目标优化函数为下式:<math> <mrow> <mi>min</mi> <munderover> <mi>Σ</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>m</mi> </munderover> <msup> <mi>f</mi> <mn>2</mn> </msup> <mrow> <mo>(</mo> <msub> <mi>x</mi> <mi>i</mi> </msub> <mo>,</mo> <msub> <mi>y</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mo>[</mo> <mn>10</mn> <mo>]</mo> </mrow> </math> ③把点Pi(xi,yi),i=1,2,...,m带入式[10],得到一个超定的关于a,b,c的非线性方程组。利用Levenberg-Marquardt非线性优化算法求出a, b,c的最优值;则点x=b,<math> <mrow> <mi>y</mi> <mo>=</mo> <mfrac> <mi>a</mi> <mn>2</mn> </mfrac> <mo>+</mo> <mi>c</mi> </mrow> </math> 即为阶跃边缘点的位置。 |