月基对地观测平台热红外传感器成像模拟方法

摘要

本发明公开了一种月基对地观测平台热红外传感器成像模拟方法，该方法能够实现大尺度甚至半球尺度热红外对地观测平台的成像模拟。测试结果表明本发明提出的月基对地观测平台热红外成像模拟方法模拟的地表温度时空分布数据与实际情况接近，能够反映真实的地表温度分布情况；最终模拟生成的月基对地观测平台热红外影像DN(Digital Number)值的空间分布与入瞳处辐射亮度影像十分相似，较好的描述了地表热红外辐射能量的空间分布特性。

说明书

技术领域

本发明涉及热红外对地观测成像模拟领域，尤其涉及一种月基对地观测平台热红外传感器成像模拟方法。

背景技术

月球作为地球的自然卫星，其轨道特性与人造卫星轨道差异明显，在月球上建立工作在热红外波段的地球观测平台(即月基对地观测平台)在获取全球及区域大尺度地表长波辐射及地表温度方面具有独特的优势，能够有效弥补现有人造卫星在对地观测领域数据获取的缺陷，为准确的获取全球尺度的地表温度及长波辐射提供了新的观测方式。但目前月基对地观测平台还处于理论论证阶段，在月基对地观测平台建成之前需要借助于模拟的方式对其整体的观测性能和影像质量进行分析。然而，月球的轨道特性及其观测特性与人造卫星存在较大差异，已有的遥感成像模拟模型无法直接用于月基对地观测平台热红外传感器的模拟。

利用计算机进行遥感成像系统模拟的物理模型能够以较低的成本对将要发射的遥感系统性能进行评估，进而为遥感系统的载荷设计等相关工作提供有力参考。因此，国内外众多研究机构和学者针对遥感成像系统模拟开展了大量研究工作，并已研究出多种遥感成像系统模拟模型且在计划发射的遥感卫星研制中进行了成功的应用(A,WiestL,Keller P,et al.SENSOR:a tool for the simulation of hyperspectral remotesensing systems[J].Isprs Journal of Photogrammetry&Remote Sensing,2001,55(5):299-312；Karl Segl,Luis Guanter,Christian Rogass,et al.EeteS—The EnMAP End-to-End Simulation Tool[J].IEEE Journal of Selected Topics in Applied EarthObservations&Remote Sensing,2012,5(2):522-530；Kaufmann H,Segl K,Richter R,etal.End-to-end sensor simulation for spectral band selection and optimizationwith application to the Sentinel-2 mission[J].Appl Opt,2012,51(4):439-449；Pandya M R,Shah D B,Trivedi H J,et al.Simulation of at-sensor radiance overland for proposed thermal channels of Imager payload onboard INSAT-3Dsatellite using MODTRAN model[J].Journal of Earth System Science,2011,120(1):19-25；Prosser C D.PICASSO:an end-to-end image simulation tool for space andairborne imaging systems II.Extension to the thermal infrared:equations andmethods[J].Proc Spie,2010,4(4):2840-2849；Schott J R,Brown S D,R V,etal.An Advanced Synthetic Image Generation Model and its Application to Multi/Hyperspectral Algorithm Development[J].Canadian Journal of Remote Sensing,1999,25(2):99-111.)。

目前已有的遥感成像模拟模型存在一些局限，主要表现在以下几个方面：

1)成像模拟模型主要是针对某一类型人造卫星或者航空平台搭载的成像系统而研发的，由于不同的遥感成像系统的观测目标、轨道特性、光谱通道、空间分辨率、成像方式等存在较大差异，在具体应用时需要有针对性的进行调整和修改；

2)热红外成像模拟模型主要是针对幅宽较小的中高分辨率遥感卫星，无法直接应用于大尺度热红外遥感卫星模拟。

发明内容

本发明所要解决的技术问题是针对月基对地观测平台的轨道特点，提出一种适合月基对地观测平台热红外传感器的大尺度成像模拟方法。

本发明的技术方案：

本发明提供了一种月基对地观测平台热红外传感器成像模拟方法，包括以下步骤：

步骤1：获取MODIS月合成地表发射率、地表温度(LST)、观测时间、NOAA提供的全球数字高程数据(DEM)；

