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

本发明专利技术公开了一种卫星遥感产品支持的单通道法地表温度反演方法，包括如下步骤：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算法得到地表温度反演结果。本发明专利技术中地表比辐射率确定精度较高，利于实现该类型传感器单通道法地表温度的反演。

【技术实现步骤摘要】
卫星遥感产品支持的单通道法地表温度反演方法
本专利技术涉及一种卫星遥感产品支持的单通道法地表温度反演方法，适用于具有单一热红外通道的卫星影像的地表温度遥感反演。
技术介绍
地表温度是进行区域和全球尺度地球表面物理过程研究的一个重要参数，是研究地气系统能量平衡、地气相互作用的基本物理量。在环境以及灾害监测中，地表温度也发挥着重要的作用，如对影响城市生活舒适度的热岛效应的监测，对影响农作物产量的旱灾的监测，对危害森林、草地的火灾监测等。因此，精确获取地表面温度具有重要的实际意义。随着定量遥感技术的发展，近年来使用具有热红外通道的卫星传感器反演地表温度的研究也得到越来越多的开展。在地表温度反演中，地表比辐射率的确定是一个关键的问题，其确定精度直接影响地表温度的反演精度。常用的地表比辐射率确定方法主要有以下两种，第一种是基于地表分类的方法。第二种是基于地表比辐射率和植被指数或植被覆盖度之间的统计关系模型，如Sobrino(2001年)提出的混合像元比辐射率估算方法。由于第二种方法获取数据容易，而且精度较高，被广泛用于地表温度反演时地表比辐射率的确定。但由于部分卫星传感器不能实时获取与热红外通道相对应的植被指数或植被覆盖信息，限制了地表温度反演工作的开展。如我国发射的HJ-1BIRS传感器与可以获取植被指数的HJ-1BCCD数据覆盖区域有较大的差异，若用HJ-1BCCD数据来确定HJ-1BIRS覆盖区域的植被指数，一方面不能实现HJ-1BIRS地表温度的实时反演，另一方面由于受两个数据上云量之和的影响，会严重降低反演的空间范围。
技术实现思路
针对现有技术中部分热红外数据在运用单通道方法开展地表温度反演时，不能精确地确定地表比辐射率的技术问题，本专利技术提出了一种卫星遥感产品支持的单通道法地表温度反演方法，地表比辐射率确定精度较高，利于实现该类型传感器单通道法地表温度的反演。为了实现上述目的，本专利技术采用如下技术方案：卫星遥感产品支持的单通道法地表温度反演方法，包括如下步骤:a、收集与待反演影像同时相的16天合成的MODIS地表植被指数产品NDVI；b、将步骤a中获得的MODIS地表植被指数产品NDVI进行投影转换，转换到与IRS数据相统一的UTM投影下，并对MODIS地表植被指数产品NDVI进行重采样运算；将MODISNDVI数值从大到小等间隔进行密度分割，并统计各个区间内的像元个数及其所占百分比；然后将IRS第1波段反射率值同样按照降序排列统计，以与MODISNDVI区间范围占相同百分比为匹配条件对图像进行分割；根据MODIS影像的NDVI区间范围，对分割后的IRS影像赋予新的NDVI值；c、基于NDVI与地表比辐射率的统计模型，计算地表比辐射率；d、获取与IRS数据同时相的MODIS水汽通道数据，用三通道比值法来实现大气水汽含量的反演，确定HJ-1BIRS的大气水汽含量；e、将步骤c和步骤d中得到的地表比辐射率ελ和大气水汽含量代入QK&B算法或JM&S算法得到地表温度反演结果。上述步骤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。上步骤c中，在NDVI确定之后，植被和土壤的发射率εv、εs分别由ASTER光谱库提供的发射率曲线与IRS热红外波段光谱响应函数进行卷积运算获得：εv＝0.9833，εs＝0.9732；地表比辐射率ελ的计算公式为：ελ＝0.9732+0.0101*NDVI。本专利技术具有如下优点：本专利技术充分考虑到HJ-1BIRS与HJ-1BCCD数据的空间覆盖差异，提出用MODIS的植被指数产品协同确定地表比辐射率的方法，利用NDVI和近红外波段反射率的关联关系，通过直方图匹配的方法，完成了IRS数据比辐射率的提取。由于HJ-1B传感器自身不具有水汽反演能力，本专利技术采用MODIS水汽通道协助提取水汽数据，保证了水汽参数的高精度获取。本专利技术使用卫星遥感产品协同确定地表比辐射率以及大气水汽含量，从而提高了单通道法地表温度的反演精度，解决了部分卫星传感器不能实时获取与热红外通道相对应的植被指数或植被覆盖信息，限制了地表温度反演工作的开展的问题，能够更准确地实现地表温度反演，利于提升遥感产品在环境以及灾害监测中的应用水平。附图说明图1为本专利技术中卫星遥感产品支持的单通道法地表温度反演方法的流程示意图；图2为2008年12月31日统计的各个温度下像元个数的灰度直方图；图3为2009年1月4日统计的各个温度下像元个数的灰度直方图。具体实施方式下面结合附图以及具体实施方式对本专利技术作进一步详细说明：结合图1所示，一种卫星遥感产品支持的单通道法地表温度反演方法，具体步骤如下：第一步、收集与待反演影像同时相的16天合成的MODIS地表植被指数产品NDVI。第二步、基于合成的MODIS地表植被指数产品对IRS热红外波段的NDVI指数赋值。首先，由于HJ-1BIRS与MODIS数据投影系统和空间分辨率不完全一致，IRS采用UTM投影，近红外和热红外波段分辨率分别为150m和300m，MODIS的NDVI产品和地表反射率数据采用Sinusoidal投影，分辨率分别为250m和1km，针对这一问题，本专利技术将各类数据在统一的UTM投影下，并对各类影像进行重采样运算；将MODISNDVI数值从大到小等间隔进行密度分割，并统计各个区间内的像元个数及其所占百分比，然后将IRS第1波段反射率值同样按照降序排列统计，以与MODISNDVI区间范围占相同百分比为匹配条件对图像进行分割；根据MODIS影像的NDVI区间范围，对分割后的IRS影像赋予新的NDVI值。具体赋值方法如下：1、根据土壤的光谱特征，当NDVI＜0.2，认为地表全部由裸土构成，对应IRS近红外波段NDVI取0.2；2、考虑当植被覆盖较高时，NDVI趋于饱和，因此当NDVI＞0.7，将对应IRS近红外波段NDVI统一取最大值0.8；3、当NDVI在0.2-0.7之间，认为地表由植被和裸土两种地表构成，以0.1的间隔等分，对应IRS近红外波段赋予各区间的上限值，NDVI依次取0.3-0.7；第三步、地表比辐射率确定在NDVI确定之后，植被和土壤的发射率εv、εs分别由ASTER光谱库提供的发射率曲线与IRS热红外波段光谱响应函数进行卷积运算获得：εv＝0.9833，εs＝0.9732；地表比辐射率通过式(1)进行估计：ελ＝0.9732+0.0101*NDVI(1)通过使用土壤与植被的地物光谱数据模拟得出，上述方法确定的植被指数最大可导致地表比辐射率的误差为0.001。第四步、获取与IRS数据同时相的MODIS水汽通道数据，用三通道比值法来实现大气水汽含量的反演。因为MODIS与环境星的过境时间一致，所以可以利用MODIS的水汽反演结果支持环境星。强吸收波段16、17、18波段根据三通道比本文档来自技高网...

