[发明专利]一种基于分布式压缩感知的星载差分层析SAR成像方法

权利要求书

1.一种基于分布式压缩感知的星载差分层析SAR成像方法，其特征在于，包括以下步骤：

步骤1、预处理SAR图像，得到差分层析SAR数据堆栈；

步骤2、基于KS检测理论，识别差分层析SAR数据堆栈中相邻同质像元，判定SAR图像像元内是否可能存在分布式散射体和强散射体；

步骤3、假设模型散射体数目为1时，估计SAR图像像元内能量最强散射体的残余高程与形变速率；

步骤4、构建有无强散射体假设检验问题模型，逐像元判别SAR像元内有无强散射体；

步骤5、对于步骤4中判别为无强散射体的SAR图像像元且步骤2中判别为可能存在分布式散射体的SAR图像像元，对SAR图像像元内分布式散射体的高程与分布式散射体的形变速率信息进行估计；

步骤6、对于步骤4中判别为有强散射体的SAR图像像元且步骤2中判别为可能存在强散射体的SAR图像像元，判别SAR图像像元内有单个还是双个散射体；

步骤7、对于步骤6中有强散射体的SAR图像像元，估计SAR图像像元内强散射体的高程与形变速率；

步骤8、根据步骤5估计的SAR图像像元内分布式散射体的高程、分布式散射体的形变速率以及散射系数，还根据步骤7中估计的SAR图像像元内强散射体的高程、形变速率以及散射系数，逐个像元判定SAR图像像元内分布式散射体和强散射体的有效数目，并获得有效的分布式散射体和强散射体的高程与形变速率信息，实现星载差分层析SAR成像。

2.根据权利要求1所述的一种基于分布式压缩感知的星载差分层析SAR成像方法，其特征在于，所述的步骤1包括以下步骤：对所有SAR图像进行幅度及相位校正，然后基于地理编码原理，对每幅SAR图像进行去斜操作，得到差分层析SAR数据堆栈。

3.根据权利要求1所述的一种基于分布式压缩感知的星载差分层析SAR成像方法，其特征在于，所述的步骤2包括以下步骤；

步骤2.1、采用两个像素统计所得的概率分布函数差的绝对值D_M来度量任意SAR图像像元P与相邻SAR图像像元P'的相似性；

步骤2.2、设定阈值c与置信度α，概率分布函数差的绝对值D_M的概率分布函数为H(t)，若D_M≤c且α≥1-H(t)时，则两个SAR图像像元识别为相邻同质像元；否则，两个SAR图像像元不为相邻同质像元；

步骤2.3、设置第一邻域窗口和第二邻域窗口，第二邻域窗口尺寸小于第一邻域窗口尺寸；若在第一邻域窗口内与SAR图像像元相邻的相似同质像元数目大于第一设定阈值，则判定该SAR图像像元内可能存在分布式散射体；若在第二邻域窗口与SAR图像像元相邻的相似同质像元数目大于第二设定阈值，则判定该SAR图像像元内可能存在强散射体，且该强散射体像元具有多个相似同质像元。

4.根据权利要求1所述的一种基于分布式压缩感知的星载差分层析SAR成像方法，其特征在于，所述的步骤4包括以下步骤；

构建有无强散射体假设检验问题模型：

其中,M为SAR图像的总数,λ为雷达波长，r为天线相位中心到参考地形的斜距，θ为入射角,b_⊥m与t_m分别为第m次观测中相对主图像的有效基线与时间差，与分别表示SAR像元内有强散射体的假设、SAR像元内无强散射体的假设，和分别表示在有强散射体的假设和无散射体假设条件下的概率密度函数，为差分层析SAR数据堆栈中任意同名像元的平均亮度，SAR时序观测值构成的矢量x＝[x₁,…,x_M]^T，T为转置符，ε表示热噪声标准差，为SAR图像像元内强散射体的高程，为SAR图像像元内强散射体的形变速率；

若成立，则SAR图像像元被判定为无强散射体；否则，SAR图像像元被判定为有强散射体，k为SAR图像像元内强散射体的个数，P₀与P_k表示SAR图像像元内没有强散射体的先验概率与有k个强散射体的先验概率。