一种基于球谐离散坐标法的三维含云大气辐射计算方法技术

本发明专利技术公开了一种基于球谐离散坐标法的三维含云大气辐射计算方法，包括如下步骤：首先利用爱丁顿近似法计算三维含云大气近似辐射亮度；然后基于三维邻云效应对三维含云大气近似辐射亮度进行处理，将三维邻云效应分为高阶散射增强、低阶散射增强、高阶散射减弱和低阶散射减弱四部分，并分别计算各部分辐射亮度；接着将三维含云大气近似辐射亮度与四部分辐射亮度结合得到三维含云大气初始辐射亮度，并利用其计算多次散射源函数项，进而结合单次散射源函数项与热源函数项，得到辐射源函数；最后基于球谐离散坐标法，将得到的辐射源函数整体作为初始条件，进行迭代计算，得到三维含云大气的辐射亮度。

A three-dimensional cloud atmosphere radiation calculation method based on spherical harmonic discrete coordinate method

【技术实现步骤摘要】
一种基于球谐离散坐标法的三维含云大气辐射计算方法
本专利技术涉及大气辐射计算领域，尤其是一种基于球谐离散坐标法的三维含云大气辐射计算方法。
技术介绍
云是地球大气中分布最为广泛的自然现象之一，在大气系统的辐射传输过程中扮演着举足轻重的角色，直接影响到地球的气候变迁、天气变化以及辐射收支平衡。由于微观结构与宏观形态的随机性和复杂性，云层辐射问题成为大气辐射传输领域中的世界性难题。同时，无论在水平方向还是垂直方向上，云层都具有显著的非均匀性结构特征，所以对于高分辨率遥感应用来讲，需要充分考虑到云的三维辐射特性。当前，计算成本是制约着三维大气辐射传输模型广泛应用的关键问题，而多次散射的计算是影响着计算成本的主要因素。一般认为，当大气光学厚度大于0.1时，便无法忽略多次散射效应带来的影响。由于多次散射源函数项的计算式内隐含了辐射亮度，所以需要进行多次迭代求解，这直接导致了三维辐射传输问题的计算成本升高。于是对多次散射源函数项的计算过程进行合理简化以降低计算成本，成为了亟待满足的需求。
技术实现思路
本专利技术所要解决的技术问题在于，提供一种基于球谐离散坐标法的三维含云大气辐射计算方法，能够提高初始辐射亮度的精确程度，从而在一定程度上减少计算过程中的迭代次数，降低计算成本。为解决上述技术问题，本专利技术提供一种基于球谐离散坐标法的三维含云大气辐射计算方法，包括如下步骤：(1)利用爱丁顿近似法计算三维含云大气近似辐射亮度；(2)将三维邻云效应分为高阶散射增强、低阶散射增强、高阶散射减弱和低阶散射减弱四部分，并分别计算各部分辐射亮度；(3)将三维含云大气近似辐射亮度与四部分辐射亮度结合得到三维含云大气初始辐射亮度，并利用其计算多次散射源函数项，进而结合单次散射源函数项与热源函数项，得到辐射源函数；(4)基于球谐离散坐标法，迭代计算三维含云大气的辐射亮度，整理并输出亮度图像。优选的，步骤(2)中，将三维邻云效应分为四部分，具体步骤如下：(1)根据大气属性文件，将三维含云场景均分为若干尺度相同的立体网格；(2)将三维邻云效应分解为高阶散射增强(EMS)、低阶散射增强(ELS)、高阶散射减弱(RMS)和低阶散射减弱(RLS)四部分：ICE＝IEMS+IELS-IRMS-IRLS其中ICE代表某一网格点处由三维邻云效应所产生的辐射亮度。优选的，步骤(2)中，计算三维邻云效应的各部分辐射亮度，具体步骤如下：(1)在三维含云场景内选取某网格点A，再选取与其水平相邻的网格点B；(2)计算网格点A与B之间的四种辐射效应系数CEMS、CELS、CRMS和CRLS；(3)计算网格点B对A产生的四种辐射亮度IEMS、IELS、IRMS和IRLS；(4)计算网格点B处由于三维邻云效应所产生的辐射亮度ICE；(5)重复步骤1-4，直到选取了三维含云场景内最后一个网格点作为A点，且与其相对应的所有B点的ICE都计算完成为止。优选的，步骤(3)中，将三维含云大气近似辐射亮度与四部分辐射亮度结合得到三维含云大气初始辐射亮度，并利用其计算多次散射源函数项，具体步骤如下：(1)在三维含云场景内选取某网格点A，再选取与其水平方向相邻的网格点B；(2)累加所有B点的ICE，令其结果为I′CE；(3)计算网格点A处的三维含云大气初始辐射亮度：IA＝IEDD+I′CE其中IEDD为通过爱丁顿近似法计算得到的三维含云大气近似辐射亮度；(4)利用三维含云大气初始辐射亮度计算网格点A处的多次散射源函数项；(5)重复步骤1-4，直到选取了三维含云场景内最后一个网格点作为A点，且该处的多次散射源函数项计算完成为止。本专利技术的有益效果为：本专利技术可提高初始辐射亮度的精确程度，从而在一定程度上减少计算过程中的迭代次数，降低计算成本。附图说明图1为本专利技术的方法流程示意图。具体实施方式现将结合附图对本专利技术的技术方案进行完整的描述。以下描述仅仅是本专利技术的一部分实施案例而已，并非全部。基于本专利技术中的实施案例，本领域技术人员在没有作出创造性劳动的前提下所获得的所有其他实施案例，都属于本专利技术的权利保护范围之内。如图1所示，一种基于球谐离散坐标法的三维含云大气辐射计算方法，包括如下步骤：(1)利用爱丁顿近似法计算三维含云大气近似辐射亮度；(2)将三维邻云效应分为高阶散射增强、低阶散射增强、高阶散射减弱和低阶散射减弱四部分，并分别计算各部分辐射亮度；(3)将三维含云大气近似辐射亮度与四部分辐射亮度结合得到三维含云大气初始辐射亮度，并利用其计算多次散射源函数项，进而结合单次散射源函数项与热源函数项，得到辐射源函数；(4)基于球谐离散坐标法，迭代计算三维含云大气的辐射亮度，整理并输出亮度图像。1.三维邻云效应基于球谐离散坐标法，在通过爱丁顿近似计算得到三维含云大气近似辐射亮度之后，再基于三维邻云效应对其进行处理，使其精确程度提高，进而减少迭代计算的次数。基于米氏散射等原理，预先利用实际观测的云层含水量数据计算三维含云大气的反照率、消光系数和相函数，且上述三类参量的空间分布情况与云层含水量数据保持一致，并将反照率、消光系数和相函数以及它们的三维空间坐标储存在大气属性文件内。首先根据大气属性文件，将三维含云场景均分为若干尺度相同的立体网格，于是在每个网格点处都储存着与其相对应的辐射亮度，然后根据三维邻云效应可能对辐射亮度值造成的影响，将其分为高阶散射增强、低阶散射增强、高阶散射减弱和低阶散射减弱四部分(如图1)：ICE＝IEMS+IELS-IRMS-IRLS其中ICE代表某一网格点处由三维邻云效应所产生的辐射亮度。接下来，开始在三维含云场景内选取某网格点A，再选取与其水平相邻的网格点B。2四部分辐射亮度2.1辐射效应系数计算在确定了网格点A与B之后，计算A与B之间的四种辐射效应系数CEMS、CELS、CRMS和CRLS：其中H为网格点A和B的间距；kA和kB分别表示A和B处的消光系数；ωA和ωB分别表示A和B处的反照率；R表示光束在A和B的交界处发生反射的比例，其计算方法类似于一种大气反照率的计算：其中P为网格点B处的散射相函数，其以勒让德级数内子项系数的形式储存在大气属性文件中；μ为散射相函数内散射角的余弦。2.2辐射亮度计算在辐射效应系数计算完成后，再利用其计算网格点B对A产生的四种辐射亮度IEMS、IELS、IRMS和IRLS：其中IA→B表示从A射向B的辐射亮度；IB→A表示从B射向A的辐射亮度；IB表示B的上行方向辐射亮度。于是，便可以计算得到网格点B处由于三维邻云效应所产生的辐射亮度ICE：ICE＝IEMS+IELS-IRMS-IRLS接着选取新的A点，重复2.1与2.2本文档来自技高网...

