黎 晓，李增元，陈尔学，程甲州，姜友谊

(1.中国林业科学研究院资源信息研究所，北京 100091；2.西安科技大学测绘科学与技术学院，陕西 西安 710054)

反演森林树木高度可以更好地了解森林在全球环境变化中的作用[1-3]。现在，遥感是在区域或全球尺度上进行森林高度反演的最有效方法[4]。极化干涉SAR（PolInSAR）不但拥有微波遥感全天时、全天候的特点，也综合了极化和干涉技术对散射体高度、形状等信息敏感的优点，利用极化干涉SAR 技术反演森林树高已经成为林业研究的热点之一[5-6]。

近十年来，RVoG 模型一直是树高反演和生物量估计中最常用的模型[7]。2001 年，Cloude 利用六维非线性优化方法成功反演出植被高度[8]，但该方法受初值影响较大且难以求得全局最优解；2003年，Clode 基于RVoG (Random Volume over Ground)模型的相干系数在复平面单位圆分布的特点提出了三阶段算法[9]；但是，三阶段算法中体散射复相干估计不准确会影响反演精度[10-11]。而对于重复观测或星载极化SAR 系统，时间去相关成为主要的误差源[12]。Mette 等[13]在研究RVoG 模型反演植被高度的3 种误差源时，发现时间去相关会导致反演结果产生巨大误差。因此，研究时间去相关问题需要改进RVoG 模型或提出新的反演模型。Lee 等利用L 和P 波段SAR 数据研究不同时间间隔的时间去相关时发现，时间去相关不仅降低了相干系数，并且增大了植被区相干相位的波动性[14-16]。2003 年，Papathanassiou 和 Cloude[17]提 出 RVoG+VTD(Random Volume over Ground with Volume Temporal Decorrelation)模型用于解决时间去相关问题，但该模型对反演精度的提升有限。周勇胜等[18]在RVoG+VTD 模型的基础上，利用双基线极化干涉SAR 数据解决消光系数与植被高度模糊的问题，但只使用模拟数据验证，缺少说服力。2010 年，Lavalle 提出新的时间去相关模型——TD+RVoG(Temporal Decorrelation Random Volume over Ground)模型[19-21]，该模型将时间去相关的原因归结于植被层和地面层的随机高斯运动，但该模型未知数过多，在实际中难以应用。随后，Lavalle 提出RMoG(Random Motion over Ground)模型，该模型在TD+RVoG 模型的基础上继续发展并利用十维非线性优化方法求解[22]，然而这种方法和六维非线性优化方法有同样的缺点且耗时过长；2018 年，Nafiseh 改进了RMoG 模型，提出了RMoGL模型，该模型改用傅里叶-勒让德级数作为结构函数，但该模型的迭代过程比较繁琐、耗时过长，其反演时间大概为RMoG 模型的两倍[7]。

可以看出，消除时间去相关的办法是在反演模型中加入植被运动，然而这样就一定要引入更多的未知数，从而使问题变得复杂难解；为了改善RMoG 模型反演困难且耗时过长的问题，本研究提出简化RMoG 模型，在有效弥补时间去相关的前提下，提高反演效率。为验证该方法，利用欧空局提供的BioSAR2007P 波段数据进行试验分析，结果表明，该方法简单高效。

1 研究区及数据介绍

1.1 研究区概况

研究区位于瑞典南部的Remingstorp 地区（58°28′40″ N, 13°37′25″ E），该地区地形相对平坦，总面积超过1 500 hm2，其中大约1 200 hm2为森林，其余为湖泊，海拔高度为120～145 m，树高在10 m 到30 m 之间，主要树种为挪威云杉（Picea abies (L.) Karst.），欧洲赤松（Pinus sylvestris Linn.）和桦树（Betula spp.），主要土壤类型为耕地土壤。

1.2 机载SAR 数据获取与处理

