首页 / 专利库 / 物理 / 重子不对称 / 一种基于光子追踪的大气与三维地表耦合辐射模拟方法

一种基于光子追踪的大气与三维地表耦合辐射模拟方法

阅读:596发布:2023-03-02

专利汇可以提供一种基于光子追踪的大气与三维地表耦合辐射模拟方法专利检索,专利查询,专利分析的服务。并且本 发明 公开了一种基于 光子 追踪的大气与三维地表耦合 辐射 模拟方法,包含以下步骤:(1)读入 光源 、大气场景和地表场景的输入参数,并对大气输入参数进行计算处理;(2)从光源向场景内发射一个光子,根据步骤(1)的计算结果模拟光子行进的距离,确定碰撞点 位置 ;(3)在碰撞点更新光子位置、方向和 能量 信息,收集能量贡献;(4)根据光子位置,计算光子下一步去向,继续追踪光子,直到满足终止条件;(5)发射新光子,重复追踪过程,在达到设定的光子数量之后,得到收敛的大气顶辐 亮度 结果。该方法通过追踪光子方向和位置信息来模拟地表与大气辐射传输耦合,更加贴近实际情况,进而更精确的模拟大气对地表辐射传输 信号 的影响。,下面是一种基于光子追踪的大气与三维地表耦合辐射模拟方法专利的具体信息内容。

