CN107563012A - 基于曲线拟合的大气参量查找表生成方法 - Google Patents
基于曲线拟合的大气参量查找表生成方法 Download PDFInfo
- Publication number
- CN107563012A CN107563012A CN201710678066.XA CN201710678066A CN107563012A CN 107563012 A CN107563012 A CN 107563012A CN 201710678066 A CN201710678066 A CN 201710678066A CN 107563012 A CN107563012 A CN 107563012A
- Authority
- CN
- China
- Prior art keywords
- mrow
- path length
- msub
- mtd
- sight
- 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
Links
Landscapes
- Measurement Of Radiation (AREA)
- Testing Or Calibration Of Command Recording Devices (AREA)
Abstract
本发明公开了一种基于曲线拟合的大气参量查找表生成方法，主要用来生成光电三维场景中的大气参量查找表。其方法实现过程包括：设置曲线拟合仿真参数，划分探测器海拔高度区间，计算探测器海拔高度值，计算路径长度最大有效值，划分路径长度子区间，计算路径长度区间的路径长度值，计算大气参量，计算拟合曲线参数，生成大气参量查找表，即大气透过率查找表和大气路径辐射查找表。本发明具有保持大气参量的计算精度，减少了大气参量的计算量，缩短了大气参量查找表的生成时间，提高了大气参量查找表的生成效率的优点。
Description
技术领域
本发明属于物理技术领域，更进一步涉及物理技术领域中的一种基于曲线拟合的大气参量查找表生成方法。本发明可生成光电三维场景中的大气参量查找表，用于大规模光电三维场景仿真的实时渲染。
背景技术
随着计算机技术和仿真技术的飞速发展，大规模光电三维场景的实时仿真一直是研究的热点。光电三维场景的实时仿真包括对大气传输效应的实时仿真，大气传输效应的实时仿真需要根据仿真帧频和分辨率计算仿真图像中的每个像素所对应的辐射传输几何路径下的包括大气透过率和大气路径辐射的大气参量，计算量庞大。针对上述计算量庞大的问题，通过对大气透过率和大气路径辐射进行预计算，生成大气参量查找表是一种有效的技术方法。
中国航天科工集团第二研究院二Ο七所在其申请的专利文献“红外特性测量中大气传输快速修正方法”(专利申请号201210362216.3，申请号公布号CN103674904A)中公开了一种红外特性测量中大气传输快速修正方法。该方法的主要步骤为：第一，设置计算参数；第二，选择气象模式；第三，根据所选的气象模式调用MODTRAN软件中相应的计算模型来进行大气传输参数计算；第四，计算距离测量位置不同观测高度和水平投影距离位置点的大气透过率和背景辐射亮度；第五，调用计算出的大气透过率和背景辐射亮度网格数据表格，对数据进行查分，得出不同观测仰角和观测距离点的大气透过率和背景辐射亮度变化的曲线。该方法存在的不足之处是：该方法辐射传输几何参数计算间隔太小，计算的数据量太大，计算耗时。
孙明明在其发表论文“红外场景的大气效应高效计算与仿真研究”(西安电子科技大学硕士论文2014年)中提出了一种计算大气参量的方法。该方法的主要步骤为：第一，设置气象参数；第二，MODTRAN预计算所选辐射传输几何参数下的大气透过率和大气路径辐射；第三，选取拟合函数形式；第四，采用麦夸特算法拟合大气透过率和大气路径辐射随路径长度变化的曲线；第五，对路径长度进行密集取值，生成仿真场景的大气参量查找表。该方法存在的不足之处是：大气参量的计算精度低，计算耗时，计算量大，计算场景中所有辐射传输几何路径的大气参量时，需遍历所有探测器海拔高度值。
发明内容
本发明的目的在于克服上述现有技术存在的不足，提出一种基于曲线拟合的大气参量查找表生成方法，本发明可以显著减小生成大规模光电三维场景的大气参量查找表的计算量，提高生成大气参量查找表的速度，计算精度高，方法简单实用，能满足大规模光电三维场景仿真的实时性要求。
为实现上述目的，本发明的具体步骤如下：
(1)设置曲线拟合仿真参数：
(1a)设置曲线拟合的气象条件参数；
(1b)设置曲线拟合仿真波段；
(1c)将与大气参量关联的探测器海拔高度、视线天顶角、路径长度设置为MODTRAN软件中的辐射传输几何路径参数；
(1d)将待仿真场景的几何尺度作为探测器海拔高度的最大值、视线天顶角的最大值和路径长度的最大值；
(2)将探测器视线天顶角的初始值赋值为0；
(3)划分探测器海拔高度子区间：
在探测器海拔高度取值范围[0.2，15.4]的区间内，将探测器海拔高度的划分为[0.2，4]、(4，7.8]、(7.8，11.6]、(11.6，15.4]四个子区间；
(4)按照下式，分别计算四个高度分段子区间中的每个探测器海拔高度值：
其中，hjM表示在第M个探测器海拔高度子区间内的第j个探测器海拔高度值，j在1≤j≤12中取整数值，hmin M表示探测器海拔高度中第M个子区间的左端点；
(5)计算路径长度最大有效值：
(5a)按照下式，计算视线天顶角大于90度时，辐射在穿过大气的传输路径上不与地球表面相交的路径长度的最大有效值：
其中，l'max表示辐射在穿过大气的传输路径上不与地球表面相交的路径长度的最大有效值，h表示探测器海拔高度，-表示负操作，θ表示视线天顶角，cos表示余弦操作；
(5b)视线天顶角不大于90度时，令辐射在穿过大气的传输路径上不与地球表面相交的路径长度最大有效值等于路径长度最大值；
(6)划分路径长度子区间：
(6a)判断视线天顶角是否大于90度；若是，则执行步骤(6b)；否则，执行步骤(7)；
(6b)按照下式，计算路径长度子区间分段点处的路径长度值：
其中，lPART表示路径长度子区间分段点处的路径长度值；e表示自然常数，其值为2.71828；
(6c)在路径长度的取值范围[0,l'max]的区间内，将路径长度划分为
[0,lPART]和(lPART,l'max]两个路径长度子区间；
(7)计算路径长度值：
(7a)按照下式，计算区间[0,lPART]的路径长度的起始长度控制因子：
其中，pstart表示区间[0,lPART]的路径长度子的起始路径长度控制因子；
(7b)按照下式，计算在区间[0,lPART]中选取的路径长度控制因子的个数：
其中，Q表示在区间[0,lPART]中选取的路径长度因子个数，[]表示向下取整；
(7c)按照下式，计算视线天顶角大于90度，区间[0,lPART]中选取的每个路径长度控制因子中的第q个路径长度控制因子：
其中，pq表示视线天顶角大于90度，在路径长度子区间[0,lPART]中选取的Q个路径长度控制因子中的第q个路径长度控制因子，pq中的q在1≤q≤Q中取整数值，pq-1表示在区间[0,lPART]中总共选取的Q个路径长度控制因子中的第q-1个路径长度控制因子，pq-1中的q在1＜q≤Q中取整数值；
(7d)按照下式，计算视线天顶角大于90度时，区间(lPART,l'max]的终止路径长度控制因子：
其中，pend表示路径长度子区间(lPART,l'max]的终止路径长度控制因子；
(7e)按照下式，计算视线天顶角大于90度时，区间[0,lPART]中的每个路径长度值：
其中，lq表示在区间[0,lPART]中所选取的每个路径长度值中的第q个路径长度值，q在1≤q≤Q取整数值；
(7f)按照下式，计算视线天顶角大于90度，在区间(lPART,l'max]中选取的每个路径长度控制因子中的第r个路径长度控制因子：
其中，pr表示视线天顶角大于90度，在路径长度子区间(lPART,l'max]中选取的每个路径长度控制因子中的第r个路径长度控制因子，pr中的r在1≤r≤R中取整数值，pr-1表示视线天顶角大于90度，在区间(lPART,l'max]中选取的每个路径长度控制因子中的第r-1个路径长度控制因子，pr-1中的r在1＜r≤R中取整数值，R表示在区间(lPART,l'max]中选取的路径长度控制因子的个数，当h＜2.6km时，R等于7，当h≥2.6km时，R等于9；
(7g)按照下式，分别计算视线天顶角大于90度时区间(lPART,l'max]中的每个路径长度值：
其中，lr表示在区间(lPART,l'max]中所选取的每个路径长度值中的第r个路径长度值，r在1≤r≤R取整数值；
(7h)按照下式，计算视线天顶角不大于90度时路径长度的取值范围[0,l'max]中的每个路径的长度：
其中，ln表示路径长度的取值范围[0,l'max]中所选取的每个路径长度值中的第n个路径的长度值；
(8)计算大气参量：
(8a)将辐射传输几何路径参数(θ,hjM,lq)和(θ,hjM,lr)输入MODTRAN软件中，分别计算视线天顶角大于90度时大气参量E(θ,hjM,lq)和大气参量E(θ,hjM,lr)；
(8b)将辐射传输几何路径参数(θ,hjM,ln)输入MODTRAN软件中，计算视线天顶角不大于90度时的大气参量E(θ,hjM,ln)；
(9)计算拟合曲线参数；
(9a)根据最小二乘方法，将视线天顶角不大于90度的ln和E(θ,hjM,ln)作为原始数据，计算使得拟合数据与原始数据的误差平方和最小时的θ和hjM取确定值时的大气参量随路径长度变化的曲线的八个曲线参数，其中，表示的第s+1项参数，l表示路径长度；
(9b)根据最小二乘方法，将视线天顶角大于90度的路径长度子区间[0,lPART]中的lq和E(θ,hjM,lq)作为原始数据，计算使得拟合数据与原始数据的误差平方和最小时的θ和hjM取确定值时的大气参量随路径长度变化的曲线的七个曲线参数，其中，表示的第s+1项参数；
(9c)根据最小二乘方法，将视线天顶角大于90度的路径长度子区间(lPART,l'max]中的lr和E(θ,hjM,lr)作为原始数据，计算使得拟合数据与原始数据的误差平方和最小时的θ和hjM取确定值时的大气参量随路径长度变化的曲线的七个参数，其中，表示的第s+1项参数；
(9d)根据最小二乘方法，将视线天顶角不大于90度的hjM和作为原始数据，计算s+1条使得拟合数据与原始数据的误差平方和最小时的θ取确定值时的大气参量随路径长度变化的曲线系数随探测器海拔高度变化的曲线A(θ)s(h)的七个曲线参数，d(θ)st表示A(θ)s(h)的第t+1项参数，h表示探测器海拔高度；
(9e)按照下式，计算视线天顶角大于90度，θ和hjM取确定值，路径子区间[0,lPART]的任意高度的大气参量随路径长度变化的曲线的各项系数：
其中，b(θ)s(h)表示视线天顶角大于90度，θ和hjM取确定值，路径子区间[0,lPART]任意探测器海拔高度的大气参量随路径长度变化的曲线的第s+1项参数，表示视线天顶角为θ，探测器海拔高度值为h(j-1)M时大气参量随路径长度变化的曲线的第s+1项系数，其中，h(j-1)M表示第M个探测器海拔高度子区间中第j-1个探测器海拔高度值；
(9f)按照下式，计算视线天顶角大于90度，θ和hjM取确定值时路径子区间(lPART,l'max]的任意高度的大气参量随路径长度变化的曲线的各项系数：
其中，c(θ)s(h)表示视线天顶角大于90度，θ和hjM取确定值，探测器海拔高度等于h时路径子区间(lPART,l'max]的大气参量随路径长度变化的曲线的第s+1项参数，表示θ和h(j-1)M取确定值时大气参量随路径长度变化的曲线的第s+1项系数；
(10)对视线天顶角进行累加5操作，得到累加后视线天顶角；
(11)判断累加后视线天顶角是否大于180度，若是，则执行步骤(12)；否则，执行步骤(3)；
(12)生成大气参量查找表：
(12a)按步骤(2)和步骤(10)中的方法对视线天顶角取值，在零至探测器海拔高度最大值范围内对探测器海拔高度随机取值，在零至路径长度最大有效值范围内对路径长度随机取值；
(12b)将选取的视线天顶角值、探测器海拔高度值和路径长度值代入步骤(9)确定的曲线中，计算大气参量，生成以视线天顶角、探测器海拔高度和路径长度为索引的大气参量查找表。
本发明与现有技术相比具有以下优点：
第一，由于本发明划分了路径长度子区间，按照不同的计算方法计算了不同路径长度子区间的路径长度值，克服了现有技术中辐射传输几何参数计算间隔太小，计算的数据量太大，计算耗时的缺陷。使得本发明具有能用更少的计算量，生成网格间隔更密的大气参量查找表，计算效率高的优点。
第二，由于本发明计算了路径长度最大有效值，对路径长度在零至路径长度最大有效值范围内取值，并且计算了大气参量随路径长度变化的曲线参数和大气参量随路径长度变化的曲线参数随探测器海拔高度变化的曲线参数，克服了现有技术中大气参量的计算精度低，计算耗时，计算大，计算场景中所有辐射传输几何路径的大气参量时，需遍历所有探测器海拔高度值的缺陷，使得本发明具有计算精度高，计算量少的优点。
附图说明
图1为本发明的流程图；
图2为本发明生成的大气透过率查找表中的大气透过率值与MODTRAN计算的标准大气透过率值的对比图；
图3为本发明生成的大气路径辐射查找表中的大气路径辐射值与MODTRAN计算的标准大气路径辐射值的对比图。
具体实施方式
下面结合附图对本发明做进一步的详细描述。
参照附图1，本发明的具体步骤如下。
步骤1，设置曲线拟合仿真参数。
第1步，设置仿真的气象条件参数。
本发明的实施例中，气象条件参数设置如下：
大气模式为1976年美国标准大气，温度和气压的高度廓线使用模式默认值，水汽的高度廓线使用模式默认值，臭氧的高度廓线使用模式默认值，甲烷的高度廓线使用模式默认值，一氧化二氮的高度廓线使用模式默认值，一氧化碳的高度廓线使用模式默认值，二氧化碳混合率330，气溶胶模式为乡村型能见度23km，气溶胶的季节性修正为由模式决定的修正，高层大气气溶胶为同温层背景气溶胶。
第2步，设置曲线拟合的仿真波段。
本发明的实施例中，仿真波段为3μm～5μm。
第3步，将与大气传输效应物理量关联的探测器海拔高度、视线天顶角、路径长度设置为MODTRAN软件中的辐射传输几何路径参数。
第4步，将待仿真场景的几何尺度作为探测器海拔高度的最大值、视线天顶角的最大值和路径长度的最大值。
本发明的实施例中，仿真场景的几何尺度为70km*70km，探测器海拔高度的最大值为15.4km，视线天顶角的最大值为180度，路径长度最大值为89.32km。
步骤2，将探测器视线天顶角的初始值赋值为0。
步骤3，划分探测器海拔高度子区间。
在探测器海拔高度取值范围[0.2，15.4]的区间内，将探测器海拔高度的划分为[0.2，4]、(4，7.8]、(7.8，11.6]、(11.6，15.4]四个子区间；
步骤4，确定探测器海拔高度值。
按照下式，分别对四个高度分段子区间中的每个探测器海拔高度取值：
其中，hjM表示在第M个探测器海拔高度子区间内的第j个探测器海拔高度值；hminM表示探测器海拔高度中第M个子区间的左端点；j在1≤j≤12中取整数值。
步骤5，获得路径长度最大有效值。
按照下式，计算视线天顶角大于90度时辐射在传输路径上不与地球表面相交的路径长度最大有效值：
其中，l'max表示辐射在穿过大气的传输路径上不与地球表面相交的路径长度的最大有效值，θ表示视线天顶角，cos表示余弦操作；h表示探测器海拔高度。
视线天顶角不大于90度时，令辐射在传输路径上不与地球表面相交的路径长度最大有效值等于路径长度最大值。
步骤6，划分路径长度子区间。
第1步，判断视线天顶角是否大于90度，若是，则执行步骤6中的第2步；否则，执行步骤7。
第2步，按照下式，确定路径长度子区间分段点处的路径长度值，得到划分路径长度后的路径长度子区间：
其中，lPART表示路径长度子区间分段点处的路径长度值；e表示自然常数，其值为2.71828。
第3步，在路径长度的取值范围[0,l'max]的区间内，将路径长度划分为[0,lPART]和(lPART,l'max]两个路径长度子区间。
步骤7，计算路径长度值。
按照下式，计算区间[0,lPART]的路径长度的起始长度控制因子：
其中，pstart表示区间[0,lPART]的路径长度子的起始路径长度控制因子；
按照下式，计算在区间[0,lPART]中选取的路径长度控制因子的个数：
其中，Q表示在区间[0,lPART]中选取的路径长度因子个数，[]表示向下取整；
按照下式，计算视线天顶角大于90度，区间[0,lPART]中选取的每个路径长度控制因子中的第q个路径长度控制因子：
其中，pq表示视线天顶角大于90度，在路径长度子区间[0,lPART]中选取的Q个路径长度控制因子中的第q个路径长度控制因子，pq中的q在1≤q≤Q中取整数值，pq-1表示在区间[0,lPART]中总共选取的Q个路径长度控制因子中的第q-1个路径长度控制因子，pq-1中的q在1＜q≤Q中取整数值；
按照下式，计算视线天顶角大于90度时，区间(lPART,l'max]的终止路径长度控制因子：
其中，pend表示路径长度子区间(lPART,l'max]的终止路径长度控制因子
按照下式，计算视线天顶角大于90度时，区间[0,lPART]中的每个路径长度值：
其中，lq表示在区间[0,lPART]中所选取的每个路径长度值中的第q个路径长度值，q在1≤q≤Q取整数值；
按照下式，计算视线天顶角大于90度，在区间(lPART,l'max]中选取的每个路径长度控制因子中的第r个路径长度控制因子：
其中，pr表示视线天顶角大于90度，在路径长度子区间(lPART,l'max]中选取的每个路径长度控制因子中的第r个路径长度控制因子，pr中的r在1≤r≤R中取整数值，pr-1表示视线天顶角大于90度，在区间(lPART,l'max]中选取的每个路径长度控制因子中的第r-1个路径长度控制因子，pr-1中的r在1＜r≤R中取整数值，R表示在区间(lPART,l'max]中选取的路径长度控制因子的个数，当h＜2.6km时，R等于7，当h≥2.6km时，R等于9；
按照下式，分别计算视线天顶角大于90度时区间(lPART,l'max]中的每个路径长度值：
其中，lr表示在区间(lPART,l'max]中所选取的每个路径长度值中的第r个路径长度值，r在1≤r≤R取整数值；
按照下式，计算视线天顶角不大于90度时路径长度的取值范围[0,l'max]中的每个路径的长度：
其中，ln表示路径长度的取值范围[0,l'max]中所选取的每个路径长度值中的第n个路径的长度值；
步骤8，计算大气参量。
将辐射传输几何路径参数(θ,hjM,lq)和(θ,hjM,lr)输入MODTRAN软件中，分别计算视线天顶角大于90度时大气参量E(θ,hjM,lq)和大气参量E(θ,hjM,lr)；
将辐射传输几何路径参数(θ,hjM,ln)输入MODTRAN软件中，计算视线天顶角不大于90度时的大气传输效应物理量E(θ,hjM,ln)；
步骤9，计算拟合曲线参数。
根据最小二乘方法，将视线天顶角不大于90度的ln和E(θ,hjM,ln)作为原始数据，计算使得拟合数据与原始数据的误差平方和最小时的θ和hjM取确定值时的大气参量随路径长度变化的曲线的八个曲线参数，其中，表示的第s+1项参数，l表示路径长度。
根据最小二乘方法，将视线天顶角大于90度的路径长度子区间[0,lPART]中的lq和E(θ,hjM,lq)作为原始数据，计算使得拟合数据与原始数据的误差平方和最小时的θ和hjM取确定值时的大气参量随路径长度变化的曲线的七个曲线参数，其中，表示的第s+1项参数。
根据最小二乘方法，将视线天顶角大于90度的路径长度子区间(lPART,l'max]中的lr和E(θ,hjM,lr)作为原始数据，计算使得拟合数据与原始数据的误差平方和最小时的θ和hjM取确定值时的大气参量随路径长度变化的曲线的七个参数，其中，表示的第s+1项参数。
根据最小二乘方法，将视线天顶角不大于90度的hjM和作为原始数据，计算s+1条使得拟合数据与原始数据的误差平方和最小时的θ取确定值时的大气参量随路径长度变化的曲线系数随探测器海拔高度变化的曲线A(θ)s(h)的七个曲线参数，d(θ)st表示A(θ)s(h)的第t+1项参数，h表示探测器海拔高度。
按照下式，计算视线天顶角大于90度，θ和hjM取确定值，路径子区间
[0,lPART]的任意高度的大气参量随路径长度变化的曲线的各项系数：
其中，b(θ)s(h)表示视线天顶角大于90度，θ和hjM取确定值，路径子区间[0,lPART]任意探测器海拔高度的大气参量随路径长度变化的曲线的第s+1项参数，表示视线天顶角为θ，探测器海拔高度值为h(j-1)M时大气参量随路径长度变化的曲线的第s+1项系数，其中，h(j-1)M表示第M个探测器海拔高度子区间中第j-1个探测器海拔高度值。
按照下式，计算视线天顶角大于90度，θ和hjM取确定值时路径子区间
(lPART,l'max]的任意高度的大气参量随路径长度变化的曲线的各项系数：
其中，c(θ)s(h)表示视线天顶角大于90度，θ和hjM取确定值，探测器海拔高度等于h时路径子区间(lPART,l'max]的大气参量随路径长度变化的曲线的第s+1项参数，表示θ和h(j-1)M取确定值时大气参量随路径长度变化的曲线的第s+1项系数。
步骤10，对视线天顶角进行累加5操作，得到累加后视线天顶角。
步骤11，判断累加5后视线天顶角θ是否大于180度，若是，则执行步骤12；否则，执行步骤3。
步骤12，生成大气参量查找表。
按步骤2和步骤10中的方法对视线天顶角取值，在零至探测器海拔高度最大值范围内对探测器海拔高度随机取值，在零至路径长度最大有效值范围内对路径长度随机取值；
将选取的视线天顶角值、探测器海拔高度值和路径长度值代入步骤9确定的曲线中，计算大气参量，生成以视线天顶角、探测器海拔高度和路径长度为索引的大气参量查找表。
本发明的实施例中，对视线天顶角在[0，180]范围内，每间隔5度进行取值；对探测器海拔高度在[0.2，15.4km]范围内，每间隔0.2km进行取值；对路径长度在[0.01218，l'max]范围内，每间隔进行取值，25≤k≤114。
下面结合仿真实验对本发明的效果做进一步的描述。
1.仿真内容
为了验证基于曲线拟合的大气参量查找表的有效性，生成了以视线天顶角、探测器海拔高度和路径长度为索引的大气透过率查找表和大气路径辐射查找表。
气象条件参数设置如下：大气模式为1975年美国标准大气；天气状况为晴天；近地面大气温度为0；边界层气溶胶模式中的消光系数和气象视距类型为乡村消光系数，缺省气象视距为23km；地面能见度为23km；光谱波段为3um-5um。大气透过率和大气路径辐射查找表中视线天顶角的范围为[0，180]，间隔为5度，进行取值探测器海拔高度范围为[0.2，15.4km]，间隔为0.2km，路径长度最大值为[0.01218，l'max]，间隔为25≤k≤114。
2.仿真效果分析
图2是本发明生成的大气透过率查找表中的大气透过率值与MODTRAN计算的标准大气透过率值的对比图。其中图2中的横坐标表示路径长度，纵坐标表示大气透过率。图2中的所有的圆形点表示视线天顶角为30度，探测器海拔高度为13.6km时MODTRAN软件计算的标准大气透过率值。所有的十字点表示视线天顶角为30度，探测器海拔高度为13.6km时本发明计算的大气透过率值。所有的正方形点表示视线天顶角为30度，探测器海拔高度为6.4km时MODTRAN软件计算的标准大气透过率值。所有的星形点表示视线天顶角为30度，探测器海拔高度为6.4km时本发明计算的大气透过率值。
比较图2中的数据，可以发现所有的圆形点与所有十字点几乎重合，所有的正方形点与所有的星形点几乎重合，可见本发明生成的大气参量查找表中大气透过率值与MODTRAN软件计算的标准大气透过率值具有很高的吻合度，计算精度高。
图3是本发明生成的大气路径辐射查找表中的大气路径辐射值和MODTRAN软件计算的标准大气路径辐射值的对比图，其中图2的横坐标表示路径长度，纵坐标表示大气路径辐射。图3中的所有的圆形点表示视线天顶角为30度，探测器海拔高度为13.6km时MODTRAN软件计算的标准大气路径辐射值。所有的十字点表示视线天顶角为30度，探测器海拔高度为13.6km时本发明计算的大气路径辐射值。所有的正方形点表示视线天顶角为30度，探测器海拔高度为6.4km时MODTRAN软件计算的标准大气路径辐射值。所有的星形点表示视线天顶角为30度，探测器海拔高度为6.4km时本发明计算的大气路径辐射值。
比较图3中的数据，可以发现所有的圆形点与所有十字点几乎重合，所有的正方形点与所有的星形点几乎重合，可见本发明生成的大气参量查找表中大气路径辐射值与MODTRAN软件计算的标准大气路径辐射值具有很高的吻合度，计算精度高。
Claims (2)
1.一种基于曲线拟合的大气参量查找表生成方法，包括如下步骤：
(1)设置曲线拟合仿真参数：
(1a)设置曲线拟合的气象条件参数；
(1b)设置曲线拟合仿真波段；
(1c)将与大气参量关联的探测器海拔高度、视线天顶角、路径长度设置为MODTRAN软件中的辐射传输几何路径参数；
(1d)将待仿真场景的几何尺度作为探测器海拔高度的最大值、视线天顶角的最大值和路径长度的最大值；
(2)将探测器视线天顶角的初始值赋值为0；
(3)划分探测器海拔高度子区间：
在探测器海拔高度取值范围[0.2，15.4]的区间内，将探测器海拔高度的划分为[0.2，4]、(4，7.8]、(7.8，11.6]、(11.6，15.4]四个子区间；
(4)按照下式，分别计算四个高度分段子区间中的每个探测器海拔高度值：
其中，hjM表示在第M个探测器海拔高度子区间内的第j个探测器海拔高度值，j在1≤j≤12中取整数值，hminM表示探测器海拔高度中第M个子区间的左端点；
(5)计算路径长度最大有效值：
(5a)按照下式，计算视线天顶角大于90度时，辐射在穿过大气的传输路径上不与地球表面相交的路径长度的最大有效值：
<mrow>
<msubsup>
<mi>l</mi>
<mi>max</mi>
<mo>&prime;</mo>
</msubsup>
<mo>=</mo>
<mfrac>
<mrow>
<mo>-</mo>
<mi>h</mi>
</mrow>
<mrow>
<mi>cos</mi>
<mi>&theta;</mi>
</mrow>
</mfrac>
</mrow>
其中，l'max表示辐射在穿过大气的传输路径上不与地球表面相交的路径长度的最大有效值，h表示探测器海拔高度，-表示负操作，θ表示视线天顶角，cos表示余弦操作；
(5b)视线天顶角不大于90度时，令辐射在穿过大气的传输路径上不与地球表面相交的路径长度最大有效值等于路径长度最大值；
(6)划分路径长度子区间：
(6a)判断视线天顶角是否大于90度；若是，则执行步骤(6b)；否则，执行步骤(7)；
(6b)按照下式，计算路径长度子区间分段点处的路径长度值：
<mrow>
<msub>
<mi>l</mi>
<mrow>
<mi>P</mi>
<mi>A</mi>
<mi>R</mi>
<mi>T</mi>
</mrow>
</msub>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mfrac>
<msubsup>
<mi>l</mi>
<mi>max</mi>
<mo>&prime;</mo>
</msubsup>
<msup>
<mi>e</mi>
<mn>3</mn>
</msup>
</mfrac>
</mtd>
<mtd>
<mrow>
<mi>h</mi>
<mo>&GreaterEqual;</mo>
<mn>2.6</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mfrac>
<msubsup>
<mi>l</mi>
<mi>max</mi>
<mo>&prime;</mo>
</msubsup>
<msup>
<mi>e</mi>
<mn>2</mn>
</msup>
</mfrac>
</mtd>
<mtd>
<mrow>
<mi>h</mi>
<mo><</mo>
<mn>2.6</mn>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
其中，lPART表示路径长度子区间分段点处的路径长度值；e表示自然常数，其值为2.71828；
(6c)在路径长度的取值范围[0,l'max]的区间内，将路径长度划分为[0,lPART]和(lPART,l'max]两个路径长度子区间；
(7)计算路径长度值：
(7a)按照下式，计算区间[0,lPART]的路径长度的起始长度控制因子：
<mrow>
<msub>
<mi>p</mi>
<mrow>
<mi>s</mi>
<mi>t</mi>
<mi>a</mi>
<mi>r</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mn>10</mn>
<mi>l</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<mn>1000</mn>
<msubsup>
<mi>l</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
<mo>&prime;</mo>
</msubsup>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mn>20</mn>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>h</mi>
<mo><</mo>
<mn>2.6</mn>
<mi>k</mi>
<mi>m</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mn>10</mn>
<mi>l</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<mn>1000</mn>
<msubsup>
<mi>l</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
<mo>&prime;</mo>
</msubsup>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mn>30</mn>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>h</mi>
<mo>&GreaterEqual;</mo>
<mn>2.6</mn>
<mi>k</mi>
<mi>m</mi>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
其中，pstart表示区间[0,lPART]的路径长度子的起始路径长度控制因子；
(7b)按照下式，计算在区间[0,lPART]中选取的路径长度控制因子的个数：
<mrow>
<mi>Q</mi>
<mo>=</mo>
<mo>&lsqb;</mo>
<mfrac>
<mrow>
<msub>
<mi>p</mi>
<mrow>
<mi>s</mi>
<mi>t</mi>
<mi>a</mi>
<mi>r</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>-</mo>
<mn>15</mn>
</mrow>
<mn>7</mn>
</mfrac>
<mo>&rsqb;</mo>
<mo>+</mo>
<mn>6</mn>
</mrow>
其中，Q表示在区间[0,lPART]中选取的路径长度因子个数，[]表示向下取整；
(7c)按照下式，计算视线天顶角大于90度，区间[0,lPART]中选取的每个路径长度控制因子中的第q个路径长度控制因子：
其中，pq表示视线天顶角大于90度，在路径长度子区间[0,lPART]中选取的Q个路径长度控制因子中的第q个路径长度控制因子，pq中的q在1≤q≤Q中取整数值，pq-1表示在区间[0,lPART]中总共选取的Q个路径长度控制因子中的第q-1个路径长度控制因子，pq-1中的q在1＜q≤Q中取整数值；
(7d)按照下式，计算视线天顶角大于90度时，区间(lPART,l'max]的终止路径长度控制因子：
<mrow>
<msub>
<mi>p</mi>
<mrow>
<mi>e</mi>
<mi>n</mi>
<mi>d</mi>
</mrow>
</msub>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mn>10</mn>
<mi>l</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<mn>1000</mn>
<msubsup>
<mi>l</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
<mo>&prime;</mo>
</msubsup>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mn>20</mn>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>h</mi>
<mo><</mo>
<mn>2.6</mn>
<mi>k</mi>
<mi>m</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mn>10</mn>
<mi>l</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<mn>1000</mn>
<msubsup>
<mi>l</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
<mo>&prime;</mo>
</msubsup>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mn>30</mn>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>h</mi>
<mo>&GreaterEqual;</mo>
<mn>2.6</mn>
<mi>k</mi>
<mi>m</mi>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
其中，pend表示路径长度子区间(lPART,l'max]的终止路径长度控制因子；
(7e)按照下式，计算视线天顶角大于90度时，区间[0,lPART]中的每个路径长度值：
<mrow>
<msub>
<mi>l</mi>
<mi>q</mi>
</msub>
<mo>=</mo>
<mfrac>
<msup>
<mi>e</mi>
<mrow>
<msub>
<mi>p</mi>
<mi>q</mi>
</msub>
<mo>/</mo>
<mn>10</mn>
</mrow>
</msup>
<mn>1000</mn>
</mfrac>
</mrow>
其中，lq表示在区间[0,lPART]中所选取的每个路径长度值中的第q个路径长度值，q在1≤q≤Q取整数值；
(7f)按照下式，计算视线天顶角大于90度，在区间(lPART,l'max]中选取的每个路径长度控制因子中的第r个路径长度控制因子：
<mrow>
<msub>
<mi>p</mi>
<mi>r</mi>
</msub>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mo>&lsqb;</mo>
<mn>10</mn>
<mi>l</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<mn>1000</mn>
<msubsup>
<mi>l</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
<mo>&prime;</mo>
</msubsup>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>r</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>p</mi>
<mrow>
<mi>r</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>-</mo>
<mi>r</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>p</mi>
<mrow>
<mi>r</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>-</mo>
<mi>r</mi>
<mo>+</mo>
<mn>1</mn>
<mo>></mo>
<msub>
<mi>p</mi>
<mi>q</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>p</mi>
<mrow>
<mi>e</mi>
<mi>n</mi>
<mi>d</mi>
</mrow>
</msub>
</mtd>
<mtd>
<mrow>
<msub>
<mi>p</mi>
<mrow>
<mi>r</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>-</mo>
<mi>r</mi>
<mo>+</mo>
<mn>1</mn>
<mo>></mo>
<msub>
<mi>p</mi>
<mi>q</mi>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
其中，pr表示视线天顶角大于90度，在路径长度子区间(lPART,l'max]中选取的每个路径长度控制因子中的第r个路径长度控制因子，pr中的r在1≤r≤R中取整数值，pr-1表示视线天顶角大于90度，在区间(lPART,l'max]中选取的每个路径长度控制因子中的第r-1个路径长度控制因子，pr-1中的r在1＜r≤R中取整数值，R表示在区间(lPART,l'max]中选取的路径长度控制因子的个数，当h＜2.6km时，R等于7，当h≥2.6km时，R等于9；
(7g)按照下式，分别计算视线天顶角大于90度时区间(lPART,l'max]中的每个路径长度值：
<mrow>
<msub>
<mi>l</mi>
<mi>r</mi>
</msub>
<mo>=</mo>
<mfrac>
<msup>
<mi>e</mi>
<mrow>
<msub>
<mi>p</mi>
<mi>r</mi>
</msub>
<mo>/</mo>
<mn>10</mn>
</mrow>
</msup>
<mn>1000</mn>
</mfrac>
</mrow>
其中，lr表示在区间(lPART,l'max]中所选取的每个路径长度值中的第r个路径长度值，r在1≤r≤R取整数值；
(7h)按照下式，计算视线天顶角不大于90度时路径长度的取值范围[0,l'max]中的每个路径的长度：
<mrow>
<msub>
<mi>l</mi>
<mi>n</mi>
</msub>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mn>0.001</mn>
<mo>*</mo>
<msup>
<mi>e</mi>
<mrow>
<mo>(</mo>
<mn>0.6</mn>
<mo>*</mo>
<mi>n</mi>
<mo>+</mo>
<mn>1.9</mn>
<mo>)</mo>
</mrow>
</msup>
</mrow>
</mtd>
<mtd>
<mrow>
<mn>1</mn>
<mo>&le;</mo>
<mi>n</mi>
<mo>&le;</mo>
<mn>3</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mn>0.001</mn>
<mo>*</mo>
<msup>
<mi>e</mi>
<mrow>
<mo>(</mo>
<mn>1.1</mn>
<mo>*</mo>
<mi>n</mi>
<mo>+</mo>
<mn>0.4</mn>
<mo>)</mo>
</mrow>
</msup>
</mrow>
</mtd>
<mtd>
<mrow>
<mn>4</mn>
<mo>&le;</mo>
<mi>n</mi>
<mo>&le;</mo>
<mn>9</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mn>0.001</mn>
<mo>*</mo>
<msup>
<mi>e</mi>
<mrow>
<mo>(</mo>
<mn>0.2</mn>
<mo>*</mo>
<mi>n</mi>
<mo>+</mo>
<mn>10.3</mn>
<mo>)</mo>
</mrow>
</msup>
</mrow>
</mtd>
<mtd>
<mrow>
<mn>10</mn>
<mo>&le;</mo>
<mi>n</mi>
<mo>&le;</mo>
<mn>13</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mn>0.001</mn>
<mo>*</mo>
<msup>
<mi>e</mi>
<mrow>
<mo>(</mo>
<mn>0.6</mn>
<mo>*</mo>
<mi>n</mi>
<mo>+</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
</msup>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>n</mi>
<mo>=</mo>
<mn>14</mn>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
其中，ln表示路径长度的取值范围[0,l'max]中所选取的每个路径长度值中的第n个路径的长度值；
(8)计算大气参量：
(8a)将辐射传输几何路径参数(θ,hjM,lq)和(θ,hjM,lr)输入MODTRAN软件中，分别计算视线天顶角大于90度时大气参量E(θ,hjM,lq)和大气参量E(θ,hjM,lr)；
(8b)将辐射传输几何路径参数(θ,hjM,ln)输入MODTRAN软件中，计算视线天顶角不大于90度时的大气参量E(θ,hjM,ln)；
(9)计算拟合曲线参数；
(9a)根据最小二乘方法，将视线天顶角不大于90度的ln和E(θ,hjM,ln)作为原始数据，计算使得拟合数据与原始数据的误差平方和最小时的θ和hjM取确定值时的大气参量随路径长度变化的曲线的八个曲线参数，其中，表示的第s+1项参数，l表示路径长度；
(9b)根据最小二乘方法，将视线天顶角大于90度的路径长度子区间[0,lPART]中的lq和E(θ,hjM,lq)作为原始数据，计算使得拟合数据与原始数据的误差平方和最小时的θ和hjM取确定值时的大气参量随路径长度变化的曲线的七个曲线参数，其中，表示的第s+1项参数；
(9c)根据最小二乘方法，将视线天顶角大于90度的路径长度子区间(lPART,l'max]中的lr和E(θ,hjM,lr)作为原始数据，计算使得拟合数据与原始数据的误差平方和最小时的θ和hjM取确定值时的大气参量随路径长度变化的曲线的七个参数，其中，表示的第s+1项参数；
(9d)根据最小二乘方法，将视线天顶角不大于90度的hjM和作为原始数据，计算s+1条使得拟合数据与原始数据的误差平方和最小时的θ取确定值时的大气参量随路径长度变化的曲线系数随探测器海拔高度变化的曲线A(θ)s(h)的七个曲线参数，d(θ)st表示A(θ)s(h)的第t+1项参数，h表示探测器海拔高度；
(9e)按照下式，计算视线天顶角大于90度，θ和hjM取确定值，路径子区间[0,lPART]的任意高度的大气参量随路径长度变化的曲线的各项系数：
<mrow>
<msub>
<mi>b</mi>
<mrow>
<mo>(</mo>
<mi>&theta;</mi>
<mo>)</mo>
<mi>s</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>h</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<msub>
<mi>b</mi>
<mrow>
<mo>(</mo>
<mi>&theta;</mi>
<mo>,</mo>
<msub>
<mi>h</mi>
<mrow>
<mi>j</mi>
<mi>M</mi>
</mrow>
</msub>
<mo>)</mo>
<mi>s</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>b</mi>
<mrow>
<mo>(</mo>
<mi>&theta;</mi>
<mo>,</mo>
<msub>
<mi>h</mi>
<mrow>
<mo>(</mo>
<mi>j</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
<mi>M</mi>
</mrow>
</msub>
<mo>)</mo>
<mi>s</mi>
</mrow>
</msub>
</mrow>
<mrow>
<msub>
<mi>h</mi>
<mrow>
<mi>j</mi>
<mi>M</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>h</mi>
<mrow>
<mo>(</mo>
<mi>j</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
<mi>M</mi>
</mrow>
</msub>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mo>*</mo>
<mi>h</mi>
</mrow>
其中，b(θ)s(h)表示视线天顶角大于90度，θ和hjM取确定值，路径子区间[0,lPART]任意探测器海拔高度的大气参量随路径长度变化的曲线的第s+1项参数，表示视线天顶角为θ，探测器海拔高度值为h(j-1)M时大气参量随路径长度变化的曲线的第s+1项系数，其中，h(j-1)M表示第M个探测器海拔高度子区间中第j-1个探测器海拔高度值；
(9f)按照下式，计算视线天顶角大于90度，θ和hjM取确定值时路径子区间(lPART,l'max]的任意高度的大气参量随路径长度变化的曲线的各项系数：
<mrow>
<msub>
<mi>c</mi>
<mrow>
<mo>(</mo>
<mi>&theta;</mi>
<mo>)</mo>
<mi>s</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>h</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<msub>
<mi>c</mi>
<mrow>
<mo>(</mo>
<mi>&theta;</mi>
<mo>,</mo>
<msub>
<mi>h</mi>
<mrow>
<mi>j</mi>
<mi>M</mi>
</mrow>
</msub>
<mo>)</mo>
<mi>s</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>c</mi>
<mrow>
<mo>(</mo>
<mi>&theta;</mi>
<mo>,</mo>
<msub>
<mi>h</mi>
<mrow>
<mo>(</mo>
<mi>j</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
<mi>M</mi>
</mrow>
</msub>
<mo>)</mo>
<mi>s</mi>
</mrow>
</msub>
</mrow>
<mrow>
<msub>
<mi>h</mi>
<mrow>
<mi>j</mi>
<mi>M</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>h</mi>
<mrow>
<mo>(</mo>
<mi>j</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
<mi>M</mi>
</mrow>
</msub>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mo>*</mo>
<mi>h</mi>
</mrow>
其中，c(θ)s(h)表示视线天顶角大于90度，θ和hjM取确定值，探测器海拔高度等于h时路径子区间(lPART,l'max]的大气参量随路径长度变化的曲线的第s+1项参数，表示θ和h(j-1)M取确定值时大气参量随路径长度变化的曲线的第s+1项系数；
(10)对视线天顶角进行累加5操作，得到累加后视线天顶角；
(11)判断累加后视线天顶角是否大于180度，若是，则执行步骤(12)；否则，执行步骤(3)；
(12)生成大气参量查找表：
(12a)按步骤(2)和步骤(10)中的方法对视线天顶角取值，在零至探测器海拔高度最大值范围内对探测器海拔高度随机取值，在零至路径长度最大有效值范围内对路径长度随机取值；
(12b)将选取的视线天顶角值、探测器海拔高度值和路径长度值代入步骤(9)确定的曲线中，计算大气参量，生成以视线天顶角、探测器海拔高度和路径长度为索引的大气参量查找表。
2.根据权利要求1所述的基于曲线拟合的大气参量查找表生成方法，其特征在于：步骤(1a)中所述的气象条件参数包括：大气模式，温度和气压的高度廓线，水汽的高度廓线，臭氧的高度廓线，甲烷的高度廓线，一氧化二氮的高度廓线，一氧化碳的高度廓线，二氧化碳混合率，气溶胶模式，气溶胶的季节性修正，高层大气气溶胶。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710678066.XA CN107563012B (zh) | 2017-08-10 | 2017-08-10 | 基于曲线拟合的大气参量查找表生成方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710678066.XA CN107563012B (zh) | 2017-08-10 | 2017-08-10 | 基于曲线拟合的大气参量查找表生成方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107563012A true CN107563012A (zh) | 2018-01-09 |
| CN107563012B CN107563012B (zh) | 2019-08-13 |
Family
ID=60975240
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710678066.XA Active CN107563012B (zh) | 2017-08-10 | 2017-08-10 | 基于曲线拟合的大气参量查找表生成方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107563012B (zh) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108984963A (zh) * | 2018-08-17 | 2018-12-11 | 北京电子工程总体研究所 | 一种地球红外背景辐射仿真方法及系统 |
| CN109884405A (zh) * | 2019-03-21 | 2019-06-14 | 重庆市计量质量检测研究院 | Vr眼镜综合测量方法 |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103674904A (zh) * | 2012-09-25 | 2014-03-26 | 中国航天科工集团第二研究院二〇七所 | 红外特性测量中大气传输快速修正方法 |
| CN103838968A (zh) * | 2014-02-28 | 2014-06-04 | 北京航天飞行控制中心 | 一种低仰角对流层折射修正方法 |
| US20140249753A1 (en) * | 2013-03-04 | 2014-09-04 | United States Government, As Represented By The Secretary Of The Navy | Atmospheric Density Correlation Technique |
| CN104217128A (zh) * | 2014-09-19 | 2014-12-17 | 北京空间机电研究所 | 一种起伏地形下卫星侧摆成像大气邻近效应模拟方法 |
| CN104279967A (zh) * | 2014-05-23 | 2015-01-14 | 中国科学院光电研究院 | 基于高光谱图像的气溶胶光学厚度反演方法 |
| CN104537180A (zh) * | 2015-01-04 | 2015-04-22 | 中国科学院国家天文台南京天文光学技术研究所 | 天文选址大气光学参数测量仪的数值模拟方法 |
| CN107016216A (zh) * | 2017-04-27 | 2017-08-04 | 北京环境特性研究所 | 一种计算无气象参数地点的大气透过率的方法 |
-
2017
- 2017-08-10 CN CN201710678066.XA patent/CN107563012B/zh active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103674904A (zh) * | 2012-09-25 | 2014-03-26 | 中国航天科工集团第二研究院二〇七所 | 红外特性测量中大气传输快速修正方法 |
| US20140249753A1 (en) * | 2013-03-04 | 2014-09-04 | United States Government, As Represented By The Secretary Of The Navy | Atmospheric Density Correlation Technique |
| CN103838968A (zh) * | 2014-02-28 | 2014-06-04 | 北京航天飞行控制中心 | 一种低仰角对流层折射修正方法 |
| CN104279967A (zh) * | 2014-05-23 | 2015-01-14 | 中国科学院光电研究院 | 基于高光谱图像的气溶胶光学厚度反演方法 |
| CN104217128A (zh) * | 2014-09-19 | 2014-12-17 | 北京空间机电研究所 | 一种起伏地形下卫星侧摆成像大气邻近效应模拟方法 |
| CN104537180A (zh) * | 2015-01-04 | 2015-04-22 | 中国科学院国家天文台南京天文光学技术研究所 | 天文选址大气光学参数测量仪的数值模拟方法 |
| CN107016216A (zh) * | 2017-04-27 | 2017-08-04 | 北京环境特性研究所 | 一种计算无气象参数地点的大气透过率的方法 |
Non-Patent Citations (2)
| Title |
|---|
| CALVIN K. LIANG ETC.: "《Improved VIIRS Day/Night Band Imagery With Near-Constant Contrast》", 《IEEE TRANSACTIONS ON GEOSCIENCE AND REMOTE SENSING》 * |
| 孙明明: "《红外场景的大气效应高效计算与仿真研究》", 《中国优秀硕士学位论文全文数据库信息科技辑》 * |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108984963A (zh) * | 2018-08-17 | 2018-12-11 | 北京电子工程总体研究所 | 一种地球红外背景辐射仿真方法及系统 |
| CN109884405A (zh) * | 2019-03-21 | 2019-06-14 | 重庆市计量质量检测研究院 | Vr眼镜综合测量方法 |
| CN109884405B (zh) * | 2019-03-21 | 2021-01-22 | 重庆市计量质量检测研究院 | Vr眼镜综合测量方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107563012B (zh) | 2019-08-13 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103353923B (zh) | 基于空间特征分析的自适应空间插值方法及其系统 | |
| CN105608326B (zh) | 一种山区复杂地形风场大涡模拟入口边界条件输入方法 | |
| CN111257241B (zh) | 一种deei的基于卫星观测的大气二氧化碳浓度反演算法 | |
| CN105842676A (zh) | 一种雷达盲区分析方法及装置 | |
| CN103454695B (zh) | 一种gps电离层tec层析方法 | |
| CN105224715A (zh) | 一种山区地貌下强风三维脉动风场综合模拟方法 | |
| CN104035096B (zh) | 一种基于多普勒天气雷达的垂直风廓线非线性反演方法 | |
| CN104867179B (zh) | 一种全谱段光学成像仪遥感影像仿真方法 | |
| CN105631218A (zh) | 基于idtcm的遥感地表温度时间归一化方法 | |
| CN106371627A (zh) | 一种测量绘画板位置信息的控制系统 | |
| CN106446444A (zh) | 基于贝叶斯最大熵与先验知识的土壤水分空间预测研究 | |
| CN108776340A (zh) | 一种基于遗传算法的顺轨干涉合成孔径雷达海表流场反演方法 | |
| CN104268322B (zh) | Weno差分方法的一种边界处理技术 | |
| CN103743679B (zh) | 一种加快温室气体垂直柱浓度反演速度的方法 | |
| CN106918797B (zh) | 一种基于标准数字功率源的数字化电能表检定和溯源方法 | |
| CN107563012A (zh) | 基于曲线拟合的大气参量查找表生成方法 | |
| CN114581791A (zh) | 基于modis数据大气水汽含量反演方法及系统 | |
| CN106895890A (zh) | 一种多声道超声波气体流量计声道权系数计算方法 | |
| CN102054294B (zh) | 基于曲面论和优化控制理论的曲面建模方法 | |
| CN104835203B (zh) | 一种基于OptiX的自然场景温度场计算方法及系统 | |
| CN117933095A (zh) | 一种基于机器学习的地表发射率实时反演及同化方法 | |
| CN115728842A (zh) | 一种大气可降水量获取方法、系统、装置及存储介质 | |
| CN103674904A (zh) | 红外特性测量中大气传输快速修正方法 | |
| CN101750616A (zh) | 植被风阻的测量方法及系统 | |
| CN109241212A (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 |