本研究所用的SAR 数据是德国宇航局（DLR）与瑞典国防研究局（FOI）利用机载E-SAR 系统于2007 年获取的P 波段全极化干涉数据，其水平基线为40 m，波长为0.86 m，空间基线为20 m，时间基线为50 分，距离向分辨率为0.75 m，方位向分辨率1.5 m，数据获取的具体时间为2007 年4 月2 日。数据的主辅影像已经过精确配准，相应的地平数据和有效波数也由DLR 提供。为了提高相干性和减弱噪声的影响，数据在预处理时使用了7×7 的Boxcar 滤波处理。图1 为该数据Pauli 基下RGB 彩色合成图。

图1 Pauli 基合成图Fig.1 Image in Pauli basis

1.3 Lidar 数据获取与处理

Lidar 数据于2010 年获取，其中部分区域在2009 年已进行测量。Lidar 所测的归一化数字表面模型(NDSM)，即去除地形信息的DSM，在林区相当于冠层高度模型(CHM)。CHM 为树高反演结果提供验证，其空间分辨率为0.5 m×0.5 m，如图2 所示。

2 研究方法

2.1 RMoG 模型

RMoG 模型假设植被层由随机分布在粗糙介质层上的垂直物体组成，并且植被冠层和底层都有沿垂直轴的随机运动。在这个假设的基础上，RMoG 模型的复相干表述为：

图2 Lidar CHM 图Fig.2 Lidar CHM

其中σ2（Z）是沿Z 轴的一阶近似运动方差，被定义为：

将上述式子联合并求解积分，便可得到RMoG模型的最终形式：

其中

其中hv 为植被高度，σ 为平均消光系数，θ 为雷达入射角，φg为地面相位，μ 为地体散射幅度比，θs为沿着距离方向的地面坡度角。

从RMoG 模型的公式中可以看出，要求解此方程至少需要5 种不同的极化通道；又因为μ 随着极化通道的改变而改变，所以极化通道的增加会使未知数的个数增加，当未知数个数大于极化通道个数时，此问题就转化为最优化问题。因此反演时，不可避免的存在长计算时间和对初值敏感的问题。

RMoG 模型在RVoG 模型的基础上加入了植被高斯运动来减小时间去相关的影响，但由于未知数的增加，使反演问题复杂，反演速度变慢，难以在实际中得到良好的应用。为此，本研究对RMoG 模型进行简化；为求解此模型，在借鉴三阶段算法的基础上，将其归结于3 步。

2.2 简化RMoG 模型

森林等植被区的时间去相关按照原因可分为3 类：一是植被运动造成的时间去相关，如风的影响；二是天气变化引起的介电常数的变化，如降水；三是人为或自然原因引起的场景改变造成的时间去相关。现在的研究大都集中在植被运动引起的时间去相关上[23]。由RMoG 模型可知，减弱时间去相关影响的办法为加入植被冠层和底层运动，以此来反映植被运动引起的时间去相关。根据文献[7]，为了简化RMoG 模型，在数据采集之间的时间间隔较短的情况下，可以将地面运动σg假设为0。当此假设成立时，则γtg=1，于是可以简化为：

而RVoG 模型中体散射公式为：

与之相比，γvt表达式中增加了植被冠层运动σv；对此，在进行查找表时，可以将消光系数固定，从而建立植被高度hv 和冠层运动σv的查找表；而地面相位φg可以从三阶段算法中得到。因此，简化RMoG 模型也可以分为3 步求解，其步骤如下：

1)在复平面上对多个极化干涉相干系数进行最小二乘直线拟合，该直线与单位圆有两个交点。

2)地面相位与纯体散射去相干的确定。理论上应找到地体散射比μ 为0 的点对应的相干系数作为纯体散射去相干的估计，但在实际应用中一般选择HV 极化通道；由于PD 极化相干在体散射去相干上的良好应用[24-25]，本研究选择了PDH极化方式作为纯体散射去相干。计算两交点与PDH、PDL的距离，选取距离PDH比距离PDL远的点作为地面相位点。

3)采取合适的消光系数并根据公式10 建立查找表，反演最终得到植被高度。