1.一种基于光子追踪的大气与三维地表耦合辐射模拟方法,其特征在于,它包含以下步骤:
(1)读入光源、大气场景和地表场景的输入参数,并对大气输入参数进行计算处理;
(2)从光源向场景内发射一个光子,根据步骤(1)的计算结果模拟光子行进的距离,并确定碰撞点位置
(3)在碰撞点更新光子位置、方向和能量信息,并收集光子对传感器的能量贡献;
(4)根据光子位置,计算光子下一步去向,继续追踪光子直到到达终止条件;
(5)从光源发射新光子,重复追踪过程,在达到设定的光子数量之后,得到收敛的大气层顶辐亮度结果。
2.根据权利要求1所述的一种基于光子追踪的大气与三维地表耦合辐射模拟方法,其特征在于:所述模拟方法的步骤(1)为读入光源、大气场景和地表场景的输入参数,并对输入参数进行计算处理:
三维地表场景使用几何面片组合的方法进行描述,使用三维坐标描述每个三形或圆形的大小、位置和朝向信息,同时输入每个面片的高光谱的反射率和透射率信息;使用几何面片的拼接可以描述任意复杂程度的三维地表场景;
由于地球大气的垂直异质性,三维地表场景上方的大气建模成不等间隔的12个高度分层,每一个分层内的大气平均匀并具有相应的光学参数;对于每个高度处的大气,按类型分类为大气分子,气溶胶;分别输入各个模拟波段每个分层处大气分子,气溶胶和云的消光系数和单次散射反照率;对应输入气溶胶和云的双Henyey-Greenstein非对称参数描述其散射相函数;
太阳作为大气与三维地表耦合场景的入射光源,其输入的几何参数包括太阳天顶角和方位角,光源设置在大气层顶部,其大小和下方地表三维场景的边界大小保持相同;输入的光谱参数为入射太阳光能量光谱分布;
对输入的大气光学参数进行预处理;使用瑞利散射生成大气分子散射相函数查找表;
对应不同波段和不同高度分层使用双Henyey-Greenstein参数计算云和气溶胶的散射相函数查找表;
根据大气高程调整每个大气分层实际分层高度;
对每一个波段,每一高度位置处的不同大气组分的吸收系数、单次散射反照率、散射相函数查找表进行等效合成;并分别储存作为后续模拟过程直接使用的计算参数。
3.根据权利要求1所述的一种基于光子追踪的大气与三维地表耦合辐射模拟方法,其特征在于:所述模拟方法的步骤(2)为从光源向场景内发射一个光子,根据步骤(1)的计算结果模拟光子行进的距离,并确定碰撞点位置:
光源面设置在大气层顶部,使用均匀随机数在光源面上随机确定每个光子的起始位置,光子的方向由输入的太阳天顶角和方位角确定,而每个光子的初始能量由输入光谱参数决定:
公式中N是拟追踪的光子个数,S表示模拟场景的底面积,λ为当前模拟的波长
光子进入场景后,将首先进入大气模,按设定的初始方向沿直线传播;光线在场景中行进的距离的计算方法则将根据光子所处位置确定;
后续传播过程中,若光子在大气之中,则按高度分层找到当前位置的消光系数,计算得到自由路径;光子在大气中行进的自由路径S由随机数和所在层的消光系数共同决定:
rand是0到1之间的均匀随机数,βatm是在步骤(1)中计算得到的当前位置处的大气等效消光系数,iz表示光子当前所在大气分层,λ是当前模拟的波长;在行进自由路径S之后,光子与大气组分发生碰撞;
若光子在地表场景之中,则模拟光子按当前方向传播直到与几何面片发生碰撞,根据三角形或圆片的三维坐标进行几何求交,计算得到光线碰撞的三角形或圆片编号以及碰撞点位置。
4.根据权利要求1所述的一种基于光子追踪的大气与三维地表耦合辐射模拟方法,其特征在于,其中步骤(3)在碰撞点更新光子位置、方向和能量信息,并收集光子对传感器的能量贡献:光子携带的能量将在碰撞点对应衰减:
Qout(λ)=Qin(λ)·ω(λ)
其中Qout和Qin分别代表了碰撞后和碰撞前的光子能量,ω是单次散射反照率,在大气中是当前位置的等效大气反照率,在地表场景中则为当前碰撞几何面片的反射率和透射率之和;
光子在大气中碰撞之后的新的散射方向分布满足散射相函数的分布,使用步骤(1)中计算得到的当前位置处的散射相函数查找表抽取并计算新的传输方向;光子在地表场景面片表面发生的散射假设其满足朗伯定律,根据当前面片法线朝向信息计算新的散射方向;
在每次碰撞发生的位置,使用光子扩散方法计算累计碰撞事件对设置的观测传感器的能量贡献;如果从碰撞点位置到达传感器的路径上没有受到遮挡,则根据光子入射方向和碰撞点的位置信息,计算光线向传感器方向散射的概率,并按比例收集能量:
Qcollect(Ω,λ)=Qout(λ)·P(Ω,λ)·exp(-La·βatm)
P为光线向传感器方向散射的概率,La表示从碰撞点到达传感器位置处经过的大气路径长度,公式中最后指数项描述这段距离上的上的大气消光。
5.根据权利要求1所述的一种基于光子追踪的大气与三维地表耦合辐射模拟方法,其特征在于:所述模拟方法的步骤(4)为根据光子位置,计算光子下一步去向,继续追踪光子,直到到达终止条件:
光子将按新的方向离开碰撞点,继续沿直线传播并重复碰撞过程;在整个追踪过程中,光子携带位置、方向与能量信息,使用光子的位置信息判定当前处在大气还是地表场景之中,当光子位于大气中时,后续传播过程使用大气相函数查找表随机抽取光子新的方向,使用大气中的路径长度计算方法和碰撞处理方法继续追踪光子,更新光子携带的信息并收集传感器贡献;当光子处于地表场景之中时,根据当前发生碰撞的面片的法线方向和朗伯定律随机抽取光子新的方向,使用地表场景的路径长度计算方法和碰撞处理方法继续追踪光子,更新光子携带的信息并收集传感器贡献;
当经过一定次数的碰撞之后,光子上携带的能量将变得比较小,一旦光子能量小于阈值,使用“俄罗斯轮盘”来能量守恒且无偏的终止该光子的追踪过程;
光子在从大气底部进入地表场景时具有具体的方向和位置信息,可以很好的模拟出经过大气散射之后的天空漫射光的非各向同性;
一旦光子在地表场景中通过碰撞散射离开地表场景上边界,使用光子的高度信息和具体传播方向,可以令光子再次回到大气中继续追踪,很好的再现出在大气与地表场景交界处的光子在两个部分之间的循环往复作用,做到大气和大气下的地表场景的完全耦合。
6.根据权利要求1所述的一种基于光子追踪的大气与三维地表耦合辐射模拟方法,其特征在于:所述模拟方法的步骤(5)为发射新光子,重复追踪过程,在达到设定的光子数量之后,得到收敛的模拟结果:
使用蒙特卡洛思想,在一个光子的模拟过程中,把每一个不确定的物理过程按其物理分布描述成概率过程,并使用随机数抽样的方法确定在本次模拟中发生的具体事件;每一个使用随机方法描述的位置都严格遵循其物理分布规律;
通过大数量的重复模拟过程,即使用一定数目的模拟光子数,使每一个使用随机事件描述的位置按其分布得到充分的采样,最后即可通过累计的每一次对观测传感器的能量贡献得到最终整体三维地表场景与大气耦合结果,即大气层顶传感器观测方向接收到的辐亮度信号
上式中各个符号的含义同前文。

