CN101477197A - 用于林地复杂场景高光谱遥感数据的仿真方法 - Google Patents

Abstract

用于林地复杂场景高光谱遥感数据的仿真方法。本发明涉及高光谱遥感数据仿真的技术领域。它解决了现有的高光谱遥感数据仿真方法没能考虑地物的高度因素而不适用于林地附近存在高大建筑物情况下的高光谱遥感数据仿真的问题，步骤如下：一、设定仿真参数，并输入仿真中用到的光谱库地物反射率参考曲线；二、将欲仿真得到的高光谱遥感数据分解为太阳辐射直接被目标象元反射的能量、建筑物反射的太阳辐射被目标象元再次反射的能量、天空散射辐射被目标象元反射的能量、程辐射和临近效应，并按照步骤一中仿真参数分别进行仿真；三、将步骤二中获得的仿真结果组合得到最终的仿真数据。利用该方法为高光谱遥感数据后处理的研究提供经济且高质量的仿真数据。

Description

用于林地复杂场景高光谱遥感数据的仿真方法

技术领域

本发明涉及的是高光谱遥感数据仿真的技术领域。

背景技术

高光谱遥感技术是将常规遥感成像技术与光谱技术相结合产生的新型光学遥感技术，它的出现大大拓宽了光学遥感的应用领域。高光谱遥感可以在对目标地物的空间特征成像的同时，获取每个像元近乎连续的窄波段光谱信息。这样就可以既获得地面目标的空间特征，又获得其丰富的光谱信息，进而定性定量地对被测对象进行物理分析和识别，从而使基于特征光谱信息的目标地物遥感识别成为可能。这使得高光谱遥感技术在矿质调查、水环境、生态、大气等领域都具有无法估量的作用。

正是由于高光谱遥感系统所具有的这些其他遥感手段所没有的卓越特性，其一问世就引起世界各国研究者的广泛关注。目前高光谱遥感系统核心部件——成像光谱仪的研制已经达到了相当高的水平。然而由于高光谱遥感的运行成本过高，导致现有高光谱遥感数据的质量和数量难以满足数据处理研究的要求，造成高光谱数据的处理技术显得相对滞后。因此，通过对光谱成像系统进行研究，对其建立仿真模型并获取仿真数据，是一种值得深入研究的方法。

目前已有的高光谱遥感仿真数据获取方法按照对地表情况的假设可分为四大类，基于均一、郎伯面地表的高光谱遥感数据仿真方法，基于均一、非郎伯面地表的高光谱遥感数据仿真方法和基于非均一、郎伯面地表的高光谱遥感数据仿真方法，以及基于非均一、非郎伯面地表的高光谱遥感数据仿真方法。以上所提到的均一性是指在传感器分辨率所对应的地表分辨单元内仅有一类地物，郎伯面则是指完全的漫反射面。

对于建筑物附近的人造林地这种的场景，单个象元中包含树木和土壤两种地物，所以不符合均一性假设；同时，作为树木，其不同能量入射方向和观察方向的双向反射率并不相同，也不满足郎伯面的假设。故对于建筑物附近的人造林地这种的情况只能采用基于非均一、非郎伯面地表的高光谱遥感数据仿真方法。现有该类方法非常少，2003年Verhoef等人提出了一种通过将PROSPECT模型、GeoSAIL模型和四流地表-大气辐射传输模型相耦合仿真植被地表高光谱遥感数据的方法，但该方法没能考虑地物的高度因素，亦不适用于林地附近存在高大建筑物情况下的高光谱遥感数据仿真。

发明内容

本发明为了解决现有的高光谱遥感数据仿真方法没能考虑地物的高度因素而不适用于林地附近存在高大建筑物情况下的高光谱遥感数据仿真的问题，而提出了一种用于林地复杂场景高光谱遥感数据的仿真方法。

本发明的步骤如下：

步骤一：设定仿真参数，并输入仿真中用到的光谱库地物反射率参考曲线；

