一种基于STARFM融合技术的遥感墒情监测方法技术

本发明专利技术公开了一种基于STARFM融合技术的遥感墒情监测方法，该方法主要包括：步骤一、根据TM、ETM数据对HJ星CCD数据进行预处理；步骤二、利用大气辐射传输模型MODTRAN4实现HJ星的逐像元大气校正；步骤三、对MODIS数据进行预处理；步骤四、搜索相似像素，计算HJ星和MODIS数据的相似像素权重、转化系数；步骤五、基于STARFM融合算法进行预测。本发明专利技术通过STARFM融合方法综合利用HJ星和MODIS数据进行融合，获得时空较高分辨率的干旱监测结果，克服了单一卫星MODIS数据空间分辨率低和HJ星不能实现逐日观测的难题，使得反演的时间分辨率提高到1天，墒情监测的空间分辨率提高到30米，在农业、水利、环保与生态监测领域具有重要应用价值。

【技术实现步骤摘要】
一种基于STARFM融合技术的遥感墒情监测方法
本专利技术涉及一种基于STARFM融合技术的遥感墒情监测方法，具体涉及一种综合利用HJ星和MODIS数据获得较高时空分辨率的干旱监测方法，属于环境监测

技术介绍
干旱是主要的自然灾害之一，由于其影响范围大、持续时间长，并且后续影响久等特点使得干旱造成的经济损失远远超过了其他气象灾害，因此对干旱状况进行监测，及时了解旱情信息，掌握干旱动向，为相关部门提供科学合理的依据，对减少因干旱造成的损失具有十分重要的意义。传统的干旱监测方法主要有土壤湿度计法、称重法、电阻法等，虽然具有精度较高的优点，但是由于监测花费的人力物力大、监测范围有限、监测速度较慢，难以满足大范围的实时监测。随着科学技术的进步，遥感技术因其宽广的覆盖范围、较高的时空分辨率以及获取数据快捷等优点，在干旱的监测研究工作中作为大范围内信息获取的有效手段得到了广泛的应用，较传统的干旱监测方法有不可比拟的潜在优势。目前应用较为广泛的旱情遥感监测评估方法主要有：基于土壤热惯量的旱情监测方法、基于蒸散模型的旱情监测方法和基于植被指数的旱情监测方法。然而上述利用单一卫星遥感技术的监测干旱方法虽然取得了一定成功，但由于单一卫星本身的缺陷，如分辨率低，周期长等，无法完全满足干旱监测的需要。而基于多卫星遥感的监测利用不同卫星各自的特点，具有其优势：可以有目的对各变量参数采取遥感融合模型，提高反演的时空精度，可以增大观测范围、增强观测质量、提高观测频率，弥补了单个卫星监测造成区域旱情监测不准的缺陷。因此，基于多卫星融合进行干旱监测，将能很大程度上提高干旱监测的准确度。但是，目前基于多卫星融合进行干旱监测的方案并不多见，反演的时间精度低，墒情监测的空间分辨率精度不高，干旱监测结果并不理想。
技术实现思路
本专利技术所要解决的技术问题是针对前述
技术介绍
中的缺陷，综合HJ星和MODIS数据，提出一种耦合两者数据的方法进行干旱监测，能够实现较高时空分辨率的干旱监测结果，提高干旱监测水平。本专利技术为解决上述技术问题采用以下技术方案：一种基于STARFM融合技术的遥感墒情监测方法，通过STARFM融合方法综合利用HJ星和MODIS数据进行融合，得到干旱监测结果，具体包括以下步骤：步骤1、根据TM、ETM数据对HJ星CCD数据进行预处理；步骤2、利用大气辐射传输模型MODTRAN4实现HJ星的逐像元大气校正；步骤3、MODIS数据预处理；步骤4、搜索相似像素，计算HJ星和MODIS数据的相似像素权重、转化系数；步骤5、基于STARFM融合算法进行预测：基于STARFM融合方法获取同一时间的较高分辨率的HJ星CCD数据与较低分辨率MODIS数据，通过计算影像间空间分布的差异，结合另一时间较低分辨率的MODIS数据进行相应时间较高分辨率的HJ星CCD数据的预测。进一步的，本专利技术的基于STARFM融合技术的遥感墒情监测方法，所述步骤1具体包括步骤：(1.1)对TM、ETM、HJ星CCD数据进行镶嵌处理并裁剪；(1.2)根据裁剪出的TM、ETM底图作为基准图像校正HJ星CCD数据，控制校正误差；(1.3)对HJ星CCD数据进行辐射定标，并计算表观反射率。进一步的，本专利技术的基于STARFM融合技术的遥感墒情监测方法，所述步骤2具体包括：利用MODTRAN4模型模拟辐射传输过程及对应大气参数下的表观辐亮度，建立大气校正查找表，对HJ数据进行逐像元校正，其中辐射传输方程为：其中Lm表示接收的表观辐亮度，L0表示路径辐射既程辐射，T表示地表到传感器的透过率，s为大气底层球面反照率，ρ为像元反射率，Fd为到达地表的下行辐射通量密度。进一步详细介绍步骤2：(2.1)根据输入的大气参数假设三个像元反射率ρ，MODTRAN运行三次得到三个模拟的星上幅亮度Lm；根据得到的结果ρ和Lm，利用辐射传输方程公式建立三个以L0，s及TFd/π为变量的三元一次方程组，求解方程组得出一组与输入大气参数对应的L0，s及TFd/π。(2.2)按照上述步骤输入不同大气参数建立通用查找表。(2.3)读取HJ星的CCD数据，若是DN值将其图像转化为表观幅亮度图像；读取气溶胶光学厚度(AOT)数据的每个像元，根据每个像元的气溶胶光学厚度值，找到与该大气参数对应的三个参数：L0，s及TFd/π。(2.4)逐像元校正CCD或HSI数据。进一步的，本专利技术的基于STARFM融合技术的遥感墒情监测方法，其特征在于，所述步骤3中MODIS数据的预处理包括几何校正、图像裁剪、大气校正。进一步的，本专利技术的基于STARFM融合技术的遥感墒情监测方法，所述步骤4中搜索相似像素包括搜索HJ卫星和MODIS数据成像时间相同或相近的图像，计算HJ星和MODIS数据的相似像素权重函数如下：其中，Wijk为权重系数，Cijk是由窗口中心点的预测像元与窗口中其它像元的光谱距离、时间距离、空间距离计算的结果，W是窗口的大小，n代表输入已知图像对的数目。进一步的，本专利技术的基于STARFM融合技术的遥感墒情监测方法，所述步骤5基于STARFM融合算法进行预测，其计算公式如下：其中，L(xW/2,yW/2,T2)是预测的T2时刻的高分辨像元值，W是窗口的大小，n代表输入已知图像对的数目，Wijk为权重系数，(xW/2,yW/2)是窗口中间的像元，M(xi,yi,T2)是窗口位置(xi,yi)处在T2时刻的低分辨率像元值，L(xi,yi,T1)和M(xi,yi,T1)则分别是高分辨率蒸散发ET和低分辨率蒸散发ET在T1时刻的相应像元值。本专利技术采用以上技术方案与现有技术相比，具有以下技术效果：本专利技术利用STARFM融合方法通过每天的MODIS数据得到每天较高分辨率的HJ星数据，解决了单一卫星监测得到的MODIS数据空间分辨率较低以及HJ星不能每天观测的缺陷，将HJ星与MODIS相结合，得到时空较高分辨率的合成影像，使得反演的时间精度提高到一天，墒情监测的空间分辨率精度提高到30m，从而可以得到较高时空分辨率的干旱监测结果。附图说明图1是本专利技术的一种基于STARFM融合技术的遥感墒情监测方法的流程图。具体实施方式下面结合附图对本专利技术的技术方案做进一步的详细说明：本
技术人员可以理解的是，除非另外定义，这里使用的所有术语(包括技术术语和科学术语)具有与本专利技术所属领域中的普通技术人员的一般理解相同的意义。还应该理解的是，诸如通用字典中定义的那些术语应该被理解为具有与现有技术的上下文中的意义一致的意义，并且除非像这里一样定义，不会用理想化或过于正式的含义来解释。首先，介绍一下HJ星(环境减灾卫星)的基本概念：HJ星是国家2008年9月发射的用于环境与灾害监测预报的小卫星星座。主要有HJ-1-A星和HJ-1-B星两颗星座。HJ-1-A卫星和HJ-1-B卫星的轨道完全相同，相位相差180度。HJ-1-A星搭载了2台CCD相机和1台超光谱成像仪(HSI)，HJ-1-B星搭载2台CCD相机和1台红外相机(IRS)。两卫星上装载的两台CCD相机以星下点对称放置联合完成对地观测幅宽度700公里的推扫成像，具有30米的地面象元分辨率，两台CCD相机组网后重访周期仅为2天。此外，A卫星装载有一台超光谱成像仪，完成宽度为50公里对地本文档来自技高网...

