三维辐射传输模型LESS原理及其应用

doi:10.11873/j.issn.1004-0323.2019.5.0914

三维辐射传输模型LESS原理及其应用

漆建波^,¹^,², 谢东辉², 许月², 阎广建²

1.北京林业大学省部共建森林培育与保护教育部重点实验室，北京 100083

2.北京师范大学地理科学学部，北京 100875

Principles and Applications of the 3D Radiative Transfer Model LESS

Qi Jianbo^,¹^,², Xie Donghui², Xu Yue², Yan Guangjian²

1.The Key Laboratory for Silviculture and Conservation of Ministry of Education，Beijing Forestry University，Beijing 100083，China

2.Faculty of Geographical Science，Beijing Normal University，Beijing 100875，China

收稿日期: 2019-09-11 修回日期: 2019-11-07

基金资助:

国家自然科学基金项目“植被遥感中多尺度复杂场景反射辐射的计算机模拟模型研究”. 41571341

Received: 2019-09-11 Revised: 2019-11-07

作者简介 About authors

漆建波（1989-），四川德阳人，博士，讲师，主要从事植被定量遥感研究Email:jianboqi@bjfu.edu.cn , E-mail：jianboqi@bjfu.edu.cn

摘要

三维辐射传输模型能够准确地刻画太阳辐射与异质地表之间的相互作用，近年来已成为定量遥感建模与反演研究中的重要工具。LESS模型是基于光线追踪的三维真实冠层辐射传输模型，充分利用了光线追踪的前向追踪模式模拟能量平衡问题以及后向追踪模式模拟大尺度（公里级）遥感影像，从而实现在同一模型中多种遥感数据的模拟。目前，LESS模型可以模拟多角度反射率、多/高光谱影像、鱼眼相机、复杂地形区上下行短波辐射、冠层分层FPAR等数据，可以为验证物理模型、发展参数化模型以及训练神经网络模型等提供更为可靠的模拟数据集。本文主要介绍了LESS模型的基本原理和典型的应用。LESS模型可以从www.lessrt.org网站下载。

关键词： LESS ; 光线追踪 ; 三维辐射传输 ; 真实场景

Abstract

Three-dimensional （3D） radiative transfer model can accurately describe the interactions between solar radiation and heterogeneous land surfaces. Recently， it has become an important tool for quantitative remote sensing studies. LESS is a ray-tracing based 3D radiative transfer model， which take full advantage of the forward ray-tracing techniques for simulating radiative budget and backward ray-tracing for simulating large-scale images， which makes it possible to simulate various remote sensing data in a single model. Currently， LESS can simulate multi-angle Bidirectional Reflectance Factor （BRF）， multi-spectral/high-spectral images， fish-eye cameras， upwelling/downwelling shortwave radiation in rugged terrains and layered FPAR， etc. This simulated dataset can be used for validating physical modes， developing parameterized models， as well as training neural networks. This paper presents the fundamentals of LESS and its applications. LESS can be downloaded from www.lessrt.org.

Keywords： LESS ; Ray tracing ; 3D Radiative transfer ; Realistic scene

PDF (6161KB) 元数据多维度评价相关文章导出 EndNote| Ris| Bibtex 收藏本文

本文引用格式

漆建波, 谢东辉, 许月, 阎广建. 三维辐射传输模型LESS原理及其应用. 遥感技术与应用[J], 2019, 34(5): 914-924 doi:10.11873/j.issn.1004-0323.2019.5.0914

Qi Jianbo, Xie Donghui, Xu Yue, Yan Guangjian. Principles and Applications of the 3D Radiative Transfer Model LESS. Remote Sensing Technology and Application[J], 2019, 34(5): 914-924 doi:10.11873/j.issn.1004-0323.2019.5.0914

１引言

在植被定量遥感中，辐射传输建模是理解冠层与太阳辐射之间相互作用机制以及从遥感信号反演植被参数的重要基础。然而由于地物的复杂性，难以建立起遥感信号与植被参数之间的简单数学关系，因此传统遥感模型往往对地物的结构特征进行简化和抽象描述。例如，将植被冠层简化为水平均匀分布的叶片，从而发展了一系列的一维辐射传输模型^［1-2］；或者将树冠抽象为简单的几何体而发展了几何光学模型^［3-4］。但是，这些模型往往难以描述精细的植被结构，因此近年来能够准确描述空间异质性的三维辐射传输模型受到了越来越多的重视。

