偏振反射信息在植被遥感双向反射研究中的作用[PDF全文]
孙仲秋, 赵云升, 卢珊, 吕云峰
摘要: 本文基于菲涅尔反射公式,结合现有的偏振反射模型对单叶与植被多角度偏振测量结果与理论结合进行分析;通过研究发现,偏振反射在前向散射方向可以通过菲涅尔反射公式进行解释,但是,在后向散射方向则需要考虑其他物理机理。此外,传感器获取的植被偏振信息既可以作为一种“噪声”来剔除,剔除后在可见光波段将相对误差从原来的30%降低到20%以内,又可以作为额外有效的信息源表征植被的结构特征:通过模型参数的大小判断冠层形态的平整程度。本文可以作为植被偏振探测的系统化方法,并且给出植被固有的偏振反射效应规律,即植被越平展光滑,产生的偏振信息越多;同时也将偏振光遥感在植被监测中的有效性凸显出来,即偏振信息的剥离有助于提升双向反射模型的计算精度。
关键词: 偏振遥感     多角度遥感     植被     双向反射模型     镜面反射    
DOI: 10.11834/jrs.20188173    
收稿日期: 2018-04-08
中图分类号: TP702    文献标识码: A    
作者简介: 孙仲秋,1986年生,男,副教授,研究方向为偏振光遥感理论与应用。E-mail:sunzq465@nenu.edu.cn
基金项目: 国家自然科学基金(编号:41671347,41771362)
Function of polarization on the bidirectional reflectance factor of vegetation samples
SUN Zhongqiu, ZHAO Yunsheng, LU Shan, LYU Yunfeng
1. School of Geographical Science, Northeast Normal University, Changchun 130024, China
2. College of Urban and Environmental Science, Changchun Normal University, Changchun 130032, China
Abstract: In this study, we investigate the relationship between polarization and the structural property of vegetation covers to deepen our understanding of the physical mechanism of multiangle polarized light from natural surfaces, which is basic for describing the properties of Earth’s surface using polarized remote sensing. In this study, we explain the physical mechanism of the polarization properties of vegetation on the basis of the Fresnel equation and its derivative polarized reflectance models and compare them with the polarization measurement results of a single leaf and two vegetation covers, which are measured using a goniometer system. The Fresnel equation is effective for explaining the polarization of vegetation in forward-scattering directions, however, other theories should be considered in explaining polarization in backward-scattering directions. The polarization of vegetation covers can be considered as " noise” when we use the photometric signal, such as the separation of specular portion (which is computed by bidirectional polarized reflectance factor) from the total reflectance factor, which decreases the difference between BRF model results and measured results from 30% to 20%. Moreover, it possesses potentially useful information, such as the relationship between model parameter and vegetation cover roughness, which can characterize the structure of vegetation covers. Investigating the polarization of vegetation contributes to the understanding of the optical property of natural surfaces and potentially provides an additional and effective method for remote-sensing applications of vegetation covers. Our study also potentially provides a method for measuring the polarization of vegetation covers and demonstrates the physical rules of polarization in vegetation covers, such as the flatter the vegetation samples, the more polarization can be measured, and separating polarization from the total reflectance factor contributes to the improvement of the ability of current BRF models to simulate vegetation reflection. These results highlight the efficiency of polarized remote sensing on characterizing natural surfaces such as vegetation covers.
Key Words: polarized remote sensing    multi-angle    vegetation    bidirectional reflectance distribution model    specular reflection    
1、引 言

被动光学遥感主要是利用反射光与地表参数之间的关系实现反演与监测,光强度信息在遥感应用当中被广泛关注,但是偏振信息却很少用在研究地面自然目标的特性 (Talmage和Curran,1986)。20世纪70年代末至今,研究者们已经证明了偏振测量技术在大气探测,材料鉴定以及目标识别方面具有重要的科学价值 (Tyo 等,2006)。尤其是伴随着机载与星载探测器的出现,使偏振遥感技术在表征地表参数研究中变成了可能。其中,RSP (Research Scanning Polarimeter) (Cairn 等,1999)与POLDER (Polarization and Directionality of the Earth’s Reflectances) (Deschamps 等,1994) 是非常具有代表性的偏振探测器。虽然这些研究的主要目的是为了消除地表偏振对大气偏振遥感反演应用中的贡献,但是也为地表偏振特性的研究提供了依据。众多科学机构利用这些偏振信息结合地表偏振反射模型来描述与解释自然地物的偏振特性,其中主要包括植被,土壤以及积雪 (Breon 等,1995Nadal和Breon,1999Ottaviani 等,2015)。

关于地表偏振的研究主要得益于众多研究者使用地基平台测量过程中发现了偏振信息与自然地物特征属性之间的关系。其中,关于植被偏振特性的研究备受关注,早期的研究分别从叶片与植被冠层两个层面来测量植被的偏振特性(Curran,1981Grant,1987Vanderbilt和Grant,1985),测量结果说明了偏振信息与叶片和冠层的生物特性具有一定的联系:研究发现叶片的漫反射分量来源于叶片内部结构,偏振信息是由叶片表面所产生的镜面反射引起的,同时也发现,植被冠层的几何特征也会改变其偏振信息。上述结论帮助研究者们从经验上了解到植物叶片与植被冠层是如何产生偏振信息的(Vanderbilt 等,1985)。基于偏振是由镜面反射产生的假设,研究者们利用菲涅尔公式建立偏振解析模型对植被的偏振反射进行模拟计算。结果发现偏振反射模型可以很好的评估测量结果,并且具有分析植被冠层结构特征的能力 (Vanderbilt和Grant,1985)。而非偏振部分完全是由叶片内部产生的,因此具有反映植被生理特征的潜力 (Vanderbilt 等,1985)。

然而,在随后的研究中,星载偏振探测器的出现并没有促进偏振遥感技术在植被特征参数反演方面的研究 (Deschamps 等,1994)。一方面,星载偏振探测器获取的偏振信息主要是用来反演气溶胶参数,而将地面植被背景的偏振信息作为影响因素考虑 (Tanré 等,2011);另外一方面,受到大气偏振的影响,低分辨率的星载探测器很难将地面的偏振信息与植被的特征参数建立关系。即使在这样的背景下,研究者们指出仍然需要系统地开展植被多角度偏振特性的研究 (Litvinov 等,2011Maignan 等,2009),这不仅可以将植被对大气偏振的贡献解释清楚,而且能阐明植被自身的偏振传输过程,从而确定偏振遥感在植被监测中的作用。

与此同时,基于地面偏振测量平台对植被偏振反射特性的研究也一直被广泛开展 (Peltoniemi 等,2015Sun 等,2017Suomalainen 等,2009),这些测量结果在20世纪八、九十年代研究的基础上考虑了更多的光谱信息,使研究者们对植被多角度偏振光谱反射特性有了更多的了解。无论是机载还是地面平台开展的研究,都为利用偏振信息反演地物的特征参数提供了潜在的数据源。但是,由于偏振信息与植被特征参数之间的关系尚不清晰,仍然无法利用偏振信息反映出植被的特征参数。如果要量化偏振信息与植被特征参数之间的关系,应该以地基平台开展研究,例如:野外与实验室测量。这样才能获取植被最本质的偏振反射特性,进而反演植被的结构特征参数。此外,偏振信息在众多植被双向反射模型研究中并没有被单独考虑,那么它是否会对双向反射模型的计算结果有所影响。

本研究基于多角度偏振测量,分析了信息对植被反射模型的影响以及偏振信息对植被反射模型的影响以及偏振信息在植被结构特征参数表达方面的潜力,发挥偏振测量在植被遥感研究中的作用。

2、植被多角度偏振反射原理     (2.1) 偏振测量基本原理

自然地物反射过程所产生的偏振信息可以通过斯托克斯参量表达出来,其中地表反射过程的V分量通常可以忽略(Talmage和Curran,1986):

