专利名称：基于Radon变换和极谐波变换的不变矩目标识别方法
技术领域：
本发明属于图像处理技术领域，更进一步涉及遥感应用和医学诊断领域中基于 Radon变换和极谐波变换的不变矩目标识别方法。本发明应用于遥感应用中可以提高遥感图像中目标的识别率，在医学诊断领域中可以更准确地识别医学图像中的病变部位，尤其是对于肿瘤的定位和识别效果显著，而且对于含噪图像的目标识别本发明也具有较高的识别率。
背景技术：
在遥感应用和医学诊断领域中，为了提高图像目标的识别率，采用基于图像的形状特征的目标识别方法。目前图像目标识别方法主要是基于几何不变矩作为基本特征来实现图像目标识别。北京航空航天大学的专利申请“一种基于关键点的仿射不变矩的目标识别方法”(公开号CN101464948A，申请日2009. 1. 14)中公开了一种基于关键点的仿射不变矩的目标识别方法。这种目标识别方法的具体步骤包括(1)图像预处理对原始采集图像进行处理，减少噪声等无效信息对后续处理的影响，增强有用信息，提高图像对比度；(2) 主轮廓提取对目标边缘图像进行主轮廓的提取，得到外轮廓的边缘图像，过滤掉内部的轮廓，和无用的轮廓分支；(3)关键帧提取以分割出来的目标灰度图像为基础，首先计算目标的质心，然后以质心为扩展点向周围引伸出多条射线，寻找每条射线方向上最近灰度极值点，将所有灰度极值点当作关键点集合；(4)仿射不变特征提取提取出目标主轮廓的仿射不变特征，计算多阶仿射不变矩向量；( 目标识别通过步骤中提出的特征来识别目标，输出识别结果。虽然这种目标识别方法在抗噪能力和计算复杂度方面有一定的改善效果，但是由于关键帧提取的过程难以克服高斯白噪声的干扰，这会影响到目标识别过程的适应性。同时，在步骤(5)目标识别部分所使用的隐马尔科夫模型会增加计算的复杂度。 综上所述，这些不足会影响到目标识别的稳定性、准确性以及识别效率。
发明内容
本发明针对现有的图像目标识别技术的不足，提出一种基于Radon变换和极谐波变换的不变矩目标识别方法，通过构造Radon复指数不变矩、Radon正余弦不变矩和极地复指数实虚不变矩三个新的不变矩来实现图像的目标识别，以解决现有目标识别技术中所采用的不变矩抗噪能力差、计算复杂度高和数据不稳定等缺点，提高图像目标的识别率。本发明实现的具体步骤如下(1)输入待识别图像在计算机中应用matlab软件读取待识别的彩色图像；(2)图像预处理2a)用图像颜色空间转换方法将待识别的彩色图像转换为灰度图像；2b)用Sobel边缘检测方法提取灰度图像的目标区域，将图像目标区域的像素灰度值设为1，非目标区域的像素灰度值设为0，得到预处理后的图像；
(3) Radon变换对预处理后的图像进行Radon变换，得到变换后图像P (r，θ)；⑷仿射变换4a)使用matlab软件将图像P(r，θ)进行平移、尺度归一化，得到处理后的图像 P' (r, Θ )；4b)使用图像旋转方法对图像P' (r，θ )进行Φ角度的旋转，得到图像P' (r, θ +Φ)；(5)构造不变矩5a)对图像P' (r, θ+Φ)进行极地复指数变换，得到变换后图像，根据变换前图像P' (r, θ+Φ)和变换后图像之间的旋转角度差，通过等价关系转换的方法得到旋转的相位因子，将变换后图像和相位因子相乘得到构造后的Radon复指数不变矩；5b)对图像P' (r, θ+Φ)进行极地正余弦变换，得到变换后图像，根据变换前图像P' (r, θ+Φ)和变换后图像之间的旋转角度差，通过等价关系转换的方法得到旋转的相位因子，结合变换后图像和相位因子得到构造后的Radon正余弦不变矩；5c)对图像P' (r, θ+Φ)进行极地复指数变换，得到变换后图像，根据变换后图像并应用连乘的方法得到构造后的极地复指数实虚不变矩；(6)提取不变特征对在步骤4b)中得到的图像P' (r，θ +Φ)分别提取Radon复指数不变矩、Radon正余弦不变矩和极地复指数实虚不变矩；将这三个不变矩的实部和虚部作为图像P' (r, Θ+Φ)的不变特征；(7)获得特征模型7a)将与待识别图像内容相接近的图像样本集中的图像依次输入支持向量机工具构建初始特征模型；7b)使用支持向量机工具对样本集图像进行训练学习得到特征模型的最优参数；(8)图像目标识别将特征模型的最优参数和步骤(6)得到的图像P' (r，θ +Φ) 的不变特征数据输入支持向量机工具判断图像P' (r, θ+Φ)所属于的图像类别，完成图像目标识别；(9)输出图像目标识别结果。本发明与现有技术相比具有以下优点第一，本发明使用了 Radon变换，能够很好的解决了图像中噪声干扰问题，从而能更好的体现图像的真实性，提高图像目标识别的准确率。第二，本发明由于使用了极谐波变换，很好的克服了信息冗余、核函数计算复杂度高和数据不稳定等问题，从而提高了不变矩的整体性能以及目标识别的适用性和稳定性；第三，本发明通过提取不变矩的实部和虚部作为不变特征，代替了常规的统计特征的方法，能够更好的表现出图像的不同信息，提高了目标识别的效率。
图1为本发明的流程图；图2为本发明仿真示意具体实施方式
下面结合附图1对本发明做进一步的描述步骤1，输入待识别图像在计算机中应用matlab软件读取存储在计算机硬盘空间中待识别的彩色图像。步骤2，图像预处理2a)应用图像颜色空间转换方法，将待识别彩色图像通过如下公式转换为灰度图像Gray = 0. 233R+0. 587G+0. 114B其中，Gray为转换之后的灰度图像，R，G，B分别为待识别彩色图像中像素的红，
绿，蓝三种颜色分量值。2b)用Sobel边缘检测方法提取灰度图像的目标区域，具体步骤是使用Sobel边缘检测方法对步骤2a)中得到的灰度图像进行边缘检测，得到边缘闭合的边缘图像；具体步骤是第一，分别使用水平方向和垂直方向的高斯卷积滤波器对图像进行卷积运算，得到图像的水平方向和垂直方向的偏导数；第二，计算图像的梯度和模值；第三，在图像梯度方向上比较图像模值，选择梯度方向上模值最大的图像像素作为边缘像素，并设置该像素值为1，得到边缘图像。将边缘图像中闭合边缘的外围区域像素灰度值设为0，闭合边缘的内部区域像素灰度值设为1，得到图像目标区域分割后的二值图像；将二值图像与步骤加)中得到的灰度图像进行点乘操作，得到只包含图像目标区域的灰度图像。
步骤3，I adon 变换对步骤2b!3)中得到的灰度图像按照以下公式进行Radon变换，得到Radon变换后的图像函数P(r，θ)
权利要求
1.一种基于Radon变换和极谐波变换的不变矩目标识别方法，包括以下步骤(1)输入待识别图像在计算机中应用matlab软件读取待识别的彩色图像；(2)图像预处理2a)用图像颜色空间转换方法将待识别的彩色图像转换为灰度图像；2b)用Sobel边缘检测方法提取灰度图像的目标区域，将图像目标区域的像素灰度值设为1，非目标区域的像素灰度值设为0，得到预处理后的图像；(3)Radon变换对预处理后的图像进行Radon变换，得到变换后图像P(r，θ)；(4)仿射变换4a)使用matlab软件将图像P(r，θ)进行平移、尺度归一化，得到处理后的图像 P' (r, θ )；4b)使用图像旋转方法对图像P' (r，Θ)进行Φ角度的旋转，得到图像P' (r, θ +Φ)；(5)构造不变矩5a)对图像P' (r, θ+Φ)进行极地复指数变换，得到变换后图像，根据变换前图像 P' (r, θ+Φ)和变换后图像之间的旋转角度差，通过等价关系转换的方法得到旋转的相位因子，将变换后图像和相位因子相乘得到构造后的Radon复指数不变矩；5b)对图像P' (r, θ+Φ)进行极地正余弦变换，得到变换后图像，根据变换前图像 P' (r, θ+Φ)和变换后图像之间的旋转角度差，通过等价关系转换的方法得到旋转的相位因子，结合变换后图像和相位因子得到构造后的Radon正余弦不变矩；5c)对图像P' (r, θ+Φ)进行极地复指数变换，得到变换后图像，根据变换后图像并应用连乘的方法得到构造后的极地复指数实虚不变矩；(6)提取不变特征对步骤4b)中获得的图像P'(r, θ +Φ)分别提取Radon复指数不变矩、Radon正余弦不变矩和极地复指数实虚不变矩；将这三个不变矩的实部和虚部作为图像P' (r，θ+Φ)的不变特征；(7)获得特征模型7a)将与待识别图像内容相接近的图像样本集中的图像依次输入支持向量机工具构建初始特征模型；7b)使用支持向量机工具对样本集图像进行训练学习得到特征模型的最优参数；(8)图像目标识别将特征模型的最优参数和步骤(6)得到的图像P'(r，θ +Φ)的不变特征数据输入支持向量机工具判断图像P' (r, Θ+Φ)所属的图像类别，完成图像目标识别；(9)输出图像目标识别结果。
2.根据权利要求1所述的基于Radon变换和极谐波变换的不变矩目标识别方法，其特征在于步骤2b)中所述Sobel边缘检测方法的具体步骤是第一，分别使用水平方向和垂直方向的高斯卷积滤波器对图像进行卷积运算，得到图像的水平方向和垂直方向的偏导数；第二，分别计算图像的梯度和模值；第三，在图像梯度方向上比较图像模值，选择梯度方向上模值最大的图像像素作为边缘像素，并设置该像素值为1，得到边缘图像。
3.根据权利要求1所述的基于Radon变换和极谐波变换的不变矩目标识别方法，其特征在于步骤4b)中所述的Φ角度在0° -360°之间，以1°为增量。
4.根据权利要求1所述的基于Radon变换和极谐波变换的不变矩目标识别方法，其特征在于步骤fe)中所述的极地复指数变换按以下公式得到
5.根据权利要求1所述的基于Radon变换和极谐波变换的不变矩目标识别方法，其特征在于步骤fe)中所述的Radon复指数不变矩按下式构造
6.根据权利要求1所述的基于Radon变换和极谐波变换的不变矩目标识别方法，其特征在于步骤恥)中所述的极地正余弦变换的具体步骤是第一，按以下公式对图像进行极地余弦变换
7.根据权利要求1所述的基于Radon变换和极谐波变换的不变矩目标识别方法，其特征在于步骤恥)中所述的Radon正余弦不变矩(RPCST)按下列公式构造
8.根据权利要求1所述的基于Radon变换和极谐波变换的不变矩目标识别方法，其特征在于步骤5c)中所述的极地复指数实虚不变矩(RI)按下列公式构造ZM{nJ) = f[{MPCET{nJ)),.=1其中，ZKI(n，l)为得到的η阶1重极地复指数实虚不变矩；R为极地复指数变换的实部； I为极地复指数变换的虚部；ζ为不变矩；η为不变矩的阶数；1为极地正余弦变换的重数；N 为正整数；i为1到N之间的正整数；MraT(n，1)为图像P' (r, Θ+Φ)的n阶1重极复指数变换后得到的图像函数；PCET为极地复指数变换。
全文摘要
本发明公开了基于Radon变换和极谐波变换的不变矩目标识别方法。其步骤为(1)输入待识别图像；(2)图像预处理；(3)Radon变换；(4)仿射变换；(5)构造不变矩；(6)提取不变特征；(7)特征模型；(8)图像目标识别；(9)输出图像目标识别结果。本发明成功构造出Radon复指数不变矩、Radon正余弦不变矩和极地复指数实虚不变矩三个新的不变矩，通过提取不变矩的实部和虚部作为不变特征，能够有效地解决噪声干扰的问题，更好的体现图像的真实性，提高图像目标识别的准确率。本发明具有较好的适用性和稳定性，提高了不变矩的整体性能以及图像目标识别的适用性和稳定性。
文档编号G06K9/66GK102324045SQ20111025738
公开日2012年1月18日 申请日期2011年9月1日 优先权日2011年9月1日
发明者刘娟, 史俊杰, 李伟生, 王一丁, 王煦, 苗启广, 许鹏飞, 陈为胜 申请人:西安电子科技大学