步骤二：将欲仿真得到的高光谱遥感数据分解为太阳辐射直接被目标象元反射的能量、建筑物反射的太阳辐射被目标象元再次反射的能量、天空散射辐射被目标象元反射的能量、程辐射和临近效应五个部分，并按照步骤一中设定的仿真参数分别进行仿真；

步骤三：将步骤二中获得的仿真结果组合得到最终的仿真数据。

本方法用于林地复杂场景高光谱遥感数据的仿真，它全面考虑了地表混合象元反射率的非均一性、非郎伯性，以及大气程辐射、临近效应和建筑物二次入射等作用，考虑了地物的高度影响，能够更好的近似高光谱遥感的真实情况。利用该方法可以方便的为高光谱遥感数据后处理方法的研究提供经济且高质量的仿真数据。

附图说明

图1是本发明方法的流程图；图2是场景几何关系示意图，①是天空散射辐射被目标象元反射的能量、②是临近效应、③太阳辐射直接被目标象元反射的能量、④建筑物反射的太阳辐射被目标象元再次反射的能量、⑤程辐射；图3是四流地表-大气辐射传输模型示意图，图中A层是大气顶层、B层是大气层、C层是地表，各符号下标s表示传输方向为太阳下行辐射传输方向，下标o表示传感器上行辐射传输方向，上标-表示散射传输下行方向，上标+表示散射传输上行方向，括号中t表示高度处于大气顶层，b表示高度处于地表大气底层；图4是四流交叉影响示意图，圆弧部分表明该部分接受外界的能量，矩形部分标志表明该部分向外界传输能量，各箭头表示的能量流动；图5是建筑物反射的太阳辐射能量被目标象元再次反射的几何关系示意图。

具体实施方式

具体实施方式一：结合图1说明本实施方式，本实施方式的步骤如下：

步骤一：设定仿真参数，并输入仿真中用到的光谱库地物反射率参考曲线；

步骤二：将欲仿真得到的高光谱遥感数据分解为太阳辐射直接被目标象元反射的能量、建筑物反射的太阳辐射被目标象元再次反射的能量、天空散射辐射被目标象元反射的能量、程辐射和临近效应五个部分，并按照步骤一中设定的仿真参数分别进行仿真；

步骤三：将步骤二中获得的仿真结果组合得到最终的仿真数据。

具体实施方式二：本实施方式与具体实施方式一不同点在于步骤一中的仿真参数划分为两组，分别为：地表场景参数和大气场景参数；

地表场景参数为地面场景大小、分辨单元大小(GSD)、各单元地物分类、树木类型、树木平均高度、树木冠层尺寸、植株密度和建筑物高度；其中树木冠层尺寸简化为一个椭球体，其大小有长轴和短轴确定。

大气场景参数为太阳入射方向的天顶角、太阳入射方向的方向角、传感器观测方向的天顶角、传感器观测方向的方向角、MODTRAN参数、光谱范围和光谱分辨率。

其它步骤与具体实施方式一相同。

具体实施方式三：结合图2至图5说明本实施方式，本实施方式与具体实施方式一不同点在于步骤二，步骤二是全面考虑各种因素对于高光谱遥感数据获取的影响，如图2是将高光谱遥感数据划分为太阳辐射直接被目标象元反射的能量

建筑物反射的太阳辐射被目标象元再次反射的能量

天空散射辐射被目标象元反射的能量程辐射L_λ，sky和临近效应L_{λ，adjacent}五个部分。于是总辐射也就是仿真得到的最终数据传感器处辐照能量L_λ，sensor便可由公式(1)表示：

$L_{\mathrm{\λ},\mathrm{sensor}}=L_{\mathrm{\λ},R}^I+L_{\mathrm{\λ},R}^{\mathrm{II}}+L_{\mathrm{\λ},R}^{\mathrm{III}}+L_{\mathrm{\λ},\mathrm{adjacent}}+L_{\mathrm{\λ},\mathrm{sky}}---\left(1\right)$

接下来分别仿真以上五个部分。

