发明名称 极化SAR图像舰船目标检测方法
摘要 本发明公开了一种极化SAR图像目标检测方法,主要解决现有检测方法对能量的依赖和虚警率较高的问题。其实现步骤是:选取目标和杂波两类训练样本,并进行特征提取和能量归一;定义极化字典学目标函数,监督学一个对杂波和目标具有区分能力的极化字典;对测试图像进行极化特征提取和能量归一;采用极化字典对测试特征向量进行稀疏编码;将特征向量在杂波子字典下的重构误差比例作为检验统计量;比较检验统计量与阈值的大小,将大于阈值的像素判为目标,否则判为杂波,完成目标检测。本发明仅利用极化信息进行目标检测,具有稳健性好,虚警率低,能对弱小目标进行检测的优点,适用于目标和杂波对比度较低时的极化SAR图像目标检测。
申请公布号 CN104091335B 申请公布日期 2017.01.18
申请号 CN201410320187.3 申请日期 2014.07.04
申请人 西安电子科技大学 发明人 刘宏伟;文伟;王英华;陈渤
分类号 G06T7/00(2006.01)I;G06K9/66(2006.01)I 主分类号 G06T7/00(2006.01)I
代理机构 陕西电子工业专利中心 61205 代理人 王品华;朱红星
主权项 一种极化SAR图像的舰船目标检测方法,包括:A.训练步骤:(A1)从极化SAR图像场景中,选取M个杂波像素作为杂波训练样本集合,对集合中像素i=1,...,M,采用像素i的3×3邻域窗口内像素计算杂波极化协方差矩阵C<sub>i</sub>;(A2)根据杂波极化协方差矩阵C<sub>i</sub>,计算杂波极化特征向量<img file="FDA0001030453670000013.GIF" wi="59" he="50" /><img file="FDA0001030453670000014.GIF" wi="1797" he="62" />其中(·)<sup>T</sup>表示矩阵转置;式中C<sub>i</sub>(a,b)表示极化协方差矩阵C<sub>i</sub>第a行第b列的元素,a≤b,a=1,2,3,b=1,2,3;R(·),I(·)分别表示取实部和虚部算子;(A3)对杂波极化特征向量<img file="FDA0001030453670000012.GIF" wi="41" he="57" />进行能量归一,得到归一化的杂波极化特征向量:<img file="FDA0001030453670000011.GIF" wi="253" he="76" />其中||·||<sub>2</sub>表示2范数算子;(A4)将归一化杂波特征向量x<sub>i</sub>拼接成杂波特征矩阵X<sub>1</sub>=[x<sub>1</sub>,...,x<sub>i</sub>,...,x<sub>M</sub>];(A5)从极化SAR图像场景中,选取N个目标像素作为目标训练样本集合,对集合中像素j=1,...,N,采用与步骤(A1),(A2),(A3)相同的操作,得到目标像素j归一化极化特征向量x<sub>j</sub>,将归一化目标极化特征向量拼接成目标特征矩阵X<sub>2</sub>=[x<sub>1</sub>,...,x<sub>j</sub>,...,x<sub>N</sub>];(A6)将杂波特征矩阵X<sub>1</sub>与目标特征矩阵X<sub>2</sub>进行拼接,得到特征矩阵:X=[X<sub>1</sub>,X<sub>2</sub>];(A7)定义杂波编码子字典D<sub>1</sub>=[d<sub>1</sub>,...,d<sub>p</sub>,...,d<sub>K</sub>],字典的每一列d<sub>p</sub>称为一个原子其中p=1,...,K,K为原子个数;定义目标编码子字典D<sub>2</sub>=[d<sub>1</sub>,...,d<sub>q</sub>,...,d<sub>K</sub>],q=1,...,K;将杂波编码字典与目标编码字典进行拼接,得到编码字典D=[D<sub>1</sub>,D<sub>2</sub>];(A8)利用特征矩阵X,采用稀疏编码与字典更新交替迭代的方法对编码字典D进行学习,得到输出优化后的编码字典D<sup>*</sup>=D<sup>L+1</sup>,L为字典更新迭代次数:B.测试步骤:(B1)对于一组待检测全极化SAR图像I,定义一个与水平‐水平极化图像I<sub>HH</sub>大小相同 的指示图像B,用以对检测结果中目标像素的位置进行标记;(B2)对待检测极化SAR图像I中的像素t采用与训练阶段步骤(A1)相同的操作,计算测试极化协方差矩阵C<sub>t</sub>;(B3)根据测试极化协方差矩阵C<sub>t</sub>,进行与训练阶段步骤(A2)相同的操作,提取测试极化特征向量<img file="FDA0001030453670000024.GIF" wi="59" he="52" />(B4)对测试极化特征向量<img file="FDA0001030453670000025.GIF" wi="38" he="52" />进行与训练阶段步骤(A3)相同的操作,获得归一化测试极化特征向量x<sub>t</sub>;(B5)假设测试像素点t为杂波,根据训练阶段步骤(A8)学习获得的编码字典D<sup>*</sup>,对测试极化特征向量x<sub>t</sub>进行稀疏编码,获得测试极化特征向量x<sub>t</sub>的稀疏编码系数<img file="FDA0001030453670000021.GIF" wi="218" he="149" />其中z<sub>t1</sub>,z<sub>t2</sub>分别为测试极化特征向量x<sub>t</sub>在杂波编码字典D<sub>1</sub>和目标编码字典D<sub>2</sub>的编码系数;(B6)根据稀疏编码系数z<sub>t</sub>,构建检验统计量:<img file="FDA0001030453670000022.GIF" wi="390" he="142" />(B7)根据虚警率设置检测阈值T=0.42,将检验统计量<img file="FDA0001030453670000026.GIF" wi="98" he="55" />与该阈值进行比较,若大于检测阈值T则判别为目标,在指示图像B中将像素t所在位置标记为1,否则标记为0;(B8)对待检测图像I的每一个像素,都进行与步骤(B2)~(B7)相同的操作,完成对指示图像B的赋值,指示图像B即为对应于待检测图像I的检测结果。
地址 710071 陕西省西安市太白南路2号