【技术保护点】
一种基于STARFM融合技术的遥感墒情监测方法，其特征在于，通过STARFM融合方法综合利用HJ星和MODIS数据进行融合，得到干旱监测结果，具体包括以下步骤：步骤1、根据TM、ETM数据对HJ星CCD数据进行预处理；步骤2、利用大气辐射传输模型MODTRAN4实现HJ星的逐像元大气校正；步骤3、MODIS数据预处理；步骤4、搜索相似像素，计算HJ星和MODIS数据的相似像素权重、转化系数；步骤5、基于STARFM融合算法进行预测：基于STARFM融合方法获取同一时间的较高分辨率的HJ星CCD数据与较低分辨率MODIS数据，通过计算影像间空间分布的差异，结合另一时间较低分辨率的MODIS数据进行相应时间较高分辨率的HJ星CCD数据的预测。

【技术特征摘要】
1.一种基于STARFM融合技术的遥感墒情监测方法，其特征在于，通过STARFM融合方法综合利用HJ星和MODIS数据进行融合，得到干旱监测结果，具体包括以下步骤：步骤1、根据TM、ETM数据对HJ星CCD数据进行预处理；具体包括步骤：(1.1)对TM、ETM、HJ星CCD数据进行镶嵌处理并裁剪；(1.2)根据裁剪出的TM、ETM底图作为基准图像校正HJ星CCD数据，控制校正误差；(1.3)对HJ星CCD数据进行辐射定标，并计算表观反射率；步骤2、利用大气辐射传输模型MODTRAN4实现HJ星的逐像元大气校正，具体包括：利用MODTRAN4模型模拟辐射传输过程及对应大气参数下的表观辐亮度，建立大气校正查找表，对HJ数据进行逐像元校正，其中辐射传输方程为：其中Lm表示接收的表观辐亮度，L0表示路径辐射既程辐射，T表示地表到传感器的透过率，s为大气底层球面反照率，ρ为像元反射率，Fd为到达地表的下行辐射通量密度；步骤3、MODIS数据预处理；步骤4、搜索相似像素，计算HJ星和MODIS数据的相似像素权重、转化系数；步骤5、基于STARFM融合算法进行预测：基于STARFM融合方法获取同一时间的较高分辨率的HJ星CCD数据与较低分辨率MODIS数据，通过计算影像间空间分布的差异，结合另一时间较低分辨率的MODIS数据进行相应时间较高分辨率的HJ星CCD数据的预测。2.根据权利要求1所述的基于STARF...

【专利技术属性】
技术研发人员：杨涛，何祺胜，师鹏飞，崔同，李振亚，周旭东，王晓燕，
申请(专利权)人：河海大学，
类型：发明
国别省市：江苏;32