说明书全文

一种基于光子追踪的大气与三维地表耦合辐射模拟方法

技术领域

[0001] 本发明涉及定量遥感模拟领域,尤其涉及一种基于光子追踪的大气与三维地表耦合辐射模拟方法。在定量分析大气耦合对遥感信号的影响以及地表植被生态学等研究方面具有重要意义。

背景技术

[0002] 遥感测量是通过传感器远程接收来自地物的电磁辐射信息来探究地物状态的重要信息获取工具,获取地表物体的方向性反射率等特征信息能够得到丰富的物体结构与类别信息。建立三维地表场景模型可以通过修改输入的场景状态参数直接模拟计算场景的反射光谱,对定量化的研究地表场景三维结构对遥感信号的影响具有重要意义。
[0003] 而对于遥感观测来说,地球大气对信号的影响是一个无法忽略的关键问题。地表场景在特定观测度上的反射率特性,一定程度上取决于当前的大气条件。现在主要采取的方法是使用大气顶部的遥感信号经过大气校正,得到地表场景顶部高度处的光谱信号。但这需要对大气状况做各种近似假设,而这些关于大气状况的假设在非典型情况下会带来很大的误差;另一方面,现有的大气耦合模型多采取将地表场景和大气的贡献累计在大气底部高度,再通过计算有限次相互作用的方法模拟地表场景与大气的耦合,无法还原真实的物理过程,且大多数情况下忽略了大气临近效应,得到的结果精度有限。
[0004] 因此,定量研究遥感信号与大气性质和三维地表场景之间的关系,最方便、最有发展潜的方法是建立一个能真实还原三维地表场景与大气辐射传输耦合的模型,这对更好理解分析三维地表场景的辐射传输过程以及实际的传感器指标设计研制等方面具有重要的意义及应用价值。
[0005] 需要注意的是,本部分旨在为权利要求书中陈述的本发明提供背景或上下文。此处的描述不因为包括在本部分中就承认是现有技术

发明内容

