出版日期: 2018-11-25
点击次数:
下载次数:
DOI: 10.11834/jrs.20187377
2018 | Volumn22 | Number 6
上一篇  |  下一篇


  
植被冠层立体结构与叶片倾角的偏振光效应
expand article info 杨彬1,2 , 晏磊2,3 , Yuri Knyazikhin4 , 刘思远2,3
1. 湖南大学 电气与信息工程学院,长沙 410082
2. 北京大学 空间信息集成与3S工程应用北京市重点实验室,北京 100871
3. 北京大学 地球与空间科学学院 遥感与地理信息系统研究所,北京 100871
4. 波士顿大学 地球与环境系,波士顿 MA 02215

摘要

植被偏振特性研究对于植被监测与组分定量反演具有极其重要的作用。植被冠层的反射辐射具有偏振特性,这种特性与入射辐射和植被冠层结构相关。本文分析了偏振对光子—叶片—冠层之间细微相互作用及其变化的有效探测能力,并利用研究型扫描式偏振辐射仪RSP(Research Scanning Polarimeter)数据系统对比分析了偏振对不同叶倾角分布的估测。通过上述研究得出以下结论:(1)偏振观测能够对光线在冠层立体结构中的透射反射再出射过程给出精细刻画,若不用偏振手段对这一过程进行甄别并去除,则直接测算的植被散射系数会产生高达140%的误差;(2)利用偏振手段可以为高精度大倾角、多时相遥感观测提供可能,以此可改变目前光学遥感小角度、垂直观测的较严格约束;(3)偏振辐射呈现出随波长的稳定特性(相关系数>0.96),使得利用偏振手段可以更好地研究冠层结构;(4)不同叶倾角分布对入射辐射存在不同的偏振反射,为利用多角度偏振信息进行遥感植被精细分类提供了新的途径。本文详细描述冠层结构和植被偏振特性的相互作用,通过对冠层立体结构与叶倾角的研究,刻画了植被定量遥感的方向性信息与高精度实现,为高分辨率遥感定量化的有效信息挖掘提供了新手段。

关键词

入射观测几何, 冠层立体结构, 叶倾角, 偏振光效应, RSP, 高分辨率遥感

Vegetation polarimetric effect of three-dimensional structure and leaf inclination of canopy
expand article info YANG Bin1,2 , YAN Lei2,3 , Yuri Knyazikhin4 , LIU Siyuan2,3
1.School of Electric and Information Engineering, Hunan University, Changsha 410082, China
2.Beijing Key Laboratory of Spatial Information Integration and 3S Application, Peking University, Beijing 100871, China
3.Institute of Remote Sensing and Geographic Information System, School of Earth and Space Sciences, Peking University, Beijing 100871, China
4.Earth and Environment Department, Boston University, Boston MA02215, US

Abstract

Studying the polarization characteristics of vegetation is important in vegetation monitoring and component quantitative inversion. The reflected radiation from the vegetation canopy possesses a polarization characteristic that is related to the interaction of incident radiation and vegetation canopy and the distribution of leaf inclination angles. This study analyzes the effect of photon–canopy interaction on the polarization characteristics of vegetation and Research Scanning Polarimeter (RSP) utilizes data combined with the Rondeaux and Herman model to study the polarization characteristics of vegetation under different leaf inclination distributions. A systematic comparative analysis is of more significant than a polarization study considering that leaf inclination distribution satisfies spherical distribution. This study presents the mechanism of photon–canopy interaction and some existing results related to the impact of canopy structure on the polarization characteristics of vegetation. As for the impact of leaf inclination, this study utilizes RSP data collected on September 16, 2005 in Oklahoma, USA, after atmospheric correction and dense vegetation extraction, and then combines the Rondeaux and Herman model for a systematic comparative analysis on wavelength, sun observation geometry, and leaf inclination distribution. In this study, we analyzed the impact of photon–canopy interaction on canopy structure and leaf inclination distribution on the polarization characteristics of vegetation. We found that the canopy structure is highly related with polarization characteristics, and the removal of the canopy structure makes polarization a major factor that influences vegetation reflection. We analyzed the stability of the polarimetric reflectance of the vegetation in wavelength range [400, 2500], the distribution of specular reflection in different directions, and the influence of different leaf inclinations and canopy structures on polarimetric radiation. Results of these analyses are clearly presented in the figures in this paper.(1) The effects of canopy structure and leaf polarization on the canopy reflectance of vegetation are closely related. Removing the influence of canopy structure is beneficial to the analysis of polarization characteristics and makes the polarization characteristics a dominant factor affecting canopy reflection. Neglecting the polarization property makes the scattering coefficient (obtained by removing canopy structure impact) produce up to 140% error; this part of the error is almost negligible before removing the canopy structure. (2) The polarimetric effect of multiangle and multitemporal remote sensing on the surface cannot be ignored. Multiangle polarization characteristic analysis and error rejection provide many possibilities for remote-sensing observations under low confinement conditions and help expand the space–time scale of earth observation. (3) Polarization characteristics under a given incident observation geometry is not sensitive to wavelength, and the mirror reflection of the canopy for the incident radiation exists in all directions due to the distribution of the blade inclination angle. Polarization reflectance increases with the increase in phase angle, and polarization characteristics are affected by the distribution of the leaf inclination angle, especially at large phase angles, allowing the use of polarization remote sensing to classify species and obtain canopy structural parameters. In this study, the polarization method is used to represent a high-precision and rich directional information acquisition process in high-resolution quantitative remote sensing and is the further exploration of effective information based on existing remote sensing methods. The effective information contributes to the further description and study of the greenhouse effect and provides a theoretical basis for the basic propositions in high-resolution remote sensing inversion.

Key words

sun-sensor geometry, canopy structure, leaf inclination, polarimetric characteristics, RSP, high resolution remote sensing

1 引 言

常规遥感传感器通过测量植被冠层反射或散射的入射辐射来量化植被的光学特性,其由于非接触、大面积等特点,被广泛应用于植被监测、辐射收支和全球气候变化等研究中。随着卫星传感器分辨率的提高与地面测量、模型构建等同步技术的进步,对于植被更高精度监测与组分反演的需求与日俱增。

Ollinger和Knyazikhin两位科学家团队自2008年起在美国科学院院刊PNAS(Proceedings of the National Academy of Sciences of the United States of America)上对于利用植被反射光谱进行氮含量反演的争论引起广泛关注。Ollinger等人(20082013)提出植被在近红外的反射率与氮含量之间存在正相关关系,Townsend等人(2013)指出该正相关关系并不包含氮含量的信息,而是由冠层结构与叶片光学特征造成,并提出需要偏振手段研究冠层结构与表面反射(Knyazikhin 等,2013a, 2013b, 2013c)。以Ollinger为代表的研究团队从宏观的生态学角度出发,将植被看作一个对入射辐射直接反射的整体;而Knyazikhin团队从遥感光学角度出发考虑了植被的立体特性,即光子在植被中微小的反射、透射、折射等特征,这与前者存在本质上的区别。这些特征会体现在植被的偏振反射特性上,因此偏振探测手段有利于从冠层角度更精细地刻画植被反射的光学特征,从而为甄别PNAS上引出的争议提供手段。偏振效应在植被冠层结构与高分辨率定量遥感反演中的重要性日渐凸显,对植被偏振特性的研究也成为丰富植被遥感体系的基础性工作。