如图3和图4采取了四流地表-大气辐射传输模型表示地表反射率均为0、0.5、1三种情况下传感器处辐照能量L_λ，sensor与太阳入射到大气顶层的辐照能量E_s ^o的关系，如公式(2)

$L_{\mathrm{\λ},\mathrm{sensor}}={\mathrm{\ρ}}_{\mathrm{so}}\frac{{E_s}^o}{\mathrm{\π}}\cos {\mathrm{\θ}}_s+\frac{{\mathrm{\τ}}_{\mathrm{ss}}{\stackrel{\‾}{r}}_{\mathrm{sd}}+{\mathrm{\τ}}_{\mathrm{sd}}{\stackrel{\‾}{r}}_{\mathrm{dd}}}{1-{\stackrel{\‾}{r}}_{\mathrm{dd}}{\mathrm{\ρ}}_{\mathrm{dd}}}{\mathrm{\τ}}_{\mathrm{do}}\frac{{E_s}^o}{\mathrm{\π}}\cos {\mathrm{\θ}}_s+\frac{{\mathrm{\τ}}_{\mathrm{sd}}+{\mathrm{\τ}}_{\mathrm{ss}}{\stackrel{\‾}{r}}_{\mathrm{sd}}{\mathrm{\ρ}}_{\mathrm{dd}}}{1-{\stackrel{\‾}{r}}_{\mathrm{dd}}{\mathrm{\ρ}}_{\mathrm{dd}}}{r}_{\mathrm{do}}{\mathrm{\τ}}_{\mathrm{oo}}\frac{{E_s}^o}{\mathrm{\π}}\cos {\mathrm{\θ}}_s+{\mathrm{\τ}}_{\mathrm{ss}}{r}_{\mathrm{so}}{\mathrm{\τ}}_{\mathrm{oo}}\frac{{E_s}^o}{\mathrm{\π}}\cos {\mathrm{\θ}}_s---\left(2\right)$

其中等号右边第一项表示程辐射，第二项表示临近效应，第三项表示天空散射辐射被目标象元反射的能量，第四项表示太阳辐射直接被目标象元反射的能量。其中τ_ss，τ_oo，ρ_so，ρ_dd，τ_sd，τ_do6个参数描述大气的状况；r_so，r_do，r_dd，r_sd4个参数描述地表象元的反射率。

由于大气情况不取决于地物分类，故通过MODTRAN仿真得到获取地表所有象元反射率均为0、0.5或者1三种情况下的高光谱遥感数据，再利用仿真得到高光谱遥感数据反推描述大气状况的6个参数τ_ss，τ_oo，ρ_so，ρ_dd，τ_sd，τ_do，最后再把参数τ_ss，τ_oo，ρ_so，ρ_dd，τ_sd，τ_do应用到复杂情况下林地高光谱遥感数据的仿真中。

通过公式(2)得到如下公式：

$L_{\mathrm{PATH}}=[{\mathrm{\ρ}}_{\mathrm{so}}+\frac{{\mathrm{\τ}}_{\mathrm{ss}}{\stackrel{\‾}{r}}_{\mathrm{sd}}+{\mathrm{\τ}}_{\mathrm{sd}}{\stackrel{\‾}{r}}_{\mathrm{dd}}}{1-{\stackrel{\‾}{r}}_{\mathrm{dd}}{\mathrm{\ρ}}_{\mathrm{dd}}}{\mathrm{\τ}}_{\mathrm{do}}]\frac{{E_s}^o}{\mathrm{\π}}\cos {\mathrm{\θ}}_s---\left(3\right)$

$L_{\mathrm{GSUN}}={\mathrm{\τ}}_{\mathrm{ss}}{r}_{\mathrm{so}}{\mathrm{\τ}}_{\mathrm{oo}}\frac{{E_s}^o}{\mathrm{\π}}\cos {\mathrm{\θ}}_s---\left(4\right)$