该方法也可以理解为在三阶段算法中加入了植被冠层运动，以此来抵消时间去相关的影响。

3 试验结果与分析

3.1 高度反演结果

分别用本研究方法和传统三阶段方法对研究区数据进行反演，得到的植被高度图分别如图3 和图4 所示。RVoG 模型中将时间去相关解释为体散射去相关，因此高估了树高；图4 与图2 对比，可以明显发现图像右边结果高估严重。而本研究方法考虑时间去相关的影响，在模型中加入植被冠层运动，所以对高估问题有显著的改善，并且植被区的结果与Lidar 生成的CHM 结果较好地吻合。此外，平地效应去除，干涉相干估计等预处理步骤均由PolSARpro 软件处理。

3.2 精度对比分析

为定量评价本研究算法的优劣，使用均方根误差（RMSE）和决定系数（R2）为定量评价指标。由于本研究和文献[7]使用的数据在同一地区，所以在比较精度的时候，直接利用了文献[7]中P 波段RMoG 模型和RMoGL模型的精度，如表1 所示。虽然研究区内大部分区域为植被区，但仍有部分湖泊存在；其主要植被区已由Lidar 扫描，所以在讨论精度时，只对CHM 内像素进行讨论，而对其他区域忽略。

图3 本研究方法Fig.3 Method in this paper

图4 三阶段方法Fig.4 Three-stage inversion

从表1 可以看出：在均方根误差方面，三阶段算法最差，其次是本研究方法，最好的是RMoGL模型；而在相关性方面，本研究方法和三阶段方法表现较好，最差的为RMoG 模型。通过定量的对比可知，本研究方法在精度上有一定的优势；且由于该方法避免了非线性优化过程，可以大大地提高反演速率。

4 讨论

在利用PolInSAR 数据估计森林参数时，时间去相关是反演模型的关键问题。从RMoG 模型出发，继承了其基于高斯运动的时间相干性，并通过忽略地面运动将其简化，称为简化RMoG 模型。为避免出现未知数过多现象，采用固定消光系数的方法，并结合查找表使其求解过程更加简单；在有效消除时间去相关影响的前提下，大大缩减反演时间。

Askne 等在研究时间去相关时，按照时间基线的长度将时间去相关分为三类：短时间基线内，风速等是产生时间去相关的主要影响因素；中时间基线内，植被区介电常数的改变是产生时间去相关的主要原因；而在长时间基线内，时间去相关可能是植被生长或人为破坏造成[26-28]。本研究的时间去相关主要集中在植被运动也就是风速上，并没有对介电常数或其他进行讨论；换言之，本研究方法只适用于短时间基线或机载SAR 数据。

由于消光系数并不是树高反演时的敏感参数，因此本研究采用固定消光系数的方法。但已有研究的实验都是建立树种类型差异不大或树冠形态差异不大的前提下，本研究区满足此条件。然而当多树种或树冠差异较大时，如何固定消光系数，需要更多的实验研究。

5 结论

为更好地从PolInSAR 数据估计出森林参数，时间去相关的影响成为不可不解决的问题。在时间去相关方面，RMoG 模型和RMoGL模型都有着较好的反演精度，但由于其难以求解、耗时较长的原因并不能很好地应用。本研究将RMoG 模型简化，忽略地面运动，在保证反演精度的同时，提高了反演效率。通过BioSAR2007 数据的试验，可以得到以下结论：与未加入时间去相关的RVoG 模型对比，简化RMoG 模型只加入植被冠层运动抵消时间去相关影响的方法行之有效，并且反演结果明显改善，得到的高度图更接近Lidar 生成的CHM，从中可以看出考虑时间去相关的重要性和必要性；在精度对比方面，与三阶段算法相比，本研究方法的精度有所提升，均方根误差由8.52 降到6.24；与RMoG 模型和RMoGL模型相比，本研究方法虽然略有不足，但由于忽略植被地面运动，很好地改善了RMoG 模型和RMoGL模型反演困难、耗时过长的缺点。本研究方法可以有效消除时间去相关的影响，且使反演过程更加高效。