[0006] 本发明的目的在于针对目前考虑大气的三维地表场景辐射传输模拟过程中近似较多,模拟结果不准确的问题,提供一种基于光子追踪的大气与三维地表耦合辐射模拟方法。
[0007] 本发明的技术解决方案为:通过蒙特卡洛思想,参数化描述大气部分和地表部分。在模拟过程中跟踪每一个光子方向、位置能量信息,并使用其位置信息,判断光子目前处在大气中或三维地面场景之中,对应使用大气的光学参数或地表场景的光学和结构参数计算光子的行进距离与方向。在每个光子与场景中元素发生碰撞的位置收集并累计在传感器观测方向上的能量贡献。通过一定光子数量的模拟过程,得到收敛的大气层顶辐亮度。其具体步骤如下:
[0008] (1)读入光源、大气场景和地表场景的输入参数,并对大气输入参数进行计算处理;
[0009] (2)从光源向场景内发射一个光子,根据步骤(1)的计算结果模拟光子行进的距离,并确定碰撞点位置;
[0010] (3)在碰撞点更新光子位置、方向和能量信息,并收集光子对传感器的能量贡献;
[0011] (4)根据光子位置,计算光子下一步去向,继续追踪光子,直到到达终止条件;
[0012] (5)从光源发射新光子,重复追踪过程,在达到设定的光子数量之后,得到收敛的大气层顶辐亮度结果。
[0013] 其中,步骤(1)读入光源和地表场景的三维结构参数、光谱参数以及大气的光谱参数,并对大气输入参数进行计算处理:
[0014] 三维地表场景使用几何面片组合的方法进行描述,使用三维坐标描述每个三角形或圆形的大小、位置和朝向信息,同时输入每个面片的高光谱的反射率和透射率信息;使用几何面片的拼接可以描述任意复杂程度的三维地表场景。
[0015] 由于地球大气的垂直异质性,三维地表场景上方的大气建模成不等间隔的12个高度分层,每一个分层内的大气平均匀并具有相应的光学参数;对于每个高度处的大气,按类型分类为大气分子,气溶胶;分别输入各个模拟波段每个分层处大气分子,气溶胶和云的消光系数和单次散射反照率;对应输入气溶胶和云的双Henyey-Greenstein非对称因子参数描述其散射相函数。
[0016] 太阳作为大气与三维地表耦合场景的入射光源,其输入的几何参数包括太阳天顶角和方位角,光源设置在大气层顶部,其大小和下方地表三维场景的边界大小保持相同。输入的光谱参数为入射太阳光能量光谱分布。
[0017] 对输入的大气光学参数进行预处理,用瑞利散射生成大气分子散射相函数查找表;对应不同波段和不同高度分层使用双Henyey-Greenstein参数计算云和气溶胶的散射相函数查找表;根据大气高程调整每个大气分层实际高度;对每一个波段,每一高度位置处的不同大气组分的吸收系数、单次散射反照率、散射相函数查找表进行等效合成,并将结果分别储存作为后续模拟过程直接使用的计算参数。
[0018] 其中,步骤(2)从光源向场景内发射一个光子,根据步骤(1)的计算结果模拟光子行进的距离,并确定碰撞点位置:根据设定的光源参数,确定光子从大气层顶进入场景的起始位置、方向和能量。
[0019] 光源面设置在大气层顶部,使用均匀随机数在光源面上随机确定每个光子的起始位置,光子的方向由输入的太阳天顶角和方位角确定,而每个光子的初始能量由输入光谱参数决定:
[0020]
[0021] 公式中N是拟追踪的光子个数,S表示模拟场景的底面积,λ为当前模拟的波长;光子进入场景后,将首先进入大气模,按设定的初始方向沿直线传播。光线在场景中行进的距离的计算方法则将根据光子所处位置确定。
[0022] 根据光子所处位置,确定其处在大气内部或是三维地表场景内部。若光子在大气之中,则按高度分层找到当前位置的大气等效消光系数,计算得到自由路径。光子在大气中行进的自由路径S由随机数和所在层的等效消光系数共同决定:
[0023]
[0024] rand是0到1之间的均匀随机数,βatm是在步骤(1)中计算得到的当前位置处的大气等效消光系数,iz表示光子当前所在大气分层,λ是当前模拟的波长;在行进自由路径S之后,光子与大气组分发生碰撞。
[0025] 若光子在地表场景之中,则模拟光子按当前方向传播直到与几何面片发生碰撞,根据三角形或圆片的三维坐标进行几何求交,计算得到光线碰撞的三角形或圆片编号以及碰撞点位置。
[0026] 其中,步骤(3)在碰撞点更新光子位置、方向和能量信息,并收集光子对传感器的能量贡献:在碰撞点,光子携带的能量将在碰撞点对应衰减:
[0027] Qout(λ)=Qin(λ)·ω(λ)
[0028] 其中Qout和Qin分别代表了碰撞后和碰撞前的光子能量,ω是单次散射反照率,在大气中是当前位置的等效大气反照率,在地表场景中则为当前碰撞三角形或圆片的反射率和透射率之和。
[0029] 光子在大气中碰撞之后的新的散射方向分布满足散射相函数的分布,使用步骤(1)中计算得到的当前位置处的散射相函数查找表抽取并计算新的传输方向;光子在地表场景面片表面发生的散射假设其满足朗伯定律,根据面片朝向信息计算新的散射方向。
[0030] 在每次碰撞发生的位置,使用光子扩散方法计算累计碰撞事件对设置的观测传感器的能量贡献;如果从碰撞点位置到达传感器的路径上没有受到遮挡,则根据光子入射方向和碰撞点的位置信息,计算光线向传感器方向散射的概率,并按比例收集能量:
[0031] Qcollect(Ω,λ)=Qout(λ)·P(Ω,λ)·exp(-La·βatm)
[0032] P为光线向传感器方向散射的概率,La表示从碰撞点到达传感器位置处经过的大气路径长度,公式中最后指数项描述这段距离上的大气消光。
[0033] 其中,步骤(4)根据光子位置,计算光子下一步去向,继续追踪光子,直到到达终止条件:光子将按新的方向离开碰撞点,继续沿直线传播并重复碰撞过程。在整个追踪过程中,光子携带位置、方向与能量信息,使用光子的位置信息判定当前处在大气还是地表场景之中,当光子位于大气中时,后续传播过程使用大气相函数查找表随机抽取光子新的方向,使用大气中的路径长度计算方法和碰撞处理方法继续追踪光子,更新光子携带的信息并收集传感器贡献;当光子处于地表场景之中时,根据当前发生碰撞的面片的法线方向和朗伯定律随机抽取光子新的方向,使用地表场景的路径长度计算方法和碰撞处理方法继续追踪光子,更新光子携带的信息并收集传感器贡献。
[0034] 当经过一定次数的碰撞之后,光子上携带的能量将变得比较小,一旦光子能量小于阈值,使用“俄罗斯轮盘”来能量守恒且无偏的终止该光子的追踪过程。
[0035] 光子在从大气底部进入地表场景时具有具体的方向和位置信息,可以很好的模拟出经过大气散射之后的天空漫射光的非各向同性。一旦光子在地表场景中通过碰撞散射离开地表场景上边界,使用光子的高度信息和具体传播方向,可以令光子返回大气中继续追踪,很好的再现出在大气与地表场景交界处的光子在两个部分之间的循环往复作用,做到大气和大气下的地表场景的完全耦合。
[0036] 其中,步骤(5)发射新光子,重复追踪过程,在达到设定的光子数量之后,得到收敛的大气层顶辐亮度结果:使用蒙特卡洛思想,在一个光子的模拟过程中,把每一个不确定的物理过程按其物理分布描述成概率过程,并使用随机数抽样的方法确定在本次模拟中发生的具体事件,每一个使用随机方法描述的位置都严格遵循其物理分布规律。
[0037] 通过大数量的重复模拟(即使用一定数目的模拟光子数),使每一个使用随机事件描述的位置按其分布得到充分的采样,最后即可通过累计的每一次对观测传感器的能量贡献得到最终大气层顶传感器观测方向接收到的辐亮度信号:
[0038]
[0039] 上式中各个符号的含义同前文。
[0040] 本发明与现有技术相比的优点在于:基于物理原理还原了光子在包含有大气的场景中的传输过程,真实再现出光子在大气与地表场景交界处的循环往复作用过程,做到大气和地表场景的辐射过程耦合,减少了描述该过程中的假设与近似,从而减小了大气层顶辐亮度的模拟误差。
[0041] 它具有以下的优点:(1)通过计算机模拟,在有限的计算时间下,可以计算得到不同高度处、任意观测角度、可选波长范围、不同大气和地表参数设置条件下的方向反射率因子等遥感感兴趣量。功能多样,可修改的参数丰富且具有代表性。(2)利用三维坐标几何面片的方式可以描述任意复杂程度的三维地物场景,并且三维场景的构建与辐射模拟过程互相独立,前者可引入计算机图形学的最新发展成果,具有广阔的应用场景。(3)使用该方法建立的模型为计算机模拟模型,相对于另外常见的遥感解析模型和几何光学模型的建模方法精度更高,同时在模拟过程中引入了多种加速算法,在计算时间上也得到了有效控制。随着计算机性能的提高,计算效率将随之得到进一步的提升。附图说明
[0042] 图1为本发明一种基于光子追踪的大气与三维地表耦合辐射模拟方法的流程图