此外，随着高分辨率三维数据（例如激光雷达数据）的丰富，如何充分利用这些信息并与辐射传输理论相结合来反演更高精度的地表参量成为定量遥感研究中的一个新问题，并且已得到国内外相关学者的关注^［5-6］。由于三维辐射传输模型可以描述任意复杂的场景结构，能够更大限度地利用多源遥感数据所提供的信息量，因此三维辐射传输模型正是解决这一问题的理想工具。

在过去20多年中，遥感领域已经发展了一系列的三维辐射传输模型，例如基于辐射度算法的RGM模型^［7］和RAPID模型^［8］、基于光线追踪的Raytran模型^［9］以及基于离散坐标法的DART模型^［10］等。从现有模型来看，光线追踪算法由于其良好的扩展性，在三维辐射传输模型中使用的频率更高，并且已经得到了较多的应用^［11-12］，但光线追踪算法本身又有多种实现方式，从而也导致了各个模型之间存在显著的差异。

根据发射光线起点的不同，光线追踪算法可分为前向光线追踪和后向光线追踪。前向光线追踪从光源发射光线，并跟踪光线在场景中的传播过程，从而模拟所需的遥感信号，典型的模型有Raytran、FLIGHT^［13］、FLiES^［14］、Rayspread^［15］以及FluorWPS^［16］等。Raytran模型从光源发射单色光，根据场景元素的光学属性决定光线是否被吸收或者被散射，最终通过统计逃逸出场景的光线数量以及光线所携带的能量来估算反射率。虽然Raytran模型这种实现方式简单，但是计算效率较低，且难以同时模拟多波段的反射信号。为了解决Raytran模型的效率问题，Widlowski等^［15］在Raytran的基础上发展了Rayspread模型，通过在每个光线碰撞点处向观测方向发射“子光线”（secondary rays）来估算反射率，提高了计算效率和精度。利用同样的原理，Zhao等^［17］提出了进一步考虑权重变化、可进行高光谱和荧光信号辐射模拟的Weighted Photon Spread （WPS）模型和FluorWPS模型。后向光线追踪从传感器发射光线，并计算辐射源沿该光线传播到传感器的辐亮度值。后向光线追踪的主要优势在于只计算了恰好进入传感器部分的能量，因此在模拟影像方面具有较大的优势，但是难以计算能量平衡问题。典型的后向光线追踪模型是DIRSIG模型^［18］，由于后向光线追踪的特点，该模型尤其适合于研究传感器参数对遥感信号的影响，是辅助传感器设计的有力工具，但是该模型的主要问题是难以共享使用。因此，能够同时进行前向光线追踪和后向光线追踪以应对不同的模拟需求则显得尤其重要，同时也是现有模型所欠缺的功能。

现有三维辐射传输计算机模拟模型所面临的另一个问题是计算效率低以及模拟场景小。为了提高计算效率，常用的做法是简化场景结构，例如FLIGHT模型将树冠表示为简单几何体（椭球体等）以简化光线与场景的求交运算；RAPID模型将树冠表示为孔隙面元，从而减少场景中的面元数量，以扩大模拟范围以及提高计算效率。因此目前遥感领域急需发展一个既能表达真实场景结构又能高效率模拟大场景的三维辐射传输模型。

近年来，光线追踪技术在计算机图形学领域内得到了长足的发展，涌现了一批高效率的物理渲染引擎，这些引擎能够从虚拟的三维场景中渲染逼真的图像，例如LuxCoreRender（https：//luxcorerender.org/）、PBRT^［19］以及POV-Ray^［20］等。这些渲染引擎充分地使用了最新的加速计算技术（例如GPU），极大地提高了场景的渲染效率。然而直接将这些渲染引擎应用于遥感信号的模拟并不是一件容易的工作，其主要原因在于这类渲染引擎主要关注图像的逼真程度，而不是物理量的准确性。部分模型只能模拟RGB 3个宽波段，而无法模拟遥感场景的多光谱或者高光谱影像。因此，在这些引擎基础上进行修改或者利用这些引擎提供的加速结构重新实现遥感信号的模拟是一个更为可行的方案。例如，Auer等^［21］基于POV-Ray开发了三维场景的SAR信号模拟器；DIRSIG模型最近对其核心部分进行了重构，以利用最新的图形学技术加速其计算过程^［22］；DART模型在其核心部分集成了Intel公司开发的Embree光线追踪库，其计算效率得到了较大的提升。

