
Citation: | |
森林高度是估算森林蓄积量及生物量的重要基础数据,对于研究森林资源状况以及分析全球生态环境、气候变化具有重要意义。极化合成孔径雷达干涉测量技术(Polarimetric SAR Interferometry, PolInSAR)采用微波监测模式,其回波信号不仅记录垂直结构及其属性信息,且可以区分同一分辨单元内不同散射体高度的能力,已被视为大范围、高分辨率、高精度反演森林高度的有效手段之一[1]。
为了实现利用PolInSAR观测量准确地提取森林高度,Thrauhft等人[2]建立了随机地体二层散射模型(Random Volume over Ground, RVoG),该模型将森林散射场景抽象为两层,即由随机均匀分布的散射体组成的植被层,以及微波信号不可穿透的地表层。随后,Papathanassiou等人[3,4]进一步分析PolInSAR复相干性与RVoG模型的关联,建立了利用PolInSAR反演森林高度的框架。实质上,该框架是基于体散射去相干的模型表达来反演森林高度等参数的,并且利用不同PolInSAR数据都获得了较高的反演精度[5-7]。
由于森林场景具有显著的时变性,具有长时间间隔的星载重轨干涉SAR(如ALOS-1至少为46天,ALOS-2为14天)散射场景内介电常数变化(如降雨)和风动都会产生严重的时间去相干。因此,除了体去相干的影响,时间去相干也是星载重轨极化干涉SAR数据中不可忽略的去相干因素,决定了森林参数反演的精度,甚至是反演成败的关键。为此,Yang等人[8]在随机移动散射模型(Random Motion over Ground, RMoG)[9]和体时去相干散射模型 (Volume Temporal Decorrelation, VTD)模型[10]基础上,提出一种时间去相干半经验森林高度反演方法。该方法结合少量机载LiDAR森林高度数据辅助时间去相干半经验模型解算,利用ALOS-1 PARSAR-1 HV极化相干幅度成功实现了大尺度森林高度反演。
然而,该方法需假设HV极化不包含地表散射回波能量贡献,事实上,L波段SAR信号具有较强的穿透性,尤其当森林高度较低或密度较小时,HV极化方式会记录显著地表回波信号。此外,该方法只适用于单基线干涉数据,尚未考虑多基线条件下,如何充分利用观测几何的多样性提升反演结果的可靠性。因此本文的目的是针对上述反演方法的限制,利用ALOS-2 PARSAR-2多基线PolInSAR数据更为准确地提取森林高度。主要思路如下:首先利用相干最大分离算法(Maximum Coherence Difference, MCD)在极化空间内寻求具有最少地面散射能量贡献的极化方式,以获得更为纯净的森林冠层散射贡献。然后利用该极化方式的相干幅度,在少量森林高度地面调查数据辅助下基于时间去相干半经验模型进行森林高度反演。在此基础之上,结合多基线数据根据PolInSAR相干集在复数平面内的几何表达,甄选最优观测干涉数据的反演结果作为森林高度反演最终结果。
综合顾及垂直方向上散射体分布产生的体去相干、散射场景内介电特性改变和植被风动引起的时间去相干同时占主导地位,星载重轨PolInSAR复相干系数一般形式表示为[8]
γ(ω)=eiφ0⋅γvd⋅γv/m+γgd⋅μ(ω)1+μ(ω) | (1) |
式中,
γv/m=∫h0exp[−12(4πλ)2σ2r(z)]⋅f(z)⋅exp(ikzz)dz∫h0f(z)dz | (2) |
其中,
σr(z)=σrhrz | (3) |
式中,
为了解决上述模型过参数化问题,Yang等人[8]对散射场景做如下假设:(1) 散射场景内时间去相干与消光系数在空间上具有一致性;(2) HV极化方式具有较小地面散射能量贡献,可假定其地体幅度比
γ=Sscene⋅exp[−12(4πσrαλhr)2h2]≈Sscene⋅sinc(hCscene) | (4) |
其中,
Sscene=|γvd|; Cscene=λhr2π2σrα | (5) |
已有方法主要选用对森林冠层较为敏感的HV极化方式进行模型求解,但是ALOS-2 PALSAR-2发射具有较强穿透能力的L波段电磁波,HV极化方式回波信号中同样会记录显著地表回波信号。鉴于此,本文利用ALOS-2 PALSAR-2全极化数据结合极化相干最优理论,尽可能抑制地表回波信号的干扰。具体方法如下:
在主辅极化SAR影像散射机制相同的情况下,极化干涉SAR的复相干系数表示为[12]
γ=⟨ωHΩ12ω⟩√⟨ωHT11ω⟩⟨ωHT22ω⟩ | (6) |
式中,自相关矩阵
Re(γeiϕ)=ωHAωωHTωA=eiϕΩ12+e−iϕΩH122,T=T11+T222} | (7) |
式中,
γ1=ωH1Ω12ω1ωH1Tω1, γ2=ωH2Ω12ω2ωH2Tω2 | (8) |
相比传统InSAR技术只能获取HH, HV或VV极化方式对应的复相干系数,相干集中包含了特定极化散射机理对应的复相干系数,为寻求极化空间内具有更为纯净森林冠层散射贡献的极化方式提供了可能。相干区域范围示意如图1所示,其中在相干区域成对边界点中距离最远的一对相干系数点
kz>0:ifarg(γAγ∗B)>0thenγ(μmin)=γA,γ(μmax)=γBifarg(γAγ∗B)<0thenγ(μmin)=γB,γ(μmax)=γAkz<0:ifarg(γAγ∗B)<0thenγ(μmin)=γA,γ(μmax)=γBifarg(γAγ∗B)>0thenγ(μmin)=γB,γ(μmax)=γA} | (9) |
其中,
kz=4πB⊥λRsinθ | (10) |
即便简化了模型参数和采用多基线PolInSAR数据增加了观测量,利用传统多维非线性迭代求解时间去相干半经验模型仍存在秩亏问题。因此本文采用一种外部数据辅助反演法[8],即先利用小范围真实森林高度数据辅助解算出模型参数
主成分分析思想[17]为实现上述思路提供了契机,即通过对训练数据中
X=[Var(hreal)Cov(hreal,hinvert)Cov(hinvert,hreal)Var(hinvert)]=[P11P12P21P22][λ100λ2][P11P12P21P22]−1k=P21P11} | (11) |
其中
b=M(hreal)−M(hinvert)[M(hreal)+M(hinvert)]/2 | (12) |
式中,M 表示取平均运算。
散点点阵椭圆主轴确定后,显然可以通过建立使逼近参数
(k−1)2+(b−0)2=min | (13) |
该目标函数可以利用高斯-牛顿迭代算法进行非线性最小二乘求解,如式(14)所示
[S∗sceneC∗scene]=(JT0J0)−1JT0[1−k00−b0]+[Sscene0Cscene0] | (14) |
式中,
J0=[∂k∂Sscene∂k∂Cscene∂b∂Sscene∂b∂Cscene]|Sscene0Cscene0 | (15) |
然后将得到修正后的模型参数作为新的初始点进行下一次迭代,经过多次迭代后即可获得最佳模型参数
|[S∗sceneC∗scene]−[S(∗−1)sceneC(∗−1)scene]|<ε | (16) |
在利用训练数据求得时间去相干半经验模型参数后,对每个像元求解一元非线性方程得到整个散射场景内的森林高度结果。
时间去相干、体去相干以及其他噪声等因素会共同影响PolInSAR复相干性在复平面单位圆上的几何表达[18]。在多基线配置下,不同干涉对在同一分辨单元内往往呈现出不同的相干区域结构(如图1所示)。而相干特性
P=|γ(μmin)−γ(μmax)||γ(μmin)+γ(μmax)| | (17) |
式中,
通过相干特性指标P 甄选出不同干涉对在同一分辨单元内反演出的最优森林高度值作为多基线PolInSAR森林高度融合结果,多基线PolInSAR融合反演框架可表示为
max‖P1(γ1(μmin),γ1(μmax))P2(γ2(μmin),γ2(μmax))⋮PN(γN(μmin),γN(μmax))‖ | (18) |
式中,N 为极化干涉SAR观测基线数。
研究区域黄丰桥国有林场(27°
地面实测数据由中南林业科技大学于2016年6~7月采集得到,通过在林区范围内选取60个相互独立的林分样地以确保避免空间自相关,每个林分样地规格为30×30 m。树高则基于单木测高原理利用激光测高仪测得,林分高度范围为4.60~20.20 m,平均高度为13.24 m。本文通过随机采样,将60个林分样地数据随机分为45个训练数据(图3黄点所示)和15个验证数据(图3红点所示)两组。
多基线星载重轨PolInSAR数据是利用日本宇航局(JAXA)提供的5景覆盖研究区域的ALOS-2 PALSAR-2 L波段全极化数据。该SAR影像范围如图3蓝色虚线所示,获取时间为2016年6月至8月,获取模式为StripMap2(SM2),影像主要参数信息如表1所示。将5景SAR影像组成3个时间基线为14天的干涉影像对(BL1, BL2, BL3),然后各自进行配准处理,并进行公共带通滤波以确保去除几何去相干。相干性以11×11窗口进行估计,并应用Boxcar滤波进行平滑处理以消除斑点效应。最后利用SRTM DEM对SAR影像进行地理编码,并将其重采样至与DEM空间分辨率一致(30×30 m)。图4为3个干涉对的HV极化和
日期(2016年) | 垂直有效波数(rad/m) | 时间基线(天) | 距离向/方位向分辨率(m) | 中心入射角 (°) | 极化方式 |
0616—0630 (BL1) | 0.013~0.015 | ||||
0630—0714 (BL2) | 0.010~0.011 | 14 | 2.86/2.97 | 38.99 | Full |
0811—0825 (BL3) | 0.009~0.010 |
以选取的15个验证林分的实测森林高度(H-field)对反演结果(H-invert)进行分析评价,图5为3个干涉对利用HV极化反演得到的散点图结果,均方根误差RMSE分别为:4.20 m, 4.03 m和3.42 m。利用
从上述单基线森林高度反演结果看,不同干涉对反演整体精度较为接近,但是对于同一林分在利用不同干涉对反演的结果却存在明显差异。因此,当多基线数据可用时,我们进一步在单基线PolInSAR森林高度反演结果的基础上挖掘PolInSAR数据本身特性并对其森林高度反演能力进行评判。与时间去相干相关的参数
模型参数 | BL1 | BL2 | BL3 |
Sscene | 0.69 | 0.78 | 0.78 |
Cscene | 9.88 | 10.08 | 11.14 |
3个干涉对在验证林分的相干特性P值、反演森林高度值以及多基线融合森林高度值如表3所示,从整体看,根据相干特性P值大小从3个单基线PolInSAR反演结果中甄选出的森林高度结果更接近于实测真实森林高度。整个实验区的多基线PolInSAR融合反演结果以及精度评定如图7所示,均方根误差RMSE为2.05 m,相比于已有的方法,本文提出的多基线PolInSAR融合反演策略精度至少提高了40%(与图5中BL3基线结果对比),同时,相关系数也提升至0.81。
林分样地编号 | BL1 P值 / 森林高度(m) | BL2 P值 / 森林高度(m) | BL3 P值 / 森林高度(m) | 多基线融合结果(m) | 实测森林高度(m) |
1 | 0.130 / 17.82 | 0.113 / 17.02 | 0.081 / 16.89 | 17.82 | 14.43 |
2 | 0.116 / 14.38 | 0.104 / 15.30 | 0.091 / 16.52 | 14.38 | 14.20 |
3 | 0.092 / 12.46 | 0.075 / 15.83 | 0.135 / 11.34 | 11.34 | 9.80 |
4 | 0.103 / 15.21 | 0.111 / 15.34 | 0.119 / 14.19 | 14.19 | 16.00 |
5 | 0.106 / 6.86 | 0.106 / 7.24 | 0.131 / 8.31 | 8.31 | 10.70 |
6 | 0.110 / 12.98 | 0.083 / 14.67 | 0.118 / 11.89 | 11.89 | 13.50 |
7 | 0.114 / 13.35 | 0.096 / 15.30 | 0.101 / 16.10 | 13.35 | 13.43 |
8 | 0.079 / 14.29 | 0.106 / 16.15 | 0.117 / 16.22 | 16.22 | 16.95 |
9 | 0.069 / 12.12 | 0.090 / 17.63 | 0.060 / 12.30 | 17.63 | 20.10 |
10 | 0.104 / 12.33 | 0.089 / 13.67 | 0.102 / 11.72 | 12.33 | 15.60 |
11 | 0.075 / 18.34 | 0.103 / 16.75 | 0.154 / 10.16 | 10.16 | 13.30 |
12 | 0.113 / 9.08 | 0.134 / 9.46 | 0.106 / 12.69 | 9.46 | 11.00 |
13 | 0.086 / 13.76 | 0.096 / 9.07 | 0.109 / 16.00 | 16.00 | 16.40 |
14 | 0.197 / 10.17 | 0.230 / 8.71 | 0.186 / 9.51 | 8.71 | 6.00 |
15 | 0.103 / 14.59 | 0.064 / 19.17 | 0.128 / 15.40 | 15.40 | 14.70 |
在多基线全极化数据可用条件下,弥补单基线InSAR观测信息不足以及几何结构单一的问题,对于反演结果整体精度提升具有重要作用。本文提出了一种星载重轨多基线PolInSAR反演森林高度的策略,对InSAR极化空间和观测几何空间进行扩展,主要结论如下:
(1) 该方法利用MCD相干优化算法获得对体散射最为敏感的极化方式,并基于时间去相干半经验模型进行森林高度反演,使每条单基线反演精度在一定程度上都有所提高。
(2) 利用由相干特性指标P确定的相干区域最优准则可以优选出同一分辨单元内最优的单基线森林高度反演结果。因此,相比仅利用单基线单一极化反演方法,多基线PolInSAR融合策略具有更好的稳定性,精度也更高。
[1] |
LIU Y C, WANG W, PAN X Y, et al. Inverse omega-K algorithm for the electromagnetic deception of synthetic aperture radar[J]. IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing, 2016, 9(7): 3037–3049. doi: 10.1109/JSTARS.2016.2543961
|
[2] |
ZHAO B, ZHOU F, SHI X R, et al. Multiple targets deception jamming against ISAR using electromagnetic properties[J]. IEEE Sensors Journal, 2015, 15(4): 2031–2038. doi: 10.1109/JSEN.2014.2368985
|
[3] |
SHI X R, ZHOU F, ZHAO B, et al. Deception jamming method based on micro-Doppler effect for vehicle target[J]. IET Radar, Sonar & Navigation, 2016, 10(6): 1071–1079. doi: 10.1049/iet-rsn.2015.0371
|
[4] |
ZHAO B, ZHOU F, and BAO Z. Deception jamming for squint SAR based on multiple receivers[J]. IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing, 2015, 8(8): 3988–3998. doi: 10.1109/JSTARS.2014.2322612
|
[5] |
ZHAO B, HUANG L, ZHOU F, et al. Performance improvement of deception jamming against SAR based on minimum condition number[J]. IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing, 2017, 10(3): 1039–1055. doi: 10.1109/JSTARS.2016.2614957
|
[6] |
徐少坤. SAR欺骗式干扰技术研究[D]. [硕士论文], 国防科学技术大学, 2007.
XU Shao-kun. Research on SAR deception jamming methods[D]. [Master dissertation], National University of Defense Technology, 2007.
|
[7] |
陈志勇, 付耀文, 黎湘. 实时SAR欺骗干扰系统研究[J]. 信号处理, 2009, 25(8A): 546–548.
CHEN Zhi-yong, FU Yao-wen, and LI Xiang. A study on the realtime deceptive jamming system to SAR[J]. Signal Processing, 2009, 25(8A): 546–548.
|
[8] |
汪连栋, 柏仲干, 王国良. 考虑失配的SAR欺骗干扰分析与仿真[J]. 航天电子对抗, 2006, 22(4): 30–32. doi: 10.3969/j.issn.1673-2421.2006.04.010
WANG Lian-dong, BAI Zhong-gan, and WANG guo-liang. Simulation and analysis of deception jamming effects against SAR considering mis-matched filter effection[J]. Aerospace Electronic Warfare, 2006, 22(4): 30–32. doi: 10.3969/j.issn.1673-2421.2006.04.010
|
[9] |
LI C and ZHU D Y. The detection of deception jamming against SAR based on dual-aperture antenna cross-track interferometry[C]. Proceedings of 2006 CIE International Conference on Radar, Shanghai, China, 2006: 1–4. doi: 10.1109/ICR.2006.343274.
|
[10] |
QIN J M, YANG J, HE Z H, et al. Analysis of target loss due to suppressing SAR jamming using dual-channel cancellation[C]. Proceedings of 2006 CIE International Conference on Radar, Shanghai, China, 2006: 1–4. doi: 10.1109/ICR.2006.343324.
|
[11] |
张双喜, 孙光才, 刘艳阳, 等. 基于四通道SAR欺骗式干扰抑制算法[J]. 现代雷达, 2011, 33(2): 22–26, 30. doi: 10.3969/j.issn.1004-7859.2011.02.006
ZHANG Shuang-xi, SUN Guang-cai, LIU Yan-yang, et al. Deception interference supression algorithm based on four-channel for SAR[J]. Modern Radar, 2011, 33(2): 22–26, 30. doi: 10.3969/j.issn.1004-7859.2011.02.006
|
[12] |
SUN B Z and LI J W. A new interference elimination method for multi-satellite SAR system[C]. Proceedings of 2008 IEEE International Geoscience and Remote Sensing Symposium, Boston, MA, USA, 2008: IV-1316–IV-1319. doi: 10.1109/IGARSS.2008.4779973.
|
[13] |
XIONG W, ZHANG G, WEN F Q, et al. Trilinear decomposition-based spatial-polarisational filter method for deception jamming suppression of radar[J]. IET Radar, Sonar & Navigation, 2016, 10(4): 765–773. doi: 10.1049/iet-rsn.2015.0348
|
[14] |
WANG B, GUI G L, ZHANG S, et al. Deceptive jamming suppression based on coherent cancelling in multistatic radar system[C]. Proceedings of 2016 IEEE Radar Conference (RadarConf), Philadelphia, PA, USA, 2016: 1–5. doi: 10.1109/RADAR.2016.7485304.
|
[15] |
SOUMEKH M. SAR-ECCM using phase-perturbed LFM chirp signals and DRFM repeat jammer penalization[J]. IEEE Transactions on Aerospace and Electronic Systems, 2006, 42(1): 191–205. doi: 10.1109/TAES.2006.1603414
|
[16] |
赵宗锋. 基于信号波形设计的SAR抗干扰技术研究[D]. [硕士论文], 国防科学技术大学, 2007.
ZHAO Zong-feng. SAR anti-jamming based on waveform design[D]. [Master dissertation], National University of Defense Technology, 2007.
|
[17] |
李江源, 王建国. 利用复杂调制LFM信号的SAR抗欺骗干扰技术[J]. 电子与信息学报, 2008, 30(9): 2111–2114. doi: 10.3724/SP.J.1146.2007.00231
LI Jiang-yuan and WANG Jian-guo. The anti-jamming technology for SAR cheat jamming using complex modulated LFM signal[J]. Journal of Electronics &Information Technology, 2008, 30(9): 2111–2114. doi: 10.3724/SP.J.1146.2007.00231
|
[18] |
冯祥芝, 许小剑. 随机线性调频斜率SAR抗欺骗干扰方法研究[J]. 系统工程与电子技术, 2009, 31(1): 69–73. doi: 10.3321/j.issn:1001-506X.2009.01.018
FENG Xiang-zhi and XU Xiao-jian. Study of countermeasures to deceptive jamming using random linear modulation frequency ratio SAR[J]. Systems Engineering and Electronics, 2009, 31(1): 69–73. doi: 10.3321/j.issn:1001-506X.2009.01.018
|
[19] |
LI W, LU X Q, DA X Y, et al. Anti-jamming method based on orthogonal codes jittered and random initial phase for SAR[C]. Proceedings of 2007 IET International Conference on Radar Systems, Edinburgh, U.K., 2007: 1–5. doi: 10.1049/cp: 20070651.
|
[20] |
FENG Q Q, XU H P, WU Z F, et al. Deceptive jamming suppression for SAR based on time-varying initial phase[C]. Proceedings of 2016 IEEE International Geoscience and Remote Sensing Symposium (IGARSS), Beijing, China, 2016: 4996–4999. doi: 10.1109/IGARSS.2016.7730303.
|
[21] |
HOSSAIN M A, ELSHAFIEY I, ALKANHAL M A, et al. Anti-jamming capabilities of UWB-OFDM SAR[C]. Proceedings of the 8th European Radar Conference, Manchester, U.K., 2011: 313–316.
|
[22] |
罗双才, 唐斌. 一种基于盲分离的欺骗干扰抑制算法[J]. 电子与信息学报, 2011, 33(12): 2801–2806. doi: 10.3724/SP.J.1146.2011.00380
LUO Shuang-cai and TANG Bin. An algorithm of deception jamming suppression based on blind signal separation[J]. Journal of Electronics &Information Technology, 2011, 33(12): 2801–2806. doi: 10.3724/SP.J.1146.2011.00380
|
[23] |
ZHAO B, HUANG L, HE C L, et al. SAR deception jamming identification via differential feature enhancement[C]. Proceedings of 2016 CIE International Conference on Radar, Guangzhou, China, 2016: 1–5. doi: 10.1109/RADAR.2016.8059174.
|
[24] |
ZHAO B, HUANG L, and ZHANG J H. Single channel SAR deception jamming suppression via dynamic aperture processing[J]. IEEE Sensors Journal, 2017, 17(13): 4225–4230. doi: 10.1109/JSEN.2017.2695001
|
[25] |
ZHAO B, HUANG L, LI J, et al. Target reconstruction from deceptively jammed single-channel SAR[J]. IEEE Transactions on Geoscience and Remote Sensing, 2018, 56(1): 152–167. doi: 10.1109/TGRS.2017.2744178
|
[26] |
ELAD M, STARCK J L, QUERRE P, et al. Simultaneous cartoon and texture image inpainting using morphological component analysis (MCA)[J]. Applied and Computational Harmonic Analysis, 2005, 19(3): 340–358. doi: 10.1016/J.ACHA.2005.03.005
|
1. | 汪思源,曲毅,陈怡君. 基于U-Net的涡旋电磁波雷达成像方法. 空军工程大学学报. 2024(03): 77-85 . ![]() | |
2. | 潘浩然,马晖,胡敦法,刘宏伟. 基于涡旋电磁波新体制的雷达前视三维成像. 雷达学报. 2024(05): 1109-1122 . ![]() | |
3. | 毛德庆,杨建宇,杨明杰,张永超,张寅,黄钰林. IAA-Net:一种实孔径扫描雷达迭代自适应角超分辨成像方法. 雷达学报. 2024(05): 1073-1091 . ![]() | |
4. | 马晖,胡敦法,师竹雨,刘宏伟. 基于涡旋电磁波的雷达应用研究进展. 现代雷达. 2023(05): 27-41 . ![]() | |
5. | 袁航,罗迎,陈怡君,苏令华. 基于反正弦圆环天线阵列的二维成像. 北京航空航天大学学报. 2023(06): 1487-1494 . ![]() |
日期(2016年) | 垂直有效波数(rad/m) | 时间基线(天) | 距离向/方位向分辨率(m) | 中心入射角 (°) | 极化方式 |
0616—0630 (BL1) | 0.013~0.015 | ||||
0630—0714 (BL2) | 0.010~0.011 | 14 | 2.86/2.97 | 38.99 | Full |
0811—0825 (BL3) | 0.009~0.010 |
模型参数 | BL1 | BL2 | BL3 |
Sscene | 0.69 | 0.78 | 0.78 |
Cscene | 9.88 | 10.08 | 11.14 |
林分样地编号 | BL1 P值 / 森林高度(m) | BL2 P值 / 森林高度(m) | BL3 P值 / 森林高度(m) | 多基线融合结果(m) | 实测森林高度(m) |
1 | 0.130 / 17.82 | 0.113 / 17.02 | 0.081 / 16.89 | 17.82 | 14.43 |
2 | 0.116 / 14.38 | 0.104 / 15.30 | 0.091 / 16.52 | 14.38 | 14.20 |
3 | 0.092 / 12.46 | 0.075 / 15.83 | 0.135 / 11.34 | 11.34 | 9.80 |
4 | 0.103 / 15.21 | 0.111 / 15.34 | 0.119 / 14.19 | 14.19 | 16.00 |
5 | 0.106 / 6.86 | 0.106 / 7.24 | 0.131 / 8.31 | 8.31 | 10.70 |
6 | 0.110 / 12.98 | 0.083 / 14.67 | 0.118 / 11.89 | 11.89 | 13.50 |
7 | 0.114 / 13.35 | 0.096 / 15.30 | 0.101 / 16.10 | 13.35 | 13.43 |
8 | 0.079 / 14.29 | 0.106 / 16.15 | 0.117 / 16.22 | 16.22 | 16.95 |
9 | 0.069 / 12.12 | 0.090 / 17.63 | 0.060 / 12.30 | 17.63 | 20.10 |
10 | 0.104 / 12.33 | 0.089 / 13.67 | 0.102 / 11.72 | 12.33 | 15.60 |
11 | 0.075 / 18.34 | 0.103 / 16.75 | 0.154 / 10.16 | 10.16 | 13.30 |
12 | 0.113 / 9.08 | 0.134 / 9.46 | 0.106 / 12.69 | 9.46 | 11.00 |
13 | 0.086 / 13.76 | 0.096 / 9.07 | 0.109 / 16.00 | 16.00 | 16.40 |
14 | 0.197 / 10.17 | 0.230 / 8.71 | 0.186 / 9.51 | 8.71 | 6.00 |
15 | 0.103 / 14.59 | 0.064 / 19.17 | 0.128 / 15.40 | 15.40 | 14.70 |
日期(2016年) | 垂直有效波数(rad/m) | 时间基线(天) | 距离向/方位向分辨率(m) | 中心入射角 (°) | 极化方式 |
0616—0630 (BL1) | 0.013~0.015 | ||||
0630—0714 (BL2) | 0.010~0.011 | 14 | 2.86/2.97 | 38.99 | Full |
0811—0825 (BL3) | 0.009~0.010 |
模型参数 | BL1 | BL2 | BL3 |
Sscene | 0.69 | 0.78 | 0.78 |
Cscene | 9.88 | 10.08 | 11.14 |
林分样地编号 | BL1 P值 / 森林高度(m) | BL2 P值 / 森林高度(m) | BL3 P值 / 森林高度(m) | 多基线融合结果(m) | 实测森林高度(m) |
1 | 0.130 / 17.82 | 0.113 / 17.02 | 0.081 / 16.89 | 17.82 | 14.43 |
2 | 0.116 / 14.38 | 0.104 / 15.30 | 0.091 / 16.52 | 14.38 | 14.20 |
3 | 0.092 / 12.46 | 0.075 / 15.83 | 0.135 / 11.34 | 11.34 | 9.80 |
4 | 0.103 / 15.21 | 0.111 / 15.34 | 0.119 / 14.19 | 14.19 | 16.00 |
5 | 0.106 / 6.86 | 0.106 / 7.24 | 0.131 / 8.31 | 8.31 | 10.70 |
6 | 0.110 / 12.98 | 0.083 / 14.67 | 0.118 / 11.89 | 11.89 | 13.50 |
7 | 0.114 / 13.35 | 0.096 / 15.30 | 0.101 / 16.10 | 13.35 | 13.43 |
8 | 0.079 / 14.29 | 0.106 / 16.15 | 0.117 / 16.22 | 16.22 | 16.95 |
9 | 0.069 / 12.12 | 0.090 / 17.63 | 0.060 / 12.30 | 17.63 | 20.10 |
10 | 0.104 / 12.33 | 0.089 / 13.67 | 0.102 / 11.72 | 12.33 | 15.60 |
11 | 0.075 / 18.34 | 0.103 / 16.75 | 0.154 / 10.16 | 10.16 | 13.30 |
12 | 0.113 / 9.08 | 0.134 / 9.46 | 0.106 / 12.69 | 9.46 | 11.00 |
13 | 0.086 / 13.76 | 0.096 / 9.07 | 0.109 / 16.00 | 16.00 | 16.40 |
14 | 0.197 / 10.17 | 0.230 / 8.71 | 0.186 / 9.51 | 8.71 | 6.00 |
15 | 0.103 / 14.59 | 0.064 / 19.17 | 0.128 / 15.40 | 15.40 | 14.70 |