本发明涉及到邮政分拣设备系统
技术领域：
，具体地说涉及一种邮政分拣系统所获取图像的处理方法。
背景技术：
：邮政分拣是国家公共安全的重点行业，百姓有意无意邮寄违禁品包裹投入到邮寄流程将对本人及他人生命和国家财产构成一种潜在的威胁，极易造成严重事故。更值得注意的是一些犯罪分子和恐怖分子为达到破坏社会安定的目的，有意邮寄违禁品来“考核”国家公共安全的建设水平更是一种严重的威胁，发生在国内外的恐怖事件已证明了这一点。随着国际反恐形势的不断深化，恐怖分子的破坏手段变得越来越高科技、越来越复杂化，使得国家职能机构应对不暇。采用高科技手段来提高安全检测的准确性、查堵破坏手段的复杂性变得越来越重要。多年来广泛使用的安全检查设备，如磁力针、金属武器检测门、X射线检测仪等，能发现武器和普通炸药等危险品，在安全检查工作中发挥了重要作用。但是，这些设备受原技术条件所限，不尽如人意，漏检失误较多。再加之当今随着科学技术的发展，犯罪分子和恐怖分子也利用高新技术，制造新的武器、爆炸物等。如利用集成电路技术制造的高精密炸弹、塑性炸药及毒品等，上述传统检测手段就无能为力。技术实现要素：本发明的目的在于针对上述现有技术的缺陷，提供一种邮政分拣系统所获取图像的处理方法；采用图像特征的提取、图像增强、图像复原、图像识别及图像编码等处理方法，对邮政分拣系统所获取的图像进行有效处理。为了实现上述目的，本发明的技术方案是：一种邮政分拣系统所获取图像的处理方法，该图像处理方法的步骤是：S1、将邮政分拣系统所获取图像采用正交变换的方法，得到变换域系数阵列；A1、用灰度直方图对图像进行处理，选择图像某一灰度级D，然后定义一条轮廓线，该轮廓线连接了图像上所有具有灰度级D的点；所得到的轮廓线形成了包围灰度级大于等于D的区域的封闭曲线；将一幅连续图像中被具有灰度级D的所有轮廓线包围的面积称为该图像的阈值面积函数，直方图定义为：H(D)=limΔD→0A(D)-A(D+ΔD)ΔD=-ddDA(D)---(1)]]>将图像看成是一个二维的随机变量，这面积函数相当于其累积分布函数，而灰度直方图相当于其概率密度函数；对于离散函数，我们固定ΔD为1，则上式变为：H(D)＝A(D)-A(D+1)(2)A2、灰度的线性变换将图像中所有的点的灰度按照线性灰度函数进行变换；该线性灰度变换函数f(x)是一个一维线性函数：f(x)＝fA·x+fB(3)灰度变换方程为：DB＝f(DA)＝fA·DA+fB(5)式中参数fA为线性函数的斜率，fB是线性函数在y轴的截距，DA表示输入图像的灰度，DB表示输出图像的灰度；A3、灰度的阈值变换采用灰度的阈值变换将一幅灰度图像转换成黑白二值图像；灰度的阈值变换的变换函数表达式：f(x)=0x<T255x≥T---(4)]]>T是指定的阈值；A4、采用灰度拉伸的方法控制输出灰度参数的分布，有选择地拉伸某段灰度区间以改善输出图像；灰度拉伸是分段进行线性变换的函数，其表达式为：f(x)=y1x1xx<x1y2-y1x2-x1(x-x1)+y1x1≤x≤x2255-y2255-x2(x-x2)+y2x>x2---(5)]]>A5、灰度均衡，通过点运算使输入图像转换为在每一灰度级上都有相同的像素点数的输入图像(即输出的灰度参数是平的)；按照图像的概率密度函数(归一化到单位面积的直方图)的定义：p(x)=1A0H(x)---(6)]]>其中H(x)为直方图，A0为图像的面积；设转换前图像的概率密度函数为pr(r)，转换后的图像的概率密度函数为ps(s)，转换函数为s＝f(r)；由概率知识我们可以得到：ps(s)=pr(r)·drds---(7)]]>这样，如果想使转换后图像的概率密度函数为1(即直方图为平的)，则必须满足：pr(r)=dsdr---(8)]]>等式两边对r积分，可得：s=f(r)=∫0rpr(μ)dμ=1A0∫0rH(μ)dμ---(9)]]>该转换公式被称为图像的累积分布函数；上面的公式是被归一化后推导的，对于没有归一化的情况，只要乘以最大灰度值(Dmax，对于灰度图就是255)即可；灰度均衡的转换公式为：DB=f(DA)=DmaxA0∫0DAH(μ)dμ---(10)]]>对于离散图像，转换公式为：DB=f(DA)=DmaxA0Σi=0DAHi---(11)]]>式中Hi为第i级灰度的像素的个数；S2、图像的几何变换处理B1、图像的平移，将图像中所有的点都按照指定的平移量水平垂直移动；如图像中两个任意点x1和x2的关系用矩阵表示如下：x1y11=10tx01ty001x0y01---(12)]]>对该矩阵求逆，可以得到逆变换：即这样，平移后的图像上的每一点都可以在原图像中找到对应点；B2、图像的镜像变换设图像高度为lHeight，宽度为lWidth。原图中(x0，y0)经过水平镜向后坐标将变为(lWidth-x0，y0)，其矩阵表达式为：x1y11=-10lWidth010001x0y01---(14)]]>逆运算矩阵表达式为：即同样，(x0，y0)经过垂直镜向后坐标将变为(x0，lHeight-y0)，其矩阵表达式为：x1y11=1000-1lHeight001x0y01---(16)]]>其逆矩阵表达式为：即B3、图像的转置将图像的x坐标和y坐标互换。此变换将是图像的高度和宽度互换。其变换矩阵表达式：x1y11=010100001x0y01---(18)]]>它的逆矩阵变换表达式为：即B4、图像的缩放假设图像X轴方向缩放比率为fx，Y轴方向缩放比率是fy，那么，原图中点(x0，y0)对应与新图中的点(x1，y1)的转换矩阵为：x1y11=fx000fy0001x0y01---(20)]]>它的逆运算表达式为：x0=x1fxxy=y1fy---(21)]]>如果fx＝fy＝0.5时，图像将被缩小一半；此时，缩小后图像中的(0，0)像素将对应原图中的(0，0)像素；图像中的(0，1)像素对应于原图中的(0，2)像素；(1，0)像素对应于原图中的(2，0)像素，依此类推；在原图基础上，每行隔一个像素取一点，每个以行进行操作；同理，当fx＝fy＝2时，图像将放大2倍；实际是将原图每行中的像素重复取值一遍，然后每行重复一次。依此类推；B5、图像的旋转在旋转前：x0=rcos(α)y0=rsin(α)---(22)]]>在旋转后写成矩阵表达式：x1y11=cos(θ)sin(θ)0-sin(θ)cos(θ)0001x0y01---(23)]]>其逆运算如下：x0y01=cos(θ)-sin(θ)0sin(θ)cos(θ)0001x1y11---(24)]]>上述旋转是绕坐标轴原点(0，0)进行的，如果是绕一个指定点(a，b)旋转，则先要将坐标系平移到该点，再进行旋转，然后平移回新的坐标原点；两种坐标系的坐标变换表达式为：xIIyII1=10-a0-1b001xIyI1---(25)]]>其逆变换矩阵表达式为：xIyI1=10a0-1b001xIIyII1---(26)]]>假设图像未旋转时中心坐标为(a，b)，旋转后中心坐标为(c，d)，(在新的坐标系下，以旋转后新图像左上角为原点)，则旋转后变换矩阵表达式为：x1y11=10c0-1d001x1II′y1II′1=10c0-1d001cos(θ)sin(θ)0-sin(θ)cos(θ)0001x1IIy1II1=10c0-1d001cos(θ)sin(θ)0-sin(θ)cos(θ)000110-a0-1b001x0y01---(27)]]>其逆变换矩阵表达式：x0y01=10a0-1b001cos(θ)-sin(θ)0sin(θ)cos(θ)000110c0-1d001x1y11---(28)]]>即x0y00=cos(θ)sin(θ)-ccos(θ)-dsin(θ)+a-sin(θ)cos(θ)csin(θ)-dcos(θ)+b001---(29)]]>因此{x0=x1cos(θ)+y1sin(θ)-csos(θ)-dsin(θ)+ay0=-x1sin(θ)+y1cos(θ)+csin(θ)-dcos(θ)+b---(30)]]>根据上面的转换公式，我们可以方便的编写出实现图像旋转的函数。首先计算出公式中需要的几个参数：a，b，c，d和旋转后新图像的高、宽度。现在已知图像的原始高度为lHeight，宽度为lWidth，以图像中心为坐标系原点，则原始图像四个角的坐标分别为(-lWidth-12,lHeight-12),(lWidth-12,lHeight-12)]]>(-lWidth-12,-lHeight-12),(-lWidth-12,-lHeight-12)]]>按照旋转公式，在旋转后的新图中，这四个点的坐标(根据上面的公式)为：(fDstX1,fDstY1)=(-lWidth-12cos(θ)+lHeight-12sin(θ),lWidth-12sin(θ)+lHeight-12cos(θ))]]>(fDstX1,fDstY1)=(lWidth-12cos(θ)+lHeight-12sin(θ),-lWidth-12sin(θ)+lHeight-12cos(θ))]]>(fDstX1,fDstY1)=(-lWidth-12cos(θ)-lHeight-12sin(θ),-lWidth-12sin(θ)-lHeight-12cos(θ))]]>(fDstX1,fDstY1)=(-lWidth-12cos(θ)-lHeight-12sin(θ),lWidth-12sin(θ)-lHeight-12cos(θ))]]>则新图像的宽度lNewWidth和高度lNewHeight为：lNewWidth＝max(|fDstX4-fDstx1|，|fDstX3-fDstX2|)lNewHeight＝max(|fDstY4-fDstY1|，|fDstY3-fDstY2|)令因为a=lWidth-12,b=lHeight-12,]]>c=lNewWidth-12,d=lNewHeight-12]]>所以则x0=xlcos(θ)+ylsin(θ)+f1y0=xlsin(θ)+ylcos(θ)+f2---(31)]]>S3、图像的正交变换处理采用沃尔什变换法对图像进行正交变换；S4、图像增强处理空间域的图像增强用下面的数学运算表达式来表示：g(x，y)＝f(x，y)·h(x，y)(36)式中f(x，y)是处理前的图像，g(x，y)表示处理后的图像，h(x，y)表示空间处理函数；在图像的某种变换中(通常是频率域中)对图像的变换值进行某种运算处理，然后变换回空间域；数学描述如下：G(μ，v)＝H(μ，v)·F(μ，v)(33)在式中，表示某种频域的正变换，表示该种频域的反变换。C1、图像的平滑处理采用低通滤波法来减少噪声；C2、图像的锐化处理采用图像消除模糊算法进行图像的锐化处理，使变模糊的图像变得更加清晰；设图像为f(x，y)，f(x，y)在点(x，y)处的梯度矢量G→[f(x,y)]=∂f∂x∂f∂y---(35)]]>我们知道梯度的方向是在函数f(x，y)最大变化率方向上，且其幅值为：G[f(x,y)]=(∂f∂x)2+(∂f∂y)2---(36)]]>由此，我们得出下列结论：梯度的数值就是f(x，y)在其变化率方向上的单位距离所增加的量；对于离散数字图像，上式可以写为：G[f(i,j)]=[f(i,j)-f(i+1,j)]2+[f(i,j)-f(i,j+1)]2---(37)]]>为了方便起见，采用近似计算公式：在图像锐化过程中，图像变化缓慢的地方梯度很小，图像显得很暗，给出一个阈值Δ，如果G[f(x，y)]小于该阈值，则图像保持原灰度值不变；如果大于或等于阈值，则赋值为G[f(x，y)]；用公式表示如下：g(x,y)=G[f(x,y)](G[f(x,y)]≥Δ)f(x,y)(G[f(x,y)]<Δ)---(39)]]>如果以La表示一个固定的灰度值，上式可以写成：g(x,y)=La(G[f(x,y)]≥Δ)f(x,y)(G[f(x,y)]<Δ)---(40)]]>在进行图像增强时，我们可以使图像有一个固定背景的灰度值Lb，以突出边缘灰度的变换，变换公式如下：g(x,y)=G[f(x,y)](G[f(x,y)]≥Δ)Lb(G[f(x,y)]<Δ)---(41)]]>只保留两个灰度值时，公式就变为：g(x,y)=La(G[f(x,y)]≥Δ)Lb(G[f(x,y)]<Δ)---(42)]]>S5、图像边缘检测、提取及轮廓跟踪D1、边缘检测采用Roberts边缘检测算子、Sobel边缘检测算子、Prewitt边缘检测算子、Krisch边缘检测算子、高斯-拉普拉斯算子中的一种方法考察图像的每个像素在某个邻域内灰度的变化，利用边缘临近一阶或二阶方向导数变化规律；D2、轮廓的提取和跟踪对于二值图像，进行轮廓提取，将图像内部点掏空，得到图像的轮廓；先根据“探测准则”找出目标物体轮廓上的像素，再根据这些像素的某些特征用“跟踪准则”找出目标物体上的其他像素，从而对图像进行跟踪。作为对上述方法的改进，所述Roberts边缘检测算子的公式为：g(x,y)={[f(x,y)-f(x+1,y+1)]2+[f(x+1,y)-f(x,y+1)]2}12---(43)]]>其中f(x，y)是具有整数像素坐标的输入图像，平方根运算使该处理类似于在人类视觉系统中发生的过程。作为对上述方法的改进，所述“探测准则”的算法是：首先找到第一个边界像素的“探测准则”是：按照从左到右、从下到上的顺序搜索，找到的第一个黑点一定是最左下方的边界点，记为A；它的右、右上、左、左上四个邻点中至少有一个是边界点，记为B；从B开始找起，按右、右上、上、左上、左、左下、下、右下的顺序查找相邻点中的边界点C；如果C就是A点，则表示已经转了一圈，程序结束；否则从C点继续找，直到找到A为止；如果它的上下左右四个邻点都不是黑点则它既为边界点。作为对上述方法的改进，所述“探测准则”的另一种算法是：首先按照上面所说的“探测准则”找到最下方的边界点；以这个边界点起始，假设已经沿顺时针方向环绕整个图像一圈找到了所有的边界点；由于边界是连续的，所以每一个边界点都可以用这个边界点对前一个边界点所形成的夹角来表示；使用下面的跟踪准则：从第一个边界点开始，定义初始的搜索方向为沿左上方；如果左上方的点是黑点，则为边界点，否则搜索方向顺时针旋转45°；这样一直到找到第一个黑点为止；然后把这黑点作为新的边界点，在当前搜索方向的基础上逆时针旋转90°，继续用同样的方法搜索下一个黑点，直到返回最初的边界点为止。与现有技术相比，本发明具有的优点和积极效果是：本发明的邮政分拣系统所获取图像的处理方法，实现了旅客行包的图形图像处理软件系统。可在双显示器上显示出旅客行包的黑白图像和彩色图像，并包括图形的边缘增强、图形放大、图形局部放大、有机物剔出、无机物剔出、黑白图像反转、伪彩色图像、图像增暗、图像增亮、危险品报警等功能；实现细化分析、区域分析及特殊图像处理等数学模型；提高物质识别的准确性。具体实施方式下面结合具体实施方式对本发明技术作进一步详细说明。本发明的一种邮政分拣系统所获取图像的处理方法，该图像处理方法的步骤是：S1、将邮政分拣系统所获取图像采用正交变换的方法，得到变换域系数阵列；A1、用灰度直方图对图像进行处理，选择图像某一灰度级D，然后定义一条轮廓线，该轮廓线连接了图像上所有具有灰度级D的点；所得到的轮廓线形成了包围灰度级大于等于D的区域的封闭曲线；将一幅连续图像中被具有灰度级D的所有轮廓线包围的面积称为该图像的阈值面积函数，直方图定义为：H(D)=limΔD→0A(D)-A(D+ΔD)ΔD=-ddDA(D)---(1)]]>将图像看成是一个二维的随机变量，这面积函数相当于其累积分布函数，而灰度直方图相当于其概率密度函数；对于离散函数，我们固定ΔD为1，则上式变为：H(D)＝A(D)-A(D+1)(2)A2、灰度的线性变换将图像中所有的点的灰度按照线性灰度函数进行变换；该线性灰度变换函数f(x)是一个一维线性函数：f(x)＝fA·x+fB(3)灰度变换方程为：DB＝f(DA)＝fA·DA+fB(5)式中参数fA为线性函数的斜率，fB是线性函数在y轴的截距，DA表示输入图像的灰度，DB表示输出图像的灰度；A3、灰度的阈值变换采用灰度的阈值变换将一幅灰度图像转换成黑白二值图像；灰度的阈值变换的变换函数表达式：f(x)=0x<T255x≥T---(4)]]>T是指定的阈值；A4、采用灰度拉伸的方法控制输出灰度参数的分布，有选择地拉伸某段灰度区间以改善输出图像；灰度拉伸是分段进行线性变换的函数，其表达式为：f(x)=y1x1xx<x1y2-y1x2-x1(x-x1)+y1x1≤x≤x2255-y2255-x2(x-x2)+y2x>x2---(5)]]>A5、灰度均衡，通过点运算使输入图像转换为在每一灰度级上都有相同的像素点数的输入图像(即输出的灰度参数是平的)；按照图像的概率密度函数(归一化到单位面积的直方图)的定义：p(x)=1A0H(x)---(6)]]>其中H(x)为直方图，A0为图像的面积；设转换前图像的概率密度函数为pr(r)，转换后的图像的概率密度函数为ps(s)，转换函数为s＝f(r)；由概率知识我们可以得到：ps(s)=pr(r)·drds---(7)]]>这样，如果想使转换后图像的概率密度函数为1(即直方图为平的)，则必须满足：pr(r)=dsdr---(8)]]>等式两边对r积分，可得：s=f(r)=∫0rpr(μ)dμ=1A0∫0rH(μ)dμ---(9)]]>该转换公式被称为图像的累积分布函数；上面的公式是被归一化后推导的，对于没有归一化的情况，只要乘以最大灰度值(Dmax，对于灰度图就是255)即可；灰度均衡的转换公式为：DB=f(DA)=DmaxA0∫0DAH(μ)dμ---(10)]]>对于离散图像，转换公式为：DB=f(DA)=DmaxA0Σi=0DAHi---(11)]]>式中Hi为第i级灰度的像素的个数；S2、图像的几何变换处理B1、图像的平移，将图像中所有的点都按照指定的平移量水平垂直移动；如图像中两个任意点x1和x2的关系用矩阵表示如下：x1y11=10tx01ty001x0y01---(12)]]>对该矩阵求逆，可以得到逆变换：即这样，平移后的图像上的每一点都可以在原图像中找到对应点；B2、图像的镜像变换设图像高度为lHeight，宽度为lWidth。原图中(x0，y0)经过水平镜向后坐标将变为(lWidth-x0，y0)，其矩阵表达式为：x1y11=-10lWidth010001x0y01---(14)]]>逆运算矩阵表达式为：即同样，(x0，y0)经过垂直镜向后坐标将变为(x0，lHeight-y0)，其矩阵表达式为：x1y11=1000-1lHeight001x0y01---(16)]]>其逆矩阵表达式为：即B3、图像的转置将图像的x坐标和y坐标互换。此变换将是图像的高度和宽度互换。其变换矩阵表达式：x1y11=010100001x0y01---(18)]]>它的逆矩阵变换表达式为：即B4、图像的缩放假设图像X轴方向缩放比率为fx，Y轴方向缩放比率是fy，那么，原图中点(x0，y0)对应与新图中的点(x1，y1)的转换矩阵为：x1y11=fx000fy0001x0y01---(20)]]>它的逆运算表达式为：x0=x1fxxy=y1fy---(21)]]>如果fx＝fy＝0.5时，图像将被缩小一半；此时，缩小后图像中的(0，0)像素将对应原图中的(0，0)像素；图像中的(0，1)像素对应于原图中的(0，2)像素；(1，0)像素对应于原图中的(2，0)像素，依此类推；在原图基础上，每行隔一个像素取一点，每个以行进行操作；同理，当fx＝fy＝2时，图像将放大2倍；实际是将原图每行中的像素重复取值一遍，然后每行重复一次。依此类推；B5、图像的旋转在旋转前：x0=rcos(α)y0=rsin(α)---(22)]]>在旋转后写成矩阵表达式：x1y11=cos(θ)sin(θ)0-sin(θ)cos(θ)0001x0y01---(23)]]>其逆运算如下：x0y01=cos(θ)-sin(θ)0sin(θ)cos(θ)0001x1y11---(24)]]>上述旋转是绕坐标轴原点(0，0)进行的，如果是绕一个指定点(a，b)旋转，则先要将坐标系平移到该点，再进行旋转，然后平移回新的坐标原点；两种坐标系的坐标变换表达式为：xIIyII1=10-a0-1b001xIyI1---(25)]]>其逆变换矩阵表达式为：xIyI1=10a0-1b001xIIyII1---(26)]]>假设图像未旋转时中心坐标为(a，b)，旋转后中心坐标为(c，d)，(在新的坐标系下，以旋转后新图像左上角为原点)，则旋转后变换矩阵表达式为：x1y11=10c0-1d001x1II′y1II′1=10c0-1d001cos(θ)sin(θ)0-sin(θ)cos(θ)0001x1IIy1II1=10c0-1d001cos(θ)sin(θ)0-sin(θ)cos(θ)000110-a0-1b001x0y01---(27)]]>其逆变换矩阵表达式：x0y01=10a0-1b001cos(θ)-sin(θ)0sin(θ)cos(θ)000110c0-1d001x1y11---(28)]]>即x0y00=cos(θ)sin(θ)-ccos(θ)-dsin(θ)+a-sin(θ)cos(θ)csin(θ)-dcos(θ)+b001---(29)]]>因此{x0=x1cos(θ)+y1sin(θ)-csos(θ)-dsin(θ)+ay0=-x1sin(θ)+y1cos(θ)+csin(θ)-dcos(θ)+b---(30)]]>根据上面的转换公式，我们可以方便的编写出实现图像旋转的函数。首先计算出公式中需要的几个参数：a，b，c，d和旋转后新图像的高、宽度。现在已知图像的原始高度为lHeight，宽度为lWidth，以图像中心为坐标系原点，则原始图像四个角的坐标分别为(-lWidth-12,lHeight-12),(lWidth-12,lHeight-12)]]>(-lWidth-12,-lHeight-12),(-lWidth-12,-lHeight-12)]]>按照旋转公式，在旋转后的新图中，这四个点的坐标(根据上面的公式)为：(fDstX1,fDstY1)=(-lWidth-12cos(θ)+lHeight-12sin(θ),lWidth-12sin(θ)+lHeight-12cos(θ))]]>(fDstX1,fDstY1)=(lWidth-12cos(θ)+lHeight-12sin(θ),-lWidth-12sin(θ)+lHeight-12cos(θ))]]>(fDstX1,fDstY1)=(-lWidth-12cos(θ)-lHeight-12sin(θ),-lWidth-12sin(θ)-lHeight-12cos(θ))]]>(fDstX1,fDstY1)=(-lWidth-12cos(θ)-lHeight-12sin(θ),lWidth-12sin(θ)-lHeight-12cos(θ))]]>则新图像的宽度lNewWidth和高度lNewHeight为：lNewWidth＝max(|fDstX4-fDstx1|，|fDstX3-fDstX2|)lNewHeight＝max(|fDstY4-fDstY1|，|fDstY3-fDstY2|)令因为a=lWidth-12,b=lHeight-12,]]>c=lNewWidth-12,d=lNewHeight-12]]>所以则x0=xlcos(θ)+ylsin(θ)+f1y0=xlsin(θ)+ylcos(θ)+f2---(31)]]>S3、图像的正交变换处理采用沃尔什变换法对图像进行正交变换；S4、图像增强处理空间域的图像增强用下面的数学运算表达式来表示：g(x，y)＝f(x，y)·h(x，y)(36)式中f(x，y)是处理前的图像，g(x，y)表示处理后的图像，h(x，y)表示空间处理函数；在图像的某种变换中(通常是频率域中)对图像的变换值进行某种运算处理，然后变换回空间域；数学描述如下：G(μ，v)＝H(μ，v)·F(μ，v)(33)在式中，表示某种频域的正变换，表示该种频域的反变换。C1、图像的平滑处理采用低通滤波法来减少噪声；C2、图像的锐化处理采用图像消除模糊算法进行图像的锐化处理，使变模糊的图像变得更加清晰；设图像为f(x，y)，f(x，y)在点(x，y)处的梯度矢量G→[f(x,y)]=∂f∂x∂f∂y---(35)]]>我们知道梯度的方向是在函数f(x，y)最大变化率方向上，且其幅值为：G[f(x,y)]=(∂f∂x)2+(∂f∂y)2---(36)]]>由此，我们得出下列结论：梯度的数值就是f(x，y)在其变化率方向上的单位距离所增加的量；对于离散数字图像，上式可以写为：G[f(i,j)]=[f(i,j)-f(i+1,j)]2+[f(i,j)-f(i,j+1)]2---(37)]]>为了方便起见，采用近似计算公式：在图像锐化过程中，图像变化缓慢的地方梯度很小，图像显得很暗，给出一个阈值Δ，如果G[f(x，y)]小于该阈值，则图像保持原灰度值不变；如果大于或等于阈值，则赋值为G[f(x，y)]；用公式表示如下：g(x,y)=G[f(x,y)](G[f(x,y)]≥Δ)f(x,y)(G[f(x,y)]<Δ)---(39)]]>如果以La表示一个固定的灰度值，上式可以写成：g(x,y)=La(G[f(x,y)]≥Δ)f(x,y)(G[f(x,y)]<Δ)---(40)]]>在进行图像增强时，我们可以使图像有一个固定背景的灰度值Lb，以突出边缘灰度的变换，变换公式如下：g(x,y)=G[f(x,y)](G[f(x,y)]≥Δ)Lb(G[f(x,y)]<Δ)---(41)]]>只保留两个灰度值时，公式就变为：g(x,y)=La(G[f(x,y)]≥Δ)Lb(G[f(x,y)]<Δ)---(42)]]>S5、图像边缘检测、提取及轮廓跟踪D1、边缘检测采用Roberts边缘检测算子、Sobel边缘检测算子、Prewitt边缘检测算子、Krisch边缘检测算子、高斯-拉普拉斯算子中的一种方法考察图像的每个像素在某个邻域内灰度的变化，利用边缘临近一阶或二阶方向导数变化规律；D2、轮廓的提取和跟踪对于二值图像，进行轮廓提取，将图像内部点掏空，得到图像的轮廓；先根据“探测准则”找出目标物体轮廓上的像素，再根据这些像素的某些特征用“跟踪准则”找出目标物体上的其他像素，从而对图像进行跟踪。所述Roberts边缘检测算子的公式为：g(x,y)={[f(x,y)-f(x+1,y+1)]2+[f(x+1,y)-f(x,y+1)]2}12---(43)]]>其中f(x，y)是具有整数像素坐标的输入图像，平方根运算使该处理类似于在人类视觉系统中发生的过程。所述“探测准则”的算法是：首先找到第一个边界像素的“探测准则”是：按照从左到右、从下到上的顺序搜索，找到的第一个黑点一定是最左下方的边界点，记为A；它的右、右上、左、左上四个邻点中至少有一个是边界点，记为B；从B开始找起，按右、右上、上、左上、左、左下、下、右下的顺序查找相邻点中的边界点C；如果C就是A点，则表示已经转了一圈，程序结束；否则从C点继续找，直到找到A为止；如果它的上下左右四个邻点都不是黑点则它既为边界点。所述“探测准则”的另一种算法是：首先按照上面所说的“探测准则”找到最下方的边界点；以这个边界点起始，假设已经沿顺时针方向环绕整个图像一圈找到了所有的边界点；由于边界是连续的，所以每一个边界点都可以用这个边界点对前一个边界点所形成的夹角来表示；使用下面的跟踪准则：从第一个边界点开始，定义初始的搜索方向为沿左上方；如果左上方的点是黑点，则为边界点，否则搜索方向顺时针旋转45°；这样一直到找到第一个黑点为止；然后把这黑点作为新的边界点，在当前搜索方向的基础上逆时针旋转90°，继续用同样的方法搜索下一个黑点，直到返回最初的边界点为止。与现有技术相比，本发明具有的优点和积极效果是：本发明的邮政分拣系统所获取图像的处理方法，实现了旅客行包的图形图像处理软件系统。可在双显示器上显示出旅客行包的黑白图像和彩色图像，并包括图形的边缘增强、图形放大、图形局部放大、有机物剔出、无机物剔出、黑白图像反转、伪彩色图像、图像增暗、图像增亮、危险品报警等功能；实现细化分析、区域分析及特殊图像处理等数学模型；提高物质识别的准确性。以上所述仅为本发明的较佳实施案例而已，并不用以限制本发明，凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。当前第1页1 2 3