基于以上考虑，利用Mitsuba^［23］渲染软件所提供的加速框架发展了基于光线追踪算法的三维真实冠层辐射传输模型LESS^［24］（LargE-Scale Remote Sensing Data and Image Simulation Framework），该模型能够同时进行前向光线追踪以及后向光线追踪，同时充分利用了最新图形学技术以实现更准确和更高效的大场景遥感信号的模拟。本文主要介绍了LESS的基本原理、主要功能以及典型的应用案例，为利用LESS模型进行定量遥感建模研究提供参考。

2 LESS基本原理

在LESS模型中，根据光线出发的位置，分别实现了前向光子追踪算法和后向光线追踪算法。前向光子追踪从光源发射光子，并模拟光子在整个场景中的传播过程，通过对光子能量的记录，可以模拟场景的反射率、FPAR（Fraction of Absorbed Photosynthetically Active Radiation）以及上下行辐射等遥感数据。后向光线追踪从传感器发射光线，因此适合模拟遥感多/高光谱影像以及各种投影模式的相机。

2.1　前向光子追踪算法

在LESS前向光子追踪模式中，光子从光源出发进入场景并与场景中的元素进行求交运算。典型的光源有太阳直射光和天空散射光，一般太阳直射光作为平行光处理，而天空散射光作为各向同性入射光处理。对于直射光而言，LESS将其表示为一个投影后能够覆盖整个场景的圆盘，该圆盘的大小和位置由场景的外接球确定，这种方式能够保证在任意太阳入射角下，入射辐射能够覆盖整个场景。因此，整个直射光源的发射功率为 $L_{d i r} = π r^{2} E_{d i r}$ ，其中 $r$ 为圆盘半径， $E_{d i r}$ （w/m²）为沿太阳入射方向的辐照度。对于散射光而言，光子的起始位置位于包围整个场景的一个半球之上，其半径通常为外接球的2~3倍，因此天空散射光的发射功率为 $L_{d i f f} = 2 π^{2} R^{2} I_{d i f f}$ ，其中 $R$ 为天空半球的半径， $I_{d i f f}$ 为天空散射光的辐亮度。

当产生入射光子时，每个光子的入射能量（瓦特）可表示为：

P^{0} (λ) = \frac{L_{d i r} + L_{d i f f}}{N}

(1)

其中： $N$ 为一次模拟中产生的总光子数量，在这 $N$ 个光子中，直射光线和散射光线的数量按照各自的发射功率进行分配，即 $w_{k} = \frac{L_{k} (λ)}{\sum_{k = 1}^{K} L_{k} (λ)}$ ，式中 $K$ 为光源数量。需要注意的是，此处的光子为能量包的概念，即该光子包含了不同波长的能量，且该能量在传输过程中可能根据场景元素表面散射特性进行改变，从而达到在同一个模拟中同时模拟多波段遥感信号的目的。

2.1.1　真实光子跟踪算法

当光子进入场景后，可能与场景中的元素发生碰撞，在每一个碰撞点，光子的能量将变为：

P^{Q} (λ) = P^{0} (λ) ∙ \prod_{q = 1}^{Q} [π f (q, ω_{i,} ω_{o,} λ) / p_{q}]

(2)

其中： $Q$ 为光子传输过程中发生碰撞的次数， $f (q, ω_{i,} ω_{o,} λ)$ 表示在碰撞点 $q$ 处的BSDF（Bidirectional Scattering Distribution Function）函数； $p_{q}$ 为光子发生反射或者透射的概率，可以根据场景元素（例如叶片）的反射率和透过率进行设定，但是在多波段情况下每个波段的反射率和透射率有较大的差异，因此通常情况下 $p_{q}$ 可以简单地设置为0.5，由于计算光子能量时已经考虑了此概率带来的影响，在光子数量较大时，同样能够得到较为准确的结果。与传统模型（例如DART）将光线拆为多个子光线做法不同的是，LESS模型中这种基于概率的方式能够保持光子数量不变，提高了计算效率。

当光子从场景顶部逃离场景后，能够被放置于场景顶部的虚拟半球所搜集，该半球通常被划分为多个等面积的面元，利用每个面元的立体角可以估算该方向的场景反射率。

f_{B R F_{i}} = \frac{π P_{i}^{A}}{∆ Ω_{i} ∙ c o s θ_{i}^{c} ∙ P_{s c e n e}}

(3)

