发明名称 基于SRC-DP与LDA相结合的图像特征提取方法
摘要 本发明公开了一种基于SRC-DP与LDA相结合的图像特征提取方法,主要解决现有特征提取方法因为仅考虑重构关系或判别关系,导致无法更加准确描述样本信息从而使得图片分类结果不理想的问题。其实现步骤是:1.输入训练样本,计算样本的类内、类间离散度矩阵,初始化投影矩阵,2.对训练样本进行投影,并依次求解投影样本的稀疏表示系数;3.分别计算投影样本的类内、类间重构离散度矩阵;4.构造目标函数求解新的投影矩阵;5.对步骤2-4进行迭代,直至循环次数大于给定初始值,输出最终投影矩阵。本发明增强了图像分类的准确性,提高了分类识别率,可用于警务系统中人物身份的判别或影像拍摄中物体的查找。
申请公布号 CN104715266A 申请公布日期 2015.06.17
申请号 CN201510109065.4 申请日期 2015.03.12
申请人 西安电子科技大学 发明人 刘阳;高全学;王勇;王前前
分类号 G06K9/66(2006.01)I;G06K9/46(2006.01)I 主分类号 G06K9/66(2006.01)I
代理机构 陕西电子工业专利中心 61205 代理人 王品华;张问芬
主权项 一种基于SRC‑DP与LDA相结合的图像特征提取方法,包括如下步骤:(1)将所选图像库中的每一幅图片转化为向量进行存储,并分别抽取图像库中每类图像的一部分构成训练样本集<img file="FDA0000681150680000011.GIF" wi="718" he="92" />其中R表示实数域,d表示原始空间中训练样本的维度,C表示训练样本的类别数,N<sub>i</sub>表示第i类的训练样本数,,<img file="FDA0000681150680000012.GIF" wi="234" he="85" />表示所有训练样本的总数;(2)利用线性判别分析方法LDA分别计算原始训练样本的类间离散度矩阵S<sub>b</sub>和类内离散度矩阵S<sub>w</sub>;(3)给定权重系数β,用来调节判别结构与重构结构的比例关系,β取值范围为[0,1],设定最大迭代步数m,最大迭代步数m由收敛效果决定,为了求解投影到低维空间中训练样本的稀疏系数,初始化投影矩阵P=P<sub>1</sub>∈R<sup>d×p</sup>,其中p为低维空间中训练样本的维度,p<<d,假设当前迭代次数k=1;(4)将训练样本集X中的每一个训练样本投影到p维空间上,即<img file="FDA0000681150680000016.GIF" wi="235" he="80" />i=1,...,C,j=1...,N<sub>i</sub>,<img file="FDA0000681150680000017.GIF" wi="58" he="77" />表示原始空间中第i类的第j个训练样本,<img file="FDA0000681150680000018.GIF" wi="51" he="73" />表示p维空间中第i类的第j个训练样本;(5)通过计算下式求解<img file="FDA0000681150680000019.GIF" wi="57" he="78" />对应的稀疏表示系数向量<img file="FDA00006811506800000110.GIF" wi="70" he="79" />且保证所求的<img file="FDA00006811506800000111.GIF" wi="44" he="79" />的l<sub>1</sub>范数最小:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>min</mi><msub><mrow><mo>|</mo><mo>|</mo><msubsup><mi>s</mi><mi>j</mi><mi>i</mi></msubsup><mo>|</mo><mo>|</mo></mrow><mn>1</mn></msub><mi>s</mi><mo>.</mo><mi>t</mi><mo>.</mo><msubsup><mi>y</mi><mi>j</mi><mi>i</mi></msubsup><mo>=</mo><msubsup><mi>As</mi><mi>j</mi><mi>i</mi></msubsup><mo>,</mo></mrow>]]></math><img file="FDA0000681150680000013.GIF" wi="499" he="101" /></maths>其中<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>A</mi><mo>=</mo><mo>[</mo><msubsup><mi>y</mi><mn>1</mn><mn>1</mn></msubsup><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><msubsup><mi>y</mi><mrow><mi>j</mi><mo>-</mo><mn>1</mn></mrow><mi>i</mi></msubsup><mo>,</mo><msubsup><mi>y</mi><mrow><mi>j</mi><mo>+</mo><mn>1</mn></mrow><mi>i</mi></msubsup><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><msubsup><mi>y</mi><msub><mi>N</mi><mi>c</mi></msub><mi>C</mi></msubsup><mo>]</mo><mo>,</mo></mrow>]]></math><img file="FDA0000681150680000014.GIF" wi="568" he="87" /></maths>i=1,...,C,j=1...,N<sub>i</sub>;(6)由第(5)步中求解所得的训练样本的稀疏表示向量<img file="FDA0000681150680000015.GIF" wi="76" he="82" />并利用基于判别投影的稀疏表示分类器方法SRC‑DP分别计算原始训练样本的类间重构离散度矩阵E<sub>b</sub>和类内重构离散度矩阵E<sub>w</sub>;(7)根据(2)和(7)计算出的离散度矩阵S<sub>b</sub>,S<sub>w</sub>,E<sub>b</sub>,E<sub>w</sub>构建如下广义特征方程(E<sub>b</sub>+βS<sub>b</sub>)α<sub>i</sub>=λ<sub>i</sub>(E<sub>w</sub>+βS<sub>w</sub>)α<sub>i</sub>    i=1,...,p;(8)求解(7)中广义特征方程中λ<sub>i</sub>对应的特征向量α<sub>i</sub>,且满足λ<sub>1</sub>≥...≥λ<sub>i</sub>...≥λ<sub>p</sub>,得到λ<sub>i</sub>对应特征向量α<sub>i</sub>的投影矩阵P=[α<sub>1</sub>,...α<sub>i</sub>,...α<sub>p</sub>];(9)令k=k+1,当k≤m时,返回步骤(4),直至k>m时,输出最终的投影矩阵P'。
地址 710071 陕西省西安市太白南路2号