具体实施方式

[0043] 为了更好的说明本发明涉及的基于光子追踪的三维地表与大气辐射耦合模拟方法,利用可见光到短波红外波长范围的高光谱三维植被冠层数据和大气数据进行了辐射传输模拟,计算得到了垂直观测方向的大气层顶反射率。具体实现步骤如下:
[0044] (1)读入光源、大气场景和地表场景的输入参数,并对大气输入参数进行计算处理:三维地表场景的输入参数包括其结构参数和光谱参数,读入处于拔节期的玉米冠层场景(使用不同大小的三角形面片拼接描述)的结构参数文件和反射率和透射率光谱文件。读入描述土壤的土壤结构参数文件和反射光谱文件。
[0045] 读入对应波长范围的大气光谱参数,使用标准中纬度夏季大气类型作为大气输入参数,按高度区间划分为下端较为稠密,上端较为稀疏的12个分层高度分层,波段区间可以是为400-2500纳米,但不限于此;
[0046] 太阳作为大气与三维地表耦合场景的入射光源,入射太阳天顶角设置为30度,方位角设置为0度。光源设置在大气层顶部,其大小和下方地表三维场景的边界大小保持相同。输入的光谱参数为入射太阳光能量光谱分布。
[0047] 对输入的大气光学参数进行预处理,用瑞利散射生成大气分子散射相函数;对应不同波段和不同高度分层使用双Henyey-Greenstein参数计算云和气溶胶的散射相函数查找表;根据大气高程调整每个大气分层实际高度;对每一个波段,每一高度位置处的不同大气组分的吸收系数、单次散射反照率、散射相函数查找表进行等效合成;并分别储存作为后续模拟过程直接使用的计算参数。
[0048] (2)从光源向场景内发射一个光子,根据步骤(1)的计算结果模拟光子行进的距离,并确定碰撞点位置:光源面设置在大气层顶部,使用均匀随机数在光源面上随机确定一个随机位置作为光子起点,按30度天顶角、0度方位角方向进入场景。每个光子的初始能量由输入光谱参数决定:
[0049]
[0050] 公式中N是拟追踪的光子个数,S表示模拟场景的底面积,λ为当前模拟的波长;光子进入场景后,将首先进入大气模块,按设定的初始方向沿直线传播。光线在场景中行进的距离的计算方法则将根据光子所处位置确定。
[0051] 根据光子所处位置,确定其处在大气内部或是玉米冠层内部。若光子在大气之中,则按高度分层找到当前位置的大气等效消光系数,计算得到自由路径。
[0052] 光子在大气中行进的自由路径S由随机数和所在层的消光系数共同决定:
[0053]
[0054] rand是0到1之间的均匀随机数,βatm是在步骤(1)中计算得到的当前位置处的大气等效消光系数,iz表示光子当前所在大气分层,λ是当前模拟的波长;在行进自由路径S之后,光子与大气组分发生碰撞。
[0055] 若光子在地表场景之中,则模拟光子按当前方向传播直到与几何面片发生碰撞,根据三角形或圆片的三维坐标进行几何求交,计算得到光线碰撞的三角形或圆片编号以及碰撞点位置。
[0056] (3)在碰撞点更新光子位置、方向和能量信息,并收集光子对传感器的能量贡献:在碰撞点,光子携带的能量将在碰撞点对应衰减:
[0057] Qout(λ)=Qin(λ)·ω(λ)
[0058] 其中Qout和Qin分别代表了碰撞后和碰撞前的光子能量,ω是单次散射反照率,在大气中是当前位置的等效大气反照率,在地表场景中则为当前碰撞三角形或圆片的反射率和透射率之和。
[0059] 光子在大气中碰撞之后的新的散射方向分布满足散射相函数的分布,使用步骤(1)中计算得到的当前位置处的散射相函数查找表抽取并计算新的传输方向;光子在地表场景面片表面发生的散射假设其满足朗伯定律,根据面片朝向信息计算新的散射方向。
[0060] 在每次碰撞发生的位置,使用光子扩散方法计算累计碰撞事件对设置的观测传感器的能量贡献;如果从碰撞点位置到达传感器的路径上没有受到遮挡,则根据光子入射方向和碰撞点的位置信息,计算光线向传感器方向散射的概率,并按比例收集能量:
[0061] Qcollect(Ω,λ)=Qout(λ)·P(Ω,λ)·exp(-La·βatm)
[0062] P为光线向传感器方向散射的概率,La表示从碰撞点到达传感器位置处经过的大气路径长度,公式中最后指数项描述这段距离上的上的大气消光。
[0063] (4)根据光子位置,计算光子下一步去向,继续追踪光子,直到到达终止条件:光子将按新的方向离开碰撞点,继续沿直线传播并重复碰撞过程;在整个追踪过程中,光子携带位置、方向与能量信息,使用光子的位置信息判定当前处在大气还是地表场景之中,当光子位于大气中时,后续传播过程使用大气相函数查找表随机抽取光子新的方向,使用大气中的路径长度计算方法和碰撞处理方法继续追踪光子,更新光子携带的信息并收集传感器贡献;当光子处于地表场景之中时,根据当前发生碰撞的面片的法线方向和朗伯定律随机抽取光子新的方向,使用地表场景的路径长度计算方法和碰撞处理方法继续追踪光子,更新光子携带的信息并收集传感器贡献。
[0064] 当经过一定次数的碰撞之后,光子上携带的能量将变得比较小,一旦光子能量小于阈值,使用“俄罗斯轮盘”来能量守恒且无偏的终止该光子的追踪过程。
[0065] 光子在从大气底部进入地表场景时具有具体的方向和位置信息,可以很好的模拟出经过大气散射之后的天空漫射光的非各向同性。一旦光子在地表场景中通过碰撞散射离开地表场景上边界,使用光子的高度信息和具体传播方向,可以令光子返回大气中继续追踪,很好的再现出在大气与地表场景交界处的光子在两个部分之间的循环往复作用,做到大气和大气下的地表场景的完全耦合。
[0066] (5)发射新光子,重复追踪过程,在达到设定的光子数量之后,得到收敛的大气层顶辐亮度结果:使用蒙特卡洛思想,在一个光子的模拟过程中,把每一个不确定的物理过程(光子行进的距离,碰撞后散射的方向)按其物理分布描述成概率过程,并使用随机数抽样的方法确定在本次模拟中发生的具体事件,每一个使用随机方法描述的位置都严格遵循其物理分布规律。
[0067] 通过N=10万光子数的模拟,使每一个使用随机事件描述的位置按其分布得到充分的采样,最后即可通过累计的每一次对观测传感器的能量贡献得到最终整体三维地表场景与大气耦合结果,即大气层顶传感器观测方向接收到的辐亮度信号:
[0068]
[0069] 上式中各个符号的含义同前文。在未进行多线程优化的情况下,在个人笔记本电脑(CPU Intel Core i7-7500U@2.70GHz)的上计算用时约为5.2分钟。
高效检索全球专利

专利汇是专利免费检索,专利查询,专利分析-国家发明专利查询检索分析平台,是提供专利分析,专利查询,专利检索等数据服务功能的知识产权数据服务商。

我们的产品包含105个国家的1.26亿组数据,免费查、免费专利分析。

申请试用

分析报告

专利汇分析报告产品可以对行业情报数据进行梳理分析,涉及维度包括行业专利基本状况分析、地域分析、技术分析、发明人分析、申请人分析、专利权人分析、失效分析、核心专利分析、法律分析、研发重点分析、企业专利处境分析、技术处境分析、专利寿命分析、企业定位分析、引证分析等超过60个分析角度,系统通过AI智能系统对图表进行解读,只需1分钟,一键生成行业专利分析报告。

申请试用

QQ群二维码
意见反馈