其中： $P_{i}^{A}$ 是面元 $i$ 搜集到的总的光子能量，即 $P_{i}^{A} = \sum_{P^{Q} \in ∆ Ω_{i}} P^{Q}$ ， $∆ Ω_{i}$ 为每一个面元对应的立体角； $θ_{i}^{c}$ 是面元 $i$ 的中心天顶角； $P_{s c e n e}$ 是场景顶部的参考平面（与场景大小一致）所搜集到的能量。如果光子从场景边缘逃离场景，往往将其重新移动到场景的对立面以相同的方向再进入场景，从而模拟水平无限延伸的场景模式。当光子离开场景或者散射次数达到一个预设值（例如5次），光子将被终止。为了避免能量损失，每个光子按照一定的概率 $p$ 被终止传播（通常 $p = 0.05$ ），如果光子没有成功地被终止，其能量将会乘于概率值 $\frac{1}{1 - p}$ ，以保持能量平衡。

2.1.2　虚拟光子算法

在真实光子算法中，反射率的估算要借助于一个微小立体角，但是该立体角内的光子方向仍然有一定的差异，从而导致在热点附近产生较大的低估。此外，当立体角内搜集的光子数量不足时，估算的反射率会有较大的误差。为解决该问题，LESS模型提出了虚拟光子算法，即当光子与场景元素发生碰撞且未被完全吸收的情况下，将会沿观测方向（虚拟方向）发射一个虚拟光子，用于计算入射光经过碰撞点可能沿观测方向所散射的能量（ $W ∙ s r^{- 1}$ ），即：

I = V ∙ P^{q - 1} ∙ f (q, ω_{i}, ω_{v}) ∙ c o s < ω_{v}, ω_{n} >

(4)

其中： $P^{q - 1}$ 表示在第 $q$ 个碰撞点处光子入射的能量； $ω_{v}$ 是需要观测方向； $V$ 是可视因子，如果虚拟光子沿虚拟方向传播路径被场景元素遮挡则 $V = 0$ ，否则 $V = 1$ 。最终反射率的计算公式为：

f_{B R F_{v}} = \frac{π I_{v}^{A}}{c o s θ_{v} ∙ P_{s c e n e}}

(5)

式中： $I_{v}^{A}$ 为虚拟方向 $v$ 的辐射强度（ $W ∙ s r^{- 1}$ ）； $θ_{v}$ 是虚拟方向的天顶角。由于虚拟光子算法得到的能量是基于无限小的立体角，因此在热点方向能够得到较为准确的反射率值。同时，在光子传播路径的每一个碰撞点都发射虚拟光子，而真实光子算法只在最终光子离开场景后才有能量贡献，因此虚拟光子算法相较于真实光子算法的效率有显著的提升，能利用较少的入射光子数量得到较为准确的结果^［25］。

2.2　后向光线追踪算法

2.2.1　一次散射

后向光线追踪从传感器发射光线，经过与场景元素求交后找到与光源之间的一条路径，并计算光源沿该路径散射到传感器的辐亮度，因此后向光线追踪也称为后向路径跟踪。如图1所示，当光线找到与场景元素的交点 $q$ 时，需要计算太阳或者大气经过 $q$ 点散射到传感器的辐亮度。由于场景中可能同时存在太阳直射光和天空散射光，因此将根据一定的概率在这两个光源中随机选择一个，通常它们的选择概率可以设置为0.5，即 $p_{d i r} = p_{d i f f} = 0.5$ ，根据选择的光源的不同，计算方式也有一定的差异：

图1

新窗口打开| 下载原图ZIP| 生成PPT

图1 后向光线追踪表面散射计算过程

Fig.1 Surface scattering calculation in backward path tracing

（1）如果选择了太阳，那么将从 $q$ 点出发沿太阳入射的反方向在太阳圆盘上找到一点 $q^{'}$ 。如果 $q$ 与 $q^{'}$ 之间无遮挡，则太阳的辐照度为 $E_{d i r} / p_{d i r e c t i o n}^{s u n}$ ，其中 $p_{d i r e c t i o n}^{s u n}$ 表示太阳产生该方向的光线的概率，由于太阳入射方向固定，因此 $p_{d i r e c t i o n}^{s u n} = p_{d i r} ∙ 1$ 。

（2）如果选择了大气，将在 $q$ 点上半球空间随机产生一个方向，并计算 $q$ 点沿该方向与大气半球的交点 $q^{'}$ ，若 $q$ 与 $q^{'}$ 之间无遮挡，则大气的辐亮度为 $I_{d i f f} / p_{d i r e c t i o n}^{a t m}$ ，其中 $p_{d i r e c t i o n}^{a t m}$ 表示在 $q$ 处产生该方向光线的概率，其计算公式为 $p_{d i r e c t i o n}^{a t m} = p_{d i f f} \frac{c o s θ_{a t m}}{π}$ ， $θ_{a t m}$ 为光线与 $q$ 点平面法线之间的夹角。