目前,国内外利用偏振手段进行的遥感反演研究主要集中于大气方面,实现了气溶胶和多种云光学参数的遥感反演(Guo 等,2014Dubovik 等,2011Sano,2004)。利用偏振信息进行参数反演的研究较少,主要存在于实验室条件或利用卫星影像进行大区域研究。例如以单叶片偏振BRDF模型模拟其偏振反射特性,并利用实测值有效反演叶片漫反射率和表面粗糙度等结构参数(Xie 等,2010);利用POLDER卫星数据结合偏振信息和辐射强度信息构建模型,可以同时反演大区域大气气溶胶和地表反射率的区域分布(Duan和Lü,2008)。对于植被的偏振遥感研究则多集中于其偏振反射特性与对叶片和冠层结构的模拟和经验、物理模型的构建方面。较早期的植被偏振测量和物理特性分析(Egan,1970Curran, 1981, 1982Fitch 等,1984)以及后来分别对叶片和植被冠层的偏振信息提取(Grant 等,1993Vanderbilt 等,1985Woessner和Hapke,1987)集中地得出了一些一致性结论,即植被的散射特性由辐射在植被结构中的多次散射造成,叶表的镜面反射是导致反射光具有偏振态的直接原因(Grant,1987Grant 等,1987),并且偏振特性会随入射观测几何的改变而改变(Rondeaux和Herman,1991Vanderbilt and Grant, 1985Yang 等,2016)。这些结论为光在叶片和植被内部的传输、散射和偏振特性研究提供了理论基础,利用菲涅尔反射折射定律所构建的诸多模型也广泛地用来模拟计算植被的偏振二向反射值(Rondeaux和Herman,1991Waquet 等,2009Nadal和Breon,1999Maignan 等,2009Xie 等,2017),与实验室和野外测量的多种植被类型偏振反射实测值相比较均产生了相近的结果(Sun 等,2016)。随着仪器和测量手段的发展,更高精度、更高的角度和光谱分辨率的实测数据和现存模型的验证和改进工作越来越多(Suomalainen 等,2009bWu和Zhao,2005Puttonen 等,2009),从而扩展了利用偏振信息进行植被遥感研究的应用范围,并不断证明了其巨大潜能。此外,利用偏振信息还可以对岩石、水体、雪等地表典型地物开展多种研究(Zhao 等,2004Shen 等,2016Suomalainen 等,2009aPeltoniemi 等,1994)。

多种偏振反射率二向分布模型精细刻画了光在植被冠层中的辐射传输过程。从理论上来说,冠层反射的偏振特性主要由冠层结构决定(Myneni和Ross,1991),即叶片在冠层内部的分布对光子—冠层作用的影响及冠层内的叶倾角状况。Vanderbilt和Grant(1985)提出了一种模型,利用叶面积指数LAI(Leaf Area Index),叶倾角分布函数LIDF(Leaf Inclination Distribution Function),粗糙度因子和逃逸概率来描述植被在2π空间内的偏振反射分布,并指出偏振信息有助于区分树种。之后,Rondeaux和Herman(1991)基于均匀冠层提出了一个更简单的偏振反射率模型。基于Rondeaux-Herman模型,Breon等人(1995)推导了裸土和植被的偏振反射模型。随着POLDER数据的获取,Nadal和Breon(1999)建立了地表双参数经验双向偏振分布模型。一些研究考虑阴影效应,提出了具有阴影函数的偏振反射率模型(Waquet 等,2009Litvinov 等,2011)。Maignan等人(2009)还引入了结合归一化差分植被指数(NDVI)和地表偏振的单参数模型。基于以上研究成果,Xie等人(2017)提出了城市地区的偏振反射率模型,取得了较高精度结果,这个模型对于植被也具有适用性。实际上,Rondeaux-Herman模型之后的大多数模型均基于叶片球形分布的假设构建,并没有严格考虑不同叶倾角分布函数的影响。

基于上述研究,本文从理论推导出发,介绍光子—叶片—冠层的相互作用,即植被的叶片光学特征与冠层立体结构,分析其对于植被偏振光效应的影响,以期甄别PNAS上关于植株氮含量反演的争议,并在Rondeaux-Herman模型的基础上,通过引入不同叶倾角分布来研究叶倾角对植被偏振光效应的影响,以完善偏振反射模型,实现植被叶倾角分布的反演,并探索由叶倾角产生的多角度偏振信息的更多应用价值。

2 叶片偏振光学特性与植被冠层结构

2.1 反射过程中的叶片偏振光学特性

植被反射辐射具有偏振特性的本质原因是光子—叶片的单次散射,多次散射具有褪偏作用而在探测到的偏振反射率中忽略不计。光子—叶片的单次散射可描述如下:被叶片散射的入射辐射包括两部分,第1部分是在空气和叶面分界处发生镜面反射,这部分反射辐射没有进入叶片内部,是部分偏振的(Grant 等,1993)。第2部分进入叶片,与叶片内部的生化物质相互作用后被传感器接收,这部分通常是非偏振的(图1)。基于强度的常规遥感可以探测到这两部分的总和,但不能对二者进行有效区分(Vanderbilt和Grant,1985)。利用偏振探测能够提取其中的线偏振部分,实现二者有效区分。

图 1 传感器接收到的来自叶片的两部分辐射
Fig. 1 Two-part radiation from the leaf received by the sensor

图1(Yang 等,2016),可以假设叶表的散射反照率为 $\omega $ ,入射辐射在叶表直接被反射的概率为 ${r_{\rm{s}}}$ ,在叶片内部向上散射的概率为 $\beta $ ,向下散射的概率为 $\tau $ 。根据菲涅尔原理,可以知道在叶表直接被反射的辐射是部分偏振的,也就是说 ${r_{\rm{s}}}$ 可以进一步分解为

${r_{\rm{s}}} = {r_{{\rm{sp}}}} + {\rm{ r\_snp}}$ (1)

式中, ${r_{{\rm{sp}}}}$ 代表线偏振部分,而 ${\rm{r\_snp}}$ 代表非偏振部分。

进入到叶片内部,与其相互作用后,辐射被向上散射的总概率为

${r_{\rm{i}}} = (1 - {r_{\rm{s}}})\omega \beta $ (2)

则辐射被吸收的概率 ${a_{\rm{i}}}$ 以及向下透射的概率 ${t_{\rm{i}}}$

${a_{\rm{i}}} = (1 - {r_{\rm{s}}})(1 - \omega)$ (3)
${t_{\rm{i}}} = (1 - {r_{\rm{s}}})\omega \tau $ (4)

