[发明专利]基于压缩多信号分类的层析合成孔径雷达测量树高的方法

权利要求书

1.一种基于压缩多信号分类的层析合成孔径雷达测量树高的方法，其特征在于，包括：

步骤A：对观测区域在不同基线位置发射脉冲信号进行重复观测，接收观测区域的反射回波信号Y_m，其中，m＝1、2、……，M，M为航过数；

步骤B：根据接收到的反射回波信号Y_m，提取其HH通道、HV通道和VV通道的像素点数据，全部反射回波信号对应的像素点数据构成方位-距离分辨单元的多通道像素点矩阵Y；

步骤C：根据SAR层析成像几何模型得到观测矩阵A；

步骤D：根据回波信号构建的像素点矩阵Y以及观测矩阵A，采用正则化算法求出k-3个场景支撑集的索引值，得到支撑集I_k-3，其中，支撑集索引值是指场景中散射目标点在高程向出现的位置，k为场景稀疏度，即高程方向上散射目标点的个数；

步骤E：利用经典谱估计算法求出所述k-3个支撑集索引值之外的其余三个支撑集元素，将其并入支撑集中，得到完整支撑集I_k；以及

步骤F，由完整支撑集I_k中的索引值获知地表和树冠的位置，由两者的位置差获得树高信息，其中，地表和树冠分别对应支撑集I_k中的索引值集中的两个高度位置。

2.根据权利要求1所述的方法，其特征在于，所述步骤D中，所述正则化算法为同步的正交基追踪算法或IST算法。

3.根据权利要求2所述的方法，其特征在于，采用同步的正交基追踪算法求出k-3个场景支撑集的索引值包括：

子步骤D1：初始化像素点矩阵Y₀和存放索引值的支撑集S₀，令Y₀＝Y，S₀＝φ，φ为空集，设置当前迭代次数j＝1；

子步骤D2：在第j次迭代中，寻找满足与残差信号Y_j乘积的2-范数最大时，所对应的观测矩阵A的列序号l_j，即为本次迭代所要求得的支撑集索引值：

||alj*Yj||2=max1≤l≤N||al*Yj||2]]>

其中，为观测矩阵A的第l_j列，Y_j为第j次迭代时的像素点矩阵的残差信号：为S_j中索引值对应的观测矩阵A的相应列构成的支撑矩阵Asupp=al1al2······alj-1]]>的正交基；

子步骤D3：由支撑集的索引值l_j更新支撑集：

S_j＝S_j-1∪{l_j}

子步骤D4：判断若j＞k-3，则终止迭代，并令I_k-3＝S_j，像素点矩阵的残差信号为Y_k-3，执行步骤E，否则，j递增1，继续子步骤D2。