一种静轨面阵遥感相机分时动态成像模拟方法与流程

本发明涉及一种静轨面阵遥感相机分时动态成像模拟方法，属于工程科学技术中的遥感成像模拟技术领域，在遥感成像仿真方面能够发挥重要作用。

背景技术：

地球静止轨道卫星运行高度约为36000公里，其轨道平面与地球的赤道平面重合，地球静止轨道卫星从地球上看是静止于赤道上空的某点，因此静轨对地观测系统具有随时观测的优点，时间分辨率的优势十分明显，采用面阵探测器的凝视型遥感相机则能够充分利用静止轨道对地观测的优势，实现多光谱谱段的分时动态成像，即能够实现单谱段单次成像、全谱段单次成像、单谱段连续成像、全谱段连续成像等不同模式。

针对静轨面阵遥感相机进行成像建模与仿真，可以深化人们对于此类型遥感相机成像过程和成像数据特性的理解，也可用于预测此类型载荷的成像质量并提供仿真数据，为新载荷的系统设计和地面数据处理系统建立提供支持。

目前还没有专门针对静轨面阵遥感相机的成像仿真模型，传统的遥感相机成像仿真模型一般依赖于载荷的基本特性参数，对载荷本身的物理成像过程描述不够全面，这导致传统相机仿真模型无法在成像模式上对新相机进行成像模拟，限制了相机仿真模型的应用能力。

技术实现要素：

本发明的目的是：提供一种静轨面阵遥感相机分时动态成像模拟方法，建立了针对静轨面阵遥感相机的动态成像模型，克服了现有相机仿真模型无法在成像模式上进行仿真及现有成像模拟方法对高时频动态成像适用性不高的问题，实现了面阵遥感相机在对地观测时分时动态成像模拟，有效地提高了遥感相机仿真的真实性，极大地提高了成像仿真模型的应用能力。

本发明的技术解决方案是：基于真实的静轨面阵遥感相机成像过程，结合凝视相机成像模式，构建出充分考虑了时空变化特性的分时成像仿真机制，并在遥感相机结构参数层面上计算不同时间、不同成像模式情况下的相机输出DN值，实现分时动态成像的真实模拟。

本发明一种静轨面阵遥感相机分时动态成像模拟方法，其步骤如下：

步骤一：输入指定成像区域对应的高时频入瞳辐亮度图像；

步骤二：设置面阵遥感相机的成像指令；

步骤三：设置面阵遥感相机系统参数；

步骤四：根据步骤二确定的成像指令，利用步骤一输入的入瞳辐亮度、步骤三输入的相机系统参数进行模拟计算，得到DN值，并生成DN值图像；

步骤五：根据步骤四得到的DN值图像判断成像指令是否全部执行完成，是则输出0级仿真数据，否则重复步骤四。

所述步骤二中，设置面阵遥感相机的成像指令如下：

任意选择成像模式中的1种或多种进行组合，形成成像指令，其中所述成像模式如下：

a)可见光至红外0.4～14μm光谱范围中的任意1个波段进行1次仿真成像，定义该模式为单谱段单次成像；

b)可见光至红外0.4～14μm光谱范围中的任意M(M≥2)个波段进行1次仿真成像，定义该模式为多谱段单次成像；M≥2

c)可见光至红外0.4～14μm光谱范围中的任意1个波段进行T₁、T₂、…T_N时刻的连续仿真成像，定义该模式为单谱段连续成像；N≥2

d)可见光至红外0.4～14μm光谱范围中的任意M个波段进行T₁、T₂、…T_N时刻的连续仿真成像，定义该模式为多谱段连续成像。

所述步骤四中，生成DN值图像的方法如下：

a)单谱段单次成像模式时，假设对T₁时刻的λ₁波段入瞳辐亮度L(T₁,λ₁)进行成像模拟，从而得到DN值，DN值计算公式如下：

其中，C是转换因子，G、B分别是运算放大器的增益和偏置，N_full位探测器满阱电子数，bit为量化位数，N(T₁,λ₁)是探测器在T₁时刻λ₁波段产生的光生电子数，

其中，A为探测器探元面积；t(T₁,λ₁)为T₁时刻λ₁波段时相机的积分时间；ε为光学孔径面积遮拦比；η(λ₁)为量子效率；F为相机的F数；h为普朗克常数；c为光速；τ(λ₁)为光学系统透过率；R(λ₁)光谱响应函数；