传感器接收到的辐射与叶表及内部相互作用后的总能量则可以简单表述为

${r} = {r_{\rm{s}}} + (1 - {r_{\rm{s}}})\omega \beta $ (5)

则传感器探测到的偏振度p

$p = \frac{{{r_{{\rm{sp}}}}}}{{{r_{\rm{s}}} + (1 - {r_{\rm{s}}})\omega \beta }}$ (6)

而一般情况下, ${r_{{\rm{sp}}}}$ 可以认为是常量,即是与波长无关的量。当叶片的反射率较低时,即 $\omega \beta $ 表现为较小值,从视觉效果来看,目标会显得很暗,则偏振度 $p$ 表现为较大值,则实现了“弱光强化”的过程;反之,当地表的反射率较大时,即 $\omega \beta $ 表现为较大值,从视觉效果来看,目标会显得很亮,则偏振度 $p$ 表现为较小值,则实现了“强光弱化”的过程。这是偏振光能够敏感植被冠层立体结构精细的透射反射再出射的物理基础,也是偏振具有“强光弱化,弱光强化”的光学遥感特性的通用物理表征。

2.2 冠层结构与偏振反射特性解耦

上述光子—叶片的单次散射造成的偏振特性在叶片尺度可较为明显的探测与分辨。而在冠层尺度,这种直接探测得到的偏振特性与冠层结构是相互耦合的。

在可见光—近红外范围内,尽管冠层的吸收过程与波长相关,但是光子与冠层相互作用的概率却是仅仅由冠层结构决定,与波长无关,称之为光谱不变量理论(Knyazikhin 等,1998)。利用光谱不变量理论可以精确描述冠层结构及光子与冠层的相互作用。

图2,假设冠层底部为不反射任何辐射的黑体,辐射入射到冠层后,被冠层截获的概率为截获率i0。没有被截获的那部分辐射将直接透过冠层,到达冠层底部,这一部分的概率为透过率t0=1–i0。被截获的辐射与叶片第一次相互作用后,有a1=(1–ωλ)i0部分辐射被吸收,其中ωλ为叶片单次散射反照率。光子与冠层内部的散射体相互作用后再次与其他散射体相互作用的概率定义为再碰撞概率p,则有i1=ωλi0p部分辐射会与冠层内部的其他叶片相互作用,另外有s1=(1–p)ωi0的辐射会被散射出冠层。i1部分会激发与叶片的第二次相互作用,分解为吸收、散射出冠层和继续与其他叶片相互作用3部分,以此类推,发生多次散射。从辐射多次散射后的分布可以发现,每次散射的吸收和逃逸能量呈现等比数列关系。

假设冠层底部是不反射辐射的黑体,且再碰撞概率在各散射级数保持不变,则冠层与辐射相互作用的吸收a( $ \lambda$ )和散射s( $\lambda$ )分别为

$a\left( \lambda \right) = {i_0}\left( {1 - {\omega _\lambda }} \right)\sum\limits_{k = 0}^\infty {{{\left( {{\omega _\lambda }p} \right)}^k}} = \frac{{1 - {\omega _\lambda }}}{{1 - p{\omega _\lambda }}}{i_0}$ (7)
$s\left( \lambda \right) = {i_0}\left( {1 - p} \right){\omega _\lambda }\sum\limits_{k = 0}^\infty {{{\left( {{\omega _\lambda }p} \right)}^k}} = \frac{{\left( {1 - p} \right){\omega _\lambda }}}{{1 - p{\omega _\lambda }}}{i_0}$ (8)
图 2 光子在冠层内的多次散射与多次吸收
Fig. 2 Multiple scattering and multiple absorption of photons in the canopy

若将逃逸概率(1–p)分解到空间各个方向,并用ρ表示(即ρ( $ \varOmega $ )表示光子与叶片碰撞后从方向 $ \varOmega $ 逃逸出冠层的概率),则可以得到如图2所示的光子多次散射与逃逸关系,基于此可以获得冠层双向反射率因子BRF(Bidirectional Reflectance Factor)

${\rm{BRF}}\left({{{\rm{\varOmega }}_0}, {\rm{\varOmega }}} \right) = \frac{{\rho \left({\rm{\varOmega }} \right){i_0}\left({{{\rm{\varOmega }}_0}} \right)}}{{1 - p{\omega _\lambda }}}{\omega _\lambda }$ (9)

这样,冠层的散射和吸收以及光子逃逸关系可以通过光谱不变量精确描述,利用光谱不变量理论构建结构参数并进行冠层结构—叶片光学特性解耦至关重要。利用光谱不变量理论可以建立起方向散射面积因子DASF和散射系数W $ _\lambda$ 的概念(Knyazikhin 等,2013c)。DASF是冠层沿某方向的可视叶面积比例,是一个光谱不变量。对于浓密植被,DASF可直接通过710—790 nm的反射率获得,并在此基础上获得冠层的光学因子,即散射系数W $ _\lambda$ W $_\lambda$ 可由BRF通过与DASF的比值直接获得,相比于BRF,散射系数与叶片内部的信息(叶绿素含量、化学成分等)更相关,因而使得从遥感高精度反演叶片内部信息成为可能。

2.3 基于叶片与冠层尺度的偏振光效应定量化探测

由上述分析可以看出,遥感传感器观测到的BRF同时受叶片表面反射与冠层结构影响,这两种影响是紧密耦合的。冠层的立体结构使得植被在各角度上的BRF出现不同值,而引入冠层结构的考虑是为了描述光子在冠层内部的二次碰撞与逃逸过程。考虑这些因素计算得到的散射系数(即上述散射系数的计算过程)可以帮助恢复各个观测角度下植被冠层的反射率光谱,使其归化到同一值上。因此冠层结构的解耦为研究植被的多角度(尤其是大角度)偏振特性提供了必要条件,而偏振特性的精细研究也为实现传感器大角度观测、给出精确解算结果提供了可能。

针对引言中提到的争议,一方面从植物冠层出发(Yang 等,2017a)刻画反射特征,另一方面从2009年开始进行大量地面和实验室试验,于2016年在对针叶叶簇样本的研究中发现,植被冠层与叶片表面反射(部分偏振)特性对探测到的叶簇反射率均会产生误差,而只有偏振探测能够敏感出上述作用带来的微弱光线变化导致的误差;在40°天顶角进行观测时,若不考虑此误差,则测算得到的叶簇样本的散射系数在蓝光波段可产生140%的误差(Yang 等,2016)。现阶段,遥感对地观测对于多角度(尤其是大倾角)、多时相观测需求巨大,越来越多的观测平台上开始搭载多角度传感器,因此偏振效应在多角度观测上带来的误差不可忽略。由此看出,考虑光子在冠层内多次碰撞与逃逸的细微过程进行冠层结构的解耦,使得表面反射对植被光谱的影响占据了主导地位,而利用偏振手段可敏感出表面反射产生的变化,从而将表面反射与内部散射部分有效分离,实现对遥感反演参量的精确刻画。

