卫星遥感产品支持的单通道法地表温度反演方法专利检索-小卫星人造卫星太空飞行器与操作专利检索查询-专利查询网

卫星遥感产品支持的单通道法地表 温度反演方法

阅读：990发布：2024-02-25

专利汇可以提供卫星遥感产品支持的单通道法地表温度反演方法专利检索，专利查询，专利分析的服务。并且本发明公开了一种卫星遥感产品支持的单通道法地表温度反演方法，包括如下步骤：a、收集与待反演影像同时相的16天合成的MODIS地表植被指数产品NDVI；b、通过MODIS地表植被指数产品NDVI与HJ-1B？IRS 近红外波段反射率在降序排列下各个数值区间的像元个数百分比匹配，实现对HJ-1B？IRS热红外波段对应像元的NDVI植被指数赋值；c、基于NDVI与地表比辐射率的统计关系模型，计算地表比辐射率；d、获取与IRS数据同时相的MODIS 水汽通道数据，用三通道比值法来实现大气水汽含量的反演，确定HJ-1B？IRS的水汽含量；e、将地表比辐射率和水汽含量代入QK&B 算法或JM&S算法得到地表温度反演结果。本发明中地表比辐射率确定精度较高，利于实现该类型传感器单通道法地表温度的反演。，下面是卫星遥感产品支持的单通道法地表温度反演方法专利的具体信息内容。

权利要求

1.卫星遥感产品支持的单通道法地表温度反演方法，其特征在于包括如下步骤:
a、收集与待反演影像同时相的16天合成的MODIS地表植被指数产品NDVI；
b、将步骤a中获得的MODIS地表植被指数产品NDVI进行投影转换，转换到与IRS数据相统一的UTM投影下，并对MODIS地表植被指数产品NDVI进行重采样运算；
将MODIS NDVI数值从大到小等间隔进行密度分割，并统计各个区间内的像元个数及其所占百分比；然后将IRS第1波段反射率值同样按照降序排列统计，以与MODIS NDVI区间范围占相同百分比为匹配条件对图像进行分割；
根据MODIS影像的NDVI区间范围，对分割后的IRS影像赋予新的NDVI值；
c、基于NDVI与地表比辐射率的统计模型，计算地表比辐射率；
d、获取与IRS数据同时相的MODIS水汽通道数据，用三通道比值法来实现大气水汽含量的反演，确定HJ-1B IRS的水汽含量；
e、将步骤c和步骤d中得到的地表比辐射率ελ和水汽含量代入QK&B 算法或JM&S算法得到地表温度反演结果。
2.根据权利要求1所述的卫星遥感产品支持的单通道法地表温度反演方法，其特征在于，所述步骤b中，具体赋值方法如下：
b1、根据土壤的光谱特征，当NDVI＜0.2，认为地表全部由裸土构成，对应IRS近红外波段NDVI取0.2；
b2、考虑当植被覆盖较高时，NDVI趋于饱和，因此当NDVI＞0.7，将对应IRS近红外波段NDVI统一取最大值0.8；
b3、当NDVI在0.2-0.7之间，认为地表由植被和裸土两种地表构成，以0.1的间隔等分，对应IRS近红外波段赋予各区间的上限值，NDVI依次取0.3-0.7。
3.根据权利要求1所述的卫星遥感产品支持的单通道法地表温度反演方法，其特征在于，所述步骤c中，在NDVI确定之后，植被和土壤的发射率εv、εs分别由ASTER光谱库提供的发射率曲线与IRS热红外波段光谱响应函数进行卷积运算获得：εv=0.9833，εs=0.9732；
地表比辐射率ελ的计算公式为：ελ＝0.9732+0.0101*NDVI。

说明书全文

卫星遥感产品支持的单通道法地表温度反演方法

技术领域

[0001] 本发明涉及一种卫星遥感产品支持的单通道法地表温度反演方法，适用于具有单一热红外通道的卫星影像的地表温度遥感反演。

背景技术

[0002] 地表温度是进行区域和全球尺度地球表面物理过程研究的一个重要参数，是研究地气系统能量平衡、地气相互作用的基本物理量。在环境以及灾害监测中，地表温度也发挥着重要的作用，如对影响城市生活舒适度的热岛效应的监测，对影响农作物产量的旱灾的监测，对危害森林、草地的火灾监测等。因此，精确获取地表面温度具有重要的实际意义。随着定量遥感技术的发展，近年来使用具有热红外通道的卫星传感器反演地表温度的研究也得到越来越多的开展。

