一种顾及时空特征的云覆盖像元地表温度重建方法技术

本发明专利技术公开了一种顾及时空特征的云覆盖像元地表温度重建方法，包括：1、对目标区域在时间段T

【技术实现步骤摘要】
一种顾及时空特征的云覆盖像元地表温度重建方法
本专利技术属于遥感
，具体涉及一种对遥感影像数据中云覆盖像元地表温度进行重建的方法。
技术介绍
地表温度(landsurfacetemperature，LST)是研究地表能量平衡、水文变化和气候变化的重要变量。随着航天遥感技术的发展，很多传感器可提供空间连续的地表温度产品，如美国的静止轨道环境气象卫星GOES系列(GeostationaryOperationalEnvironmentalStatellites,GOES)、极轨卫星EOS/MODIS，中国的风云2号气象卫星系列、静止气象卫星FY-4A/AGRI等。但现有的地表温度产品皆是对其遥感影像进行云检测和云覆盖像元标识，而未估算云覆盖像元地表温度。云覆盖像元地表温度信息缺失在一定程度上破坏了地表温度时空分布的全局性和连续性，严重阻碍和限制了地表温度产品的应用和发展。因此，如何估算地表温度产品中云覆盖像元地表温度，形成无缝的卫星地表温度产品是热红外遥感的前沿研究难题。目前，LST数据重建方法主要有基于空间相关性的空间插值重建和以时间本文档来自技高网...

【技术保护点】
1.一种顾及时空特征的云覆盖像元地表温度重建方法，其特征在于，包括：/nS1、对目标区域在时间段T

【技术特征摘要】
1.一种顾及时空特征的云覆盖像元地表温度重建方法，其特征在于，包括：
S1、对目标区域在时间段Tlast＝[dbegin,tbegin,dstop,tstop]内的LST图像序列进行像元判断，将云覆盖像元处的LST值设为无效，其他像元处的LST值设为有效；所述LST图像序列由静止气象卫星采集；Tlast的起始时刻为dbegin日tbegin时，终止时刻为dstop日tstop时；
根据目标区域的快速大气订正云图产品计算目标区域的NDVI数据，并对NDVI数据进行处理，得到时间段Tlast内目标区域的NDVI图像序列；
S2、遍历时间段T′last＝[dbegin+nt,tbegin,dstop-nt,tstop]内每一个LST值无效的像元，对其进行时间域LST重建，nt为时间窗口尺寸值；
S3、如目标区域在时间段T′last内仍有无效LST值，遍历每一个LST值无效的像元，对其进行空间域LST重建；
S4、如目标区域在时间段T′last内仍有无效LST值，跳转至步骤S2进行下一次时空重建，直到目标区域在时间段T′last内的LST值均有效；
S5、对步骤S4得到的目标区域在时间段T′last内的LST图像序列进行S-G滤波。


2.根据权利要求1所述的云覆盖像元地表温度重建方法，其特征在于，所述静止气象卫星为风云二号、风云四号静止气象卫星、美国的静止轨道环境气象卫星GOES系列之一。


3.根据权利要求1所述的云覆盖像元地表温度重建方法，其特征在于，对NDVI数据采用最大值合成法得到NDVI图像序列，相邻NDVI图像之间的时间间隔为一日、一旬或一月；
对所述NDVI图像序列进行降采样，使其空间分辨率与LST图像的空间分辨率一致。


4.根据权利要求1所述的云覆盖像元地表温度重建方法，其特征在于，所述步骤S2中时间域LST重建具体步骤为：
S21、设待重建(xp,yp)处像元在dp日tp时的LST值；根据LST图像序列获取(xp,yp)处像元在时间段Tlast内的LST值序列，截取其中LST值无效的子序列，获取子序列两端LST值有效的时刻，分别记为LST(xp,yp,d0,t0)和LST(xp,yp,dend,tend)；如果d0＝dend＝dp，按照步骤S22对(xp,yp)处像元在dp日tp时的LST值进行重建，否则结束(xp,yp)处的时间域LST重建，跳转至S21进行下一个LST无效像元的时间域重建；
S22、构建(xp,yp)处像元以dp...

【专利技术属性】
技术研发人员：王爱辉，杨英宝，潘鑫，章勇，胡解君德，朱晶晶，姜松秀，
申请(专利权)人：河海大学，
类型：发明
国别省市：江苏;32