光子—叶片与光子—冠层的相互作用效应最直观地体现在偏振反射率上。偏振反射率二向分布模型BPDF(Bidirectional Polarized Reflectance Distribution Function)描述了理想情况下经地物反射后的偏振辐射在2π空间中的分布,目前已存在包括物理模型和半物理—半经验模型在内的十余种BPDF模型,这些模型起初均以植被中光子—叶片—冠层的相互作用为物理基础进行构建,并将其应用于其他地物分析。作者基于偏振传感器POLDER的全球尺度地表偏振反射率数据,针对5种IGBP分类地物,对6种常用模型进行了拟合精度对比(Yang 等,2017b),表1归纳了整体结果。显然,6种模型的精度均处于较高水平(Cor普遍大于0.86,RMSE普遍小于0.3%),且其中模型对植被(表1中第一、二列)的拟合精度最高。这说明,现存的BPDF模型由光子—叶片—冠层的相互作用精确刻画,这种相互作用可以作为植被遥感定量化的依据,为植被辐射传输过程以及高分辨率定量反演提供理论基础。

表 1 6种BPDF模型拟合结果比较
Table 1 Comparison of six BPDF model results

下载CSV 
模型 参数 IGBP
混合森林
IGBP
农田
IGBP
建成区
IGBP
冰雹
IGBP
稀疏植被
Nadal-Breon RMSE/% 0.168 0.147 0.231 0.367 0.106
Cor 0.922 0.955 0.934 0.948 0.987
Waquet RMSE/% 0.209 0.213 0.253 0.461 0.154
Cor 0.886 0.917 0.920 0.925 0.974
Maignan RMSE/% 0.175 0.156 0.269 0.423 0.130
Cor 0.915 0.949 0.910 0.934 0.981
Litvinov RMSE/% 0.155 0.150 0.201 0.408 0.109
Cor 0.934 0.953 0.950 0.936 0.986
Diner RMSE/% 0.233 0.277 0.342 0.666 0.296
Cor 0.871 0.881 0.863 0.863 0.918
Xie-Cheng RMSE/% 0.174 0.155 0.227 0.386 0.104
Cor 0.916 0.950 0.935 0.942 0.987
 注:RMSE为均方根误差,Cor为模型值与实测值的拟合相关系数。

3 植被叶倾角分布的偏振特性

由上述分析可知,光子与叶片和冠层的相互作用是产生植被偏振光效应的本质原因。偏振具有角度敏感特征,因而对植被叶倾角分布可以实现精密计量,因此相对于传统光学遥感,偏振对叶倾角具备更有效的探测能力。因此基于上节偏振对叶片光学特性与冠层结构的精密探测模型,本节将采用偏振手段进一步实现不同种类植被叶倾角的反演。

3.1 数据与方法

3.1.1 RSP数据

偏振传感器能够将传感器接收的辐射分为偏振部分和非偏振部分(Rondeaux和Herman,1991),使得更精细化的对地观测成为可能。研究型扫描式偏振辐射仪RSP(Research Scanning Polarimeter)可以提供机载和地面的高精度偏振测量。RSP具有9个波段410 nm(30),470 nm(20),550 nm(20),670 nm(20),865 nm(20),960 nm(20),1590 nm(60),1880 nm(90)和2250 nm(120)(Cairns 等,1999)。括号内的数字是对应波段的半峰全宽(FWHM)。RSP数据空间分辨率高,覆盖了具有不同叶倾角分布的多种植被类型,是研究植被偏振特性对叶倾角分布灵敏性的最合适的数据库之一。

本文利用2005年9月16日俄克拉荷马州实验的RSP低空数据(Knobelspiesse 等,2008)。表2展示了本文中用到的航带信息。RSP的观测天顶角(VZA)具有[–40°, 80°]大范围,由于较大VZA可能引入更多的不确定性,本文只利用[–40°, 60°]的VZA范围。θs, θv, φs, φv可以确定入射观测几何,分别代表太阳天顶角、观测天顶角、太阳方位角和观测方位角。由于RSP数据包含了大气和土壤偏振信息,为有效实现叶倾角分布对植被偏振特性研究,需要进行大气校正及浓密植被提取。

表 2 测量航带信息
Table 2 Measurement parameters

下载CSV 
类别 参数
日期 2005年9月16日
获取时间,UTC 16:32:25
飞行高度(海平面)/m 511.3 (均值),15.3 (标准差)
飞行高度(地面)/m 318
太阳—传感器几何 θs: 42.6° (均值),0.1° (标准差)
φs: 136.8° (均值),0.3° (标准差)
θv: –40°—60°
φv: 91.4° (均值),1.7° (标准差)
AERONET τa(λ=670 nm) 0.04 (quite clear)

3.1.2 大气校正

本文选取的RSP数据的获取高度为318 m,传感器获取的上行辐射数据包含了地表和大气的信息。假设传感器在距地高程为z处探测到的上行偏振反射率为 $R_{{\rm{p}}{\rm{ }}}^{{\rm{meas}}}$ ,则(Waquet 等,2009, Wang 等,2014):

$\begin{array}{l}R_{\rm{p}}^{{\rm{meas}}}\left( {z,{\theta _{\rm{s}}},{\theta _{\rm{v}}},{\varphi _{\rm{r}}}} \right) = \\ R_{\rm{p}}^{{\rm{atm}}}\left( {z,{\theta _{\rm{s}}},{\theta _{\rm{v}}},{\varphi _{\rm{r}}}} \right) + \\ R_{\rm{p}}^{{\rm{surf}}}\left( {z,{\theta _{\rm{s}}},{\theta _{\rm{v}}},{\varphi _{\rm{r}}}} \right){T^ \downarrow }\left( {{\theta _{\rm{s}}}} \right){T^ \uparrow }\left( {{\theta _{\rm{v}}}} \right)\end{array}$ (10)

式中, $R_{\rm{p}}^{{\rm{atm}}}$ $R_{\rm{p}}^{{\rm{surf}}}$ 分别代表大气偏振反射率和地表偏振反射率,φr代表相对方位角。T(θs)和T(θv)分别代表大气下行和上行透过率。偏振反射率 $R_{\rm{p}}^{{\rm{atm}}}$ 可由6SV模型精确估算(Kotchenova和Vermote,2007)。下行、上行透过率T(θs)和T(θv)可由Waquet(2009)提出的方法计算,地表偏振反射率 $R_{\rm{p}}^{{\rm{surf}}}$ 则可通过式(1)获得。

3.1.3 浓密植被提取

本文研究对象为植被地表,为有效降低土壤背景的影响,需要选取浓密植被作为研究对象。本文选用抗气溶胶植被指数(ARVI)作为浓密植被选取参数,ARVI对大气效应不敏感(Kaufman和Tanré,1992),可以通过对其设置阈值来提取浓密植被像元。ARVI的表达式为