$L_{\mathrm{GTOT}}=\frac{{\mathrm{\τ}}_{\mathrm{sd}}+{\mathrm{\τ}}_{\mathrm{ss}}{\stackrel{\‾}{r}}_{\mathrm{sd}}{\mathrm{\ρ}}_{\mathrm{dd}}}{1-{\stackrel{\‾}{r}}_{\mathrm{dd}}{\mathrm{\ρ}}_{\mathrm{dd}}}{r}_{\mathrm{do}}{\mathrm{\τ}}_{\mathrm{oo}}\frac{{E_s}^o}{\mathrm{\π}}\cos {\mathrm{\θ}}_s---\left(5\right)$

其中L_PATH表示程辐射与邻近效应能量的和，L_GSUN表示被目标反射的太阳直接入射的能量，L_GTOT表示被目标反射的天空的散射入射的能量；通过推导得到公式(6)

τ_ss＝exp(-b/cosθ_s)

τ_oo＝exp(-b/cosθ_o)

τ_ssτ_oo＝exp[-b(1/cosθ_s+1/cosθ_o)]＝GSUN₁₀₀×π/(E_s ^ocosθ_s)

ρ_so＝PATH₀×π/(E_s ^ocosθ_s)(6)

${\mathrm{\ρ}}_{\mathrm{do}}=\frac{{\mathrm{GTOT}}_{100}-2\×{\mathrm{GTOT}}_{50}}{{\mathrm{GTOT}}_{100}-{\mathrm{GTOT}}_{50}}$

${\mathrm{\τ}}_{\mathrm{sd}}=[\frac{{\mathrm{GTOT}}_{100}(1-{\mathrm{\ρ}}_{\mathrm{dd}})}{{\mathrm{GSUN}}_{100}}-{\mathrm{\ρ}}_{\mathrm{dd}}]{\mathrm{\τ}}_{\mathrm{ss}}$

${\mathrm{\τ}}_{\mathrm{do}}=\frac{{\mathrm{PATH}}_{100}-{\mathrm{PATH}}_0}{{\mathrm{GTOT}}_{100}}{\mathrm{\τ}}_{\mathrm{oo}}$

公式中，b是消光系数；通过MODTRAN仿真获得的PATH₀，PATH₅₀，PATH₁₀₀分别表示地表所有象元反射率均为0，0.5，1三种情况下对应的程辐射与邻近效应能量的和；GSUN₀，GSUN₅₀，GSUN₁₀₀分别表示地表所有象元反射率均为0，0.5，1三种情况下对应的被目标反射的太阳直接入射的能量；GTOT₀，GTOT₅₀，GTOT₁₀₀分别表示地表所有象元反射率均为0，0.5，1三种情况下对应的被目标反射的天空的散射入射的能量；

所述的地表象元反射率的4个参数r_so，r_do，r_dd，r_sd计算方法如下：

公式中，μ_o为传感器观测方向的天顶角的余弦值，μ_i为散射入射方向的天顶角的余弦值，μ_s为太阳入射方向的天顶角的余弦值，μ_v为散射观察方向的天顶角的余弦值，

为传感器观测方向的方向角，为散射入射方向的方向角，

为太阳入射方向的方向角，

为散射观察方向的方向角，ρ′为双向反射率由光谱库直接提供。

求得的6个描述大气的状况τ_ss，τ_oo，ρ_so，ρ_dd，τ_sd，τ_do的参数和4个描述地表象元的反射率的r_so，r_do，r_dd，r_sd参数后便通过公式(8)计算求得

L_{λ，adjacent}，L_λ，sky

$L_{\mathrm{\λ},R}^I={\mathrm{\τ}}_{\mathrm{ss}}{r}_{\mathrm{so}}{\mathrm{\τ}}_{\mathrm{oo}}\frac{{E_s}^o}{\mathrm{\π}}\cos {\mathrm{\θ}}_s$