上述过程相当于对光源进行了采样，要计算入射光线经过 $q$ 散射到传感器的能量，还需要考虑对 $q$ 点处的BSDF函数进行采样。由于此处物体表面为朗伯面，因此在 $q$ 产生特定入射角度光线的概率为 $p_{b s d f}^{a t m} = \frac{c o s θ_{a t m}}{π}$ （或者 $p_{b s d f}^{s u n} = \frac{c o s θ_{s u n}}{π}$ ），该过程同时对光源和BSDF函数进行了采样，这种采样方式称之为多重重要性采样（Multiple importance sampling）^［26］，它们的组合概率为：

p_{w e i g h t}^{s u n} = \frac{{p_{d i r e c t i o n}^{s u n}}^{2}}{{p_{d i r e c t i o n}^{s u n}}^{2} + {p_{b s d f}^{s u n}}^{2}}, p_{w e i g h t}^{a t m} = \frac{{p_{d i r e c t i o n}^{a t s}}^{2}}{{p_{d i r e c t i o n}^{a t m}}^{2} + {p_{b s d f}^{a t m}}^{2}}

(6)

最终太阳对于传感器辐亮度贡献的计算公式为：

L_{s u n} = \frac{E_{d i r}}{p_{d i r e c t i o n}^{s u n}} ∙ p_{w e i g h t}^{s u n} ∙ f (q, ω_{i}, ω_{o}) ∙ c o s θ_{s u n} = \frac{E_{d i r} c o s θ_{s u n}}{p_{d i r}} ∙ p_{w e i g h t}^{s u n} ∙ f (q, ω_{i}, ω_{o})

(7)

大气对传感器辐亮度贡献的计算公式为：

L_{a t m} = \frac{I_{d i f f}}{p_{d i r e c t i o n}^{a t m}} ∙ p_{w e i g h t}^{a t m} ∙ f (q, ω_{i}, ω_{o}) ∙ c o s θ_{a t m}

(8)

= \frac{π I_{d i f f}}{p_{d i f f}} ∙ p_{w e i g h t}^{a t m} ∙ f (q, ω_{i}, ω_{o})

在此计算过程中所涉及到的多重重要性采样技术能够同时考虑到光源和BSDF的特性，尤其适合于非朗伯表面散射的计算。

2.2.2　多次散射

对于多次散射的计算，将从第一个碰撞点 $q$ 根据BSDF随机产生一条新的光线，并与场景进行求交测试，如果找到交点 $q_{1}$ （如图2），将在 $q_{1}$ 重复2.2.1节中所描述的一次散射计算过程，从 $q_{1}$ 散射的辐亮度沿着 $ω_{1}$ 方向到达 $q$ ，最终经过 $q$ 返回到传感器，这一部分能量就是二次散射项。该过程一直重复直到光线离开场景或者达到用户设定的最大散射次数。同样地，为了避免能量损失，光线仍以一定的概率 $p$ 被终止传播，对于未被终止的光线，其能量将乘于概率 $\frac{1}{1 - p}$ 。对于多波段的模拟，所有波段均采用相同的路径以节省计算时间，从而实现多光谱以及高光谱影像的模拟。

图2

新窗口打开| 下载原图ZIP| 生成PPT

图2 后向光线追踪多次散射计算过程

Fig.2 Multiple scattering calculation in backward path tracing

3 LESS功能模块

造成三维模型难以使用的原因之一是三维模型输入参数多且复杂，因此开发一个操作方便的软件是使用三维模型的重要前提。为此，基于JavaFX开发了可视化界面，并且提供了一套构建三维场景、管理输入参数以及处理输出参数的工具集。图3（a）展示了LESS模型的主界面，主要包含场景显示区域、运行记录输出区域以及参数设置区域，其中参数设置区域又包含了传感器、观测几何、入射光照、光学属性、地形以及场景元素等若干子区域。

图3

新窗口打开| 下载原图ZIP| 生成PPT

图3 LESS模型主界面及典型功能

Fig.3 The main screen and major functionalities of LESS

3.1　三维场景构建

一个典型的三维场景包含了地形和一些场景元素（例如树和房子）。LESS模型支持3种类型地形文件格式，分别为平面（PLANE）、栅格影像（RASTER）和三角网（MESH）。其中栅格影像支持ENVI标准格式和GeoTIFF格式，而三角网则支持Wavefront OBJ格式。