${\rm{ARVI}} = \frac{{{R_{{\rm{NIR}}}} - {R_{{\rm{rb}}}}}}{{{R_{{\rm{NIR}}}} + {R_{{\rm{rb}}}}}}$ (11)
${R_{{\rm{rb}}}} = {R_{{\rm{red}}}} - \gamma \left( {{R_{{\rm{blue}}}} - {R_{{\rm{red}}}}} \right)$ (12)

式中,RNIRRredRblue是近红外与红、蓝波段的反射率。对于RSP,波段5(865 nm)、4(670 nm)和2(470 nm)分别被用来作为RNIRRredRblueγ取决于气溶胶类型,对于大陆型气溶胶,其值为0.9。在本文中,ARVI值介于0.600和0.775之间的像元被标记为候选浓密植被像元。在沿传感器扫描方向上超过50%的像元被标记为候选浓密植被像元时,该地面点被标记为浓密植被,如图3所示。其中水平轴线表示传感器扫描方向,垂直轴线表示观测天顶角方向。图3(b)中的每一个绿色点代表该方向的ARVI在[0.600, 0.775]范围内。

图 3 RSP数据的假彩色影像与筛选得到的浓密植被点
Fig. 3 False color image of RSP data (R (670 nm), G (865 nm), B (470 nm)) and the selected dense vegetation

3.1.4 Rondeaux-Herman模型

在去除了大气影响和提取浓密植被区域之后,可以获取浓密植被的偏振反射率分布。Rondeaux与Herman指出,对于单位体积内的叶面积S(z)均匀并且土壤贡献可忽略的冠层而言,其偏振反射率分布表达式为

${\rho _{{\rm{sp}}}}\left( {{\theta _{\rm{i}}},{\varphi _{\rm{i}}};{\theta _{\rm{r}}},{\varphi _{\rm{r}}}} \right) = \frac{\text{π}}{4}R\left( \omega \right)\frac{{f\left( {{\theta _{\rm{n}}}} \right)}}{{G\left( {{\theta _{\rm{i}}}} \right){\mu _{\rm{r}}} + G\left( {{\theta _{\rm{r}}}} \right){\mu _{\rm{i}}}}}$ (13)

式中,坐标系和太阳—传感器几何定义如图4μ=cosθ(θi是入射天顶角,θr是反射天顶角),R(ω)是菲涅尔反射系数。f(θn)是LIDF,G(θ)是G函数(Ross,1981)。

图 4 坐标系与太阳—传感器几何定义
Fig. 4 The coordinate system and sun—sensor geometry definition

LIDF比较复杂,不容易精确估算。Goel和Strebel(1984)提出了双参数Beta分布模型来描述叶倾角分布。然而Beta分布需要叶倾角的平均值和方差,不易估算。Rondeaux和Herman(1991)则采用高斯函数替代之,但是这种估算方法精度低。在本文中,采用SAIL模型的“图形法”(Verhoef和Bunnik, 1981, Zhao 等,2010)确定两个参数去描述冠层结构,运用的两个参数为ab,其中a控制平均叶倾角,b反映分布形状。ab之间满足:

$\left| {{a}} \right| + \left| b \right| < 1$ (14)

平均叶倾角 $\overline {{\theta _n}} $

$\overline {{\theta _n}} = \frac{\text{π}}{4} - \frac{2}{\text{π}}a$ (15)

通过设置叶倾角分布函数中的参数,即可实现不同叶倾角分布情况下的植被偏振反射率模拟,在已知偏振反射率的情况下,也可通过数据拟合,实现叶倾角分布函数的反演。本节在分析RSP数据的基础上,进一步探讨不同叶倾角情况下的植被偏振反射率分布及叶倾角反演,将在3.2节中具体讨论。

3.2 结果与讨论

3.2.1 植被偏振反射特性

将NIR(865 nm),红(670 nm)和绿(555 nm)波段与2250 nm的偏振反射率对比分析(图5表3),可以发现不同波长偏振反射率之间高度相关,决定系数很接近于1.0,截距很小(在10–4量级)。对于给定的冠层结构和入射观测几何,折射率在 (400 nm,2500 nm)波段范围内有略微降低(Jacquemoud和Baret,1990),说明偏振特性和波长之间仅有微弱相关性。另外,由于大气校正在555 nm的不确定性大于865 nm和670 nm,斜率随着波长的减小而略微减小。

图 5 浓密植被在NIR(865 nm),红(670 nm)和绿(555 nm)波段的偏振反射率分别与2250 nm偏振反射率相关性分析
Fig. 5 Surface polarimetric reflectances over dense vegetation at NIR band(865 nm), red band(670 nm) and green band(555 nm) are regressed as a function of that of 2250 nm respectively

表 3 线性回归参数
Table 3 Parameters of the liner regression

下载CSV 
两个波段线性拟合
0.865 mm/2.25 μm 0.67 mm/2.25 μm 0.555 mm/2.25 μm
斜率 1.08 1.07 1.01
截距 3×10–4 3×10–4 –1×10–4
r2 0.96 0.97 0.96

偏振反射率具有的光谱不变特性,使得利用偏振手段研究植被冠层结构相对于常规光学遥感而言具有更大优势。常规光学传感器获取的BRF同时包含光谱与结构信息,即BRF为波长与结构的函数,二者耦合在一起,为研究冠层结构带来了困难;而利用偏振探测时,由于偏振反射率与光谱无关,则获取到的值仅为结构信息的函数,从而可以更容易地开展植被冠层结构研究。偏振的光谱不变特性在遥感反演中具有诸多优势,利用其反射率发掘其更有效的价值有待深化。

3.2.2 不同叶倾角分布情况下的偏振反射率模拟

为了模拟不同叶倾角分布情况下的偏振反射率,选择3种类型的植被:平展型,直立型和球形。表4所示即为描述这些分布的参数。

表 4 叶倾角分布模拟
Table 4 Simulations of leaf inclination distribution

下载CSV 
分布型 参数
a b ${\overline {{\theta }} }_{\rm n}$
平展型 0.8 0.1 74°
直立性 –0.8 0.1 16°
球形 –0.35 –0.15 58°
 注:ab为拟合参数。

基于式(4)的模拟结果表明,对于直立型,平展型和球型3种叶倾角分布类型,偏振反射率分布具有一定的相似性,但是均随入射—观测几何有显著的变化。图6(a)(c)为3种LIDF的植被偏振反射率模拟结果,白色虚线代表RSP的观测。可以看出:(1)植被偏振反射存在于上半球空间各个方向,叶片在冠层内部的不同取向具有不同的概率分布,从而将入射辐射以对应概率向各方向镜面反射。(2)偏振反射率随着散射相位角(太阳和传感器之间的角度)的增大而增加。(3)当太阳入射方向和观测方向一致时,偏振反射率达到最小值,即在热点方向没有偏振现象。图7表示NIR波段RSP的偏振反射率数据,其中的相对密度是与最大密度相比时的归一化密度。结果表明,在传感器观测平面上,偏振特性随观测天顶角变化,在观测平面各个方向都大于0,在接近最小相位角的过程中降低,然后随相位角的增大而增大。