步骤2：计算本初子午线整点时刻全球各地对应的当地太阳时；

步骤3：计算一天中任意时刻LST，根据某一时刻全球各地对应的当地太阳时，通过日变化模型计算出对应时刻的LST，即可得到该时刻全球LST空间分布数据；

步骤4：根据星历数据通过计算月球高度角获得月基对地观测平台热红外传感器观测区域二值影像；

步骤5：将步骤3中计算得到的某一时刻全球LST空间分布数据与步骤4中计算出的对应时刻的月基对地观测平台热红外传感器观测区域二值影像相乘获取观测范围内的LST空间分布；

步骤6：通过大气辐射传输模型MODTRAN模型计算出月基对地观测平台热红外传感器入瞳处辐射亮度；

步骤7：计算到达月基对地观测平台热红外传感器面元的热红外辐射通量E；

步骤8：计算得到月基对地观测平台热红外波段灰度影像。

优选地，在步骤1之前还包括对MODIS数据中LST和地表发射率中的低质量或者无效数据进行筛选过滤，高质量或有效LST数据的判定条件为(式1)-(式4)：

GP1＝((QC≤1)or(16≤QC≤17)or(64≤QC≤65)or(80≤QC≤81))(式1)

GP2＝(LST1＜LST2＞LST3＞LST4)(式2)

GP3＝(T4＜T1＜T2＜T3)(式3)

GP＝((GP1＝1)and(GP2＝1)and(GP3＝1))(式4)。

优选地，上述步骤中还包括通过MODIS LST的筛选结果对MODIS LST影像、发射率影像、观测时间影像及NOAA提供的全球数字高程数据(DEM)进行掩膜，MODIS LST筛选结果为二值影像，即有效像元用1表示，无效象元用0表示，以此作为掩膜基准影像，分别与MODISLST影像、发射率影像、观测时间影像及NOAA提供的全球数字高程数据(DEM)相乘，最终得到有效LST像元所对应的LST值、地表发射率、观测时间及其高程数值。

优选地，步骤7利用公式15计算到达传感器面元的热红外辐射通量

式中F＝f₀/D，为传感器光学系统F数，D为传感器光学系统的有效孔径，f₀是传感器光学系统的有效焦距，τ是光学系统的总透过率，L_s入瞳处辐射亮度，cos⁴α是光学系统导致的成像阴影效应，α是红外焦平面上探测象元与中心轴线之间的夹角，可通过公式16计算得到

式中d_x,d_y分别表示探测元之间的垂直和水平间距，m,n分别为探测器阵列行列数，即探测器成像尺寸，i,j分别为探测元在探测器阵列中的行列号，α_ij为探测器阵列中第行第列探测元所对应的α角。

优选地，步骤8利用公式17计算得到月基对地观测平台热红外波段灰度影像；

NearestInteger[]为取整运算函数，E代表某个象元的辐射强度，E_min观测范围内单个象元辐射强度最小值，E_max观测范围内单个象元辐射强度最大值，r代表红外成像系统A/D采样位数。

本发明的有益技术效果：本发明所述方法主要是针对中低分辨率大尺度对地观测平台的热红外成像模拟，与现有小尺度、中高分辨率星载热红外对地观测成像模拟方法相比，本发明所述方法能够实现大尺度甚至半球尺度热红外对地观测平台的成像模拟。测试结果表明本发明提出的月基对地观测平台热红外成像模拟方法模拟的地表温度时空分布数据与实际情况接近，能够反映真实的地表温度分布情况；最终模拟生成的月基对地观测平台热红外影像DN(Digital Number)值的空间分布与入瞳处辐射亮度影像十分相似，较好的描述了地表热红外辐射能量的空间分布特性。

附图说明

为了更清楚地说明本发明具体实施方式或现有技术中的技术方案，下面将对具体实施方式或现有技术描述中所需要使用的附图作简单地介绍。

图1是月基对地观测平台热红外传感器成像模拟流程图。

图2是月基对地观测平台视角下地表温度空间分布场景。

图3是月基对地观测平台热红外传感器各波段光谱响应函数。

图4是月基对地观测平台热红外传感器入瞳处辐射亮度影像模拟结果。