在LESS中，场景元素均以三角网格式进行描述（Wavefront OBJ），这些元素的结构一般由外部软件产生，例如通过Onyx Tree（http：//www.onyxtree.com/）可以产生较为逼真的单木结构。当这些元素被导入到LESS模型中时，需要对其中包含的每一个组分（例如枝干和叶片）设定对应的光谱属性，但是此时的场景元素并没有立即出现在场景中，而是需要对其设定相应的地理位置。为了提高模拟效率，LESS模型采用几何变换的形式将每一个物体放置在场景中，因此一个物体可以被“复制”多次，但在此过程中只保存了一系列的变换矩阵，避免了将整个元素的所有三角网进行完全复制。通过这种方式可以以较小的计算资源模拟更大的场景（例如单木组成的森林），节省计算内存。

在三维场景构建中，LESS也提供了一些便利的工具来辅助场景的构建，例如多边形编辑工具可以在任意绘制的多边形内放置随机分布的物体，或者删除任意绘制的多边形内的物体（图3（b））。

3.2　三维场景展示

由于三维模型的场景往往比较复杂且计算时间较长，因此在计算之前有必要检查所创建的场景是否正确或者符合预定的目标。LESS模型基于JavaFX 3D开发了三维场景展示工具，该工具能够将栅格地形文件自动转换为三角网，并自动计算场景元素对应的高程信息，从而展示正确的三维场景。此外，三维场景展示工具也能够显示当前所设定的传感器的视场范围，可以辅助操作人员检查所设置传感器的参数是否正确。图3（c）展示了典型的三维场景以及对应的传感器视场范围。

3.3　遥感数据模拟

根据所选择的传感器的不同，LESS模型可以模拟并输出不同的遥感数据，例如多光谱影像、鱼眼相机、多角度反射率、四分量影像等。以FPAR模拟为例，LESS模型利用前向光子跟踪算法在每一个光子与场景的碰撞点计算场景元素所吸收的能量 $A (λ) = σ_{λ} ∙ P^{Q - 1} (λ)$ ，其中 $σ_{λ}$ 为元素吸收率，然后对可见光波段进行积分即可得到 $P A R = \int_{400}^{700} A (λ) d λ$ ，该值与冠层顶部总入射的比值即为FPAR。此外，对光子进行分层统计，可以输出分层分组分的吸收能量，这对应农作物光截获的研究具有十分重要的作用。

3.4　脚本模块

LESS的核心部分由C++编写而成，且提供相应的Python接口，因此可以通过Python调用LESS所提供的接口以实现其他模拟功能或进行二次开发。在LESS中，已经完整地集成了一套Python解释器，并且已经安装了常用的库。当模拟工程以及三维场景建立之后，可以通过主界面的Python命令行窗口运行Python脚本程序。图4展示了一段Python脚本，该脚本通过朝不同方向发射光线计算交点，可以模拟场景的点云数据。

图4

新窗口打开| 下载原图ZIP| 生成PPT

图4 LESS中运行Python脚本

Fig.4 Run Python script in LESS

4 LESS应用案例

本节主要介绍了LESS各种模拟功能的应用案例，包括遥感影像模拟、下行辐射模拟、FPAR模拟以及分析木质元素对冠层反射率的影响，并且对部分模拟结果（4.1~4.3节）与实测数据进行了对比验证。

4.1　遥感图像模拟

从激光雷达获取三维结构信息并结合光谱信息模拟遥感信号是目前三维辐射传输模型的应用方向之一，因为激光雷达提供了前所未有的精细结构信息，这些信息正是三维辐射传输模型所必需的。通常而言，从激光雷达数据重建三维虚拟森林场景有两种方式，一是基于体元的方法^［27］，二是基于单木的方法^［28］。本节主要介绍基于单木重建场景并进行遥感影像模拟的方法。

首先利用Onyx Tree生成了6棵典型的单木，其中3棵白桦，3棵落叶松（如图5所示）。单木的位置和树冠由机载激光雷达数据得到，对于每一个检测到的单木，为了从这6棵单木中找到最合适的单木，定义了一个相似因子，即 $S F = \frac{0.5 \cdot |H_{d} - H_{o}|}{H_{o}} + \frac{0.5 \cdot |C_{d} - C_{o}|}{C_{o}}$ ，其中 $H_{d}$ （ $C_{d}$ ）为激光雷达数据中检测到的树高（冠幅）， $H_{o}$ （ $C_{o}$ ）为生成的单木的树高（冠幅）。在这6棵单木中找到 $S F$ 值最小的单木。

图5

