第6卷第1期 智能系统学报 Vo1.6 No.1 2011年2月 CAAI Transactions on Intelligent Systems Feb.2011 doi:10.3969/j.issn.1673-4785.201 1.01.006 采用独立分量分析Zernike矩的 遥感图像飞机目标识别 刘富 ,于鹏 ,刘坤 (1.吉林大学通信工程学院,吉林长春130022;2,哈尔滨工业大学深空探测基础研究中心,黑龙江哈尔滨150001) 摘要:为了提高遥感图像目标自动识别系统的准确性,提出了一种新的采用独立分量分析(ICA)Zernike矩的飞机 目标识别方法.首先对分割后的目标区域进行独立分量分析处理,将待识别目标形状转换到标准形式,然后对标准 化后的图像目标提取Zenfike矩作为特征向量进行识别.通过实验表明此方法具有鲁棒性,能有效地消除遥感图像 目标尺度、旋转、平移、反转和扭曲影响,能够有效地识别遥感图像飞机目标. 关键词:独立分量分析;Zernike矩;遥感图像;飞机目标识别 中图分类号:TP391.41 文献标识码:A文章编号:1673-4785(2011)01 051-06 Research concerning aircraft recognition of remote sensing images based on ICA Zernike invariant moments LIU Fu ,YU Peng ,LIU Kun (1.College of Communications Engineering,Jilin University,Changchun 130022,China;2.Deep Space Exploration Research Center, Harbin Institute of Technology,Harbin 150001,China) Abstract:To promote the accuracy of an automatic target recognition(ATR)system for remote sensing images,a novel feature recognition method was proposed for an airplane target based on the independent component algorithm (ICA)with Zernike invariant moments.First,the divided region of interest(ROI)was analyzed with the ICA method.Also,the shape of the target was changed to canonical form.Then,the invariant moments of normalized shapes could be extracted.They would potentially be used as a feature vector to do further recognition.The experi— ment demonstrates that the performance of this method is robust.It can eliminate the effects of scaling,rotation, translation,reflection,and skewing of the remote sensing image target.Furthermore,this method can recognize an airplane’S target from the remote sensing images effectively. Keywords:independent component algorithm(ICA);Zernike invariant moments;remote sensing images;aircraft target recognition 随着遥感技术的发展,遥感影像广泛应用于环 在物体识别的过程中,物体的形状是非常重要 境监测、资源调查、军事侦查等众多领域.遥感图像 的特征之一.形状特征是图像的核心特征,也是人类 目标自动识别成为研究的热点,特别是军事目标的 视觉系统进行物体识别时所识别的关键信息之一. 自动识别领域更是研究的重点.遥感影像在拍摄过 具有仿射不变关系的目标匹配、识别一直是物体形 程中,摄像机的拍照距离要远远大于景物本身的尺 状特征研究的重点 4。. 寸.由于遥感影像是远距离成像,因此在高空的遥感 遥感图像在提取目标物体的形状特征前,首先 感应元件的主光轴与地面的关系可以看成是垂直关 要对遥感图像进行预处理,消除噪声对图像的影响. 系,由摄像机透视成像原理可知同一景物目标在不 然后通过阈值分割的方法将目标区域提取出来,对 同遥感图像中可以用仿射变换来近似表示…. 感兴趣区域进行处理并提取不变特征.这里提出了 一种基于独立分量分析Zernike矩的目标形状特征 收稿日期:2010 ̄34-20. 基金项目:军队科研基金资助项目(9140A03040508HT0178) 识别方法.首先将分割结果进行区域连通和形态学 通信作者:刘富.E・mail:liufu@jlu.edu.cn. 操作,提取目标的轮廓.然后对轮廓进行独立分量分 -52・ 智能系统学报 第6卷 析处理,并利用高阶中心矩将目标区域轮廓转换成 标准形式.对标准化后的图像提取Zernike不变矩作 为特征向量进行识别.实验表明具有仿射变换关系 的图像经过ICA处理后具有一致的标准形式,本方 法能够有效地实现遥感图像目标的自动识别. X =∑aijSj(t),i=1,2,…, =1 写成向量形式为 X(t)W--As( ). ICA用于图像特征提取的方法 引通常是将k 个类别的m幅图像(共k×m幅)看作是由 ×m个 统计独立的源图像经过线性混合而得到的,通过寻 找分离矩阵w,得到分离的各独立分量s,并以独立 1 仿射变换 由遥感图像的成像原理可知,同一景物在不同 遥感图像中表现形式不同,这种不同可由仿射变换 关系表示.1个仿射变换可以表示成1个非奇异的 分量为基向量构造特征空间,将每一幅待识别目标 图像投影到特征空间,即用独立分量的线性组合表 线性变换加上1个平移变换: =Ax+ 式中:A是仿射矩阵,满足非奇异性,即det(A)≠0. 该仿射矩阵用来表示图像的尺度、旋转和视角改变 的程度,A可以分解表示成: A=[ 吕一…sin 0 ]. 式中:s、 、0分别代表了变化尺度、扭曲和旋转角 度; 和 分别是沿着 和Y轴方向的尺度变化因 子.2×1阶矩阵 表示相对坐标平移,X是坐标X经 过仿射变换后的结果.这样,仿射变换可以看成是一 系列变换产生的组合. 这里假设有图形 和图形y,两者具有仿射变 换关系,可以表示为 m1 x i ]・ ㈩ 式中:(互 , )、( ,Y)分别为图形1y和图形l,上点的 坐标;tx和t 表示平移,将它们移到等式的左边可 以得到: m1 l2- 22㈣Yl由式(2)可知通过坐标平移,将坐标原点移到图形 的中心 ],就可将式(1)简化为 :Al,. 2 基于ICA的目标形状标准化 独立分量分析(independent component analysis, ]CA)是盲信号分离(blintt souTce separation,zss)的 一个重要分支 ],它在不知道盲信号的混合矩阵和 数量的情况下就能得到各个输人信号.假设存在Ⅳ 个独立的源信号,表示成矢量形式:S(t)= [ I(f)S2( )…SⅣ(£)] ,t=0,1,2,…,以及 个观 测信号X(t)=[X (t)X (t)…X (t) 个观测 信号是由Ⅳ个源信号线性瞬时混合而成的,即在每 个时刻t都有如下关系式: 示一个目标图像.此方法多用于降低图像噪声和提 取图像特征. 设 是待识别的目标,则 =aISl+a2¥2+…+ak ̄mS × . 式中: 是独立分量基,(o ,o ,…, )则是待识 别目标 在特征空间的投影系数,用式(3)求得: (。l,n2,…,0 )= ×( )一 . (3) 提取投影系数作为目标提取的特征值,构造合适的 分类器实现分类. 本文提出的算法是将分割后的目标形状数据转 化为标准形式,使用独立分量分析和高阶中心矩来实 现标准化.基于ICA的目标形状标准化分2步完成. 1)是先对目标形状区域进行PCA白化处理,由 上文知 和y是具有仿射变换关系的目标图形.由 文献[8]可知对具有仿射变换关系的图形区域进行 PCA白化处理,可将它们之问的仿射变化关系转换 为正交等式关系: Yh:Qyb. 式中:y’h和yh分别是'y和y经过白化处理所得到的 结果,Q是正交矩阵可以由看出图形yh可以由yb经 过尺度、平移等变换获得.图1中(a1)、(a2)和(a3)是 具有仿射变换关系的图形,(b1)、(b2)和(b3)分别是 (a1)、(a2)和(a3)经过PCA白化处理后的结果. (a1) 目 (a2) (a3) b,-A (b2) (b3) 图1 形状图片和经过PCA白化处理结果 Fig.1 Shape pictures and processing results by PCA whitening 第1期 刘富,等:采用独立分量分析Zernike矩的遥感图像飞机目标识别 ・53・ 2)采用快速独立分量分析算法(FastICA) ], 求得矩阵Q.这里如果假设’yh为标准型,根据求得 的矩阵Q可将y 转化为'y .同理可知任意具有仿 射变换关系的图形均可由ICA转化为统一的标准 形式. 假设信号源含有2个独立分量,有2个观测信 号的独立分量分析模式为 X[k]=M・f[k]. : 三 y、。。, l2J>J 1J; 其他. 这样就能确保 是具有仿射关系的图形的标 准形式.图2中(aI)、(b1)和(c1)是具有仿射关系 的图像,(a2)、(b2)和(c2)是经过ICA处理后的图 像,(a3)、(b3)和(c3)是经过三阶中心矩旋转之后 的图像.可以看到经过ICA处理和高阶矩旋转后的 图形具有一致的标准形式. 式中: k]=[i [k]i [k]r表示是2个独立信号 源的矢量表示形式; [k]=[X [k]X [k]r是2个 观测信号的矢量表示形式; 是索引数据;M是一个 2×2的混合矩阵.如果观测信号 [k]已知,就可以 通过ICA算法计算出2个独立状态分量f [ ]、 i [k].如果观测信号X [k]、X [k]是待识别目标形 状区域坐标,而 [k]和 :[k]是与目标形状有仿射 变换关系的坐标,通过公式,有 『k]: .f『k]. (4) 式中: 是一个由A和 构成的新的混合矩阵,对 于式(4)通过ICA能够得相同的独立分量i[k]= [f [ ]f [ ]r,无论是X还是 的独立分量都是相 同的.将计算出的独立分量f [k]和 [k]作为目标 形状标准形式的像素坐标,得到新的目标标准形状 ,即上文中的 . 从图2可以看到具有仿射关系的图形经过独立 分量分析处理后得到的图像之间有旋转,采用三阶 中心距将经过ICA处理后的形状坐标转化到一致 A㈤V㈦的标准形式 l .计算公式如下: = 式中: 是经过ICA处理后的图形像素坐标;v㈦V㈦V 是 经过三阶中心距转化后的目标标准形式.定义三阶 中心距11,】I2、“2.1为: 1.2=,nl2—2 玎L11一otmo2+2卢 玎}1o; H2.1=m2I一2am11一卢m20+2a m01. 式中:m。。是目标形状的二维矩,定义为 1 N -I。口 mpq :y ; = =m——oo , ——, 。 mo0 和 即是给定目标形状的重心坐标值.因为采用 的是FastICA算法,ICA处理过程是从初始矩阵 Q(0)开始的,由于Q(0)的不确定性使得l,。 可能 存在关于直线Y= 对称的2种可能,利用三阶中心 矩的性质执行以下操作: 图2Fig.2 Shape pi 形状图片和经过ICA处理的标准型 ctures and normalized results by ICA 3 基于标准化形状的不变矩特征提取 Zemike矩特征具有旋转不变性等特征,目前已 在目标识别领域中得到较为广泛的研究应用¨ . Zernike多项式是{( ( ,Y)}=R (r)e ,其中尺 是实值径向多项式设.设_厂( ,Y)是图像亮度函数, 阶的Zernike矩定义为 z = n 233[ ( ,y)]Ef( ,y)]dxdy. +v2≤1 Zemike矩具有旋转不变性、对噪声的健壮性、 没有信息冗余和表示方法有效等优点.这里提取标 准形式的前36个Zernike矩作为特征向量进行识 别.图3中(a)、(b)和(C)是对图2中标准形式 (a3)、(b3)和(c3)提取的Zernike矩.实验结果表明 经过ICA标准化后,得到具有仿射关系的图形的 Zernike矩非常相似,可以大大减少仿射变换对识别 的影响. ・54・ 智能系统学报 第6卷 6 5 4 3 2 1 0 0 O 图3标准型提取的Zernike不变矩 Fig.3 Zernike invariant moments extraction based On normalized model 4实验结果及分析 4.1实验1 实验采用基于图像特征的匹配方法对飞机目标 进行识别.实验共采用了80幅卫星图像数据,对 180个飞机目标进行识别,所用数据为QuickBird卫 星图像,图像分辨率为0.64 m.飞机形状提取采用 的方法是改进区域分割方法¨引,使用快速ICA算法 完成目标形状的标准化,采用最小欧氏距离准则的 特征匹配实现分类. 图4和图5是含有相同型号飞机的不同卫星遥 感图片处理结果,图6是含有不同型号飞机的卫星 图片处理结果.图4(b)、图5(b)和图6(b)分别是 飞机目标形状经过ICA处理后的标准形式;图4 (c)、图5(c)和图6(c)是对标准形式图像提取的 Zernike不变矩特征. fa)某型号 机遥感罔像 1O 2O 30 40 (b)标准化结果 fc)提取的Zemike矩 图4实验结果1 Fig.4 The results of experiment 1 机 图5实验结果2 Fig.5 The results of experiment 2 (a)不同型号毪机遥感图像 一~ ㈣标准化结果 (c)提取的Zemike矩 图6实验结果3 Fig・6 The result of experiment 3 实验结果分析:可以看到在不同遥感图像中的 同一型号飞机的ICA标准形式相似,同一类型飞机 的Zernike矩几乎相同的;同时也可以看到不同型号 飞机的Zernike矩有明显不同,提取的特征具有良好 的可分性.本实验正确识别飞机174个,识别率为 96.7%,实验结果表明本方法是一种有效的遥感图 像飞机目标识别算法. 4.2实验2 本实验的图片是从卫星图片中分割出的5种型 号飞机的形状,如图7.对飞机形状进行仿射变换, 分别取仿射变换的参数:0∈{0。,30。,60。,108。}、 (s =S =1)、(S =S =2)、(s =1,5 =2)、 ∈{0, 第l期 刘富,等:采用独立分量分析Zernike矩的遥感图像飞机目标识别 .55. 0.25,0.5,0.75},每种飞机得到48种经过尺度变换 和角度变换的测试样本图像,共得到240幅测试样 5 结束语 本图像. 图7从卫星图片中分割出的5种型号飞机形状 Fig.7 The 5 kinds of plane contour divided from sens- ing images 图8为遥感图像飞机形状仿射变化结果,(a) 的变换参数为0=0。,S =S =l,OL=0;(b)的变换参 数为0=30。,s 一 =1,O/=0.25;(C)的变换参数为 0=60。,S =S =2, =0.75;(d)的变换参数为0= 108。,s =s =2,OL=0.5;(e)的变换参数为0= 108。,S =1, =2; =0. 口圜囵 口 ra) (b) (C) (d) (e) 图8 飞机形状仿射变换结果 Fig.8 The results of afifne transform 分别使用本文提出的基于ICAZernike不变矩、 Zernike不变矩” 、ICA不变矩 和仿射不变矩 。 作为特征向量,设计神经网络分类器进行识别,验证 各种方法的分类效果,识别结果如表1所示. 表1实验结果 Table 1 The results of experiment % 实验结果分析:通过实验可以看到1CAZernike 矩比其他不变矩表现好,经典Zernike矩对图像尺度 变化和旋转是不变的,但无法消除扭曲对图像的影 响.本文算法相对于经典Zemike矩利用ICA标准 化消除了仿射变换关系图像中尺度变换和图像扭曲 的影响,所以分类效果好. 我们也可以看到本文算法识别率高于ICA不 变矩,这是因为在提取ICA不变矩时,提取的是待 识别图像在全部独立分量为基向量的特征空间中的 投影系数,将投影系数作为目标识别的特征量,这就 使得这种方法易受到噪声的干扰.而在本文提出的 算法中,使用的是前36个Zernike矩,从而消除了高 阶噪声的影响.实验结果表明ICAZernike矩优于其 他不变矩. 本文提出了一种采用独立分量分析Zernike矩 的飞机目标识别方法,该方法基于目标分割,对目标 区域进行独立分量分析处理并利用高阶中心矩将目 标形状转换成标准形式,提取Zernike不变矩作为特 征向量进行飞机目标识别.实验结果证明,该方法性 能鲁棒、抗干扰性能好、识别率高,能够准确识别高 分辨率遥感图像中的飞机目标.随着遥感技术的不 断发展,遥感图像分辨率越来越高,能够获得更准 确、更完善的飞机轮廓,提出的方法将具有更高的实 际应用价值.下一步的工作是利用本文方法运算速 度快、准确率高的特点,将算法应用于实时遥感图像 的飞机目标自动识别,也可扩展到其他目标的自动 识别. 参考文献: [1]安玮,李宏,徐晖,等.模式识别巾的透射变换与仿射变 换[J].系统工程与电子技术,1999,21(1):57—62. AN Wei,LI Hong,XU Hui,et a1.Projective tansformation and afine tanstbrmation in pattern recognition[J].Systems Engineering and Electronics,1999,21(1):57-62. [2]ZULIANI M,BHAGAVATHY S,MANJUNATH S.Afifne— invariant eurve matching[C]//IEEE International Confer— ence on Image Processing.Singapore,2004,5:3041— 3044. [3]AVRITRIS Y,XIROUHAKIS Y,KOI I I AS S.Affine—invar— iant curve normalization f0r object shape representation. classiifcation and retrieval_J].Machine Vision and Appli— cations,2001,13(2):80-94. [4]WOLOVICH W A,UNEI M.The determination of implicit polynomial canonical cnl3/es l J].IEEE Transactions on Pat— tern Analysis and Machine Intelligence,1 998,20(1 0): 1080—1 090. [5]LEE T W.Independent conlponent analysis theory and ap— plications[M].Boston,USA:Kluwer Academic Publish— eps,1998:87—109. [6]LEE J J,UDDIN M Z,KIM T S.Spatiotemporal human fa— cial expression recognition using fisher independent compo— nent analysis and hidden Markov model[C]//The 30th An— nua]International 1EEE EMBS Conference.Vancouver, Canada。2008:2546—2549. [7]LI Yunxia,FAN Changyuan.Face recognition by nonnega— tive independent component analysis[C]//2009 Fifth Inter- national Conference on Natural Computation._rianjin.Chi— na,2009,2:555—558. [8]FUKUNAGA K.Introduction to statistical pattern recogni— tion[M].2nd ed.San Diego,USA:A(:ademic Press, ・56・ 1990:417-425. 智能系统学报 第6卷 Journal of Inner Mongolia Agricultural University,2009, [9]HYVRINEN A.Fast and robust fixed—point algorithms for independent component analysis[J].IEEE Transactions on Neural Networks,1999,10(3):626-634. 30(2):222—225. 作者简介: 刘富,男,1968年生,教授、博士生 导师,主要研究方向为计算机视觉、模 式识别等.承担国家“863”计划项目等 [1O]王晓红,赵荣椿.矩技术在计算机图像学中的应用综述 [J].中国体视学与图像分析,2002,7(1):53-57. WANG Xiaohong,ZHAO Rongchun.Moment technique 30余项,获得发明专利2项,发表学术 论文40余篇. applications in computer image science:a review[J].Chi— nese Journal of Stereology and Image Analysis,2002,7 (1):53-57. [1 1]KHOTANZAD A,HONG Y H.Invariant image recognition by Zernike moments『J].IEEE Transactions on Pattern A— 于鹏,男,1978年生,硕士研究生,主 要研究方向为计算机视觉、模式识别. nalysis and Machine Intelligence,1990,12(5):489-497. [12]郑楠.遥感图像中飞机目标的识别分类算法研究[D]. 长春:吉林大学电子科学与工程学院,2007. ZHENG Nan.Algorithm research for aircraft recognition and classification in remote sensing images[D].Chang— chun:College of Electronic Science and Engineering,Jilin University,2007. 刘坤,男,1985年生,硕士研究生,主 要研究方向为计算机视觉、模式识别. [13]郝敏,麻硕士,侯振杰.Zernike矩的不变性与计算实现 [J].内蒙古农业大学学报,2009,30(2):222—225. HA0 Min,MA Shuoshi,HOU Zhenjie.Invarance analysis and computation implementation of Zernike moments[J]. ■■■ 2011中国智能(产业)博览会 201 1.9.16一l8北京全国农业展览馆 China International Intelligent Industries Expo 20 1 1 为贯彻落实温家宝总理提出“感知中国”的科技纲领,推进国家“十二五”建设的新举措,激发广大科技工作 者投身建设创新型国家的积极性和创造热情,进一步推动智能科学技术担当起信息化向智能化高端发展进程中 的历史重任,欣逢纪念中国人工智能学会成立30周年和世界人工智能科学诞生55周年,经中国科学技术协会批 准,得到国家科学技术部、工业和信息化部、教育部、住房和城乡建设部、中国科学院、中国工程院、国家自然科学 基金委员会悉心指导,由中国人工智能学会主办,中国电子学会、中国计算机学会、中国通信学会、中国自动化学 会、中国系统仿真学会、中国仪器仪表学会、中国图象图形学学会、中国中文信息学会、中国计算机用户协会支持 合作,北京东方炬峰展览有限公司承办的“首届中国智能博览会”(简称智博会)将于2011年9月16—18日在北 京全国农业展览馆盛大开幕.届时,本届“智博会”组委会将诚邀科技、教育、农业、机械、轻工、通信、电力、交通、建 筑、金融、卫生、安全、国防等专业领域的政府官员、专家学者、国际友人和企业精英参会,共同见证与分享“智能体 验・智慧生活”的核心主题. 本届智博会将秉承和弘扬理解、交流、创新与合作的智博理念,向社会各界传达智能科技和智慧生活的最新 信息,展现丰富多彩的智能新型能源、创新成果、新兴产业以及智慧生活体验,搭建智能产业与科技民生的桥梁, 激发广大热爱智能科学技术的有志之士相互交流,以及实现更美好生活的强烈愿望,深入探讨智能化进程中的多 元时代现象,广泛推进学术、科技、工程、产业以及资本合作的多边共赢之路,关注全球重大问题给国际社会和科 技民生带来的挑战和影响,共同展望物联网、云计算、智能电网、智能交通、智能通信、智能计算机、智能机器人、智 能建筑等智能化领域的美好未来. 网址:www.ciiie.con.cn.