$L_{\mathrm{\λ},R}^{\mathrm{III}}=\frac{{\mathrm{\τ}}_{\mathrm{sd}}+{\mathrm{\τ}}_{\mathrm{ss}}{\stackrel{\‾}{r}}_{\mathrm{sd}}{\mathrm{\ρ}}_{\mathrm{dd}}}{1-{\stackrel{\‾}{r}}_{\mathrm{dd}}{\mathrm{\ρ}}_{\mathrm{dd}}}{r}_{\mathrm{do}}{\mathrm{\τ}}_{\mathrm{oo}}\frac{{E_s}^o}{\mathrm{\π}}\cos {\mathrm{\θ}}_s$

$L_{\mathrm{\λ},\mathrm{adjacent}}=\frac{{\mathrm{\τ}}_{\mathrm{ss}}{\stackrel{\‾}{r}}_{\mathrm{sd}}{+\mathrm{\τ}}_{\mathrm{sd}}{\stackrel{\‾}{r}}_{\mathrm{dd}}}{1-{\stackrel{\‾}{r}}_{\mathrm{dd}}{\mathrm{\ρ}}_{\mathrm{dd}}}{\mathrm{\τ}}_{\mathrm{oo}}\frac{{E_s}^o}{\mathrm{\π}}\cos {\mathrm{\θ}}_s$

$L_{\mathrm{\λ},\mathrm{sky}}={\mathrm{\ρ}}_{\mathrm{so}}\frac{{E_s}^o}{\mathrm{\π}}\cos {\mathrm{\θ}}_s$ (8)

公式(2)中为建筑物反射的太阳辐射被目标象元再次反射的能量。

为计算

将建筑简化为一个如图5所示的一个侧面垂直于地面且表面为郎伯面的立方体。这样对于目标象元(x，y)，侧面上的微小的面元ds就可以看成一个点光源，整个侧面对于目标象元(x，y)的作用可以看成所有面元作用的累加，公式如下：

(10)

公式中，是目标像元的方向性反射率，其入射方向为墙壁的二次入射方向，观察方向是传感器方向，目标像元的方向性反射率由光谱库直接查询获得；

是墙壁的方向性反射率；x，x₀，n，y₀，y，m是空间坐标参数，见图5。其它步骤与具体实施方式一相同。

本发明内容不仅限于上述各实施方式的内容，其中一个或几个具体实施方式的组合同样也可以实现发明的目的。

Claims (7)

1、用于林地复杂场景高光谱遥感数据的仿真方法，其特征在于它的步骤如下：

步骤一：设定仿真参数，并输入仿真中用到的光谱库地物反射率参考曲线；

步骤二：将欲仿真得到的高光谱遥感数据分解为太阳辐射直接被目标象元反射的能量、建筑物反射的太阳辐射被目标象元再次反射的能量、天空散射辐射被目标象元反射的能量、程辐射和临近效应五个部分，并按照步骤一中设定的仿真参数分别进行仿真；

步骤三：将步骤二中获得的仿真结果组合得到最终的仿真数据。

2、根据权利要求1所述的用于林地复杂场景高光谱遥感数据的仿真方法，其特征在于步骤一中的仿真参数划分为两组，分别为：地表场景参数和大气场景参数；

地表场景参数为地面场景大小、分辨单元大小、各单元地物分类、树木类型、树木平均高度、树木冠层尺寸、植株密度和建筑物高度；其中树木冠层尺寸简化为一个椭球体，其大小有长轴和短轴确定；

大气场景参数为太阳入射方向的天顶角、太阳入射方向的方向角、传感器观测方向的天顶角、传感器观测方向的方向角、MODTRAN参数、光谱范围和光谱分辨率。

3、根据权利要求2所述的用于林地复杂场景高光谱遥感数据的仿真方法，其特征在于步骤二将高光谱遥感数据划分为太阳辐射直接被目标象元反射的能量

建筑物反射的太阳辐射被目标象元再次反射的能量

天空散射辐射被目标象元反射的能量

程辐射L_λ，sky和临近效应L_{λ，adjacent}五个部分；于是总辐射也就是仿真得到的最终数据传感器处辐照能量L_λ，sensor便可由公式(1)表示：