$\left\{ \begin{aligned} I = &({L_{0}} + {L_{90}} + {L_{45}} + {L_{135}})/2 \\ &Q = {L_{0}} - {L_{90}} \\ &U = {L_{45}} - {L_{135}} \\ &V = {L_{\rm{r}}} - {L_{\rm{l}}} \\ \end{aligned} \right.$ (1)

式中,Lx是偏振片透光轴在不同方向(Lr为右方向,Ll为左方向)获得的偏振反射辐射亮度值。本研究中主要考虑测量的偏振反射比BPRF(Bidirectional Polarized Reflectance Factor),偏振反射比值可以表示为在相同入射与探测环境下,目标偏振反射光与理想朗伯体反射光之间的比值 (Sun 等,2017b):

$\left\{ \begin{array}{l}{\rm{BPR}}{F_Q} =\displaystyle \frac{{{L_0} - {L_{90}}}}{L} =\displaystyle \frac{{{L_Q}}}{L}\\{\rm{BPR}}F =\displaystyle \frac{{{L_{\rm{p}}}}}{L} =\displaystyle \frac{{\sqrt {{Q^2} + {U^2}} }}{L}\end{array} \right.$ (2)

式中,L为理想漫反射白板的反射辐射亮度值,Lp为测量目标的偏振反射辐射亮度值,两者是在完全相同的环境中获取的。此外,双向反射系数BRF(Bidirectional Reflectance Factor)的定义可以表示为

${\rm{BRF}} = \frac{{{L_{{\rm{sample}}}}({\theta _{\rm{s}}}, {\theta _{\rm{v}}};{\phi _{\rm{s}}}, {\phi _{\rm{v}}})}}{{L({\theta _{\rm{s}}}, {\theta _{\rm{v}}};{\phi _{\rm{s}}}, {\phi _{\rm{v}}})}}$ (3)

式中,θs表示入射天顶角,θv为探测天顶角,φs为入射方位角,φv为探测方位角,Lsample为测量目标的反射辐射亮度值,这里面需要注意的是,双向反射系数与偏振反射系数都是与入射测量几何相关的参量。

    (2.2) 植被偏振反射理论

本研究中,将使用偏振反射模型对单叶与植被覆被的偏振信息进行计算。根据研究者们得出的结论,菲涅尔偏振反射系数(Fp(α, n))可以表征植被的偏振反射特性(Litvinov 等,2011Maignan 等,2009谢东海 等, 2012):

${F_{\rm{p}}}(\alpha, n) = 0.5\left({\frac{{{{\sin }^2}({\alpha _t} - \alpha)}}{{{{\sin }^2}({\alpha _t} + \alpha)}} - \frac{{{{\tan }^2}({\alpha _t} - \alpha)}}{{{{\tan }^2}({\alpha _t} + \alpha)}}} \right)$ (4)
$\sin {\alpha _t} = \sin \alpha /n$ (5)
$\begin{array}{l}\cos 2\alpha {\rm{ = }}\cos \xi = \cos {\theta _{\rm{s}}}\cos {\theta _{\rm{v}}} + \\\sin {\theta _{\rm{s}}}\sin {\theta _{\rm{v}}}\cos ({\phi _{\rm{s}}} - {\phi _{\rm{v}}})\end{array}$ (6)

式中,α为反射角,是相角的一半,ξ 表示相角,αt表示折射角,n代表植被的折射率,取值1.5。本文使用由Diner等人(2012)提出的偏振反射模型来计算单个叶片以及植被覆被的偏振反射比值,选择这个模型是由于该模型可以用来模拟平整表面所产生的偏振信息:

$\begin{array}{l}BPR{F_\lambda }({\theta _{\rm{s}}}, {\theta _{\rm{v}}};{\phi _{\rm{s}}}, {\phi _{\rm{v}}}) = {f_{\lambda, {\rm{volumetric}}}}({\theta _{\rm{s}}}, {\theta _{\rm{v}}};{\phi _{\rm{s}}}, {\phi _{\rm{v}}})+\\\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\; {f_{\lambda, {\rm{microfacet}}}}({\theta _{\rm{s}}}, {\theta _{\rm{v}}};{\phi _{\rm{s}}}, {\phi _{\rm v}})\end{array}$ (7)

式中,体散射部分表示为

$\begin{array}{l}{f_{\lambda, {\rm{volumetric}}}}({\theta _{\rm{s}}}, {\theta _{\rm{v}}};{\phi _{\rm{s}}}, {\phi _{\rm{v}}}) =\displaystyle \frac{{{a_\lambda }}}{{\text{π}}}[(\cos {\theta _{\rm{v}}} + \cos {\theta _{\rm{s}}})\\\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\cos {\theta _{\rm{v}}}\cos {\theta _{\rm{s}}}{]^{k - 1}}\exp (b\cos \gamma){D}\end{array}$ (8)

式中,aλbk为拟合参数,D为消偏矩阵,本文中等于1。第二部分可以表示为

$\begin{split}{f_{\lambda, {\rm{microfacet}}}}({\theta _{\rm{s}}}, - {\theta _{\rm{v}}}, {\phi _{\rm{v}}}) =& \zeta \frac{{f(\sigma, \beta)}}{{4\cos ({\theta _{\rm{s}}})\cos ({\theta _{\rm{v}}})\cos \beta }}\\&{M}(- m){PM}({m_0})\end{split}$ (9)
$\text{式中,}f(\sigma, \beta) = \frac{1}{{{\text{π }}{{\cos }^3}\beta 2{\sigma ^2}}}\exp \left({ - \frac{{{{\tan }^2}\beta }}{{2{\sigma ^2}}}} \right)$ (10)

式中,f(σ, β)表示小平面分布,式(9)中ζσ为需要拟合的参数。

$\cos \beta = \frac{{\cos {\theta _{\rm{v}}} + \cos {\theta _{\rm{s}}}}}{{2\cos \alpha }}$ (11)

式(9)中的菲涅尔反射矩阵P以及旋转矩阵M的详细介绍请参照Diner等人(2012)的文献。

    (2.3) 植被样品与测量过程

本研究中包括3个样本,一个单叶片和两个植株样本,如图1所示。其中为了模拟不同植被样本结构特征,图1(b)表示平展型植被样本,图1(c)表示直立性植被样本。其中,单叶为橡皮树叶片样本a(图1(a)),叶片平整,表面无毛,厚蜡质;样本b为豆叶九里香(图1(b)),冠层平整,叶片平展分布,无毛;样本c为菜豆树(图1(c)),叶片近似直立分布,无毛。选择这3种样品是由于它们为阔叶树,而且它们的叶片无毛,有较厚的蜡质层覆盖,在研究过程中不用考虑叶片绒毛引起的反射差异。

图 1 本研究中所使用的植被样品 Figure 1 Samples measured in this study

本研究在实验室内进行测量,所使用的多角度偏振测量系统为东北师范大学偏振光遥感实验室自行研发设计的装置(Sun 等,2014),测量当中入射天顶角为40与50°,方位角为15°间隔,探测天顶角间隔为10°,每个方位从0°开始,60°探测天顶角结束。偏振测量过程每个方向进行10次测量,随后求平均,用于表示样本的偏振信息。由于本研究中假设测量目标的反射以入射主平面呈对称分布,因此,测量方向的总数量为78个。测量示意图如图2所示,图2中180°方位表示前向散射探测方向,0°表示后向散射方向,实线为测量方向。进行测量时,为了保证测量面积相同,通过调整载物台,确保所有植被样品的表面与仪器的测量基准面相同。本研究将选择560 nm,670 nm以及865 nm的偏振与反射信息进行分析,这是由于RSP与POLDER传感器可以在这3个波段对地进行观测。

图 2 测量示意图 Figure 2 Measurement schematic diagram
3、偏振信息对植被双向反射模型结果的影响     (3.1) 植被反射模型

本部分主要研究双向反射模型在计算植被覆被反射特性中的效力,这是由于研究过程中发现植被b在可见光波段的反射具有明显的前向散射峰值。这归结于平展的叶片分布特性产生了镜面反射,从而将植被具有的后向热点效应所掩盖(Sandmeier 等,1998)。因此,本研究选择了两个常用的线性核驱动BRF反射模型来说明这些模型是否可以有效的计算出植被的反射信息。因为Ross-Roujean 模型和Ross-Li dense 模型分别被用在了POLDER和MODIS星载探测器上,更具有实际应用价值。其中,线性模型表示如下

${\rm{BRF}} = {k_0} + {k_1}{F_1}({\theta _{\rm{s}}}, {\theta _{\rm{v}}}, \phi) + \cdots + {k_n}{F_n}({\theta _{\rm{s}}}, {\theta _{\rm{v}}}, \phi)$ (12)

式中,kn为反演参数,Fn表示核,一般研究中,线性模型中包含3个参数,即n=2。本研究所使用的第一个模型为Ross-Roujean模型(Roujean 等,1992),表达式为

$\begin{array}{l}{F_1} =\displaystyle \frac{1}{{2{\text{π }}}}\left[ {({\text{π }} - \phi)\cos \phi + \sin \phi } \right]\tan {\theta _{\rm{s}}}\tan {\theta _{\rm{v}}} - \\\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\displaystyle\frac{1}{{\text{π }}}\left[ {\tan {\theta _{\rm{s}}} + \tan {\theta _{\rm{v}}} + \Delta ({\theta _{\rm{s}}}, {\theta _{\rm{v}}}, \phi)} \right]\end{array}$ (13)
$\Delta ({\theta _{\rm{s}}}, {\theta _{\rm{v}}}, \phi) = \sqrt {{{\tan }^2}{\theta _{\rm{s}}} + {{\tan }^2}{\theta _{\rm{v}}} - 2\tan {\theta _{\rm{s}}}\tan {\theta _{\rm{v}}}\cos \phi } $ (14)
${F_2} = \frac{{(0.5{\text{π }} - \xi)\cos \xi + \sin \xi }}{{\cos {\theta _{\rm{s}}} + \cos {\theta _{\rm{v}}}}} - \frac{{\text{π}}}{4}$ (15)

第二个模型为Ross-Li致密核模型(Li和Strahler,1992),与Ross-Roujean不同的是,两个模型的F1核有所差异,Li致密核可以表示为

${F_1} = \frac{{(1 + \cos \theta ')\sec \theta _{\rm{s}}'\sec \theta _{\rm{v}}'}}{{\sec \theta _{\rm{s}}' + \sec \theta _{\rm{v}}' - O({\theta _{\rm{s}}}, {\theta _{\rm{v}}}, \phi)}}$ (16)
$O({\theta _{\rm{s}}}, {\theta _{\rm{v}}}, \phi) = \frac{1}{{\text{π }}}(t - \sin t\cos t)(\sec \theta _{\rm{s}}' + \sec \theta _{\rm{v}}')$ (17)
$\cos t = \frac{h}{b}\frac{{\sqrt {{D^2} + {{(\tan \theta _{\rm{s}}'\tan \theta _{\rm{v}}'\sin \phi)}^2}} }}{{\sec \theta _{\rm{s}}' + \sec \theta _{\rm{v}}'}}$ (18)
$D = \sqrt {{{\tan }^2}\theta _{\rm{s}}' + {{\tan }^2}\theta _{\rm{v}}' - 2\tan \theta _{\rm{s}}'\tan \theta _{\rm{v}}'\cos \phi } $ (19)
$\cos {\theta '} = \cos \theta _{\rm{s}}'\cos \theta _{\rm{v}}' + \sin \theta _{\rm{s}}'\sin \theta _{\rm{v}}'\cos \phi $ (20)
$\theta _{\rm{s}}' = {\tan ^{ - 1}}\left({\frac{b}{r}\tan {\theta _{\rm{s}}}} \right), \theta _{\rm{v}}' = {\tan ^{ - 1}}\left({\frac{b}{r}\tan {\theta _{\rm{v}}}} \right)$ (21)

致密核中包含两个参数h/bb/r,通常赋值为2和1。

    (3.2) 植被覆被反射特性

图3给出了不同植被覆被在670 nm处BRF的分布图,入射天顶角为50°。从图3中可以发现样本b与样本c的反射分布特性具有显著差异。其中,具有平展型的样本b具有明显的前向散射峰值,如图3(a)所示;而直立型冠层的样本c则具有典型的后向反射峰值,如图3(b)所示。以往的研究已经表明(Sun 等,2017a, b),由于镜面反射在670 nm波段的贡献大,平展型植被在该波段处的反射出现前向反射峰值,因此会出现图3(a)的结果。随后,本研究利用两个BRF模型来模拟平展型与直立型植被覆被的双向反射信息。表1表2分别给出了样本b与样本c在不同波段所反演的模型参数,以及均方根误差(RMSE),同时包括模拟值与测量之间的平均差异的绝对值(ARD) (Sun 等,2017a)。其中,ARD是用来评价模型有效性的一个重要指标,根据Goel(1988)的研究表明,如果有一个模型的模拟结果与测量值之间的差异在20%以内,那么这个模型是有效的。从表1中的结果可以发现,两个常用的BRF模型在670 nm波段的计算结果与测量结果之间的差异达到了30%,而在其他波段的差异均在20%之内;表2中样本c在各个波段的测量结果与模拟值之间的差异均小于20%。这些结果表明,在670 nm波段处,BRF模型不能有效地计算样本b的反射值。出现这种结果是由于样本b在670 nm是以镜面反射按为主,而本研究所使用的BRF模型在计算植被反射过程中没有考虑镜面反射 (Sun 等,2017a)。因此,如果能将镜面反射从总反射当中分离,则可以利用现有的BRF模型对漫反射部分信息进行计算。

图 3 不同植被覆被在670 nm的反射分布图 Figure 3 Distribution of BRF at 670 nm for different vegetation covers

表 1 植被b在不同波段反演的最优模型参数,入射天顶角为50° Table 1 Best fitted parameters of vegetation b at selected wavelengths while the incident zenith angle was 50°

表 2 植被c在不同波段反演的最优模型参数,入射天顶角为50° Table 2 Best fitted parameters of vegetation c at selected wavelengths while the incident zenith angle was 50°
    (3.3) 偏振信息计算镜面反射

由于偏振反射主要是由镜面反射引起的,因此,研究者们提出利用偏振反射信息计算镜面反射。首先,将总的双向反射分成两部分,镜面部分BRFspecular与漫反射部分BRFdiffuse

${\rm BRF} = {\rm BR{F}_{{\rm{specular}}}} + {\rm BR{F}_{{\rm{diffuse}}}}$ (22)

镜面反射部分则可以根据测量的偏振信息与式(23)计算得到(Vanderbilt 等,1985)

$F(\alpha, n){\rm{ = }}\frac{{{{\sin }^2}({\alpha _t} - \alpha){{\tan }^2}({\alpha _t} + \alpha) - {{\sin }^2}({\alpha _t} + \alpha){{\tan }^2}({\alpha _t} - \alpha)}}{{{{\sin }^2}({\alpha _t} - \alpha){{\tan }^2}({\alpha _t} + \alpha) + {{\sin }^2}({\alpha _t} + \alpha){{\tan }^2}({\alpha _t} - \alpha)}}$ (23)

随后,将计算得到的镜面反射从总反射中减去,得到漫反射部分;最后,利用BRF模型对漫反射部分进行模拟计算,得到的结果如表3所示 (Sun 等,2017a)。通过表3可以发现,利用现有的BRF模型可以很好的模拟植被b的漫反射部分。其中在670 nm处的ARD已经从原来的30%降低到20%以内,而在其他波段,模型计算值与测量之间的差异也有不同程度的降低。这说明,将镜面反射从总反射中排出,有助于提高现有双向反射模型对植被覆被反射特性计算的精度。

表 3 植被b在不同波段漫反射部分反演的最优模型参数,入射天顶角为50° Table 3 Best fitted parameters of vegetation b inverted by the diffuse portion of BRF at selected wavelengths while the incident zenith angle was 50°
4、植被结构特征与偏振反射比之间关系     (4.1) 不同植被的偏振反射特性

大部分研究者用线偏振度作为指标来描述植被的偏振反射特性,并提出利用偏振度可以不需要对测量结果进行归一化处理。这是由于偏振度是偏振反射比与双向反射比的比值,而两个分量的定义都是用标准白板进行归一化计算的,所以在直接计算线偏振度时可以忽略归一化过程。为了说明本论文中偏振信息的准确性,图4给出了不同植被样品的线偏振度光谱曲线,该结果是在前向散射方向主平面测量的,入射天顶角为40°。从图4中可以发现,植被样品的偏振度曲线与植被反射光谱曲线呈现反比例关系(Breon 等,1995Maignan 等,2009),这一现象被研究者们解释为偏振反射比与波长无关(王涵 等,2014相坤生 等,2015谢东海 等,2012),所以低反射波段对应高的偏振度,高反射波段对应低偏振度值。为了进一步说明不同植被样品的偏振差异,在图5中给出了它们的偏振度分布特性,入射天顶角为40°。在图5中表示方位角过程中,180表示前向散射方向,0表示后向散射方向。从图5中可以发现,所有样品的最大偏振均出现在前向散射方向,这一规律可以通过式(4)的菲涅尔偏振反射系数来解释。然而,在图5中发现,最大偏振度值不是出现在单叶的测量结果中,而是样本b的偏振度最大。这进一步说明,偏振度是两个反射系数的比值,它不能完全代表植被的偏振特性。因此,这里用BPRF来表示植被样本的偏振特性,一方面由于BPRF是直接获取的偏振信息;另一方面现有的偏振反射模型可以直接用来计算植被的BPRF (Diner 等,2012)。

图 4 不同植被样品的偏振度光谱曲线(−Q/I) Figure 4 Spetrum of dolp (degree of linear polarization) of vegetation samples (−Q/I)

图 5 不同植被样品在670 nm处的偏振度分布图 Figure 5 Polar plot of dolp of vegetation samples at 670 nm
    (4.2) 植被结构与偏振反射比的关系

本研究中将植被样品的粗糙度按照逐渐增加的程度依次分为样本a、样本b和样本c,进而分析植被结构特征与偏振反射比之间的关系。同时利用Diner等人(2012)的模型对3个样本进行模拟计算,进一步说明菲涅尔反射系数可以要来表征植被样品的偏振特性。反演得到的模型参数如表4所示。在这里本研究没有给出偏振反射模型计算结果与测量结果的相对差异,这是由于偏振反射比值本身较小,一部分BPRF值在10–3左右。非常小的差异值就会引起较大的相对差异。因此也就意味着RMSE越小,说明偏振反射模型与测量的BPRF值之间的拟合程度越高。在偏振反射模型当中,包含一个粗糙度参数σ,通过分析表4的结果可以发现,在波段相同的情况下,粗糙度参数σ增加的顺序为单叶,平展植被(b),直立型植被(c)。这表明,粗糙度参数可以在一定程度上反映出植被样品的结构特征,即σ值越大,植被冠层越粗糙,或者说叶片越接近直立状态。这是由于叶片的分布状态会直接影响其产生镜面反射的多少,从而影响偏振信号的大小。为了进一步说明偏振反射比的产生于植被样品的结构特征有关,图6给出了不同植被样品在865 nm处BPRFQ(式(2))的分布特性,与图5的方位相同,这里面考虑了Q参量自身的符号,在前向为负,后向为正。选择BPRFQ是由于Q参量可以表示植被反射过程中所产生的大部分偏振信息,另一方面Q参量正负值的变化可以体现出植被偏振信息或者以垂直分量为主,或者以平行分量为主。从图6中可以发现,无论是单叶还是植被覆被,均在所有方向都产生了偏振反射比值,其中,图6(b)(c)共用一个图例。根据以往的研究结果表明,平整的单叶除了前向散射方向之外,在其他方向不会有镜面反射发生(Bousquet 等,2005)。那么没有镜面反射的区域所产生的偏振信息则是由其他物理机制所产生的。近期的研究表明,多重散射也可以产生偏振信息 (Kallel和Gastellu-Etchegorry,2017)。这里需要注意的是,图6(a)是植被单叶在865 nm处的偏振反射比结果,在非前向散射方向的反射值都是由叶片内部所产生的多重散射所引起的,本研究虽然没有利用模型计算多重散射结果,但是也从实验结果角度说明:除了镜面反射之外,多重散射也是产生偏振的一个原因,而多重散射产生的偏振是以平行分量为主的。

表 4 不同植被拟合Diner偏振反射模型的参数结果,入射天顶角为40° Table 4 Best fitted parameters of Diner model for vegetation samples while the incident zenith angle was 40°

图 6 不同植被样品在865 nm处BPRFQ的分布图,入射天顶角为40° Figure 6 Polar plot of BPRFQ for different vegetation samples at 865 nm while the incident zenith angle was 40°
5、结 论

本文系统地分析了地表偏振光多角度反射物理效应,并且研究了不同植被样品的多角度偏振反射信息对双向反射的影响,以及利用偏振信息反映植被的结构特征参数:植被覆被的粗糙度。研究发现,多角度偏振反射可以根据菲涅尔偏振反射理论进行解释,并且可以直接利用偏振反射信息或者将其从总反射中移除。通过本研究总结出以下结论,可以作为植被偏振测量方法以及植被偏振特性规律研究所参考。

首先,地表测量是详细了解植被偏振特性的最有效方法,因为可以获得植被最本质的偏振反射信息;这不仅对建立偏振信息与植被参数之间的关系至关重要,同时也可以为下垫面对大气偏振研究提供支持。在测量过程中应在尽可能多的方向均匀获取植被偏振信息,这样才能完整地表达偏振分布特性。

其次,植被偏振反射比可以用来计算镜面反射部分,如果将这部分镜面反射部分从总反射中减去,那么现有的双向反射模型可以以更高的精度来模拟漫反射部分。这种提升对镜面反射为主的植被覆被具有明显的效果。同时,也为未来的研究提供了一个参考,即,把偏振反射比当成一种“噪声”排出来提升植被遥感的精度,可以将BRF模型计算误差从原来的30%减少到20%,从而使模型在可见光波段有效的计算植被反射信息。

最后,研究中发现,偏振反射比还可以直接应用在对植被的研究当中,因为偏振反射模型可以有效的计算植被的偏振反射比,同时,不同粗糙度的植被覆被可以通过模型参数的大小来表示出来。但是这一结论需要后续对更多类型植被偏振特性的研究来证明。这些结果表明,偏振反射比本身也具有研究植被的结构特征参数的潜能。

本研究不仅证明了偏振信息自身具有的独特能力,而且还可以把偏振信息排出,提升基于反射强度信息可为主的研究。这进一步说明偏振遥感作为一种非常有效的探测手段,可以用于地表植被监测研究。从地表研究来看,偏振探测技术已经发挥了其特有的作用,但是在未来机载与星载应用中,还需要研究者们开展大量的研究来实现偏振遥感在地表监测中的应用,因为像元尺度上植被异质性以及与地面测量范围的差异会对反射信息产生影响。本研究可以为后续关于植被偏振遥感测量以及规律的研究提供一些参考,同时也说明了偏振遥感技术是遥感领域中非常有效的一种反演与估算植被特征参数的手段。

参考文献
[1] Bousquet L, Lacherade S, Jacquemoud S and Moya I. Leaf BRDF measurements and model for specular and diffuse components differentiation[J]. Remote Sensing of Environment, 2005, 98 (2/3) : 201 –211. DOI: 10.1016/j.rse.2005.07.005
[2] Breon F M, Tanre D, Lecomte P and Herman M. Polarized reflectance of bare soils and vegetation: measurements and models[J]. IEEE Transactions on Geoscience and Remote Sensing, 1995, 33 (2) : 487 –499. DOI: 10.1109/36.377949
[3] Cairns B, Russell E E and Travis L D. 1999. Research Scanning Polarimeter: calibration and ground-based measurements//Proceeding of SPIE, Polarization: Measurement, Analysis, and Remote Sensing II. Denver, CO, United States: SPIE: 186–196 [DOI: 10.1117/12.366329]
[4] Curran P. The relationship between polarized visible light and vegetation amount[J]. Remote Sensing of Environment, 1981, 11 : 87 –92. DOI: 10.1016/0034-4257(81)90011-0
[5] Deschamps P Y, Bréon F M, Leroy M, Podaire A, Bricaud A, Buriez J C and Seze G. The POLDER mission: instrument characteristics and scientific objectives[J]. IEEE Transactions on Geoscience and Remote Sensing, 1994, 32 (3) : 598 –615. DOI: 10.1109/36.297978
[6] Diner D J, Xu F, Martonchik J V, Rheingans B E, Geier S, Jovanovic V M, Davis A, Chipman R A and McClain S C. Exploration of a polarized surface bidirectional reflectance model using the ground-based multiangle spectropolarimetric imager[J]. Atmosphere, 2012, 3 (4) : 591 –619. DOI: 10.3390/atmos3040591
[7] Goel N S. Models of vegetation canopy reflectance and their use in estimation of biophysical parameters from reflectance data[J]. Remote Sensing Reviews, 1988, 4 (1) : 1 –212. DOI: 10.1080/02757258809532105
[8] Grant L. Diffuse and specular characteristics of leaf reflectance[J]. Remote Sensing of Environment, 1987, 22 (2) : 309 –322. DOI: 10.1016/0034-4257(87)90064-2
[9] Kallel A and Gastellu-Etchegorry J P. Canopy polarized BRDF simulation based on non-stationary Monte Carlo 3-D vector RT modeling[J]. Journal of Quantitative Spectroscopy and Radiative Transfer, 2017, 189 : 149 –167. DOI: 10.1016/j.jqsrt.2016.11.006
[10] Li X and Strahler A H. Geometric-optical bidirectional reflectance modeling of the discrete crown vegetation canopy: effect of crown shape and mutual shadowing[J]. IEEE Transactions on Geoscience and Remote Sensing, 1992, 30 (2) : 276 –292. DOI: 10.1109/36.134078
[11] Litvinov P, Hasekamp O and Cairns B. Models for surface reflection of radiance and polarized radiance: comparison with airborne multi-angle photopolarimetric measurements and implications for modeling top-of-atmosphere measurements[J]. Remote Sensing of Environment, 2011, 115 (2) : 781 –792. DOI: 10.1016/j.rse.2010.11.005
[12] Maignan F, Bréon F M, Fédèle E and Bouvier M. Polarized reflectances of natural surfaces: spaceborne measurements and analytical modeling[J]. Remote Sensing of Environment, 2009, 113 (12) : 2642 –2650. DOI: 10.1016/j.rse.2009.07.022
[13] Nadal F and Breon F M. Parameterization of surface polarized reflectance derived from POLDER spaceborne measurements[J]. IEEE Transactions on Geoscience and Remote Sensing, 1999, 37 (3) : 1709 –1718. DOI: 10.1109/36.763292
[14] Ottaviani M, van Diedenhoven B and Cairns B. Photopolarimetric retrievals of snow properties[J]. The Cryosphere, 2015, 9 (5) : 1933 –1942. DOI: 10.5194/tc-9-1933-2015
[15] Peltoniemi J I, Gritsevich M and Puttonen E. 2015. Reflectance and polarization characteristics of various vegetation types//Kokhanovsky A A, ed. Light Scattering Reviews 9: Light Scattering and Radiative Transfer. Berlin, Heidelberg: Springer: 257–294 [DOI: 10.1007/978-3-642-37985-7_7]
[16] Roujean J L, Leroy M and Deschamps P Y. A bidirectional reflectance model of the Earth's surface for the correction of remote sensing data[J]. Journal of Geophysical Research, 1992, 97 (D18) : 20455 –20468. DOI: 10.1029/92JD01411
[17] Sandmeier S, Müller C, Hosgood B and Andreoli G. Physical mechanisms in hyperspectral BRDF data of grass and watercress[J]. Remote Sensing of Environment, 1998, 66 (2) : 222 –233. DOI: 10.1016/S0034-4257(98)00060-1
[18] Sun Z Q, Huang Y H, Bao Y L and Wu D. Polarized remote sensing: a note on the Stokes parameters measurements from natural and man-made targets using a spectrometer[J]. IEEE Transactions on Geoscience and Remote Sensing, 2017, 55 (7) : 4008 –4021. DOI: 10.1109/TGRS.2017.2686485
[19] Sun Z Q, Wu D, Lv Y F and Zhao Y S. Bidirectional polarized reflectance factors of vegetation covers: influence on the BRF models results[J]. IEEE Transactions on Geoscience and Remote Sensing, 2017a, 55 (10) : 5687 –5701. DOI: 10.1109/TGRS.2017.2712137
[20] Sun Z Q, Wu D, Lv Y F and Zhao Y S. Polarized reflectance factors of vegetation covers from laboratory and field: a comparison with modeled results[J]. Journal of Geophysical Research: Atmospheres, 2017b, 122 (2) : 1042 –1065. DOI: 10.1002/2016JD025892
[21] Sun Z Q, Wu Z F and Zhao Y S. Semi-automatic laboratory goniospectrometer system for performing multi-angular reflectance and polarization measurements for natural surfaces[J]. Review of Scientific Instruments, 2014, 85 : 014503 . DOI: 10.1063/1.4862814
[22] Suomalainen J, Hakala T, Puttonen E and Peltoniemi J. Polarised bidirectional reflectance factor measurements from vegetated land surfaces[J]. Journal of Quantitative Spectroscopy and Radiative Transfer, 2009, 110 (12) : 1044 –1056. DOI: 10.1016/j.jqsrt.2009.02.017
[23] Talmage D A and Curran P J. Remote sensing using partially polarized light[J]. International Journal of Remote Sensing, 1986, 7 (1) : 47 –64. DOI: 10.1080/01431168608954660
[24] Tanré D, Bréon F M, Deuzé J L, Dubovik O, Ducos F, François P, Goloub P, Herman M, Lifermann A and Waquet F. Remote sensing of aerosols by using polarized, directional and spectral measurements within the A-Train: the PARASOL mission[J]. Atmospheric Measurement Techniques, 2011, 4 (7) : 1383 –1395. DOI: 10.5194/amt-4-1383-2011
[25] Tyo J S, Goldstein D L, Chenault D B and Shaw J A. Review of passive imaging polarimetry for remote sensing applications[J]. Applied Optics, 2006, 45 (22) : 5453 –5469. DOI: 10.1364/AO.45.005453
[26] Vanderbilt V C and Grant L. Plant canopy specular reflectance model[J]. IEEE Transactions on Geoscience and Remote Sensing, 1985, GE-23 (5) : 722 –730. DOI: 10.1109/TGRS.1985.289390
[27] Vanderbilt V C, Grant L, Biehl L L and Robinson B F. Specular, diffuse, and polarized light scattered by two wheat canopies[J]. Applied Optics, 1985, 24 (15) : 2408 –2418. DOI: 10.1364/AO.24.002408
[28] 王涵, 孙晓兵, 孙斌, 洪津. 基于航空多角度偏振辐射计遥感数据评估陆地表面偏振反射模型[J]. 光学学报, 2014, 34 (1) : 0128002 . Wang H, Sun X B, Sun B and Hong J. Evaluation of land surface polarization models based on airborne advanced atmosphere multi-angle polarization radiometer measurements[J]. Acta Optica Sinica, 2014, 34 (1) : 0128002 .
[29] 相坤生, 程天海, 顾行发, 郭红, 陈好, 王颖, 魏曦, 包方闻. 基于多角度偏振载荷数据的中国典型地物偏振特性研究[J]. 物理学报, 2015, 64 (22) : 0227801 . Xiang K S, Cheng T H, Gu X F, Guo H, Chen H, Wang Y, Wei X and Bao F W. Polarized properties of typical surface types over China based on the multi-angular polarized remote sensing measurements[J]. Acta Physica Sinica, 2015, 64 (22) : 0227801 . DOI: 10.7498/aps.64.227801
[30] 谢东海, 顾行发, 程天海, 余涛, 李正强, 陈兴峰, 陈好, 郭婧. 基于多角度偏振相机的城市典型地物双向反射特性研究[J]. 物理学报, 2012, 61 (7) : 077801 . Xie D H, Gu X F, Cheng T H, Yu T, Li Z Q, Chen X F, Chen H and Guo J. Research on the bidirectional reflectance of typical urban surface types measured by the directional polarimetric camera[J]. Acta Physica Sinica, 2012, 61 (7) : 077801 . DOI: 10.7498/aps.61.077801