新窗口打开| 下载原图ZIP| 生成PPT

图5 Onyx Tree生成的单木

Fig.5 Individual Trees generated by Onyx Tree

整个场景范围为1 km×1 km，约有50 000棵单木，模拟影像的分辨率为1 m，波段数量为64，模拟时间约为30 min（16 GB内存，8核心笔记本）。实测数据为图5展示了模拟影像与实测影像之间的对比（假彩色），其中实测影像为机载高光谱数据，高光谱传感器为AISA Eagle II，该仪器能够记录从400~970 nm之间的64个波段，空间分辨率约为1 m。总体来看，模拟影像与实测影像在纹理上比较相似，尤其是在树木比较稀疏的区域，因为在此区域单木分割的精度更高。在样方尺度，树高树冠幅的分布也与真实影像比较接近，因为在构建场景时主要根据单木的树高的冠幅进行单木的选择。但是，两幅影像最大的区别在于颜色（光谱），模拟影像的颜色更为单一，而实际影像表现较强的空间异质性，这主要是因为模拟场景中所有的树木都采用了相同的组分光谱。尽管如此，模拟数据可以为提供场景真实的结构参数（比如LAI），以及对应的反射率，可以为验证其他物理模型提供数据基础。此外，基于单木的方法具有每棵单木的细节（到叶片层次），在模拟一些高分辨率影像（例如无人机影像）时具有较大的优势。

图6

新窗口打开| 下载原图ZIP| 生成PPT

图6 模拟影像与实测影像的对比

Fig.6 Comparison of the simulated image and actual image

4.2　下行辐射的模拟

地形起伏对于太阳短波下行辐射的估算具有较大的影响，具体表现在起伏地形对太阳和天空光的遮挡以及周围地形的相互作用。通过模拟复杂地形条件下的短波下行辐射，可以定量评估不同地形对下行辐射的影响，也可以为实际测量中的布站优化提供参考依据，以最少的测量站点获得更具有区域代表性的观测值。

LESS模型可以模拟具有真实结构复杂地形条件每个像元的短波下行辐射值。以承德月亮山试验区（42.05778N， 117.0656E）为例，利用该地区真实的地形数据，开展短波下行辐射的模拟分析并与实际测量值进行了比较。该试验区一共有4个观测站点，分别位于山顶（Station1）、东坡（Station2）、南坡（Station3）以及西坡（Station4）。由于山顶站点受地形影响最小，因此在该模拟中将山顶站点测量的直散射数据作为大气底层入射辐射输入到LESS模型中，并模拟其余3个站点的下行辐射值。在模拟过程中，用到的DEM数据为从资源三号（ZY-3）立体像对得到的12.5 m分辨率的地形数据。

图7（b）展示了2015年8月15日（晴空）的模拟结果以及与实测数据的对比，模拟的时间间隔为9 min。从对比结果可以看出，模拟数据与实测数据在各个时刻都具有较高的一致性。东坡、南坡以及西坡的相关系数都为0.99，RMSE分别为33.7 W·m^-2·μm^-1、14.6 W·m^-2·μm^-1和22.9 W·m^-2·μm^-1。对于东坡，模拟值在12：00之后比实测值偏高，而南坡和西坡则在中午时出现较大的偏差，这些偏差主要是由于实测数据以及地形数据的不确定性引起的，例如辐射表安装时可能与实际坡面不平行，DEM分辨率较粗而与实际地形有一定的偏差。从3个站点之间的对比结果来看，东坡首先到达最大值（10：30左右），其次是南坡（12：00左右），最后是西坡（14：00左右），这种最大值的变化规律也符合太阳一天的方位角变化规律。此外，由于地形遮蔽的影响，东坡在18：00附近出现辐射值的突变，而西坡和南坡分别在6：30和7：30附近出现突变。以上分析说明，LESS模拟的辐射值与实测值具有较好的一致性，正确地反映地形对于下行辐射的影响。

图7

新窗口打开| 下载原图ZIP| 生成PPT

图7 复杂地形区LESS短波下行辐射模拟（Obsv表示观测值，LESS表示模拟值）

Fig.7 Solar shorwatve radiation simulation over rugged terrains （Obsv stands for actual measurement and LESS stands for simulated data）

4.3　FPAR模拟