DN值图像为数据立方体，XY维空间为T₁时刻λ₁波段的图像信息，、Z维对应成像时刻T₁和波段λ₁组合(T₁,λ₁)，其中Z维的维度为1。

b)多谱段单次成像模式时，假设对T₁时刻的λ₁、λ₂、…λ_M波段入瞳辐亮度L(T₁,λ₁,λ₂...λ_M)进行成像模拟，从而生成DN值，DN值计算公式如下：

其中，N(T₁,λ₁,λ₂...λ_M)是探测器在T₁时刻λ₁、λ₂、…λ_M波段产生的光生电子数；

DN值图像为数据立方体，XY维空间为T₁时刻λ₁、λ₂、…λ_M波段的图像信息，Z维对应成像时刻T₁和波段λ₁、λ₂、…λ_M组合(T₁,λ₁)、(T₁,λ₂)…(T₁,λ_M)，其中Z维的维度为M；

c)单谱段连续成像模式时，假设对T₁、T₂、…T_N时刻的λ₁波段入瞳辐亮度L(T₁,T₂...T_N,λ₁)进行连续成像模拟，从而生成DN值，DN值计算公式如下：

其中，N(T₁,T₂...T_N,λ₁)是探测器在T₁、T₂、…T_N时刻λ₁波段产生的光生电子数。

DN值图像为数据立方体，XY维空间为T₁、T₂、…T_N时刻λ₁波段的图像信息，Z维对应成像时刻T₁、T₂、…T_N和波段λ₁组合(T₁,λ₁)、(T₂,λ₁)…(T_N,λ₁)，其中Z维的维度为N；

d)多谱段连续成像模式时，假设对T₁、T₂、…T_N时刻的λ₁、λ₂、…λ_M波段入瞳辐亮度L(T₁,T₂...T_N,λ₁,λ₂...λ_M)进行成像模拟，从而生成DN值，DN值计算公式如下：

其中，N(T₁,T₂...T_N,λ₁,λ₂...λ_M)是探测器在T₁、T₂、…T_N时刻λ₁、λ₂、…λ_M波段产生的光生电子数。

DN值图像为数据立方体，XY维空间为T₁、T₂、…T_N时刻下λ₁、λ₂、…λ_M波段的图像信息，Z维对应成像时刻T₁、T₂、…T_N和波段λ₁、λ₂、…λ_M组合(T₁,λ₁)、(T₁,λ₂)…(T₁,λ_M)、(T₂,λ₁)、(T₂,λ₂)…(T₂,λ_M)、(T_N,λ₁)、(T_N,λ₂)…(T_N,λ_M)，其中Z维的维度为N*M。

本发明一种静轨面阵遥感相机分时动态成像模拟方法，与现有技术相比的优点在于：本发明通过设置成像指令，即利用单谱段单次成像、多谱段单次成像、单谱段连续成像、多谱段连续成像等成像模式模拟遥感相机的分时动态成像过程，真实反映了静轨面阵遥感相机获取数据的过程，克服了现有成像模拟方法对高时频动态成像适用性不高的问题，有效地提高了遥感相机仿真的真实性。

附图说明

图1为本发明流程框图。

具体实施方式

如图1所示，本发明的具体实施方法如下：

步骤一：输入指定成像区域对应的高时频入瞳辐亮度图像

输入指定成像区域对应高时频入瞳辐亮度图像，此处的高时频是指气候变化与生长季等能够引起地物表面空间形态与辐射特性发生变化的时间频率。

步骤二：面阵遥感相机成像指令设置

任意选择成像模式中的1种或多种进行组合，形成成像指令，其中成像模式的定义如下：

a)可见光至红外0.4～14μm光谱范围中的任意1个波段进行1次仿真成像，定义该模式为单谱段单次成像；

b)可见光至红外0.4～14μm光谱范围中的任意M(M≥2)个波段进行1次仿真成像，定义该模式为多谱段单次成像；

c)可见光至红外0.4～14μm光谱范围中的任意1个波段进行T₁、T₂、…T_N(N≥2)时刻的连续仿真成像，定义该模式为单谱段连续成像；

d)可见光至红外0.4～14μm光谱范围中的任意M(M≥2)个波段进行T₁、T₂、…T_N(N≥2)时刻的连续仿真成像，定义该模式为多谱段连续成像。

