发明名称 基于特征扩展和模糊支持向量机的多分组图像分类方法
摘要 基于特征扩展和模糊支持向量机的多分组图像分类方法,属于图像处理领域,本发明为解决现有多分组图像分类方法无法有效提取图像本质特征,分类精度较低的问题。本发明方法为:先去掉被噪声等严重污染而无法使用的波段,对剩下的波段依次进行二维经验模态分解,得到一些二维本征模态函数,将这些二维本征模态函数进行有机组合,扩展为多分组图像的特征,最后用模糊支持向量机分类器分类。本发明充分发挥了二维经验模态分解能自适应地提取复杂图像本质特征的优势,有效获取多分组图像特征;采用模糊支持向量机作为分类器,兼具支持向量机和模糊函数二者的优势,有利于提高分类精度。
申请公布号 CN102622611B 申请公布日期 2014.02.19
申请号 CN201210065533.9 申请日期 2012.01.13
申请人 哈尔滨工业大学 发明人 沈毅;贺智;张淼
分类号 G06K9/62(2006.01)I 主分类号 G06K9/62(2006.01)I
代理机构 哈尔滨市松花江专利商标事务所 23109 代理人 张果瑞
主权项 1.基于特征扩展和模糊支持向量机的多分组图像分类方法,其特征在于,该方法包括以下步骤:步骤一、初始化给定波段数为I<sub>0</sub>、尺寸为P×Q的多分组图像IM<sub>j</sub>(p,q),j=1,2,...,I<sub>0</sub>,p=1,2,...,P,q=1,2,...,Q,去掉多分组图像IM<sub>j</sub>(p,q)中峰值信噪比在0~8之间的被噪声严重污染而无法使用的波段,对剩下的I个有效波段重新排序,得到待扩展多分组图像IM<sub>i</sub>(p,q),i=1,2,...,I,其中I<sub>0</sub>、I、P和Q为自然数;步骤二、依次对I个待扩展多分组图像IM<sub>i</sub>(p,q)进行二维经验模态分解,得到I组二维本征模态函数BIMF<sub>i</sub><sup>ui</sup>,u<sub>i</sub>=1,2,...,U<sub>i</sub>,其中U<sub>i</sub>为第i个波段的二维本征模态函数的个数;得到I组二维本征模态函数BIMF<sub>i</sub><sup>ui</sup>的过程为:步骤A、初始化r<sub>1</sub>=IM<sub>i</sub>(p,q);u=1;v=0;SD=1000;h<sub>u</sub>,<sub>v</sub>=r<sub>1</sub>;c<sub>u</sub>=r<sub>1</sub>,其中,第i个待处理波段的待扩展多分组图像为IM<sub>i</sub>(p,q),r<sub>1</sub>为第1次二维经验模态分解后的残差,SD为终止迭代阈值,h<sub>u,v</sub>为第u次二维经验模态分解中经过第v次筛选后的剩余函数,c<sub>u</sub>为二维本征模态函数BIMF<sub>i</sub><sup>ui</sup>;步骤B、令v=v+1;h<sub>u</sub>,<sub>(v-1)</sub>=r<sub>u</sub>,并通过与相邻像素对比,找出h<sub>u</sub>,(<sub>v-1</sub>)的极大值和极小值,其中,r<sub>u</sub>为第u次二维经验模态分解后的残差,h<sub>u,(v-1)</sub>为第u次二维经验模态分解中经过第v-1次筛选后的剩余函数;步骤C、对步骤B检测出的第u次二维经验模态分解中经过第v-1次筛选后的剩余函数h<sub>u,(v-1)</sub>的极大值和极小值,构造Denaulay三角形网格,由双三次样条插值法求取图像的上包络e<sub>max</sub>和下包络e<sub>min</sub>,并计算包络均值<maths num="0001"><![CDATA[<math><mrow><mi>Ave</mi><mo>=</mo><mfrac><mrow><mo>(</mo><msub><mi>e</mi><mi>max</mi></msub><mo>+</mo><msub><mi>e</mi><mi>min</mi></msub><mo>)</mo></mrow><mn>2</mn></mfrac><mo>;</mo></mrow></math>]]></maths>步骤D、计算h<sub>uv</sub>=h<sub>u(v-1</sub>)-Ave和终止准则<maths num="0002"><![CDATA[<math><mrow><mi>SD</mi><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>p</mi><mo>=</mo><mn>1</mn></mrow><mi>P</mi></munderover><munderover><mi>&Sigma;</mi><mrow><mi>q</mi><mo>=</mo><mn>1</mn></mrow><mi>Q</mi></munderover><mo>[</mo><mfrac><msup><mrow><mo>|</mo><msub><mi>h</mi><mrow><mrow><mo>(</mo><mi>u</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mi>v</mi></mrow></msub><mrow><mo>(</mo><mi>p</mi><mo>,</mo><mi>q</mi><mo>)</mo></mrow><mo>-</mo><msub><mi>h</mi><mi>uv</mi></msub><mrow><mo>(</mo><mi>p</mi><mo>,</mo><mi>q</mi><mo>)</mo></mrow><mo>|</mo></mrow><mn>2</mn></msup><mrow><msubsup><mi>h</mi><mrow><mrow><mo>(</mo><mi>u</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mi>v</mi></mrow><mn>2</mn></msubsup><mrow><mo>(</mo><mi>p</mi><mo>,</mo><mi>q</mi><mo>)</mo></mrow></mrow></mfrac><mo>]</mo><mo>;</mo></mrow></math>]]></maths>步骤E、判断SD&lt;ε<sub>0</sub>是否成立,ε<sub>0</sub>=0.2~0.3,如果是,执行步骤F;否则,返回步骤B;步骤F、提取第u个二维本征模态函数BIMF<sub>i</sub><sup>ui</sup>:c<sub>u</sub>=h<sub>uv</sub>;和第u+1个残差r<sub>u+1</sub>=r<sub>u</sub>-c<sub>u</sub>;步骤G、判断r<sub>u+1</sub>是否单调,判断结果为否,令v=0;u=u+1,并返回步骤B;判断结果为是,终止程序,获取u个二维本征模态函数BIMF<sub>i</sub>ui和一个残差r<sub>u+1</sub>;步骤三、将所有二维本征模态函数BIMF<sub>i</sub><sup>ui</sup>有机组合,扩展为待扩展多分组图像IM<sub>i</sub>(p,q)各波段对应的特征,记为扩展后特征FBIMF;获取扩展后特征FBIMF的过程为:对第i个波段的二维本征模态函数依次进行处理,获取每个波段的二维本征模态函数对应的扩展后特征FBIMF,每个波段的二维本征模态函数的扩展后特征FBIMF的获取方法相同,为:对该波段二维本征模态函数中第u<sub>i</sub>个BIMF来说,如果u<sub>i</sub>为奇数,则将BIMF<sub>i</sub><sup>ui</sup>直接衔接到BIMF<sub>i</sub><sup>ui</sup>-<sup>1</sup>尾部,如果u<sub>i</sub>为偶数,则将BIMF<sub>i</sub><sup>ui</sup>左右翻转再衔接到BIMF<sub>i</sub><sup>ui</sup>-<sup>1</sup>尾部,依次下去,直至该波段二维本征模态函数的所有U<sub>i</sub>个BIMF全部处理完毕,获取该波段的扩展后特征FBIMF;步骤四、随机抽取扩展后特征FBIMF上的各像素形成训练样本F<sub>training</sub>BIMF<sub>k</sub>,k=1,2,...,N和测试样本F<sub>testing</sub>BIMF<sub>l</sub>,l=1,2,...,Mw其中,N为训练样本总数,M为测试样本总数;步骤五、对训练样本F<sub>training</sub>BIMF<sub>k</sub>指定模糊隶属度s<sub>k</sub>,其中,0≤s<sub>k</sub>≤1;步骤六、根据训练样本F<sub>training</sub>BIMF<sub>k</sub>及对应的模糊隶属度s<sub>k</sub>训练FSVM分类器;FSVM分类器为多分类器,由多个二分类器根据一对一策略,一对多策略或分层树策略扩展而成,每个二分类器按公式<img file="FDA00003621726200031.GIF" wi="1093" he="154" />进行构建,其中,w为权值向量,x<sub>k</sub>为训练样本F<sub>training</sub>BIMF<sub>k</sub>中的元素,<img file="FDA00003621726200032.GIF" wi="37" he="43" />(x)为从x<sub>k</sub>映射到高维特征空间的函数,b为偏置,y<sub>k</sub>为x<sub>k</sub>对应类别,且y<sub>k</sub>∈{-1,+1},α<sub>k</sub>={α<sub>1</sub>,α<sub>2</sub>,...,α<sub>N</sub>}为拉格朗日因子,且k=1,2,...,N,K(x,x<sub>k</sub>)为核函数;步骤七、通过测试样本F<sub>testing</sub>BIMF<sub>l</sub>测试训练出的模糊支持向量机,并求其分类精度,完成对多分组图像IM<sub>j</sub>(p,q)的分类。
地址 150001 黑龙江省哈尔滨市南岗区西大直街92号