$L_{\mathrm{\λ},\mathrm{sensor}}=L_{\mathrm{\λ},R}^I+L_{\mathrm{\λ},R}^{\mathrm{II}}+L_{\mathrm{\λ},R}^{\mathrm{III}}+L_{\mathrm{\λ},\mathrm{adjacent}}+L_{\mathrm{\λ},\mathrm{sky}}---\left(1\right).$

4、根据权利要求1所述的用于林地复杂场景高光谱遥感数据的仿真方法，其特征在于步骤二中的高光谱遥感数据仿真中是采用了四流地表-大气辐射传输模型表示地表反射率均为0、0.5、1三种情况下传感器处辐照能量L_λ，sensor与太阳入射到大气顶层的辐照能量E_s ^o的关系，如公式(2)

$L_{\mathrm{\λ},\mathrm{sensor}}={\mathrm{\ρ}}_{\mathrm{so}}\frac{{E_s}^o}{\mathrm{\π}}\cos {\mathrm{\θ}}_s+\frac{{\mathrm{\τ}}_{\mathrm{ss}}{\stackrel{\‾}{r}}_{\mathrm{sd}}+{\mathrm{\τ}}_{\mathrm{sd}}{\stackrel{\‾}{r}}_{\mathrm{dd}}}{1-{\stackrel{\‾}{r}}_{\mathrm{dd}}{\mathrm{\ρ}}_{\mathrm{dd}}}{\mathrm{\τ}}_{\mathrm{do}}\frac{{E_s}^o}{\mathrm{\π}}\cos {\mathrm{\θ}}_s+\frac{{\mathrm{\τ}}_{\mathrm{sd}}+{\mathrm{\τ}}_{\mathrm{ss}}{\stackrel{\‾}{r}}_{\mathrm{sd}}{\mathrm{\ρ}}_{\mathrm{dd}}}{1-{\stackrel{\‾}{r}}_{\mathrm{dd}}{\mathrm{\ρ}}_{\mathrm{dd}}}{r}_{\mathrm{do}}{\mathrm{\τ}}_{\mathrm{oo}}\frac{{E_s}^o}{\mathrm{\π}}\cos {\mathrm{\θ}}_s+{\mathrm{\τ}}_{\mathrm{ss}}{r}_{\mathrm{so}}{\mathrm{\τ}}_{\mathrm{oo}}\frac{{E_s}^o}{\mathrm{\π}}\cos {\mathrm{\θ}}_s---\left(2\right)$

其中等号右边第一项表示程辐射，第二项表示临近效应，第三项表示天空散射辐射被目标象元反射的能量，第四项表示太阳辐射直接被目标象元反射的能量；其中τ_ss，τ_oo，ρ_so，ρ_dd，τ_sd，τ_do6个参数描述大气的状况；r_so，r_do，r_dd，r_sd4个参数描述地表象元的反射率，θ_s为太阳入射方向的天顶角。

5、根据权利要求4所述的用于林地复杂场景高光谱遥感数据的仿真方法，其特征在于大气的状况的6个参数τ_ss，τ_oo，ρ_so，ρ_dd，τ_sd，τ_do，通过公式(2)得到如下公式：

$L_{\mathrm{PATH}}=[{\mathrm{\ρ}}_{\mathrm{so}}+\frac{{\mathrm{\τ}}_{\mathrm{ss}}{\stackrel{\‾}{r}}_{\mathrm{sd}}+{\mathrm{\τ}}_{\mathrm{sd}}{\stackrel{\‾}{r}}_{\mathrm{dd}}}{1-{\stackrel{\‾}{r}}_{\mathrm{dd}}{\mathrm{\ρ}}_{\mathrm{dd}}}{\mathrm{\τ}}_{\mathrm{do}}]\frac{{E_s}^o}{\mathrm{\π}}\cos {\mathrm{\θ}}_s---\left(3\right)$

$L_{\mathrm{GSUN}}={\mathrm{\τ}}_{\mathrm{ss}}{r}_{\mathrm{so}}{\mathrm{\τ}}_{\mathrm{oo}}\frac{{E_s}^o}{\mathrm{\π}}\cos {\mathrm{\θ}}_s---\left(4\right)$

