一种基于地物分布变化的地球背景红外辐射特性仿真方法技术

本发明专利技术公开了一种基于地物分布变化的地球背景红外辐射特性仿真方法，涉及地球背景红外辐射仿真技术领域

【技术实现步骤摘要】
一种基于地物分布变化的地球背景红外辐射特性仿真方法


[0001]本专利技术涉及地球背景红外辐射仿真
，具体为一种基于地物分布变化的地球背景红外辐射特性仿真方法
。

技术介绍

[0002]预警卫星对于目标的识别是天基反导体系的第一步
。
地球背景中不同地物类型的辐射特性差别较大，建立完善的地球背景红外辐射模型对于预警卫星的目标识别十分重要
。
由于卫星的探测受到自身传感器观测条件
、
大气活动等诸多因素的干扰，难以获得大量的辐射数据来满足研究的需要，因此需要对地球背景的辐射特性进行高分辨率的图像仿真
。
目前的地球背景红外辐射特性仿真技术所涉及的地物类型较少
、
研究波段较窄以及没有考虑季节变化对于地物分布变化的影响
。

技术实现思路

[0003]为了克服上述现有技术的不足，本专利技术提供了一种基于地物分布变化的地球背景红外辐射特性仿真方法，能够考虑季节变化对于地物分布的影响，实现2～
14
μ
m
任意波段，多种地物类型的地球背景红外辐射特性仿真
。
[0004]本专利技术的技术方案为：
[0005]一种基于地物分布变化的地球背景红外辐射特性仿真方法，包括以下步骤：
[0006]步骤1：以
MODIS12
全球土地覆盖类型数据为基础，根据
MODIS13
全球
NDVI
数据季节变化，基于<br/>GEE
平台和
ARCGIS
遥感软件建立两季的全球地物分布模型；
[0007]步骤2：基于全球地物分布模型和
ECOSTRESS
光谱数据库，建立地物的光谱模型；
[0008]步骤3：基于地物的光谱模型数据，利用
MODTRAN
进行辐射计算，生成地球背景红外辐射区域图像
。
[0009]作为优选，步骤1中，以两季的全球
NDVI
变化建立两季的地物分布模型，以连续若干年的
MODIS13A1
数据为研究对象，利用
GEE
平台获取
MODIS12Q1
的原始数据并进行拼接，利用
ARCGIS
进行坐标投影转化和单位转化，最终得到全球
NDVI
分布数据
。
[0010]作为优选，步骤1中，以两季的全球
NDVI
变化建立两季的地物分布模型，以连续若干年的
MODIS13A1
数据为研究对象，利用
GEE
平台获取
MODIS12Q1
的原始数据并进行拼接，利用
ARCGIS
进行坐标投影转化和单位转化，最终得到全球
NDVI
分布数据
。
[0011]作为优选，所述坐标投影转化为：将表示
MODIS
数据坐标的
SIN
正弦曲线转化为
WGS
‑
1984
坐标系；所述单位转化为：将数据的原始数值乘以比例因子得到
NDVI
值
。
[0012]作为优选，一季的全球地物分布模型建模步骤如下：将一年内间隔6个月的
NDVI
数据做差，然后对差值数据进行重分类，将重分类数据与地物类型数据关联，具体关联方法为：预设第一阈值
a、
第二阈值
b
，且
a
＜
b
，当
NDVI
小于
a
时为第一类，地物类型不做处理；当
NDVI
大于
a
小于
b
时为第二类，将绿色草地改变为枯草地，绿色森林改变为枯黄落叶林；当
NDVI
大于
b
时为第三类，将所有地物类型改变为雪地
。
本专利技术选择
a
＝
0.4
，
b
＝
0.6。
[0013]作为优选，步骤3中，所述地球背景红外辐射区域图像的生成步骤如下：
[0014]步骤
3.1
，根据所述地物的光谱模型和地表温度模型确定包括地表发射率和温度的辐射计算参数；
[0015]步骤
3.2
，将所述辐射计算参数代入到
MODTRAN
大气辐射计算软件获得不同季节和不同太阳天顶角的地球背景红外辐射模型，建立辐射计算参数
‑
季节太阳天顶角
‑
红外辐射的数据库；
[0016]步骤
3.3
，基于该数据库，输入区域经纬度和时间信息查找红外辐射数据，得到目标区域的红外辐射特性分布
。
[0017]作为优选，地表温度模型：基于
MODIS
遥感数据中的
MODIS02
辐射数据
、MODIS03
地理定位数据和
MODIS35
云掩膜数据进行拼接
、
定位
、
去云操作，最后将得到的数据通过分裂窗算法进行温度的反演，得到所述地表温度分布模型
。
[0018]作为优选，地球背景红外辐射模型：为区域个点所在的地表类型以及地表温度之间的关系，通过建立目标点经纬度与地物分布模型以及地表温度模型的关系得到
。
[0019]作为优选，所述目标点经纬度为区域目标点与经纬度之间的关系，表示如下：
[0020][0021]式中
:A
line
为行数；
L
at
为目标点所在的纬度，其中北纬为正，南纬为负；
L
为模型矩阵的总行数；
[0022][0023]式中
:A
cloumn
为列数；
L
on
为目标点所在的经度
(
西经：0～
180
°
；东经
180
～
360
°
)
；
C
为模型矩阵的总列数；
[0024]W
offset
＝
A
line
×
L+A
cloumn (6)
[0025]式中：
W
offset
为目标点对应的数据在模型文件中的偏移量
。
[0026]与现有技术相比，本专利技术的有益效果是：
[0027](1)
目前现有的地物类型数据都是一年度的，是没有考虑季节的因素，本专利技术基于
NDVI
变化实现了冬夏两季的全球地物分布模型；
[0028](2)
根据
ECOSTRESS
...