[0003] 在地表温度反演中，地表比辐射率的确定是一个关键的问题，其确定精度直接影响地表温度的反演精度。常用的地表比辐射率确定方法主要有以下两种，第一种是基于地表分类的方法。第二种是基于地表比辐射率和植被指数或植被覆盖度之间的统计关系模型，如Sobrino（2001年）提出的混合像元比辐射率估算方法。由于第二种方法获取数据容易，而且精度较高，被广泛用于地表温度反演时地表比辐射率的确定。但由于部分卫星传感器不能实时获取与热红外通道相对应的植被指数或植被覆盖信息，限制了地表温度反演工作的开展。如我国发射的HJ-1B IRS传感器与可以获取植被指数的HJ-1B CCD数据覆盖区域有较大的差异，若用HJ-1B CCD数据来确定HJ-1B IRS覆盖区域的植被指数，一方面不能实现HJ-1B IRS地表温度的实时反演，另一方面由于受两个数据上云量之和的影响，会严重降低反演的空间范围。

发明内容

[0004] 针对现有技术中部分热红外数据在运用单通道方法开展地表温度反演时，不能精确地确定地表比辐射率的技术问题，本发明提出了一种卫星遥感产品支持的单通道法地表温度反演方法，地表比辐射率确定精度较高，利于实现该类型传感器单通道法地表温度的反演。

[0005] 为了实现上述目的，本发明采用如下技术方案：

[0006] 卫星遥感产品支持的单通道法地表温度反演方法，包括如下步骤:

[0007] a、收集与待反演影像同时相的16天合成的MODIS地表植被指数产品NDVI；

[0008] b、将步骤a中获得的MODIS地表植被指数产品NDVI进行投影转换，转换到与IRS数据相统一的UTM投影下，并对MODIS地表植被指数产品NDVI进行重采样运算；

[0009] 将MODIS NDVI数值从大到小等间隔进行密度分割，并统计各个区间内的像元个数及其所占百分比；然后将IRS第1波段反射率值同样按照降序排列统计，以与MODIS NDVI区间范围占相同百分比为匹配条件对图像进行分割；

[0010] 根据MODIS影像的NDVI区间范围，对分割后的IRS影像赋予新的NDVI值；

[0011] c、基于NDVI与地表比辐射率的统计模型，计算地表比辐射率；

[0012] d、获取与IRS数据同时相的MODIS水汽通道数据，用三通道比值法来实现大气水汽含量的反演，确定HJ-1B IRS的水汽含量；

[0013] e、将步骤c和步骤d中得到的地表比辐射率ελ和水汽含量代入QK&B 算法或JM&S算法得到地表温度反演结果。

[0014] 上述步骤b中，具体赋值方法如下：

[0015] b1、根据土壤的光谱特征，当NDVI＜0.2，认为地表全部由裸土构成，对应IRS近红外波段NDVI取0.2；

[0016] b2、考虑当植被覆盖较高时，NDVI趋于饱和，因此当NDVI＞0.7，将对应IRS近红外波段NDVI统一取最大值0.8；

[0017] b3、当NDVI在0.2-0.7之间，认为地表由植被和裸土两种地表构成，以0.1的间隔等分，对应IRS近红外波段赋予各区间的上限值，NDVI依次取0.3-0.7。

[0018] 上步骤c中，在NDVI确定之后，植被和土壤的发射率εv、εs分别由ASTER光谱库提供的发射率曲线与IRS热红外波段光谱响应函数进行卷积运算获得：εv=0.9833，εs=0.9732；

[0019] 地表比辐射率ελ的计算公式为：ελ＝0.9732+0.0101*NDVI。

[0020] 本发明具有如下优点：

[0021] 本发明充分考虑到HJ-1B IRS与HJ-1B CCD数据的空间覆盖差异，提出用MODIS的植被指数产品协同确定地表比辐射率的方法，利用NDVI和近红外波段反射率的关联关系，通过直方图匹配的方法，完成了IRS数据比辐射率的提取。由于HJ-1B传感器自身不具有水汽反演能力，本发明采用MODIS水汽通道协助提取水汽数据，保证了水汽参数的高精度获取。本发明使用卫星遥感产品协同确定地表比辐射率以及水汽含量，从而提高了单通道法地表温度的反演精度，解决了部分卫星传感器不能实时获取与热红外通道相对应的植被指数或植被覆盖信息，限制了地表温度反演工作的开展的问题，能够更准确地实现地表温度反演，利于提升遥感产品在环境以及灾害监测中的应用水平。附图说明

[0022] 图1为本发明中卫星遥感产品支持的单通道法地表温度反演方法的流程示意图；

[0023] 图2为2008年12月31日统计的各个温度下像元个数的灰度直方图；

[0024] 图3为2009年1月4日统计的各个温度下像元个数的灰度直方图。

具体实施方式

[0025] 下面结合附图以及具体实施方式对本发明作进一步详细说明：

[0026] 结合图1所示，一种卫星遥感产品支持的单通道法地表温度反演方法，具体步骤如下：

[0027] 第一步、收集与待反演影像同时相的16天合成的MODIS地表植被指数产品NDVI。