【技术保护点】
1.一种基于球谐离散坐标法的三维含云大气辐射计算方法，其特征在于，包括如下步骤：/n1)利用爱丁顿近似法计算三维含云大气近似辐射亮度；/n2)将三维邻云效应分为高阶散射增强、低阶散射增强、高阶散射减弱和低阶散射减弱四部分，并分别计算各部分辐射亮度；/n3)将三维含云大气近似辐射亮度与四部分辐射亮度结合得到三维含云大气初始辐射亮度，并利用初始辐射亮度计算多次散射源函数项，进而结合单次散射源函数项与热源函数项，得到辐射源函数；/n4)基于球谐离散坐标法，迭代计算三维含云大气的辐射亮度，并输出亮度图像。/n

【技术特征摘要】
1.一种基于球谐离散坐标法的三维含云大气辐射计算方法，其特征在于，包括如下步骤：
1)利用爱丁顿近似法计算三维含云大气近似辐射亮度；
2)将三维邻云效应分为高阶散射增强、低阶散射增强、高阶散射减弱和低阶散射减弱四部分，并分别计算各部分辐射亮度；
3)将三维含云大气近似辐射亮度与四部分辐射亮度结合得到三维含云大气初始辐射亮度，并利用初始辐射亮度计算多次散射源函数项，进而结合单次散射源函数项与热源函数项，得到辐射源函数；
4)基于球谐离散坐标法，迭代计算三维含云大气的辐射亮度，并输出亮度图像。


2.如权利要求1所述的基于球谐离散坐标法的三维含云大气辐射计算方法，其特征在于，步骤2)中，将三维邻云效应分为四部分，具体为：
步骤2.1，根据大气属性文件，将三维含云场景均分为若干尺度相同的立体网格；
步骤2.2，将三维邻云效应分解为高阶散射增强EMS、低阶散射增强ELS、高阶散射减弱RMS和低阶散射减弱RLS四部分：
ICE＝IEMS+IELS-IRMS-IRLS
其中ICE代表某一网格点处由三维邻云效应所产生的辐射亮度，IEMS为高阶散射增强所产生的辐射亮度，IELS为低阶散射增强所产生的辐射亮度，IRMS为高阶散射减弱所产生的辐射亮度，IRLS为低阶散射减弱所产生的辐射亮度。


3.如权利要求1所述的基于球谐离散坐标法的三维含云大气辐射计算方法，其特征在于，步骤2)中，计算三维邻...

【专利技术属性】
技术研发人员：张寅，马俊，闫钧华，
申请(专利权)人：南京航空航天大学，
类型：发明
国别省市：江苏;32