[发明专利]一种植被参数反演方法、终端设备及存储介质

权利要求书

1.一种植被参数反演方法，其特征在于，包括以下步骤：

步骤1，对2幅SAR影像进行配准、去平地、多视以及极化干涉处理，生成单基线全极化复相干系数γ_(ω)；

步骤2，设置植被参数反演的细分区块大小为N行M列，获取区块每个像素点的复相干系数，与RVoG模型进行匹配，得到观测方程为：

其中，ω表示与基线和极化方式相关的极化状态；表示地表相位，是有关地表高程h_g的函数；μ_(ω)表示极化状态为ω所对应的地体幅度比；γ_v表示纯体去相干系数，是有关植被高度h_v和消光系数σ的函数；k_z为垂直向有效波数；θ为微波入射角；

针对每个像素的多种极化方式，均可得到至少2个复数观测方程，则区块NM个像素点可得到至少2NM个复数观测方程，即4NM个实数观测方程；

步骤3，计算区块内每个像素点的参数初值：将地体幅度比μ_(ω)的归一化参数x的初值设为0，然后利用RVoG模型和三阶段植被高度反演算法计算植被高度初值、消光系数初值和地表高程初值；

其中，归一化参数x是指：在复平面内，体散射占优通道的极化复相干系数到纯体去相干系数的距离；

步骤4，根据步骤2得到的4NM个实数观测方程以及步骤3得到的各初值，采用非线性迭代算法求解观测方程中的2NM+2个未知参数：植被高度、消光系数、地表高程、地体幅度比归一化参数；

其中，区块内NM个像素点所对应的植被高度和消光系数均恒定，地表高程和地体幅度比的归一化参数x均随像素点的变化而变化。

2.根据权利要求1所述的方法，其特征在于，利用RVoG模型和三阶段植被高度反演算法计算植被高度初值的方法为：

步骤1)，直线拟合：

对于区块内每个像素点，均利用在该基线2种极化状态下的复相干系数观测值，对RVoG模型在复平面内进行直线拟合，得到公式(3)所示直线：

步骤2)，计算地表相位和地表高程初值：

对于区块内每个像素点，计算公式(3)所示直线与复平面单位圆的2个交点A和B，将其中与体散射占优极化方式点之间的距离大于与地面散射占优极化方式点之间的距离的交点，即为地表相位点然后利用地表相位按以下公式(4)计算地表高程初值：

3.根据权利要求1所述的方法，其特征在于，利用RVoG模型和三阶段植被高度反演算法计算消光系数初值、地表高程初值的方法为：

选择区块的中心像素点计算整个区块的植被高度和消光系数初始值：

设所选择像素点的体散射占优极化通道仅包含植被层散射贡献，其地体幅度比为0，则纯体相干性的表达式为表示体散射占优通道复相干系数，是有关植被高度h_v以及消光系数σ的函数，表示为

采用查找表的方法，通过给定一系列植被高度h_v以及消光系数σ，根据公式(5)构建二维查找表，分别计算然后在表格中找到与相差最小的一组，作为植被高度初值和消光系数初值：

4.根据权利要求1所述的方法，其特征在于，区块的行数N和列数M均为大于1的奇数。

5.一种终端设备，包括存储器、处理器以及存储在所述存储器中并可在所述处理器上运行的计算机程序，其特征在于，所述处理器执行所述计算机程序时实现如权利要求1至4任一项所述的方法。

6.一种计算机可读存储介质，所述计算机可读存储介质存储有计算机程序，其特征在于，所述计算机程序被处理器执行时实现如权利要求1至4任一项所述的方法。