【技术保护点】

【技术特征摘要】
1.
一种基于地物分布变化的地球背景红外辐射特性仿真方法，其特征在于，包括以下步骤：步骤1：以
MODIS12
全球土地覆盖类型数据为基础，根据
MODIS13
全球
NDVI
数据季节变化，基于
GEE
平台和
ARCGIS
遥感软件建立两季的全球地物分布模型；步骤2：基于全球地物分布模型和
ECOSTRESS
光谱数据库，建立地物的光谱模型；步骤3：基于地物的光谱模型数据，利用
MODTRAN
进行辐射计算，生成地球背景红外辐射区域图像
。2.
根据权利要求1所述的一种基于地物分布变化的地球背景红外辐射特性仿真方法，其特征在于，步骤1中，以两季的全球
NDVI
变化建立两季的地物分布模型，以连续若干年的
MODIS13A1
数据为研究对象，利用
GEE
平台获取
MODIS12Q1
的原始数据并进行拼接，利用
ARCGIS
进行坐标投影转化和单位转化，最终得到全球
NDVI
分布数据
。3.
根据权利要求2所述的一种基于地物分布变化的地球背景红外辐射特性仿真方法，其特征在于，所述坐标投影转化为：将表示
MODIS
数据坐标的
SIN
正弦曲线转化为
WGS
‑
1984
坐标系；所述单位转化为：将数据的原始数值乘以比例因子得到
NDVI
值
。4.
根据权利要求1至3任意一项所述的一种基于地物分布变化的地球背景红外辐射特性仿真方法，其特征在于，一季的全球地物分布模型建模步骤如下：将一年内间隔6个月的
NDVI
数据做差，然后对差值数据进行重分类，将重分类数据与地物类型数据关联，具体关联方法为：预设第一阈值
a、
第二阈值
b
，且
a
＜
b
，当
NDVI
小于
a
时为第一类，地物类型不做处理；当
NDVI
大于
a
小于
b
时为第二类，将绿色草地改变为枯草地，绿色森林改变为枯黄落叶林；当
NDVI
大于
b
时为第三类，将所有地物类型改变为雪地
。5.
根据权利要求1所述的一种基于地物分布变化的地...

【专利技术属性】
技术研发人员：邢世林，史波，施杰，
申请(专利权)人：南京航空航天大学，
类型：发明
国别省市：