[0028] 第二步、基于合成的MODIS地表植被指数产品对IRS热红外波段的NDVI指数赋值。

[0029] 首先，由于HJ-1B IRS与MODIS数据投影系统和空间分辨率不完全一致，IRS采用UTM投影，近红外和热红外波段分辨率分别为150m和300m，MODIS的NDVI产品和地表反射率数据采用Sinusoidal投影，分辨率分别为250m和1km，针对这一问题，本发明将各类数据在统一的UTM投影下，并对各类影像进行重采样运算；

[0030] 将MODIS NDVI数值从大到小等间隔进行密度分割，并统计各个区间内的像元个数及其所占百分比，然后将IRS第1波段反射率值同样按照降序排列统计，以与MODIS NDVI区间范围占相同百分比为匹配条件对图像进行分割；

[0031] 根据MODIS影像的NDVI区间范围，对分割后的IRS影像赋予新的NDVI值。具体赋值方法如下：

[0032] 1、根据土壤的光谱特征，当NDVI＜0.2，认为地表全部由裸土构成，对应IRS近红外波段NDVI取0.2；

[0033] 2、考虑当植被覆盖较高时，NDVI趋于饱和，因此当NDVI＞0.7，将对应IRS近红外波段NDVI统一取最大值0.8；

[0034] 3、当NDVI在0.2-0.7之间，认为地表由植被和裸土两种地表构成，以0.1的间隔等分，对应IRS近红外波段赋予各区间的上限值，NDVI依次取0.3-0.7；

[0035] 第三步、地表比辐射率确定

[0036] 在NDVI确定之后，植被和土壤的发射率εv、εs分别由ASTER光谱库提供的发射率曲线与IRS热红外波段光谱响应函数进行卷积运算获得：εv=0.9833，εs=0.9732；

[0037] 地表比辐射率通过式(1)进行估计：

[0038] ελ＝0.9732+0.0101*NDVI (1)

[0039] 通过使用土壤与植被的地物光谱数据模拟得出，上述方法确定的植被指数最大可导致地表比辐射率的误差为0.001。

[0040] 第四步、获取与IRS数据同时相的MODIS水汽通道数据，用三通道比值法来实现大气水汽含量的反演。

[0041] 因为MODIS与环境星的过境时间一致，所以可以利用MODIS的水汽反演结果支持环境星。强吸收波段16、17、18波段根据三通道比值法分别计算出各个水汽吸收波段的透过率，因为水汽强吸收波段，透过率与水汽含量是密切相关的，利用6S辐射传输模型模拟水汽含量与透过率之间的关系，得到透过率，进而得到水汽含量。然后三个强吸收波段反演得到的水汽含量加权获得较准确的水汽反演结果。

[0042] 第五步、地表温度反演实验

[0043] 将第三步得到的地表比辐射率ελ和第四步得到的水汽含量代入QK&B算法或JM&S算法得到地表温度反演结果。

[0044] 图2和图3分别示出了2008年12月31日、以及2009年1月4日统计的各个温度下像元个数的灰度直方图，由灰度直方图可知，地表温度反演结果总体趋势上与MODIS温度标准产品具有较高的一致性，两种算法反演结果与MODIS温度产品相比，大部分像元值误差都在2k之内，可知本发明得到的地表温度产品结果稳定。

[0045] 通过对地表温度反演实验结果的精度评价，表明本发明方法反演的地表温度能够达到较高的精度，具有很强的实用性。

[0046] 当然，以上说明仅仅为本发明的较佳实施例，本发明并不限于列举上述实施例，应当说明的是，任何熟悉本领域的技术人员在本说明书的教导下，所做出的所有等同替代、或明显变形形式，均落在本说明书的实质范围之内，理应受到本发明的保护。

标题	发布/更新时间	阅读量
北斗多模授时接收装置	2020-05-08	91
一种应用于智能终端RTD定位的粗差探测方法及系统	2020-05-08	840
一种北斗和GPS卫星信号接收的正直性监测和判决方法	2020-05-08	341
一种上行同步方法、终端、以及基站	2020-05-08	132
一种卫星系统的音视频传输控制方法、装置及服务器	2020-05-08	372
卫星天线对准装置	2020-05-08	595
一种基于GPS定位与自动避障的农业机器人的实现方法	2020-05-08	756
一种基于大数据的洪水预报方法	2020-05-08	24
一种基于遥感的灌溉水有效利用系数测算方法	2020-05-08	466
一种航向角确定方法、系统、电子设备及存储介质	2020-05-08	957

卫星遥感产品支持的单通道法地表温度反演方法

卫星遥感产品支持的单通道法地表温度反演方法

技术领域

背景技术

发明内容

具体实施方式

该功能需要专业版企业版VIP权限，您可以：