发明名称 基于K‑WISHART分布的极化SAR舰船检测方法
摘要 本发明公开了一种基于K‑Wishart分布的极化SAR舰船检测方法,思路为:依次获取极化SAR雷达回波样本和极化SAR雷达回波样本的协方差矩阵C,并分别设定T为C的分类个数,依次计算T类极化SAR雷达回波样本的协方差矩阵C中n个样本的后验概率和第m类极化SAR雷达回波样本的协方差矩阵C<sub>m</sub>,进而计算F个通道的极化SAR雷达回波样本对应的实际水平值,若所述实际水平值分别大于或等于设定的显著性水平,则C<sub>m</sub>的分类结果正确,然后计算T类极化SAR雷达回波样本的协方差矩阵各自的SPAN均值,获取T个SPAN均值中的最大SPAN均值,并将所述最大SPAN均值所属类为舰船,其余T-1个类分别为海平面。
申请公布号 CN106291554A 申请公布日期 2017.01.04
申请号 CN201610616431.X 申请日期 2016.07.29
申请人 西安电子科技大学 发明人 周峰;圣铭;樊伟伟;陶明亮
分类号 G01S13/90(2006.01)I;G01S7/41(2006.01)I 主分类号 G01S13/90(2006.01)I
代理机构 西安睿通知识产权代理事务所(特殊普通合伙) 61218 代理人 惠文轩
主权项 一种基于K‑Wishart分布的极化SAR舰船检测方法,其特征在于,包括以下步骤:步骤1,获取极化SAR雷达回波样本,并根据极化SAR雷达回波样本获取极化SAR雷达回波样本的协方差矩阵C,并分别设定T为极化SAR雷达回波样本的协方差矩阵C的分类个数,设定n为极化SAR雷达回波样本的协方差矩阵C的样本个数,且T的初值为1;步骤2,对极化SAR雷达回波样本的协方差矩阵C进行分类,计算T类极化SAR雷达回波样本的协方差矩阵C中n个样本的后验概率,其中T类极化SAR雷达回波样本的协方差矩阵C中第j'个样本的后验概率γ<sub>j'</sub>={γ<sub>j'1</sub>,γ<sub>j'2</sub>,...,γ<sub>j'm</sub>,...,γ<sub>j'T</sub>},γ<sub>j'm</sub>表示第m类极化SAR雷达回波样本的协方差矩阵C中第j'个样本的后验概率;如果第m类极化SAR雷达回波样本的协方差矩阵C中第j'个样本的后验概率γ<sub>j'm</sub>最大,则将第j'个样本归为第m类极化SAR雷达回波样本的协方差矩阵C<sub>m</sub>中的样本;j'∈[1,n],进而得到第m类极化SAR雷达回波样本的协方差矩阵C<sub>m</sub>中的所有样本,m∈[1,T];步骤3,对第m类极化SAR雷达回波样本的协方差矩阵C<sub>m</sub>进行参数估计,分别计算第m类极化SAR雷达回波样本的均值协方差矩阵Σ<sub>m</sub>、第m类极化SAR雷达回波样本的协方差矩阵的多视数估计值<img file="FDA0001064536960000011.GIF" wi="83" he="70" />T类极化SAR雷达回波样本的协方差矩阵的全局有效多视数<img file="FDA0001064536960000012.GIF" wi="35" he="62" />和第m类极化SAR雷达回波样本的协方差矩阵的形状参数估计值<img file="FDA0001064536960000013.GIF" wi="83" he="63" />进而计算得到基于复K‑Wishart分布的第x个通道的第m类极化SAR雷达回波样本的协方差矩阵的概率密度函数;步骤4,根据基于复K‑Wishart分布的第x个通道的第m类极化SAR雷达回波样本的协方差矩阵的概率密度函数,对第m类极化SAR雷达回波样本的协方差矩阵C<sub>m</sub>进行适应性检测,得到F个通道的SAR雷达回波样本对应的实际水平值<img file="FDA0001064536960000014.GIF" wi="67" he="69" /><img file="FDA0001064536960000015.GIF" wi="534" he="71" /><img file="FDA0001064536960000016.GIF" wi="51" he="70" />表示第x个通道的SAR雷达回波样本的实际显著性水平值,x∈{1,2,…,F},F表示获取的极化SAR雷达回波样本包含的通道个数;如果所述F个通道的SAR雷达回波样本对应的实际显著性水平值<img file="FDA0001064536960000017.GIF" wi="42" he="68" />分别大于或等于设定的显著性水平β,则所述第m类极化SAR雷达回波样本的协方差矩阵C<sub>m</sub>的分类结果正确,然后执行步骤5;如果所述F个通道的SAR雷达回波样本对应的实际显著性水平值<img file="FDA0001064536960000021.GIF" wi="42" he="68" />中任意一个通道的SAR雷达回波的实际显著性水平值小于设定的显著性水平β,则所述第m类极化SAR雷达回波样本的协方差矩阵C<sub>m</sub>的分类结果不正确,并令T加1,返回执行步骤2;步骤5,计算第m类极化SAR雷达回波样本的协方差矩阵C<sub>m</sub>的对角线元素之和,并将所述第m类极化SAR雷达回波样本的协方差矩阵C<sub>m</sub>的对角线元素之和,作为第m类极化SAR雷达回波样本的协方差矩阵的SPAN值;步骤6,令m加1,重复执行步骤5,直到得到第T类极化SAR雷达回波样本的协方差矩阵的SPAN值,然后根据此时得到的第1类极化SAR雷达回波样本的协方差矩阵的SPAN值到第T类极化SAR雷达回波样本的协方差矩阵的SPAN值,分别计算第1类极化SAR雷达回波样本的协方差矩阵的SPAN均值到第T类极化SAR雷达回波样本的协方差矩阵的SPAN均值,获取T个SPAN均值中的最大SPAN均值,并将所述最大SPAN均值所属类作为舰船,其余T‑1个类分别作为海平面。
地址 710071 陕西省西安市太白南路2号