CN112254820A - 一种离散森林场景热红外辐射传输模拟方法 - Google Patents

一种离散森林场景热红外辐射传输模拟方法 Download PDF

Info

Publication number
CN112254820A
CN112254820A CN202011096826.4A CN202011096826A CN112254820A CN 112254820 A CN112254820 A CN 112254820A CN 202011096826 A CN202011096826 A CN 202011096826A CN 112254820 A CN112254820 A CN 112254820A
Authority
CN
China
Prior art keywords
crown
formula
canopy
visible
trunk
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202011096826.4A
Other languages
English (en)
Other versions
CN112254820B (zh
Inventor
卞尊健
杜永明
历华
曹彪
肖青
柳钦火
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Aerospace Information Research Institute of CAS
Original Assignee
Aerospace Information Research Institute of CAS
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Aerospace Information Research Institute of CAS filed Critical Aerospace Information Research Institute of CAS
Priority to CN202011096826.4A priority Critical patent/CN112254820B/zh
Publication of CN112254820A publication Critical patent/CN112254820A/zh
Application granted granted Critical
Publication of CN112254820B publication Critical patent/CN112254820B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
    • G01J5/00—Radiation pyrometry, e.g. infrared or optical thermometry
    • G01J5/007—Radiation pyrometry, e.g. infrared or optical thermometry for earth observation
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
    • G01J5/00—Radiation pyrometry, e.g. infrared or optical thermometry
    • G01J5/0066—Radiation pyrometry, e.g. infrared or optical thermometry for hot spots detection
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
    • G01J5/00—Radiation pyrometry, e.g. infrared or optical thermometry
    • G01J2005/0077—Imaging

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Environmental & Geological Engineering (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geology (AREA)
  • Photometry And Measurement Of Optical Pulse Characteristics (AREA)

Abstract

本发明公开了一种离散森林场景热红外辐射传输模拟方法，包括以下步骤：S1、考虑森林的三维结构，将整个森林场景按照类别和温度差异划分为四层，从上向下依次包括上层树冠、中层树干、林下植被和底层土壤；S2、对每层的组分进行离散，按照几何结构划分为不同的格网或者子层；S3、直接辐射模拟：按照格网内叶面积密度计算每层组分的光照和阴影可视比例，通过普朗克函数计算组分的直接发射项；S4、多次散射模拟：按照光谱不变理论计算组分间的多次散射项；S5、冠层顶发射项等于直接发射项和多次散射项之和。本发明通过建立冠层顶热辐射与冠层内各组分温度状态的联系，不仅适用性广泛且模拟精度高。

Description

一种离散森林场景热红外辐射传输模拟方法
技术领域
本发明涉及一种新的热红外波段辐射传输模拟方法，尤其涉及一种离散森林场景热红外辐射传输模拟方法。
背景技术
地表温度(LST)是森林植被生长评价、干旱和火情监测等方面的重要指标，是地球表面能量平衡和水循环过程的重要输入参量。利用遥感的手段，可以获取大范围的地表温度信息。从卫星热红外观测数据反演得到的地表温度通常是由像元内部不同组分共同作用的结果。因此像素平均温度并不能直接反应森林内部的温度分布，而需要将不同组分的作用考虑在内。由于组分的物理性质以及三维结构的异质性，地表温度反演会受到观测以及光照几何角度的显著影响。
森林的结构十分复杂，呈现出明显的层状分布，例如顶层的树冠、中间层的树干以及底层的灌木、草地以及土壤。森林场景热辐射的方向性可导致明显的温度差异，而温度的差异体现在组分间以及组分内部的温度差异。然而，在许多前向热红外模型中，仅假设了四个分量(如光照冠层、阴影冠层、光照土壤以及阴影土壤)，并未将树干以及下层植被的影响考虑在内。特别是对于落叶林等森林景观，在特定的时期，树干等成为了传感器视场内的主要组成成分，这对于亮温(亮度温度)的方向异质性有着显著的影响。除三维结构的垂直异质性外，冠层的垂直温度分布也是引起亮温角度观测差异性的重要因素。树木的叶片温度垂直变化明显，最高可达15℃，不同层次的叶片的可观测比例及其对于冠层顶亮温的贡献程度随着视角的变化而变化。在现有的许多辐射传输模型中的光照和阴影叶片平均温度的假设并不能反应这种现象。
目前，在可见光和近红外波段的辐射传输模型中，有部分模型考虑到了森林冠层结构的异质性特点，如GOSAIL模型、4SCALE模型、FRT模型等。但是在反演地表温度的热红外波段，并没有考虑林冠结构和温度的垂直异质性的相关模型。与此同时，将相关的近红外模型扩展到热红外波段的解决思路已经被许多研究所采用，如SAIL模型、TRGM模型等。
发明内容
为了解决上述技术所存在的不足之处，本发明提供了一种离散森林场景热红外辐射传输模拟方法。
为了解决以上技术问题，本发明采用的技术方案是：一种离散森林场景热红外辐射传输模拟方法，包括以下步骤：
S1、考虑森林的三维结构，将整个森林场景按照类别和温度差异划分为四层，从上向下依次包括上层树冠、中层树干、林下植被和底层土壤；通过直接辐射模拟、多次散射模拟对各组分可见比例进行计算，最后计算出热红外波段冠层顶表面热辐射量；
S2、对每层的组分进行离散，按照几何结构划分为不同的格网或者子层；
S3、直接辐射模拟：按照格网内叶面积密度计算每层组分的光照和阴影可视比例，通过普朗克函数计算组分的直接发射项；
S4、多次散射模拟：按照光谱不变理论计算组分间的多次散射项；
S5、冠层顶发射项等于直接发射项和多次散射项之和。
进一步地，步骤S1中，在热红外波段，表面辐射分为组分辐射和大气辐射两部分，特定波长的冠层顶表面热辐射量(L)由下式计算：
Figure BDA0002724032160000021
其中：θ1，θ2分别代表太阳和观测的天顶角；
Figure BDA0002724032160000022
表示观测角度和光线角度的相对方位角；Tj和εj代表j组分的温度和发射率；f代表组分的可见比例；εm代表多次散射效应，指从一个组分发射并由其他组分反射到传感器的辐射；εe代表冠层的有效发射率；
Figure 100002_DEST_PATH_IMAGE002
代表大气有效下行辐射；B代表普朗克函数。其中
Figure BDA0002724032160000034
代表j组分的直接辐射贡献；εm，j为j组分的多次散射项贡献，计算出每个组分的直接辐射和多次散射贡献，然后将其累加/积分则为冠层顶的热辐射。
进一步地，直接辐射模拟中，上层树冠可见比例的具体计算如下：
基于双向孔隙率计算出一个体元中光照和阴影树叶的可见比例，然后使用所有相应的体积元素的总和来计算光照和阴影树冠的可见比例，计算公式如下所示：
fc，s(r2)＝∑fc，k，s(r1，r2) 公式二
fc，h(r2)＝∑fc，k，h(r1，r2) 公式三
Figure BDA0002724032160000031
Figure BDA0002724032160000032
fc，k，h(r1，r2)＝fc，k(r2)-fc，k，s(r1，r2) 公式六
其中，fc，k,fc，k，h和fc，k，h分别代表第k层树冠的可见比例、光照树冠的可见比例和阴影树冠的可见比例；fc，s和fc，h为整个树冠的可见光照比例和可见阴影比例；(x，y，z)代表树冠内点M的坐标；Vk代表第k层树冠包络的空间区域； r1和r2分别代表光照和观测方向；λ代表树冠的密度；p00代表由点(x，y，z)的r1和r2计算出的双向孔隙率；uL代表叶片的面体积密度；G代表投影在观测或光照方向的叶片的比例，对于球形叶倾角分布函数(LIDF)，固定为0.5；G(θ2) 代表树叶在θ2方向的投影；
在热红外波段，假设所有组分均为朗伯体，叶片的贡献直接反映在其截面上，因此，双向孔隙率由两个独立分布概率的乘积表示，计算公式如下所示：
p00(x，y，z，r1，r2)＝p1(x，y，z，r1，r2)·p2(z，r1，r2)公式七
其中，p1表示冠层内双向孔隙率；p2表示冠层间双向孔隙率；
假设在树冠中均为均质的叶片，p1主要取决于点在树冠中的位置和叶面体积密度，p2主要取决于树密度和叶面体积密度，计算p1和p2的公式如下所示：
p1(x，y，z，r1，r2)＝ai(s1，θ1)ai(s2，θ2)CHS1(s1，s2，θ1，θ2) 公式八
ai(s，θ)＝exp[-suLG(θ)] 公式九
Figure BDA0002724032160000041
Figure BDA0002724032160000042
其中，ai(s，θ)代表方向为θ的点到冠层表面的间隙距离；CHS1为热点校正因子，描述了冠层内部叶片在光照和观测投影内的重叠效应；θ1和θ2为太阳天顶角和观测天顶角；sAB为路径长度的重叠；s1和s2分别代表从点(x，y，z)沿光照和观测方向到冠层顶面之间的路径长度；α表示r1和r2之间的立体角；b为热点因子，特定的球形叶倾角分布函数为16/(dLπ2)，水平叶倾角分布函数为 4/(dLπ)；p2用于表示树冠之间的多次阴影效应，由树冠中向上的点 M1(x1，y1，z1)和M2(x2，y2，z2)计算：
p2(z，r1，r2)＝as(z1，θ1)as(z2，θ2)CHS2(z1，z2，l12，r1，r2) 公式十二
CHS2(z1，z2，l12，r1，r2)＝exp[λcSc(z1，z2，l12，r1，r1)p0] 公式十三
其中，as(z，θ)是在方向θ、高度为z的冠层间隙的平均比例；CHS2也为热点校正因子，描述光照和观测方向上阴影冠层内部特征；参数c表征树的分布模式与泊松分布的偏差；Sc(z1，z2，l12，r1，r1)为光照和观测方向上树冠投影区域的公共部分面积，对应高度为z1，z2以及水平距离l12；p0为树冠内的联合孔隙率；孔隙率as(z，θr)在假设树为二项分布的情况下，计算如下：
as(z，θr)＝exp[-λ[b(z，θr)Scrown(z，θr)+Strunk(z，θr)]] 公式十四
b(z，θr)＝ln[1-(1-a′(z，θr)(1-c)]/(1-c) 公式十五
其中，Scrown为第z层冠层顶部的投影面积；Strunk为第z层的树干投影面积；a′(z，θr)为第z树冠在θr方向上的孔隙率，a′(z，θr)公式具体如下：
Figure BDA0002724032160000051
其中，V(z)为z高度以上树冠的体积。
进一步地，直接辐射模拟中，中层树干可见比例的具体计算如下：
将整个树干划分为不同的层次，光照树干和阴影树干的可见比例表示为：
ft(r2)＝λ∫p′00(z，r2，r2)d(z)dz 公式十七
ft，s(r1，r2)＝λ∫p′00(z，r1，r2)d(z)T(r1，r2)dz 公式十八
ft，h(r1，r2)＝ft(r2)-ft，s(r1，r2) 公式十九
其中，ft,ft，s和ft，h分别代表树干、光照树干和阴影树干的可见比例； T(r1，r2)为树干形状引起的光照-观测方位角纠正；p′00(z，r2，r2)代表从高度为z 的点的光照方向r1和观测方向r2计算的双向孔隙率，p′00由以下公式计算：
p′00(z，r1，r2)＝p′1(z，r1，r2)·p2(z′，r1，r2) 公式二十
其中，p′1代表光照和观测方向上树叶对自身冠层的双向孔隙率，p2表示相邻树冠引起的双向孔隙率，物理意义与树冠部分相同；如果从z高度的树干点发出的光在光照和观测方向均不能穿透树冠，z′＝z；如果从树干点发出的光在光照或者观测方向上接触到树冠，则相应的z′将成为其向上与树冠表面相交的高度。
进一步地，直接辐射模拟中，林下植被和底层土壤可见比例的具体计算如下：
光照叶片和阴影叶片的可见比例由下列公式求出：
fv，s(r2)＝∑fv，s，k(r1，r2) 公式二十一
fv，h(r2)＝∑fv，h，k(r1，r2) 公式二十二
fv，k(r2，r2)＝p2(zk，r1，r2)LAIkG(θv) 公式二十三
fv，k，s(r1，r2)＝p2(zk，r1，r2)LAIkG(θv) 公式二十四
fv，k，h(r1，r2)＝fv，k(r2，r2)-fv，s，k(r1，r2) 公式二十五
其中，fv，s和fv，h代表光照植被和阴影植被的可见比例，fv，k，fv，k，s和 fv，k，h分别在k层叶片可见比例，光照叶片可见比例和阴影叶片可见比例，k代表不同的层；LAIk代表第k层植被包络区域的叶面积指数；
光照土壤和阴影土壤的可见比例表示为：
fs(r2，r2)＝p2(0，r2，r2) 公式二十六
fs，s(r1，r2)＝p2(0，r1，r2) 公式二十七
fs，h(r1，r2)＝fs(r2，r2)-fs，s(r1，r2) 公式二十八
其中fs，fs，s和fs，h分别为土壤的可见比例，光照土壤的可见比例和阴影土壤的可见比例。
进一步地，多次散射模拟中，冠层内部的叶片的多次散射效应利用光谱不变理论进行模拟计算：确定冠层内部的几何结构，通过计算逐格网的半球透过率和方向透过率计算光子往上逃逸、往下方向逃逸和被格网内组分吸收的比例，基于该比例，计算组分间的散射情况。
进一步地，多次散射模拟中，对树冠体元进行光谱不变理论计算多次散射，计算树冠的多次散射εc，m；
Figure BDA0002724032160000071
Figure BDA0002724032160000072
pc＝1-eu，c-ed，c 公式三十一
εc，m＝icεcpc(1-εc)+(1-ic)(1-εu)i′cεc 公式三十二
其中，εc，m为树冠的多次散射贡献，ic为光子在观测方向的截断概率，i′G为光子在整个上半球空间的平均截断概率，为ic的积分结果；n为树冠内体体元的数目；p00为观测-太阳双向透过率，uL为树冠内叶片的体密度；λ为树冠的密度； G(θv)为叶片在θv方向的投影；au，c和ad，c为光子打到单个叶片上往上和往下散射的概率；hc为树冠的高度；x，y和z为体元的坐标；eu，c和ed，c为树冠内光子往上和往下散射的平均概率；pc为树冠内叶片发射光子打到其他叶片概率；εu为树冠下组分的等效发射率，εc为树冠叶片的发射率。
本发明通过建立冠层顶热辐射与冠层内各组分温度状态的联系，不仅适用性广泛且模拟精度高。本发明显著提高了森林冠层热辐射模拟的能力，有望用于森林冠层热辐射时空变异分析提高地表温度产品精度，有望用于森林场景组分温度反演开展干旱和火灾预警，有望用于森林生态系统的能量分布计算辅助林业管理。本发明在空间信息技术领域，尤其是定量遥感和林业领域具有重要的应用价值。
附图说明
图1为本发明离散森林场景热红外辐射传输模拟方法框架图。
图2为离散森林场景组分格网划分方法。
图3为树冠体元光子散射概念图。
图4为离散森林场景热辐射传输模型不同场景模拟结果对比图。
图5为离散森林场景热辐射传输模型与三维计算机模拟模型对比图。
图6为不同类型森林场景图。
具体实施方式
下面结合附图和具体实施方式对本发明作进一步详细的说明。
为了解决上述技术问题中的不足之处，本发明基于已有的辐射传输模型理论，将离散格网和光谱不变的理论进行有机结合，提出了一种可用于离散森林场景热红外辐射传输模拟方法。
如图1所示的一种离散森林场景热红外辐射传输模拟方法，包括以下步骤：
S1、考虑森林的三维结构，将整个森林场景按照类别和温度差异划分为四层，从上向下依次包括上层树冠、中层树干、林下植被和底层土壤；通过直接辐射模拟、多次散射模拟对各组分可见比例进行计算，最后计算出热红外波段冠层顶表面热辐射量；
S2、对每层的组分进行离散，按照几何结构划分为不同的格网或者子层；
S3、直接辐射模拟：按照格网内叶面积密度计算每层组分的光照和阴影可视比例，通过普朗克函数计算组分的直接发射项；
S4、多次散射模拟：按照光谱不变理论计算组分间的多次散射项；
S5、冠层顶发射项等于直接发射项和多次散射项之和。
步骤S1中，在热红外波段，表面辐射分为组分辐射和大气辐射两部分，特定波长的冠层顶表面热辐射量(L)由下式计算：
Figure BDA0002724032160000081
其中：θ1，θ2分别代表太阳和观测的天顶角；
Figure BDA0002724032160000092
表示观测角度和光线角度的相对方位角；Tj和εj代表j组分的温度和发射率；f代表组分的可见比例；εm代表多次散射效应，指从一个组分发射并由其他组分反射到传感器的辐射；εe代表冠层的有效发射率；
Figure 302426DEST_PATH_IMAGE002
代表大气有效下行辐射；B代表普朗克函数。其中
Figure BDA0002724032160000094
代表j组分的直接辐射贡献；εm，j为j组分的多次散射项贡献，计算出每个组分的直接辐射和多次散射贡献，然后将其累加/积分则为冠层顶的热辐射。
离散森林场景热辐射方向性模拟主要包括直接辐射模拟和多次散射模拟：
一、直接辐射模拟
树冠和树干被划分为不同的体积元，如图2所示。组分的观测、光照比例以及对于冠层顶部观测的贡献可以通过积分组建来计算相应体积元素的值。
Figure BDA0002724032160000091
在ThermalFRT 模型中，整个森林冠层被划分为四层：上层树冠、中层树干、林下植被和底层土壤。将上层树冠和林下植被的不同层次的叶片的光照和阴影部分确定为不同的组分。除去大气效应，森林的热红外方向辐射可以由直接辐射、多次散射效应的总和得到。假定组分均为朗伯体，则直接辐射部分主要受组分可见比例的影响。
直接辐射模拟中，上层树冠可见比例的具体计算如下：
在FRT模型基础上，加入组分的发射项，组分的发射项通过普朗克函数计算得到，组分的可视比例需要逐格网/层计算双向可视比例；假设各组分为朗伯体，树冠的双向孔隙率由冠层内和冠层间两个独立分布的概率相乘得到，逐格网计算，然后积分得到整个树冠的可视光照和可视阴影比例。
基于双向孔隙率计算出一个体元中光照和阴影树叶的可见比例，然后使用所有相应的体积元素的总和来计算光照和阴影树冠的可见比例，计算公式如下所示：
fc，s(r2)＝∑fc，k，s(r1，r2) 公式二
fc，h(r2)＝∑fc，k，h(r1，r2) 公式三
Figure BDA0002724032160000101
Figure BDA0002724032160000102
fc，k，h(r1，r2)＝fc，k(r2)-fc，k，s(r1，r2) 公式六
其中，fc，k,fc，k，h和fc，k，h分别代表第k层树冠的可见比例、光照树冠的可见比例和阴影树冠的可见比例；fc，s和fc，h为整个树冠的可见光照比例和可见阴影比例；(x，y，z)代表树冠内点M的坐标；Vk代表第k层树冠包络的空间区域； r1和r2分别代表光照和观测方向；λ代表树冠的密度；p00代表由点(x，y，z)的r1和r2计算出的双向孔隙率；uL代表叶片的面体积密度；G代表投影在观测或光照方向的叶片的比例，对于球形叶倾角分布函数(LIDF)，固定为0.5；G(θ2) 代表树叶在θ2方向的投影；
在近红外波段，经常使用冠层介质的散射面积相位函数。在热红外波段，假设所有组分均为朗伯体，叶片的贡献直接反映在其截面上，因此，双向孔隙率由两个独立分布概率的乘积表示，计算公式如下所示：
p00(x，y，z，r1，r2)＝p1(x，y，z，r1，r2)·p2(z，r1，r2) 公式七
其中，p1表示冠层内双向孔隙率；p2表示冠层间双向孔隙率；
假设在树冠中均为均质的叶片，p1主要取决于点在树冠中的位置和叶面体积密度，p2主要取决于树密度和叶面体积密度，计算p1和p2的公式如下所示：
p1(x，y，z，r1，r2)＝ai(s1，θ1)ai(s2，θ2)CHS1(s1，s2，θ1，θ2) 公式八
ai(s，θ)＝exp[-suLG(θ)] 公式九
Figure BDA0002724032160000111
Figure BDA0002724032160000112
其中，ai(s，θ)代表方向为θ的点到冠层表面的间隙距离；CHS1为热点校正因子，描述了冠层内部叶片在光照和观测投影内的重叠效应；θ1和θ2为太阳天顶角和观测天顶角；sAB为路径长度的重叠；s1和s2分别代表从点(x，y，z)沿光照和观测方向到冠层顶面之间的路径长度；α表示r1和r2之间的立体角；b为热点因子，特定的球形叶倾角分布函数为16/(dLπ2)，水平叶倾角分布函数为 4/(dLπ)；p2用于表示树冠之间的多次阴影效应，由树冠中向上的点 M1(x1，y1，z1)和M2(x2，y2，z2)计算：
p2(z，r1，r2)＝as(z1，θ1)as(z2，θ2)CHS2(z1，z2，l12，r1，r2) 公式十二
CHS2(z1，z2，l12，r1，r2)＝exp[λcSc(z1，z2，l12，r1，r1)p0] 公式十三
其中，as(z，θ)是在方向θ、高度为z的冠层间隙的平均比例；CHS2也为热点校正因子，描述光照和观测方向上阴影冠层内部特征；参数c表征树的分布模式与泊松分布的偏差；Sc(z1，z2，l12，r1，r1)为光照和观测方向上树冠投影区域的公共部分面积，对应高度为z1，z2以及水平距离l12；p0为树冠内的联合孔隙率；孔隙率as(z，θr)在假设树为二项分布的情况下，计算如下：
as(z，θr)＝exp[-λ[b(z，θr)Scrown(z，θr)+Strunk(z，θr)]] 公式十四
b(z，θr)＝ln[1-(1-a′(z，θr)(1-c)]/(1-c) 公式十五
其中，Scrown为第Z层冠层顶部的投影面积；Strunk为第z层的树干投影面积；a′(z，θr)为第z树冠在θr方向上的孔隙率，a′(z，θr)公式具体如下：
Figure BDA0002724032160000113
其中，V(z)为z高度以上树冠的体积。
直接辐射模拟中，中层树干可见比例的具体计算如下：
在近红外FRT模型中，由于阴影树干的辐射主要由多次散射效应引起，因此只考虑了部分树干的光照和阴影部分。但在热红外波段，阴影树干也会因为发射能量而对冠层顶辐射起到重要的作用。特别是在落叶林的落叶期，树干的每一部分都会被传感器所接收。在Thermal FRT模型中，将整个树干考虑在内，并划分为不同的层次。
将整个树干划分为不同的层次，光照树干和阴影树干的可见比例表示为：
ft(r2)＝λ∫p′00(z，r2，r2)d(z)dz 公式十七
ft，s(r1，r2)＝λ∫p′00(z，r1，r2)d(z)T(r1，r2)dz 公式十八
ft，h(r1，r2)＝ft(r2)-ft，s(r1，r2) 公式十九
其中，ft,ft，s和ft，h分别代表树干、光照树干和阴影树干的可见比例； T(r1，r2)为树干形状引起的光照-观测方位角纠正；p′00(z，r2，r2)代表从高度为z 的点的光照方向r1和观测方向r2计算的双向孔隙率，p′00由以下公式计算：
p′00(z，r1，r2)＝p′1(z，r1，r2)·p2(z′，r1，r2)公式二十
其中，p′1代表光照和观测方向上树叶对自身冠层的双向孔隙率，p2表示相邻树冠引起的双向孔隙率，物理意义与树冠部分相同；考虑到树冠与树干投影的几何关系，需要分别计算树干的三个部分：上层部分，光照和观测方向均被树冠遮挡；中层部分，单个方向被遮挡、另一方向不被遮挡；底层部分，光照和观测方向均不被树冠遮挡。计算出三部分结果后加权计算。如果从z高度的树干点发出的光在光照和观测方向均不能穿透树冠，z′＝z；如果从树干点发出的光在光照或者观测方向上接触到树冠，则相应的z′将成为其向上与树冠表面相交的高度。
直接辐射模拟中，林下植被和底层土壤可见比例的具体计算如下：
林下植被和底层土壤的可见比例也可以根据双向孔隙率来计算。由于林下植被可能存在温度的垂直异质性，因此考虑了不同层次的林下叶片。光照叶片和阴影叶片的可见比例由下列公式求出：
fv，s(r2)＝∑fv，s，k(r1，r2) 公式二十一
fv，h(r2)＝∑fv，h，k(r1，r2) 公式二十二
fv，k(r2，r2)＝p2(zk，r1，r2)LAIkG(θv) 公式二十三
fv，k，s(r1，r2)＝p2(zk，r1，r2)LAIkG(θv) 公式二十四
fv，k，h(r1，r2)＝fv，k(r2，r2)-fv，s，k(r1，r2) 公式二十五
其中，fv，s和fv，h代表光照植被和阴影植被的可见比例，fv，k，fv，k，s和 fv，k，h分别在k层叶片可见比例，光照叶片可见比例和阴影叶片可见比例，k代表不同的层；LAIk代表第k层植被包络区域的叶面积指数；
光照土壤和阴影土壤的可见比例表示为：
fs(r2，r2)＝p2(0，r2，r2) 公式二十六
fs，s(r1，r2)＝p2(0，r1，r2) 公式二十七
fs，h(r1，r2)＝fs(r2，r2)-fs，s(r1，r2) 公式二十八
其中fs，fs，s和fs，h分别为土壤的可见比例，光照土壤的可见比例和阴影土壤的可见比例。
二、多次散射计算
冠层内部的叶片的多次散射效应利用光谱不变理论进行模拟计算。冠层内部的几何结构确定，通过计算逐格网的半球透过率和方向透过率计算光子往上逃逸、往下方向逃逸和被格网内组分吸收的比例。基于该比例，计算组分间的散射情况。以树冠体元为例进行光谱不变理论计算多次散射的说明，如图3所示，为树冠内每个体元光子往上逃逸和往下逃逸的示意图，实线和虚线表示往上、逃逸路径长度。按照笛卡尔坐标系z轴的方向可以以体元为中心划分为上半球空间和下半球空间，光子的路径长度如图中箭头所示，假设冠层内叶片的分布是随机且均质的，光子不被叶片截断的概率只与路径长度有关，而该路径长度随树冠形状和散射方向所决定，因此进行整个半球的积分，随后计算树冠的多次散射εc，m。
多次散射模拟中，冠层内部的叶片的多次散射效应利用光谱不变理论进行模拟计算：确定冠层内部的几何结构，通过计算逐格网的半球透过率和方向透过率计算光子往上逃逸、往下方向逃逸和被格网内组分吸收的比例，基于该比例，计算组分间的散射情况。
多次散射模拟中，对树冠体元进行光谱不变理论计算多次散射，计算树冠的多次散射εc，m；
Figure BDA0002724032160000141
Figure BDA0002724032160000142
pc＝1-eu，c-ed，c 公式三十一
εc，m＝icεcpc(1-εc)+(1-ic)(1-εu)i′cεc 公式三十二
其中，εc，m为树冠的多次散射贡献，ic为光子在观测方向的截断概率，i′G为光子在整个上半球空间的平均截断概率，为ic的积分结果；n为树冠内体体元的数目；p00为观测-太阳双向透过率，uL为树冠内叶片的体密度；λ为树冠的密度； G(θv)为叶片在θv方向的投影；au，c和ad，c为光子打到单个叶片上往上和往下散射的概率；hc为树冠的高度；x，y和z为体元的坐标；eu，c和ed，c为树冠内光子往上和往下散射的平均概率；pc为树冠内叶片发射光子打到其他叶片概率；εu为树冠下组分的等效发射率，εc为树冠叶片的发射率。
由于考虑了森林中树冠的几何结构和垂直分层现象，相对于已有的只有树冠和底层土壤的模型，热辐射方向性模型可以模拟更加复杂的森林冠层的热辐射如图4所示。其中，SZA为太阳天顶角，tlai为树冠的叶面积指数，ulai为树冠下植被的叶面积指数，alp为树冠的密度；ht为树冠的高度；wid为树干的宽度。因此，(a)和(b)表示场景在树冠叶面积指数为0.5情况下，树冠密度分别为0.025和0.075；(c)和(d)表示场景在树冠叶面积指数为3.0情况下，树冠密度分别为0.025和0.075。Mod5crown表示场景中只有土壤和树冠；Mod6crown+uveg表示场景中只有土壤，树冠和林下植被；Mod7crown+trunk表示场景中只有土壤、树冠和树干；Mod8crown+trunk+uveg表示场景中有土壤，树冠、树干和林下植被；Brightness Temperature表示森林冠层顶亮度温度；VZA表示观测的角度。
另外，本发明提出的离散森林场景的模型模拟精度高，与三维计算机模拟模型模拟结果相比，二者的一致性好如图5所示，误差小于0.5℃。相比较图中最高和最低温度差异，即角度效应偏差，可达5℃，本发明提出的该方法是进行角度纠正和组分温度反演研究的重要工具。在实际情况下正午稀疏森林该偏差可达5-10℃。
其中，SZA为太阳天顶角，tlai为树冠的叶面积指数，ulai为树冠下植被的叶面积指数，alp为树冠的密度；ht为树冠的高度；wid为树干的宽度。因此，(a)和(b)表示场景在树冠叶面积指数为0.5情况下，树冠密度分别为 0.035和0.070；(c)和(d)表示场景在树冠叶面积指数为3.0情况下，树冠密度分别为0.035和0.070。Brightness Temperature表示森林冠层顶亮度温度；VZA表示观测的角度。TRGM表示三维模型；thermal FRT表示本专利的方法。RMSE表示均方根误差；Bias表示偏差；r2表示决定系数；σ表示离散方差；
通过本发明的研究内容，可以用于但不限于以下森林逻辑抽象场景的模拟：如图6所示，(a)均质树冠+土壤；(b)均质树冠+林下植被+土壤；(c) 均质树冠+林下树干+土壤；(d)均质树冠+林下树干+林下植被+土壤；(e) 椭球树冠+土壤；(f)椭球树冠+林下植被+土壤；(g)椭球树冠+林下树干+ 土壤；(h)椭球树冠+林下树干+林下植被+土壤。
总体而言，在热红外辐射传输理论中森林冠层的异质性可以归结为两个方面：冠层结构和组分温度分布。在Thermal FRT(forest reflectance and transmittance)模型中，森林冠层的分布不仅表现为组分之间的温度差异，还表现为组分之间的温差。由于考虑了垂直方向和水平方向的差异，Thermal FRT模型可以用于不同森林冠层下的亮温的方向性研究，且研究表明其具有良好的稳定性，可以为研究角度归一化算法和组分温度反演提供参考数据。
上述实施方式并非是对本发明的限制，本发明也并不仅限于上述举例，本技术领域的技术人员在本发明的技术方案范围内所做出的变化、改型、添加或替换，也均属于本发明的保护范围。

Claims (7)

1.一种离散森林场景热红外辐射传输模拟方法，其特征在于：所述方法包括以下步骤：
S1、考虑森林的三维结构，将整个森林场景按照类别和温度差异划分为四层，从上向下依次包括上层树冠、中层树干、林下植被和底层土壤；通过直接辐射模拟、多次散射模拟对各组分可见比例进行计算，最后计算出热红外波段冠层顶表面热辐射量；
S2、对每层的组分进行离散，按照几何结构划分为不同的格网或者子层；
S3、直接辐射模拟：按照格网内叶面积密度计算每层组分的光照和阴影可视比例，通过普朗克函数计算组分的直接发射项；
S4、多次散射模拟：按照光谱不变理论计算组分间的多次散射项；
S5、冠层顶发射项等于直接发射项和多次散射项之和。
2.根据权利要求1所述的离散森林场景热红外辐射传输模拟方法，其特征在于：所述步骤S1中，在热红外波段，表面辐射分为组分辐射和大气辐射两部分，特定波长的冠层顶表面热辐射量(L)由下式计算：
Figure FDA0002724032150000011
其中：θ1，θ2分别代表太阳和观测的天顶角；
Figure FDA0002724032150000012
表示观测角度和光线角度的相对方位角；Tj和εj代表j组分的温度和发射率；f代表组分的可见比例；εm代表多次散射效应，指从一个组分发射并由其他组分反射到传感器的辐射；εe代表冠层的有效发射率；
Figure DEST_PATH_IMAGE002
代表大气有效下行辐射；B代表普朗克函数。其中
Figure FDA0002724032150000014
代表j组分的直接辐射贡献；εm，j为j组分的多次散射项贡献，计算出每个组分的直接辐射和多次散射贡献，然后将其累加/积分则为冠层顶的热辐射。
3.根据权利要求2所述的离散森林场景热红外辐射传输模拟方法，其特征在于：所述直接辐射模拟中，上层树冠可见比例的具体计算如下：
基于双向孔隙率计算出一个体元中光照和阴影树叶的可见比例，然后使用所有相应的体积元素的总和来计算光照和阴影树冠的可见比例，计算公式如下所示：
fc，s(r2)＝∑fc，k，s(r1，r2) 公式二
fc，h(r2)＝∑fc，k，h(r1，r2) 公式三
Figure FDA0002724032150000021
Figure FDA0002724032150000022
fc，k，h(r1，r2)＝fc，k(r2)-fc，k，s(r1，r2) 公式六
其中，fc，k，fc，k，s和fc，k，h分别代表第k层树冠的可见比例、光照树冠的可见比例和阴影树冠的可见比例；fc，s和fc，h为整个树冠的可见光照比例和可见阴影比例；(x，y，z)代表树冠内点M的坐标；Vk代表第k层树冠包络的空间区域；r1和r2分别代表光照和观测方向；λ代表树冠的密度；p00代表由点(x，y，z)的r1和r2计算出的双向孔隙率；uL代表叶片的面体积密度；G代表投影在观测或光照方向的叶片的比例，对于球形叶倾角分布函数(LIDF)，固定为0.5；G(θ2)代表树叶在θ2方向的投影；
在热红外波段，假设所有组分均为朗伯体，叶片的贡献直接反映在其截面上，因此，双向孔隙率由两个独立分布概率的乘积表示，计算公式如下所示：
p00(x，y，z，r1，r2)＝p1(x，y，z，r1，r2)·p2(z，r1，r2) 公式七
其中，p1表示冠层内双向孔隙率；p2表示冠层间双向孔隙率；
假设在树冠中均为均质的叶片，p1主要取决于点在树冠中的位置和叶面体积密度，p2主要取决于树密度和叶面体积密度，计算p1和p2的公式如下所示：
p1(x，y，z，r1，r2)＝ai(s1，θ1)ai(s2，θ2)CHS1(s1，s2，θ1，θ2) 公式八
ai(s，θ)＝exp[-suLG(θ)] 公式九
Figure FDA0002724032150000031
Figure FDA0002724032150000032
其中，ai(s，θ)代表方向为θ的点到冠层表面的间隙距离；CHS1为热点校正因子，描述了冠层内部叶片在光照和观测投影内的重叠效应；θ1和θ2为太阳天顶角和观测天顶角；sAB为路径长度的重叠；s1和s2分别代表从点(x，y，z)沿光照和观测方向到冠层顶面之间的路径长度；α表示r1和r2之间的立体角；b为热点因子，特定的球形叶倾角分布函数为16/(dLπ2)，水平叶倾角分布函数为4/(dLπ)；p2用于表示树冠之间的多次阴影效应，由树冠中向上的点M1(x1，y1，z1)和M2(x2，y2，z2)计算：
p2(z，r1，r2)＝as(z1，θ1)as(z2，θ2)CHS2(z1，z2，l12，r1，r2) 公式十二
CHS2(z1，z2，l12，r1，r2)＝exp[λcSc(z1，z2，l12，r1，r1)p0] 公式十三
其中，as(z，θ)是在方向θ、高度为z的冠层间隙的平均比例；CHS2也为热点校正因子，描述光照和观测方向上阴影冠层内部特征；参数c表征树的分布模式与泊松分布的偏差；Sc(z1，z2，l12，r1，r1)为光照和观测方向上树冠投影区域的公共部分面积，对应高度为z1，z2以及水平距离l12；p0为树冠内的联合孔隙率；孔隙率as(z，θr)在假设树为二项分布的情况下，计算如下：
as(z，θr)＝exp[-λ[b(z，θr)Scrown(z，θr)+Strunk(z，θr)]] 公式十四
b(z，θr)＝ln[1-(1-a′(z，θr)(1-c)]/(1-c) 公式十五
其中，Scrown为第z层冠层顶部的投影面积；Strunk为第z层的树干投影面积；a′(z，θr)为第z树冠在θr方向上的孔隙率，a′(z，θr)公式具体如下：
Figure FDA0002724032150000041
其中，V(z)为z高度以上树冠的体积。
4.根据权利要求3所述的离散森林场景热红外辐射传输模拟方法，其特征在于：所述直接辐射模拟中，中层树干可见比例的具体计算如下：
将整个树干划分为不同的层次，光照树干和阴影树干的可见比例表示为：
ft(r2)＝λ∫p′00(z，r2，r2)d(z)dz 公式十七
ft，s(r1，r2)＝λ∫p′00(z，r1，r2)d(z)T(r1，r2)dz 公式十八
ft，h(r1，r2)＝ft(r2)-ft，s(r1，r2) 公式十九
其中，ft，ft，s和ft，h分别代表树干、光照树干和阴影树干的可见比例；T(r1，r2)为树干形状引起的光照-观测方位角纠正；p′00(z，r2，r2)代表从高度为z的点的光照方向r1和观测方向r2计算的双向孔隙率，p′00由以下公式计算：
p′00(z，r1，r2)＝p′1(z，r1，r2)·p2(z′，r1，r2) 公式二十
其中，p′1代表光照和观测方向上树叶对自身冠层的双向孔隙率，p2表示相邻树冠引起的双向孔隙率，物理意义与树冠部分相同；如果从z高度的树干点发出的光在光照和观测方向均不能穿透树冠，z′＝z；如果从树干点发出的光在光照或者观测方向上接触到树冠，则相应的z′将成为其向上与树冠表面相交的高度。
5.根据权利要求4所述的离散森林场景热红外辐射传输模拟方法，其特征在于：所述直接辐射模拟中，林下植被和底层土壤可见比例的具体计算如下：
光照叶片和阴影叶片的可见比例由下列公式求出：
fv，s(r2)＝∑fv，s，k(r1，r2) 公式二十一
fv，h(r2)＝∑fv，h，k(r1，r2) 公式二十二
fv，k(r2，r2)＝p2(zk，r1，r2)LAIkG(θv) 公式二十三
fv，k，s(r1，r2)＝p2(zk，r1，r2)LAIkG(θv) 公式二十四
fv，k，h(r1，r2)＝fv，k(r2，r2)-fv，s，k(r1，r2) 公式二十五
其中，fv，s和fv，h代表光照植被和阴影植被的可见比例，fv，k，fv，k，s和fv，k，h分别在k层叶片可见比例，光照叶片可见比例和阴影叶片可见比例，k代表不同的层；LAIk代表第k层植被包络区域的叶面积指数；
光照土壤和阴影土壤的可见比例表示为：
fs(r2，r2)＝p2(0，r2，r2) 公式二十六
fs，s(r1，r2)＝p2(0，r1，r2) 公式二十七
fs，h(r1，r2)＝fs(r2，r2)-fs，s(r1，r2) 公式二十八
其中fs，fs，s和fs，h分别为土壤的可见比例，光照土壤的可见比例和阴影土壤的可见比例。
6.根据权利要求2所述的离散森林场景热红外辐射传输模拟方法，其特征在于：所述多次散射模拟中，冠层内部的叶片的多次散射效应利用光谱不变理论进行模拟计算：确定冠层内部的几何结构，通过计算逐格网的半球透过率和方向透过率计算光子往上逃逸、往下方向逃逸和被格网内组分吸收的比例，基于该比例，计算组分间的散射情况。
7.根据权利要求6所述的离散森林场景热红外辐射传输模拟方法，其特征在于：所述多次散射模拟中，对树冠体元进行光谱不变理论计算多次散射，计算树冠的多次散射εc，m；
Figure FDA0002724032150000051
Figure FDA0002724032150000061
Figure FDA0002724032150000062
pc＝1-eu，c-ed，c 公式三十一
εc，m＝icεcpc(1-εc)+(1-ic)(1-εu)i′cεc 公式三十二
其中，εc，m为树冠的多次散射贡献，ic为光子在观测方向的截断概率，i′c为光子在整个上半球空间的平均截断概率，为ic的积分结果；n为树冠内体体元的数目；p00为观测-太阳双向透过率，uL为树冠内叶片的体密度；λ为树冠的密度；G(θv)为叶片在θv方向的投影；au，c和ad，c为光子打到单个叶片上往上和往下散射的概率；hc为树冠的高度；x，y和z为体元的坐标；eu，c和ed，c为树冠内光子往上和往下散射的平均概率；pc为树冠内叶片发射光子打到其他叶片概率；εu为树冠下组分的等效发射率，εc为树冠叶片的发射率。
CN202011096826.4A 2020-10-14 2020-10-14 一种离散森林场景热红外辐射传输模拟方法 Active CN112254820B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011096826.4A CN112254820B (zh) 2020-10-14 2020-10-14 一种离散森林场景热红外辐射传输模拟方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011096826.4A CN112254820B (zh) 2020-10-14 2020-10-14 一种离散森林场景热红外辐射传输模拟方法

Publications (2)

Publication Number Publication Date
CN112254820A true CN112254820A (zh) 2021-01-22
CN112254820B CN112254820B (zh) 2021-04-27

Family

ID=74243103

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011096826.4A Active CN112254820B (zh) 2020-10-14 2020-10-14 一种离散森林场景热红外辐射传输模拟方法

Country Status (1)

Country Link
CN (1) CN112254820B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115062260A (zh) * 2022-06-16 2022-09-16 电子科技大学 一种适用于异质森林的森林生物量PolInSAR估算方法和系统、存储介质
CN116881609A (zh) * 2023-09-06 2023-10-13 合肥工业大学 一种普适性的森林冠层间隙率的计算方法

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20020123848A1 (en) * 2000-11-20 2002-09-05 Eva Schneiderman Predictive method for polymers
CN102034027A (zh) * 2010-12-16 2011-04-27 南京大学 流域尺度土壤湿度遥感数据同化方法
CN102162844A (zh) * 2010-12-07 2011-08-24 北京理工大学 一种sar大范围森林场景遥感数据的快速模拟方法
CN102914501A (zh) * 2012-07-26 2013-02-06 南京大学 一种利用激光点云计算三维森林冠层消光系数的方法
CN103557881A (zh) * 2013-11-11 2014-02-05 北京林业大学 一种大尺度森林观测梯形网格法
US20150308896A1 (en) * 2012-06-05 2015-10-29 Hypermed Imaging, Inc. Methods and apparatus for coaxial imaging of multiple wavelengths
CN106909750A (zh) * 2017-03-10 2017-06-30 长沙理工大学 一种阔叶植被冠层反射率的计算方法及模型
CN107767364A (zh) * 2017-09-12 2018-03-06 中国林业科学研究院林业研究所 基于红外热图像精准提取树木冠层温度的方法
CN108009392A (zh) * 2017-10-19 2018-05-08 桂林航天工业学院 一种浓密植被地表的遥感反射率模型构建及标定应用方法

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20020123848A1 (en) * 2000-11-20 2002-09-05 Eva Schneiderman Predictive method for polymers
CN102162844A (zh) * 2010-12-07 2011-08-24 北京理工大学 一种sar大范围森林场景遥感数据的快速模拟方法
CN102034027A (zh) * 2010-12-16 2011-04-27 南京大学 流域尺度土壤湿度遥感数据同化方法
US20150308896A1 (en) * 2012-06-05 2015-10-29 Hypermed Imaging, Inc. Methods and apparatus for coaxial imaging of multiple wavelengths
CN102914501A (zh) * 2012-07-26 2013-02-06 南京大学 一种利用激光点云计算三维森林冠层消光系数的方法
CN103557881A (zh) * 2013-11-11 2014-02-05 北京林业大学 一种大尺度森林观测梯形网格法
CN106909750A (zh) * 2017-03-10 2017-06-30 长沙理工大学 一种阔叶植被冠层反射率的计算方法及模型
CN107767364A (zh) * 2017-09-12 2018-03-06 中国林业科学研究院林业研究所 基于红外热图像精准提取树木冠层温度的方法
CN108009392A (zh) * 2017-10-19 2018-05-08 桂林航天工业学院 一种浓密植被地表的遥感反射率模型构建及标定应用方法

Non-Patent Citations (6)

* Cited by examiner, † Cited by third party
Title
HUAGUO HUANG: "RAPID: A Radiosity Applicable to Porous IndiviDual Objects for directional reflectance over complex vegetated scenes", 《 REMOTE SENSING OF ENVIRONMENT》 *
TAO YU: "Modeling directional brightness temperature over a maize canopy in row structure", 《GEOSCIENCE AND REMOTE SENSING》 *
TIMMERMANS J: "Retrieval of canopy component temperatures through Bayesian inversion of directional thermal measurements", 《HYDROLOGY AND EARTH SYSTEM SCIENCES》 *
卞尊健: "植被冠层热辐射方向性建模及组分温度反演", 《中国博士学位论文全文数据库 基础科学辑》 *
杨贵军: "植被冠层3D辐射传输模型及热辐射方向性模拟", 《红外与毫米波学报》 *
黄华国: "利用多角度热图像提取冠层组分温度和方向亮温", 《北京师范大学学报(自然科学版)》 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115062260A (zh) * 2022-06-16 2022-09-16 电子科技大学 一种适用于异质森林的森林生物量PolInSAR估算方法和系统、存储介质
CN116881609A (zh) * 2023-09-06 2023-10-13 合肥工业大学 一种普适性的森林冠层间隙率的计算方法
CN116881609B (zh) * 2023-09-06 2023-11-21 合肥工业大学 一种普适性的森林冠层间隙率的计算方法

Also Published As

Publication number Publication date
CN112254820B (zh) 2021-04-27

Similar Documents

Publication Publication Date Title
Liu et al. Canopy occupation volume as an indicator of canopy photosynthetic capacity
CN112254820B (zh) 一种离散森林场景热红外辐射传输模拟方法
Sinoquet et al. Characterization of the light environment in canopies using 3D digitising and image processing
Pearcy et al. The functional morphology of light capture and carbon gain in the Redwood forest understorey plant Adenocaulon bicolor Hook
Dubayah et al. Topographic solar radiation models for GIS
Courbaud et al. Simulating radiation distribution in a heterogeneous Norway spruce forest on a slope
Green et al. Modeling light interception and transpiration of apple tree canopies
Malenovský et al. Discrete anisotropic radiative transfer modelling of solar-induced chlorophyll fluorescence: Structural impacts in geometrically explicit vegetation canopies
JP2003099697A (ja) 熱環境の予測方法、およびプログラム
Kükenbrink et al. Mapping the irradiance field of a single tree: Quantifying vegetation-induced adjacency effects
Mariscal et al. Light-transmission profiles in an old-growth forest canopy: simulations of photosynthetically active radiation by using spatially explicit radiative transfer models
Bian et al. A TIR forest reflectance and transmittance (FRT) model for directional temperatures with structural and thermal stratification
Béland et al. Vertical structure heterogeneity in broadleaf forests: Effects on light interception and canopy photosynthesis
Ma et al. Characterizing the three-dimensional spatiotemporal variation of forest photosynthetically active radiation using terrestrial laser scanning data
Kimes Radiative transfer in homogeneous and heterogeneous vegetation canopies
Chen et al. Unravelling the multilevel and multi-dimensional impacts of building and tree on surface urban heat islands
Stockle Canopy photosynthesis and transpiration estimates using radiation interception models with different levels of detail
Kirchner et al. Simulated directional radiances of vegetation from satellite platforms
Ogunjemiyo et al. Reflections in bumpy terrain: Implications of canopy surface variations for the radiation balance of vegetation
Mialet-Serra et al. Using plant architectural models for estimation of radiation transfer in a coconut-based agroforestry system
Stenberg et al. Modeling the spectral signature of forests: Application of remote sensing models to coniferous canopies
Rowe Incorporating landscape heterogeneity in land surface albedo models
Zeng et al. A radiative transfer model for patchy landscapes based on stochastic radiative transfer theory
Xie et al. Research on PAR and FPAR of crop canopies based on RGM
CN112198814A (zh) 一种高精度的热辐射方向性半经验半物理模拟方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant