发明名称 基于极限学机的极化SAR图像分类方法
摘要 本发明公开了一种基于极限学机的极化SAR图像分类方法,主要解决极化SAR图像分类的时间效率低的问题。其实现步骤为:输入待分类极化SAR图像标记信息和相干矩阵,进行Lee滤波;提取滤波后的相干矩阵中的特征,并将特征归一化获得数据集,并从数据集中获得标记集;将标记集划分成训练集和测试集,并用训练集训练极限学机;利用训练好的极限学机去预测数据集的类别,得到极化SAR图像分类结果。本发明应用极限学机进行极化SAR图像分类,分类精度较高,同时运行时间少,可用于地物分类与目标识别。
申请公布号 CN105205491A 申请公布日期 2015.12.30
申请号 CN201510512066.3 申请日期 2015.08.19
申请人 西安电子科技大学 发明人 焦李成;李玲玲;曾杰;马文萍;张丹;屈嵘;侯彪;王爽;马晶晶;尚荣华
分类号 G06K9/62(2006.01)I;G06K9/66(2006.01)I 主分类号 G06K9/62(2006.01)I
代理机构 陕西电子工业专利中心 61205 代理人 王品华;朱红星
主权项 一种基于极限学习机的极化SAR图像分类方法,包括如下步骤:(1)输入待分类极化SAR图像标记信息,输入一幅大小为3×3×M的待分类极化SAR图像的相干矩阵,并用Lee滤波器滤除相干噪声,得到滤波后的相干矩阵,其中,滤波后的相干矩阵中每个元素是一个3×3矩阵,M表示待分类的极化SAR图像像素点的总数;(2)将滤波后相干矩阵中每个元素对应的3×3矩阵,拉成一个9维的特征向量,获得大小为9×M的数据矩阵;(3)将数据矩阵进行归一化处理,获得数据集<img file="FDA0000784832690000011.GIF" wi="594" he="97" />其中n=9为特征的维数,x<sub>i</sub>为特征向量;(4)从数据集中获得训练集<img file="FDA0000784832690000012.GIF" wi="58" he="70" />和测试集<img file="FDA0000784832690000013.GIF" wi="96" he="70" />(4a)根据标记信息从数据集<img file="FDA0000784832690000014.GIF" wi="48" he="83" />中获得标记集<img file="FDA0000784832690000015.GIF" wi="85" he="56" /><img file="FDA0000784832690000017.GIF" wi="805" he="68" />其中m为类别总数,N为标记样本的总数,t<sub>i</sub>为类别向量;(4b)将标记集随机选取1%作为训练集<img file="FDA0000784832690000018.GIF" wi="822" he="68" />剩余99%作为测试集<img file="FDA0000784832690000019.GIF" wi="813" he="68" />其中K训练集样本个数,S为测试集的样本个数;(5)用训练集<img file="FDA00007848326900000110.GIF" wi="51" he="67" />训练极限学习机:(5a)设置极限学习机隐含层节点的激活函数为G(a,b,x)=exp(‑b||x‑a||<sup>2</sup>),隐含层单元个数为L,且K>L;(5b)随机设置输入权值a<sub>j</sub>和偏值b<sub>j</sub>,j=1,...,L;(5c)输入层为训练集<img file="FDA0000784832690000016.GIF" wi="830" he="74" />计算隐藏层输出矩阵为H,<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>H</mi><mo>=</mo><msub><mfenced open = "[" close = "]"><mtable><mtr><mtd><mrow><mi>G</mi><mrow><mo>(</mo><msub><mi>a</mi><mn>1</mn></msub><mo>,</mo><msub><mi>b</mi><mn>1</mn></msub><mo>,</mo><msub><mi>x</mi><mn>1</mn></msub><mo>)</mo></mrow></mrow></mtd><mtd><mn>...</mn></mtd><mtd><mrow><mi>G</mi><mrow><mo>(</mo><msub><mi>a</mi><mi>j</mi></msub><mo>,</mo><msub><mi>b</mi><mi>j</mi></msub><mo>,</mo><msub><mi>x</mi><mn>1</mn></msub><mo>)</mo></mrow></mrow></mtd><mtd><mn>...</mn></mtd><mtd><mrow><mi>G</mi><mrow><mo>(</mo><msub><mi>a</mi><mi>L</mi></msub><mo>,</mo><msub><mi>b</mi><mi>L</mi></msub><mo>,</mo><msub><mi>x</mi><mn>1</mn></msub><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mo>.</mo></mtd><mtd><mrow></mrow></mtd><mtd><mo>.</mo></mtd><mtd><mrow></mrow></mtd><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd><mtd><mn>...</mn></mtd><mtd><mo>.</mo></mtd><mtd><mn>...</mn></mtd><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd><mtd><mrow></mrow></mtd><mtd><mo>.</mo></mtd><mtd><mrow></mrow></mtd><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mrow><mi>G</mi><mrow><mo>(</mo><msub><mi>a</mi><mn>1</mn></msub><mo>,</mo><msub><mi>b</mi><mn>1</mn></msub><mo>,</mo><msub><mi>x</mi><mi>i</mi></msub><mo>)</mo></mrow></mrow></mtd><mtd><mn>...</mn></mtd><mtd><mrow><mi>G</mi><mrow><mo>(</mo><msub><mi>a</mi><mi>j</mi></msub><mo>,</mo><msub><mi>b</mi><mi>j</mi></msub><mo>,</mo><msub><mi>x</mi><mi>i</mi></msub><mo>)</mo></mrow></mrow></mtd><mtd><mn>...</mn></mtd><mtd><mrow><mi>G</mi><mrow><mo>(</mo><msub><mi>a</mi><mi>L</mi></msub><mo>,</mo><msub><mi>b</mi><mi>L</mi></msub><mo>,</mo><msub><mi>x</mi><mi>i</mi></msub><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mo>.</mo></mtd><mtd><mrow></mrow></mtd><mtd><mo>.</mo></mtd><mtd><mrow></mrow></mtd><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd><mtd><mn>...</mn></mtd><mtd><mo>.</mo></mtd><mtd><mn>...</mn></mtd><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd><mtd><mrow></mrow></mtd><mtd><mo>.</mo></mtd><mtd><mrow></mrow></mtd><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mrow><mi>G</mi><mrow><mo>(</mo><msub><mi>a</mi><mn>1</mn></msub><mo>,</mo><msub><mi>b</mi><mn>1</mn></msub><mo>,</mo><msub><mi>x</mi><mi>K</mi></msub><mo>)</mo></mrow></mrow></mtd><mtd><mn>...</mn></mtd><mtd><mrow><mi>G</mi><mrow><mo>(</mo><msub><mi>a</mi><mi>j</mi></msub><mo>,</mo><msub><mi>b</mi><mi>j</mi></msub><mo>,</mo><msub><mi>x</mi><mi>K</mi></msub><mo>)</mo></mrow></mrow></mtd><mtd><mn>...</mn></mtd><mtd><mrow><mi>G</mi><mrow><mo>(</mo><msub><mi>a</mi><mi>L</mi></msub><mo>,</mo><msub><mi>b</mi><mi>L</mi></msub><mo>,</mo><msub><mi>x</mi><mi>K</mi></msub><mo>)</mo></mrow></mrow></mtd></mtr></mtable></mfenced><mrow><mi>K</mi><mo>&times;</mo><mi>L</mi></mrow></msub><mo>;</mo></mrow>]]></math><img file="FDA0000784832690000021.GIF" wi="1225" he="391" /></maths>(5d)计算输出权值为<img file="FDA0000784832690000022.GIF" wi="218" he="76" />其中<img file="FDA0000784832690000023.GIF" wi="70" he="69" />为H的Moore‑Penrose广义逆,<img file="FDA0000784832690000024.GIF" wi="392" he="78" />T=[t<sub>1</sub>,....,t<sub>K</sub>]<sup>T</sup>;(6)数据集<img file="FDA0000784832690000026.GIF" wi="40" he="79" />输入到训练好的极限学习机中,得到分类结果。
地址 710071 陕西省西安市太白南路2号