$L_{\mathrm{GTOT}}=\frac{{\mathrm{\τ}}_{\mathrm{ss}}+{\mathrm{\τ}}_{\mathrm{ss}}{\stackrel{\‾}{r}}_{\mathrm{sd}}{\mathrm{\ρ}}_{\mathrm{dd}}}{1-{\stackrel{\‾}{r}}_{\mathrm{dd}}{\mathrm{\ρ}}_{\mathrm{dd}}}{r}_{\mathrm{do}}{\mathrm{\τ}}_{\mathrm{do}}\frac{{E_s}^o}{\mathrm{\π}}\cos {\mathrm{\θ}}_s---\left(5\right)$

其中L_PATH表示程辐射与邻近效应能量的和，L_GSUN表示被目标反射的太阳直接入射的能量，L_GTOT表示被目标反射的天空的散射入射的能量；通过推导得到公式(6)

τ_ss＝exp(-b/cosθ_s)

τ_oo＝exp(-b/cosθ_o)

τ_ssτ_oo＝exp[-b(1/cosθ_s+1/cosθ_o)]＝GSUN₁₀₀×π/(E_s ^ocosθ_s)

ρ_so＝PATH₀×π/(E_s ^ocosθ_s)                     (6)

${\mathrm{\ρ}}_{\mathrm{do}}=\frac{{\mathrm{GTOT}}_{100}-2\×{\mathrm{GTOT}}_{50}}{{\mathrm{GTOT}}_{100}-{\mathrm{GTOT}}_{50}}$

${\mathrm{\τ}}_{\mathrm{sd}}=[\frac{{\mathrm{GTOT}}_{100}(1-{\mathrm{\ρ}}_{\mathrm{dd}})}{{\mathrm{GSUN}}_{100}}-{\mathrm{\ρ}}_{\mathrm{dd}}]{\mathrm{\τ}}_{\mathrm{ss}}$

${\mathrm{\τ}}_{\mathrm{do}}=\frac{{\mathrm{PATH}}_{100}-{\mathrm{PATH}}_0}{{\mathrm{GTOT}}_{100}}{\mathrm{\τ}}_{\mathrm{oo}}$

公式中，b是消光系数；通过MODTRAN仿真获得的PATH₀，PATH₅₀，PATH₁₀₀分别表示地表所有象元反射率均为0，0.5，1三种情况下对应的程辐射与邻近效应能量的和；GSUN₀，GSUN₅₀，GSUN₁₀₀分别表示地表所有象元反射率均为0，0.5，1三种情况下对应的被目标反射的太阳直接入射的能量；GTOT₀，GTOT₅₀，GTOT₁₀₀分别表示地表所有象元反射率均为0，0.5，1三种情况下对应的被目标反射的天空的散射入射的能量。

6、根据权利要求1所述的用于林地复杂场景高光谱遥感数据的仿真方法，其特征在于地表象元的反射的4个参数r_so，r_do，r_dd，r_sd的计算方法如下：

公式中，μ_o为传感器观测方向的天顶角的余弦值，μ_i为散射入射方向的天顶角的余弦值，μ_s为太阳入射方向的天顶角的余弦值，μ_v为散射观察方向的天顶角的余弦值，为传感器观测方向的方向角，

为散射入射方向的方向角，

为太阳入射方向的方向角，

为散射观察方向的方向角，ρ′为双向反射率由光谱库直接提供。

7、根据权利要求1所述的用于林地复杂场景高光谱遥感数据的仿真方法，其特征在于建筑物反射的太阳辐射被目标象元再次反射的能量

是将建筑简化为一个侧面垂直于地面且为郎伯面的立方体，公式如下：

公式中，

是目标像元的方向性反射率，目标像元的方向性反射率由光谱库直接查询获得；是墙壁的方向性反射率；x，x₀，n，y₀，y，m是空间坐标参数。

Images