【技术实现步骤摘要】
本专利技术涉及雷达图像处理领域,尤其涉及一种圆迹合成孔径雷达(SAR)图像序列生成的方法。
技术介绍
圆迹SAR是一种新型遥感信息获取和处理技术,特指成像系统通过360°圆周轨迹来实现多视向大视角范围连续成像观测。在圆迹SAR成像过程中,常按方位角间隔分成若干个子孔径进行相干处理,最后再累加获得全孔径图像。若目标对所有孔径都满足方位向同散射中心假设,非相干累加处理可以获得信噪比最高的图像。然而实际地物的散射特征仅在有限的方位角内保持稳定,这对于基于极化特征参数的遥感应用,例如变化检测和地物分类,存在参数估计不准的问题。传统的窄方位角成像模式可以忽略方位向多散射中心,因为它的多视向散射特征不足以提供有价值的信息。多视向观测模式则需要进一步拓展现有方法,在保留方位向多散射中心贡献的基础上挖掘消除过程产生的信息。
技术实现思路
本专利技术提供了一种圆迹SAR图像序列生成方法,包括:步骤S1,将圆迹SAR360°全孔径平均划分成R个子孔径,子孔径序号为i,1≤i≤R,R≥k,k为允许的最小子孔径数量,R和k均为大于1的整数,并得到各个子孔径图像每点像素的极化协方差矩阵Ci; ...
【技术保护点】
一种圆迹合成孔径雷达SAR图像序列生成方法,包括:步骤S1,将圆迹SAR360°全孔径平均划分成R个子孔径,子孔径序号为i,1≤i≤R,R≥k,k为允许的最小子孔径数量,R和k均为大于1的整数,并得到各个子孔径图像每点像素的极化协方差矩阵Ci;步骤S2,如果R=k则执行步骤S5,否则执行步骤S3;步骤S3,对每个像素位置计算得到R个子孔径图像3×3统计窗口内的差异性似然比,对所述差异性似然比进行排序得到最大似然比以及与最大似然比相对应的子孔径的序号subk,作为方位向散射中心差异最大的差异最大子孔径subk;步骤S4,计算与所述差异最大子孔径subk对应的恒虚警似然比,并得 ...
【技术特征摘要】
1.一种圆迹合成孔径雷达SAR图像序列生成方法,包括:步骤S1,将圆迹SAR360°全孔径平均划分成R个子孔径,子孔径序号为i,1≤i≤R,R≥k,k为允许的最小子孔径数量,R和k均为大于1的整数,并得到各个子孔径图像每点像素的极化协方差矩阵Ci;步骤S2,如果R=k则执行步骤S5,否则执行步骤S3;步骤S3,对每个像素位置计算得到R个子孔径图像3×3统计窗口内的差异性似然比,对所述差异性似然比进行排序得到最大似然比以及与最大似然比相对应的子孔径的序号subk,作为方位向散射中心差异最大的差异最大子孔径subk;步骤S4,计算与所述差异最大子孔径subk对应的恒虚警似然比,并得到恒虚警概率,将所述恒虚警概率与门限虚警概率进行比较,若恒虚警概率大于门限虚警概率,消除所述差异最大子孔径在该像素点的散射贡献,令R=R-1,并返回步骤S2,否则执行步骤S5;步骤S5,针对在该像素点具有散射贡献的每个子孔径进行协方差矩阵的非相干累加,得到结果图像在该像素点的极化协方差矩阵,从而得到结果图像。2.根据权利要求1所述的方法,其中,k等于4。3.根据权利要求1所述的方法,其中,所述步骤S3包括:分别计算每个子孔径图像的等效视数;根据计算的等效视数分别针对每个子孔径计算差异性似然比;对计算的差异性似然比排序得到最大差异性似然比以及对应的子孔径序号subk。4.根据权利要求3所述的方法,其中,通过以下公式计算每个子孔径图像的等效视数ENL:ENL=1/τ2其中,τ为对应子孔径图像协方差矩阵的3×3统计窗口内的标准差和均值的比值。5.根据权利要求4所述的方法,其中,通过以下公式计算每个子孔径图像的差异性似然比ΛR: Λ R = Π i R | C i | ENL i | ( Σ i R C i ) / R | RENL T ]]>其中ENLi代表第i个子孔径的等效视数,ENLT代表的等效视数。6.根据权利要求5所述的方法,其中,通过以下方式计算与所述差异最大子孔径subk对应的恒虚警似然比:根据等效视数ENL,得到恒虚警似然比为ΛAB: Λ A B = | C A | N A ENL A | C B | N B ENL B | C T | N T ENL T , ]]>其中,NT=NA+NB,CT=(NACA+NBCB)/NT,NA和NB分别代表A类和B类的子孔径数量,CA和CB分别代表A类和B类的极化相干矩阵样本均值,ENL...
【专利技术属性】
技术研发人员:洪文,李洋,薛斐腾,郭小洋,王建峰,郭胜龙,
申请(专利权)人:中国科学院电子学研究所,
类型:发明
国别省市:北京;11
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。