【技术保护点】
卫星遥感产品支持的单通道法地表温度反演方法，其特征在于包括如下步骤: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算法得到地表温度反演结果。

【技术特征摘要】
1.卫星遥感产品支持的单通道法地表温度反演方法，其特征在于包括如下步骤:a、收集与待反演影像同时相的16天合成的MODIS地表植被指数产品NDVI；b、将步骤a中获得的MODIS地表植被指数产品NDVI进行投影转换，转换到与IRS数据相统一的UTM投影下，并对MODIS地表植被指数产品NDVI进行重采样运算；将MODISNDVI数值从大到小等间隔进行密度分割，并统计各个区间内的像元个数及其所占百分比；然后将IRS第1波段反射率值同样按照降序排列统计，以与MODISNDVI区间范围占相同百分比为匹配条件对图像进行分割；根据MODIS影像的NDVI区间范围，对分割后的IRS影像赋予新的NDVI值；c、基于NDVI与地表比辐射率的统计模型，计算地表比辐射率；d、获取与IRS数据同时相的MODIS水汽通道数据，用三通道比值法来实现大气水汽含量的反演，确定HJ-1BIRS的大气水汽含量；e、将步骤c和步骤d中得到的地表比辐射率ελ和大气水汽含量代入QK&B算法或JM&...

【专利技术属性】
技术研发人员：孙林，于会泳，高涛，
申请(专利权)人：山东科技大学，
类型：发明
国别省市：山东;37