步骤三：面阵遥感相机系统参数设置

设置面阵遥感相机系统参数，具体参数包括：探测器探元面积、积分时间、光学孔径面积遮拦比、量子效率、相机F数、光学系统透过率、光谱响应函数、转换因子、满阱电子数、量化位数、运算放大器的增益和偏置。

步骤四：根据步骤二确定的成像指令，利用步骤一输入的入瞳辐亮度、步骤三输入的相机系统参数进行模拟计算，得到DN值，并生成DN值图像，具体方式如下：

a)单谱段单次成像模式时，假设对T₁时刻的λ₁波段入瞳辐亮度L(T₁,λ₁)进行成像模拟，从而生成DN值，DN值计算公式如下：

其中，C是转换因子，G、B分别是运算放大器的增益和偏置，N_full位探测器满阱电子数，bit为量化位数，N(T₁,λ₁)是探测器在T₁时刻λ₁波段产生的光生电子数，

其中，A为探测器探元面积；t(T₁,λ₁)为T₁时刻λ₁波段时相机的积分时间；ε为光学孔径面积遮拦比；η(λ₁)为量子效率；F为相机的F数；h为普朗克常数；c为光速；τ(λ₁)为光学系统透过率；R(λ₁)光谱响应函数。

DN值图像为数据立方体，XY维空间为T₁时刻λ₁波段的图像信息，、Z维对应成像时刻T₁和波段λ₁组合(T₁,λ₁)，其中Z维的维度为1。

b)多谱段单次成像模式时，假设对T₁时刻的λ₁、λ₂、…λ_M波段入瞳辐亮度L(T₁,λ₁,λ₂...λ_M)进行成像模拟，从而生成DN值，DN值计算公式如下：

其中，N(T₁,λ₁,λ₂...λ_M)是探测器在T₁时刻λ₁、λ₂、…λ_M波段产生的光生电子数。

DN值图像为数据立方体，XY维空间为T₁时刻λ₁、λ₂、…λ_M波段的图像信息，Z维对应成像时刻T₁和波段λ₁、λ₂、…λ_M组合(T₁,λ₁)、(T₁,λ₂)…(T₁,λ_M)，其中Z维的维度为M。

c)单谱段连续成像模式时，假设对T₁、T₂、…T_N时刻的λ₁波段入瞳辐亮度L(T₁,T₂...T_N,λ₁)进行连续成像模拟，从而生成DN值，DN值计算公式如下：

其中，N(T₁,T₂...T_N,λ₁)是探测器在T₁、T₂、…T_N时刻λ₁波段产生的光生电子数。

DN值图像为数据立方体，XY维空间为T₁、T₂、…T_N时刻λ₁波段的图像信息，Z维对应成像时刻T₁、T₂、…T_N和波段λ₁组合(T₁,λ₁)、(T₂,λ₁)…(T_N,λ₁)，其中Z维的维度为N。

d)多谱段连续成像模式时，假设对T₁、T₂、…T_N时刻的λ₁、λ₂、…λ_M波段入瞳辐亮度L(T₁,T₂...T_N,λ₁,λ₂...λ_M)进行成像模拟，从而生成DN值，DN值计算公式如下：

其中，N(T₁,T₂...T_N,λ₁,λ₂...λ_M)是探测器在T₁、T₂、…T_N时刻λ₁、λ₂、…λ_M波段产生的光生电子数。

DN值图像为数据立方体，XY维空间为T₁、T₂、…T_N时刻下λ₁、λ₂、…λ_M波段的图像信息，Z维对应成像时刻T₁、T₂、…T_N和波段λ₁、λ₂、…λ_M组合(T₁,λ₁)、(T₁,λ₂)…(T₁,λ_M)、(T₂,λ₁)、(T₂,λ₂)…(T₂,λ_M)、(T_N,λ₁)、(T_N,λ₂)…(T_N,λ_M)，其中Z维的维度为N*M。

步骤五：根据步骤四得到的DN值图像，判断DN值图像的数量是否与本次成像指令中规定的成像模式次数相同。若相同，则成像指令全部执行完成，按照遥感数据存储标准输出0级仿真图像；若不相同，则成像指令未全部执行完成，重复步骤四。

本发明说明书中未作详细描述的内容属本领域技术人员的公知技术。