图 6 直立型,平展型和球型植被的偏振反射率模拟
Fig. 6 Polarimetric reflectance simulations of three types of distribution: erectophile, planophile and spherical

3.2.3 植被叶倾角分布反演

本文研究区域与Knobelspiesse等人(2008)的研究相同,该区域被草地和农作物覆盖。本文结果同样表明该区域存在两种不同类型的植被。根据不同的密度分布,将图7所示数据分为两组结果如图8所示,其中虚线是Rondeaux模型对数据集拟合最佳时的结果。平均叶倾角 \normalsize$\overline {{\theta }}_{\rm{n}} $ 分别是56°(a= –0.3, b= –0.1)和38°(a=0.2, b= –0.3),表明存在两个包含不同叶倾角分布函数的不同植被类型,符合该场景的真实情况(Knobelspiesse 等,2008)。当观测天顶角为–40°时,模型拟合结果表明偏振反射率可分别达到0.0078和0.0125,偏振反射率差异较大。由于植被偏振特性对叶倾角分布十分灵敏,特别是在大相位角条件下,因此采用偏振测量方法进行地表分类是可行的。且借助叶倾角分布类型和其他辅助信息共同确定植被种类,则是一个更有应用价值的工作。这也将是日后研究工作的重点之一。

图 7 NIR波段(865 nm)的偏振特性随RSP观测天顶角的变化规律
Fig. 7 Polarimetric signatures of dense vegetation at NIR (865 nm) varies with view zenith angle from RSP
图 8 RSP数据在NIR波段的两种类型浓密植被的偏振特性随观测天顶角的变化规律
Fig. 8 Polarimetric signatures of dense vegetation at NIR (865 nm) varies with view zenith angle from RSP dataset

4 结 论

利用偏振方法系统分析了植被精细反演的两个要素,即冠层结构中光子—冠层的相互作用和叶倾角分布。在对分子与冠层结构相互作用的分析中,引入方向散射面积因子即结构因子DASF,用偏振观测方法去除冠层结构对散射透射再反射光的耦合作用,以及由此在植被散射系数上产生的误差;在对叶倾角分布的分析中,选择Rondeaux-Herman模型,利用偏振RSP数据集综合分析了植被偏振特性对叶倾角分布的灵敏性,采用SAIL模型中的双参数法来描述模型中的叶倾角分布函数,利用大气校正和ARVI准则提取浓密植被,进行植被偏振特性分析,并利用真实偏振探测数据获得了模型拟合结果。

主要结论有:(1)植被冠层与叶片表面反射特性对冠层BRF均有作用并产生误差,偏振探测能够敏感出这种变化导致的误差;若不利用偏振手段对此误差进行甄别,植被冠层的散射系数最高可产生140%的误差。考虑偏振观测对冠层立体结构的灵敏探测能力和对植被生化含量的高精度反演作用,有助于甄别植被氮含量与遥感反演正负相关的国际争议。(2)在多角度、多时相遥感中,偏振探测能够灵敏感知大角度观测中光辐射传输路径和强度的微弱变化,有助于剔除常规传感器大倾角地表探测的误差,从而提升遥感反演精度,因此相比于限制在±15°天顶角内观测的常规遥感探测而言,偏振探测可将传感器有效观测天顶角扩大到±40°—60°,为低约束条件下更大角度的遥感观测提供了可能。(3)在给定的入射观测几何条件下,偏振反射率在400—2500 nm范围内具有光谱不变特性,不同波段之间的偏振反射率决定系数>0.96,偏振的这种光谱不变特性可以有效分离植被反射的光谱与结构信息,从而可以利用偏振手段更好地研究植被冠层结构;由于叶片倾角的分布,冠层对于入射辐射的镜面反射存在于各个方向上,且偏振反射率随着相位角的增大而增大。(4)偏振对叶倾角分布的变化非常灵敏,特别是在大相位角条件下,这有助于获得更丰富的多角度植被反射信息,使得利用偏振遥感对物种进行分类和获取冠层结构参数成为可能。