图5是月基对地观测平台热红外传感器红外焦平面辐照度模拟结果。

图6是月基对地观测平台热红外传感器成像模拟结果。

具体实施方式

图1是本发明的月基对地观测平台热红外传感器成像模拟方法流程示意图，包含以下步骤：

A.获取MODIS月合成地表发射率、地表温度(LST)、观测时间、NOAA提供的全球数字高程数据(DEM)。

B、MODIS高质量或有效LST数据筛选。

根据MODIS质量控制文件中8位整型二进制数据中每一位取值所代表的意义，选取无云条件下获取的质量较好的LST数据，且平均地表发射率和平均LST的误差分别不超过0.02和2K，通过二进制转十进制换算获取设定条件下MODIS质量控制文件的取值，实现高质量或有效LST和地表发射率的自动筛选。此外，由于MODIS获取的一天中四个时相的LST具有固定的时间顺序和变化规律，可通过此特点对变化异常的LST进行过滤。

高质量或有效LST和地表发射率的判定条件如下：

GP1＝((QC≤1)or(16≤QC≤17)or(64≤QC≤65)or(80≤QC≤81))(式1)

GP2＝(LST1＜LST2＞LST3＞LST4)(式2)

GP3＝(T4＜T1＜T2＜T3)(式3)

GP＝((GP1＝1)and(GP2＝1)and(GP3＝1))(式4)

其中，GP1表示MODIS LST数据质量判定结果；GP2表示MODIS四个时相LST的日变化规律判定结果；GP3表示MODIS四个时相LST获取时间的先后顺序；QC表示MODIS11C3数据的质量控制文件；LST1、LST2、LST3、LST4分别表示大约在当地时间10:30、13:30、22:30及第二天1:30左右获取的LST数据；T1、T2、T3和T4表示MODIS数据四个时相LST获取的时间，分别对应当地时间10:30、13:30、22:30及第二天1:30；GP表示最终筛选出的高质量或者有效的像元(有效像元用1表示，无效像元用0表示)。

C、通过MODIS LST的筛选结果对MODIS LST影像、发射率影像、观测时间影像及NOAA提供的全球数字高程数据(DEM)进行掩膜。

MODIS LST筛选结果为二值影像，即有效像元用1表示，无效象元用0表示，以此作为掩膜基准影像，分别与MODIS LST影像、发射率影像、观测时间影像及NOAA提供的全球数字高程数据(DEM)相乘，最终得到有效LST像元所对应的LST值、地表发射率、观测时间及其高程数值。

D、计算本初子午线整点时刻全球各地对应的当地太阳时。

本初子午线整点时刻全球各地对应的当地太阳时可由公式5计算得到

其中UTCtime表示计算得到的UTC时间，当UTCtime＜0时表示前一天，需加上24进行调整，当UTCtime＞24时表示第二天，需减去24进行调整；t₀表示本初子午线对应的时间，取值范围(0,24)；Lon表示经度，取值范围(-180°,180°)，其中东经为正，西经为负。

E、计算一天中任意时刻全球LST。

首先，以MODIS四个时相LST数据作为输入，通过(式6)-(式8)计算各像元一天中任意时刻的LST；然后通过选取t₀时刻各像元对应的当地太阳时的LST即可得到t₀时刻全球LST的空间分布数据。

其中T_day1(t)、T_day2(t)表示白天某一时刻的地表温度，T_night表示夜间地表温度，T₀表示一天中平均地表温度，t₁为上午地表温度某一时刻，T_a1、T_a2分别表示阶段一和阶段二中地表温度的变化幅度，其中T_a1(或T_a2)可由公式(6)和公式(7)在t₁处导数相等求解(见公式9)，_m表示一天中地表温度达到最大值时对应的时刻，为衰减常数。t_s为地表温度开始自由衰减的时刻，在本方法中将t_s设定为日落前一小时。ω₁、ω₂分别表示公式(6)和公式(7)的余弦函数半波宽，其中ω₁可由t₁减去日出时刻对应的时间t_rise得出，₂由一天中日照时间长度决定，是纬度和太阳赤纬δ的函数，可以通过公式(10)计算得出。δ可以表示为一年中第几天(day of the year,DOY)的函数，具体表达式见式(11)。