LESS模型可以模拟冠层分层FPAR，用于验证其他物理模型或者发展参数化反演模型。此处选择了实测的冬小麦和夏玉米FPAR数据集^［29］对LESS的FPAR模拟结果进行了对比和分析。该数据集获取于北京市昌平区小汤山国家农业精准实验基地，数据获取时间为2008年4月27日（冬小麦）和2010年8月29日（夏玉米），其中夏玉米包含两个播种期的植株：处于喇叭口期的夏玉米-I以及处于拔节期的夏玉米-II。冬小麦和夏玉米的参数（株高、垄距、光谱以及LAI等）通过实测获得。冬小麦选择了3个样点，分别沿顺垄方向（PZ8-S）和垂直垄方向（PZ8-C）各测量了一次，夏玉米选择了1个样点，同样地，在该样点处沿顺垄方向（JD28-I-S和JD28-I-S）和垂直垄方向（JD28-II-S和JD28-II-S）分别测量了1次。在每个测量样点，利用SUNSCAN冠层分析仪从冠层顶部到底部每隔16 cm（冬小麦）/20 cm（夏玉米）向上及向下各测了1次上行和下行PAR，从而得到不同高度处的FPAR。

LESS分层FPAR的模拟结果以及实测FPAR如图8所示。从该图可以看出，冬小麦实测的3个样点的FPAR值具有相似的规律，即随着观测高度的增加而降低。但是在同一高度处，不同样点之间的测量值具有一定的差异，尤其是冠层顶部差异较大，导致该差异的主要原因是不同样点之间的作物生长状态可能存在较大的不同。由于类似的原因，夏玉米在不同高度处的测量值差异也较大。不过从冬小麦和夏玉米的模拟结果来看，模拟的FPAR值与实测值比较接近，处于多次测量值的范围内，说明利用LESS进行FPAR的模拟具有较高的精度。通过模拟的方式来研究FPAR与植被结构之间的关系，可以有效地避免地面测量过程中的人为误差，对定量地分析FPAR的影响因素具有重要的意义。

图8

新窗口打开| 下载原图ZIP| 生成PPT

图8 LESS模拟FPAR模拟

Fig.8 FPAR simulations in LESS

4.4　分析木质元素对反射率的影响

辐射传输模型是利用遥感数据定量反演参数的基础。目前多数植被辐射传输模型在反演参数时，将植物结构进行简化，忽略了木质元素的影响。然而真实场景中，木质元素作为植物重要组成部分，对植物结构以及植物的理化特性都有重要影响。作为验证其他模型的验证工具之一，LESS模型不仅可以实现场景内更接近现实世界的植物结构模拟，还可以方便地改变植被结构假设，研究植被不同结构对方向反射率的影响。

图9展示了1个虚拟针叶林样方在不同叶面积情况下有枝干和无枝干的反射率差异，图中单木由Onyx Tree软件生成，通过调节叶片密度，从而实现相同枝干不同叶面积单木的生成。图中C表示场景只有土壤和枝干，F表示场景既有土壤和枝干，也有叶片。由于植被场景在红光和绿光波段的反射率值比较低，因此，枝干反射率在这两个波段引起的差异也相对较小。最大的相对差值出现在LAI为0，即没有树叶的场景下，所有差异全部来源于枝干。随着叶片的增加，LAI增大，场景在红光和绿光波段的反射率值降低，其绝对差值也逐渐减小，并在LAI大于3以后趋于稳定。与之相对应地，场景在近红外波段反射率随着LAI的增加而增大，且两类场景在该波段的反射率绝对差值也逐渐增加。

图9

新窗口打开| 下载原图ZIP| 生成PPT

图9 枝干对场景反射率的影响

Fig.9 Reflectance influenced by branches

5 结语

本文主要介绍了基于光线追踪的三维辐射传输模型LESS的基本原理以及典型的应用案例。与目前已有的模型（例如DART）相比，LESS模型的主要优势在于能够以较高的效率模拟具有真实结构的复杂大场景的各种遥感信号，并且针对不同模拟数据采用了不同的模拟策略（前向和后向光线追踪）。此外，LESS模型参数设置更为简单，与算法本身有关的参数更少，更利于普通用户的理解和使用。为简化三维模型的使用难度，LESS模型提供了简单易用的GUI程序，该程序包含了参数输入、三维场景构建、三维场景显示、模拟计算以及后处理等一系列工具。未来，LESS模型将会继续完善大气以及云辐射、热红外、激光雷达以及荧光等模拟。

参考文献

原文顺序

文献年度倒序

文中引用次数倒序

被引期刊影响因子

[1]

Verhoef

Light Scattering By Leaf Layers With Application To Canopy Reflectance Modeling： the SAIL model

［J］. Remote Sensing of Environment， 1984， 16（2）： 125-141.