参考文献(References)

  • Breon F M, Tanre D, Lecomte P and Herman M. 1995. Polarized reflectance of bare soils and vegetation: measurements and models. IEEE Transactions on Geoscience and Remote Sensing, 33 (2): 487–499. [DOI: 10.1109/36.377949]
  • Cairns B, Russell E E and Travis L D. 1999. Research scanning polarimeter: calibration and ground-based measurements. Proceedings of SPIE 3754, Polarization: Measurement, Analysis, and Remote Sensing II. Denver, CO, United States: SPIE: 186–196 [DOI: 10.1117/12.366329]
  • Curran P. 1981. The relationship between polarized visible light and vegetation amount. Remote Sensing of Environment, 11 : 87–92. [DOI: 10.1016/0034-4257(81)90011-0]
  • Curran P J. 1982. Polarized visible light as an aid to vegetation classification. Remote Sensing of Environment, 12 (6): 491–499. [DOI: 10.1016/0034-4257(82)90023-2]
  • Duan M Z and Lü D R. 2008. Simultaneously retrieving aerosol optical depth and surface albedo over land from POLDER’s multi-angle polarized measurement. Ⅱ: a case study. Chinese Journal of Atmospheric Sciences, 32 (1): 27–35. [DOI: 10.3878/j.issn.1006-9895.2008.01.03] ( 段民征, 吕达仁. 2008. 利用多角度POLDER偏振资料实现陆地上空大气气溶胶光学厚度和地表反照率的同时反演Ⅱ. 实例分析. 大气科学, 32 (1): 27–35. [DOI: 10.3878/j.issn.1006-9895.2008.01.03] )
  • Dubovik O, Herman M, Holdak A, Lapyonok T, Tanré D, Deuzé J L, Ducos F, Sinyuk A and Lopatin A. 2011. Statistically optimized inversion algorithm for enhanced retrieval of aerosol properties from spectral multi-angle polarimetric satellite observations. Atmospheric Measurement Techniques, 4 (5): 975–1018. [DOI: 10.5194/amt-4-975-2011]
  • Egan W G. 1970. Optical stokes parameters for farm cropidentification. Remote Sensing of Environment, 1 (3): 165–180. [DOI: 10.1016/S0034-4257(70)80023-2]
  • Fitch B W, Walraven R L and Bradley D E. 1984. Polarization of light reflected from grain crops during the heading growth stage. Remote Sensing of Environment, 15 (3): 263–268. [DOI: 10.1016/0034-4257(84)90036-1]
  • Goel N S and Strebel D E. 1984. Simple beta distribution representation of leaf orientation in vegetation canopies. Agronomy Journal, 76 (5): 800–802. [DOI: 10.2134/agronj1984.00021962007600050021x]
  • Grant L. 1987. Diffuse and specular characteristics of leaf reflectance. Remote Sensing of Environment, 22 (2): 309–322. [DOI: 10.1016/0034-4257(87)90064-2]
  • Grant L, Daughtry C S T and Vanderbilt V C. 1987. Variations in the polarized leaf reflectance of Sorghum bicolor . Remote Sensing of Environment, 21 (3): 333–339. [DOI: 10.1016/0034-4257(87)90016-2]
  • Grant L, Daughtry C S T and Vanderbilt V C. 1993. Polarized and specular reflectance variation with leaf surface features. Physiologia Plantarum, 88 (1): 1–9. [DOI: 10.1111/j.1399-3054.1993.tb01753.x]
  • Guo H, Gu X F, Xie D H, Yu T and Meng Q Y. 2014. A review of atmospheric aerosol research by using polarization remote sensing. Spectroscopy and Spectral Analysis, 34 (7): 1873–1880. [DOI: 10.3964/j.issn.1000-0593(2014)07-1873-08] ( 郭红, 顾行发, 谢东海, 余涛, 孟庆岩. 2014. 大气气溶胶偏振遥感研究进展. 光谱学与光谱分析, 34 (7): 1873–1880. [DOI: 10.3964/j.issn.1000-0593(2014)07-1873-08] )
  • Jacquemoud S and Baret F. 1990. PROSPECT: a model of leaf optical properties spectra. Remote Sensing of Environment, 34 (2): 75–91. [DOI: 10.1016/0034-4257(90)90100-Z]
  • Kaufman Y J and Tanré D. 1992. Atmospherically resistant vegetation index (ARVI) for EOS-MODIS. IEEE Transactions on Geoscience and Remote Sensing, 30 (2): 261–270. [DOI: 10.1109/36.134076]
  • Knobelspiesse K D, Cairns B, Schmid B, Román M O and Schaaf C B. 2008. Surface BRDF estimation from an aircraft compared to MODIS and ground estimates at the Southern Great Plains site. Journal of Geophysical Research, 113 (D20): D20105 [DOI: 10.1029/2008JD010062]
  • Knyazikhin Y, Lewis P, Disney M I, Mõttus M, Rautiainen M, Stenberg P, Kaufmann R K, Marshak A, Schull M A, Carmona P L, Vanderbilt V, Davis A B, Baret F, Jacquemoud S, Lyapustin A, Yang Y and Myneni R B. 2013a. Reply to Ollinger et al.: Remote sensing of leaf nitrogen and emergent ecosystem properties. Proceedings of the National Academy of Sciences of the United States of America, 110 (27): E2438 [DOI: 10.1073/pnas.1305930110]
  • Knyazikhin Y, Lewis P, Disney M I, Stenberg P, Mõttus M, Rautiainen M, Kaufmann R K, Marshak A, Schull M A, Carmona P L, Vanderbilt V, Davis A B, Baret F, Jacquemoud S, Lyapustin A, Yang Y and Myneni R B. 2013b. Reply to Townsend et al.: Decoupling contributions from canopy structure and leaf optics is critical for remote sensing leaf biochemistry. Proceedings of the National Academy of Sciences of the United States of America, 110 (12): E1075 [DOI: 10.1073/pnas.1301247110]
  • Knyazikhin Y, Martonchik J V, Myneni R B, Diner D J and Running S W. 1998. Synergistic algorithm for estimating vegetation canopy leaf area index and fraction of absorbed photosynthetically active radiation from MODIS and MISR data. Journal of Geophysical Research, 103 (D24): 32257–32275. [DOI: 10.1029/98JD02462]
  • Knyazikhin Y, Schull M A, Stenberg P, Mõttus M, Rautiainen M, Yang Y, Marshak A, Carmona P L, Kaufmann R K, Lewis P, Disney M I, Vanderbilt V, Davis A B, Baret F, Sacquemoud S, Lyapustin A and Myneni R B. 2013c. Hyperspectral remote sensing of foliar nitrogen content. Proceedings of the National Academy of Sciences of the United States of America, 110 (3): E185–E192. [DOI: 10.1073/pnas.1210196109]
  • Kotchenova S Y and Vermote E F. 2007. Validation of a vector version of the 6S radiative transfer code for atmospheric correction of satellite data. Part II. Homogeneous Lambertian and anisotropic surfaces. Applied Optics, 46 (20): 4455–4464. [DOI: 10.1364/AO.46.004455]
  • Litvinov P, Hasekamp O and Cairns B. 2011. Models for surface reflection of radiance and polarized radiance: comparison with airborne multi-angle photopolarimetric measurements and implications for modeling top-of-atmosphere measurements. Remote Sensing of Environment, 115 (2): 781–792. [DOI: 10.1016/j.rse.2010.11.005]
  • Maignan F, Bréon F M, Fédèle E and Bouvier M. 2009. Polarized reflectances of natural surfaces: spaceborne measurements and analytical modeling. Remote Sensing of Environment, 113 (12): 2642–2650. [DOI: 10.1016/j.rse.2009.07.022]
  • Myneni R B and Ross J. 1991. Photon-Vegetation Interactions: Applications in Optical Remote Sensing and Plant Ecology. Berlin: Springer
  • Nadal F and Breon F M. 1999. Parameterization of surface polarized reflectance derived from POLDER spaceborne measurements. IEEE Transactions on Geoscience and Remote Sensing, 37 (3): 1709–1718. [DOI: 10.1109/36.763292]
  • Ollinger S V, Reich P B, Frolking S, Lepine L C, Hollinger D Y and Richardson A D. 2013. Nitrogen cycling, forest canopy reflectance, and emergent properties of ecosystems. Proceedings of the National Academy of Sciences of the United States of America, 110 (27): E2437 [DOI: 10.1073/pnas.1304176110]
  • Ollinger S V, Richardson A D, Martin M E, Hollinger D Y, Frolking S E, Reich P B, Plourde L C, Katul G G, Munger J W, Oren R, Smith M L, Paw U K T, Bolstad P V, Cook B D, Day M C, Martin T A, Monson R K and Schmid H P. 2008. Canopy nitrogen, carbon assimilation, and albedo in temperate and boreal forests: functional relations and potential climate feedbacks. Proceedings of the National Academy of Sciences of the United States of America, 105 (49): 19336–19341. [DOI: 10.1073/pnas.0810021105]
  • Peltoniemi J I, Gritsevich M, Hakala T, Dagsson-Waldhauserová P, Arnalds Ó, Anttila K, Hannula H R, Kivekäs N, Lihavainen H, Meinander O, Svensson J, Virkkula A and de Leeuw G. 1994. Soot on snow experiment: bidirectional reflectance factor measurements of contaminated snow. The Cryosphere Discussions, 9 (3): 3075–3111.
  • Puttonen E, Suomalainen J, Hakala T and Peltoniemi J. 2009. Measurement of reflectance properties of asphalt surfaces and their usability as reference targets for aerial photos. IEEE Transactions on Geoscience and Remote Sensing, 47 (7): 2330–2339. [DOI: 10.1109/TGRS.2008.2010132]
  • Rondeaux G and Herman M. 1991. Polarization of light reflected by crop canopies. Remote Sensing of Environment, 38 (1): 63–75. [DOI: 10.1016/0034-4257(91)90072-E]
  • Ross J. 1981. The Radiation Regime and Architecture of Plant Stands. Netherlands: Springer
  • Sano I. 2004. Optical thickness and Angstrom exponent of aerosols over the land and ocean from space-borne polarimetric data. Advances in Space Research, 34 (4): 833–837. [DOI: 10.1016/j.asr.2003.06.039]
  • Shen Q, Li J S, Zhang B, Wu Y H, Zou L and Wu T X. 2016. Analyzing spectral characteristics of water involving in-situ multiangle polarized reflectance and extraction of water-leaving radiance. Spectroscopy and Spectral Analysis, 36 (10): 3269–3273. [DOI: 10.3964/j.issn.1000-0593(2016)10-3269-05] ( 申茜, 李俊生, 张兵, 吴艳红, 邹雷, 吴太夏. 2016. 水面原位多角度偏振反射率光谱特性分析与离水辐射提取. 光谱学与光谱分析, 36 (10): 3269–3273. [DOI: 10.3964/j.issn.1000-0593(2016)10-3269-05] )
  • Sun Z Q, Wu D, Lv Y F and Zhao Y S. 2016. Polarized reflectance factors of vegetation covers from laboratory and field: a comparison with modeled results. Journal of Geophysical Research, 122 (2): 1042–1065. [DOI: 10.1002/2016JD025892]
  • Suomalainen J, Hakala T, Peltoniemi J and Puttonen E. 2009a. Polarised multiangular reflectance measurements using the finnish geodetic institute field goniospectrometer. Sensors, 9 (5): 3891–3907. [DOI: 10.3390/s90503891]
  • Suomalainen J, Hakala T, Puttonen E and Peltoniemi J. 2009b. Polarised bidirectional reflectance factor measurements from vegetated land surfaces. Journal of Quantitative Spectroscopy and Radiative Transfer, 110 (12): 1044–1056. [DOI: 10.1016/j.jqsrt.2009.02.017]
  • Townsend P A, Serbin S P, Kruger E L and Gamon J A. 2013. Disentangling the contribution of biological and physical properties of leaves and canopies in imaging spectroscopy data. Proceedings of the National Academy of Sciences of the United States of America, 110 (12): E1074 [DOI: 10.1073/pnas.1300952110]
  • Vanderbilt V C and Grant L. 1985. Plant canopy specular reflectance model. IEEE Transactions on Geoscience and Remote Sensing, GE-23 (5): 722–730. [DOI: 10.1109/TGRS.1985.289390]
  • Vanderbilt V C, Grant L, Biehl L L and Robinson B F. 1985. Specular, diffuse, and polarized light scattered by two wheat canopies. Applied Optics, 24 (15): 2408–2418. [DOI: 10.1364/AO.24.002408]
  • Verhoef W and Bunnik N J J. 1981. Influence of crop geometry on multispectral reflectance determined by the use of canopy reflectance models. Proceedings of International Colloque on Signatures of Remotely Sensed Objects. Avignon, France: [s.n]
  • Wang H, Sun X B, Sun B, Liang T Q, Li C L and Hong J. 2014. Retrieval of aerosol optical properties over a vegetation surface using multi-angular, multi-spectral, and polarized data. Advances in Atmospheric Sciences, 31 (4): 879–887. [DOI: 10.1007/s00376-013-3100-5]
  • Waquet F, Cairns B, Knobelspiesse K, Chowdhary J, Travis L D, Schmid B and Mishchenko M I. 2009. Polarimetric remote sensing of aerosols over land. Journal of Geophysical Research, 114 (D1): D01206 [DOI: 10.1029/2008JD010619]
  • Waquet F, Léon J F, Cairns B, Goloub P, Deuzé J L and Auriol F. 2009. Analysis of the spectral and angular response of the vegetated surface polarization for the purpose of aerosol remote sensing over land. Applied Optics, 48 (6): 1228–1236. [DOI: 10.1364/AO.48.001228]
  • Woessner P and Hapke B. 1987. Polarization of light scattered by clover. Remote Sensing of Environment, 21 (3): 243–261. [DOI: 10.1016/0034-4257(87)90011-3]
  • Wu T X and Zhao Y S. 2005. The bidirectional polarized reflectance model of soil. IEEE Transactions on Geoscience and Remote Sensing, 43 (12): 2854–2859. [DOI: 10.1109/TGRS.2005.857905]
  • Xie D H, Cheng T H, Wu Y, Fu H, Zhong R F and Yu J. 2017. Polarized reflectances of urban areas: analysis and models. Remote Sensing of Environment, 193 : 29–37. [DOI: 10.1016/j.rse.2017.02.026]
  • Xie D H, Wang P J, Zhu Q J and Zhou H M. 2010. Modeling polarimetric BRDF of leaves surfaces. Spectroscopy and Spectral Analysis, 30 (12): 3324–3328. [DOI: 10.3964/j.issn.1000-0593(2010)12-3324-05] ( 谢东辉, 王培娟, 朱启疆, 周红敏. 2010. 单叶片偏振BRDF建模及参数反演. 光谱学与光谱分析, 30 (12): 3324–3328. [DOI: 10.3964/j.issn.1000-0593(2010)12-3324-05] )
  • Yang B, Knyazikhin Y, Lin Y, Yan K, Chen C, Park T, Choi S, Mõttus M, Rautiainen M, Myneni R B and Yan L. 2016. Analyses of impact of needle surface properties on estimation of needle absorption spectrum: case study with coniferous needle and shoot samples. Remote Sensing, 8 (7): 563 [DOI: 10.3390/rs8070563]
  • Yang B, Knyazikhin Y, Mõttus M, Rautiainen M, Stenberg P, Yan L, Chen C, Yan K, Choi S, Park T and Myneni R B. 2017a. Estimation of leaf area index and its sunlit portion from DSCOVR EPIC data: theoretical basis. Remote Sensing of Environment, 198 : 69–84. [DOI: 10.1016/j.rse.2017.05.033]
  • Yang B, Zhao H M and Chen W. 2017b. Semi-empirical models for polarized reflectance of land surfaces: intercomparison using space-borne POLDER measurements. Journal of Quantitative Spectroscopy and Radiative Transfer, 202 : 13–20. [DOI: 10.1016/j.jqsrt.2017.07.014]
  • Zhao F, Gu X F, Verhoef W, Wang Q, Yu T, Liu Q, Huang H G, Qin W H, Chen L F and Zhao H J. 2010. A spectral directional reflectance model of row crops. Remote Sensing of Environment, 114 (2): 265–285. [DOI: 10.1016/j.rse.2009.09.018]
  • Zhao H Y, Zhao H, Yan L and Zhao Y S. 2004. Model of reflection spectra of rock surface in 2π space. Acta Geologica Sinica, 78 (3): 843–847. [DOI: 10.1111/j.1755-6724.2004.tb00205.x]