日出时刻t_rise与日落时刻t_ss可通过公式

衰减常数k可通过公式(7)、公式(8)在t_s处的导数相等求解，具体表达式见公式(14)：

E、通过计算月球高度角获得月基对地观测平台热红外传感器观测区域。

通过JPL星历数据DE405计算得到某一时刻(UTC时间)地球表面各处的月球高度角(此案例中选择的时间为本初子午线对应的UTC时间2015年1月15日00:00-23:00)，本发明中将月基对地观测平台设置在月球面向地球一面的中间位置，并将月球高度角作为月基对地观测平台热红外传感器观测高度角，月球高度角大于0的区域即为月基对地观测平台热红外传感器的可观测区域。将月球高度角大于0的像元值设为1反之则设为0，即可得到月基对地观测平台热红外传感器观测区域的二值影像。

F、月基对地观测平台热红外传感器观测范围内LST的获取。

将步骤D中计算得到的某一时刻全球LST空间分布数据与步骤E中计算出的对应时刻的月基对地观测平台热红外传感器观测区域二值影像相乘获取观测范围内的LST空间分布。结果如图2所示，图中展示的为本初子午线对应的UTC时间2015年1月15日月基对地观测平台热红外传感器观测范围内的LST空间分布场景，时间间隔为6小时。

G、通过大气辐射传输模型MODTRAN模型计算出月基对地观测平台热红外传感器入瞳处辐射亮度。

将同一时刻的月球高度角、月基对地观测平台热红外传感器观测范围内的LST、地表发射率、DEM、光谱响应函数及大气廓线数据作为输入，通过大气辐射传输模型MODTRAN模型计算出月基对地观测平台热红外传感器入瞳处辐射亮度。其中所用光谱响应函数为MODIS Terra/Aqua对应的第31、32波段(band31、band32)的光谱响应函数，光谱响应函数曲线如图3所示。入瞳处辐射亮度影像计算结果如图4所示，时间为本初子午线对应的UTC时间2015年1月15日04:00。

H、计算到达月基对地观测平台热红外传感器面元的热红外辐射通量E。

利用公式15计算到达月基对地观测平台热红外传感器面元的热红外辐射通量E，结果如图5所示。

式中F＝f₀/D，为传感器光学系统F数，D为传感器光学系统的有效孔径(本案例中取值为2m)，f₀是传感器光学系统的有效焦距(本案例中取值为2000mm)，τ是光学系统的总透过率(本案例中取值为0.5)，_s入瞳处辐射亮度，cos⁴α是光学系统导致的成像阴影效应，α是红外焦平面上探测象元与中心轴线之间的夹角，可通过公式16计算得到。

式中d_x,d_y分别表示探测元之间的垂直和水平间距，m,n分别为探测器阵列行列数，即探测器成像尺寸，i,j分别为探测元在探测器阵列中的行列号，α_ij为探测器阵列中第行第列探测元所对应的α角。

I、获取月基对地观测平台热红外传感器成像结果

利用公式17计算得到月基对地观测平台热红外波段灰度影像，结果如图6所示。

NearestInteger[]为取整运算函数，E代表某个象元的辐射强度，E_min观测范围内单个象元辐射强度最小值，E_max观测范围内单个象元辐射强度最大值，r代表红外成像系统A/D采样位数(本案例中取值为12bit)。

测试结果表明本发明提出的月基对地观测平台热红外成像模拟方法模拟的地表温度时空分布数据与实际情况接近，能够反映真实的地表温度分布情况；最终模拟生成的月基对地观测平台热红外影像DN(Digital Number)值的空间分布与入瞳处辐射亮度影像十分相似，较好的描述了地表热红外辐射能量的空间分布特性。

以上所述实施例，仅为本发明的具体实施方式，用以说明本发明的技术方案，而非对其限制，本发明的保护范围并不局限于此，尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，其依然可以对前述实施例所记载的技术方案进行修改或可轻易想到变化，或者对其中部分技术特征进行等同替换；而这些修改、变化或者替换，并不使相应技术方案的本质脱离本发明实施例技术方案的精神和范围，都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应所述以权利要求的保护范围为准。