一种遥感地表温度数据的时空一体化融合方法技术

本发明专利技术涉及定量遥感数据处理技术领域，公开了一种遥感地表温度数据的时空一体化融合方法，首先对遥感地表温度数据进行预处理，其次，利用多时相温度数据和对应的时间求解地表温度昼夜变化模型参数，实现空间无缝，最后在极轨卫星地表温度影像数据上逐窗口筛选相似像元，结合由多时相静止卫星与上述模型构成的时间差异、尺度差异、相似度差异以及欧氏距离等权重因子，共同构建时空一体化融合模型并求解生成空间无缝的高时空分辨率地表温度数据集。本发明专利技术充分利用极轨卫星与多时相静止卫星地表温度数据的时空互补信息和静止卫星地表温度的昼夜变化规律，不需要任何辅助数据，易于实现，实用价值高，可用于时空无缝的遥感地表温度数据业务化运行。

【技术实现步骤摘要】

本专利技术涉及定量遥感数据处理
，尤其是涉及一种遥感地表温度数据的时 空一体化融合方法。
技术介绍
地表温度（LST，Land Surface Temperature)是研究地-气间物质与能量交换的 重要参数，广泛应用于城市热岛效应、土壤湿度估算、森林火灾监测、气候变化等领域。卫星 热红外遥感是获取地表温度的重要手段，然而，遥感观测系统需要在空间分辨率和时间分 辨率之间进行权衡与取舍，难以同时获得高时间和高空间分辨率的温度数据，极大限制了 其应用的广度与深度。例如：静止卫星地表温度数据的时间分辨率最高达到15分钟，却具 有至少3千米的低空间分辨率；而极地卫星提供的地表温度数据虽然具有较高的空间分辨 率（约100米-1000米），但其时间分辨率一般都较低。另一方面，现有发布的地表温度数 据由于云、气溶胶等不利天气的影响往往出现空值现象，导致空间不连续。因此，如何在不 改变现有观测系统的条件下，实现遥感地表温度数据空间无缝的同时提高遥感LST的空间 分辨率和时间分辨率，成为热红外遥感应用中的关键问题之一。 为了提高遥感LST数据的时空分辨率，国内外学者已经开展了一些研究。早期的 研究主要有遥感LST空间降尺度/时间扩展：在空间降尺度方面，李小文院士提出将传统地 学趋势面用于遥感尺度转换研究，其主要原理都是引入更高空间分辨率的地理要素构造趋 势面，通过建立趋势面与LST间的关系完成降尺度转换。而对于遥感LST时间尺度扩展，一 般是对多时相数据进行时相滤波、谐波分析等方法重建出中间时间的数据。以上针对单一 遥感参量来提高时空分辨率的方法，对辅助数据质量或长时序性要求较高，且提高的分辨 率程度有限，在某些应用也显得无能为力。因此，需发展一种新的LST时空降尺度方法。 多传感器遥感LST时空融合技术正是为了满足这种需求而发展起来的一种新技 术。与上述针对单一遥感参量的空间降尺度、时间扩展方法不同的是，遥感LST时空融合是 通过对具有不同空间尺度和时间尺度的LST进行融合处理，生成时间相对密集的高空间分 辨率LST数据；其优势在于无需其它辅助数据条件，同时生成的数据具有输入数据中最高 空间分辨率和最高时间分辨率。 然而，无论是传统时空降尺度还是时空融合，在算法实施过程中都假设所输入数 据空间无缝，这个假设对于地表温度数据（由于受云等影响）经常不能满足；此外，在基 于多时相、多尺度温度数据时空融合时，很少顾及参量昼夜变化规律、遥感观测条件等的差 另IJ，直接影响了融合影像的求解精度。因此如何利用多时相、多尺度时空互补信息实现遥感 地表温度数据空间无缝同时提高其时空分辨率，需要研究新的时空一体化融合方法。
技术实现思路
本专利技术所要解决的问题是：针对现有观测系统的技术限制和遥感地表温度数据时 空不连续现象，提供。 为解决上述技术问题，本专利技术提出的技术方案为：一种遥感地表温度数据的时空 一体化融合方法，设有一个空间无缝的极轨卫星地表温度数据M和同一地区的多时相静止 卫星地表温度数据Gtp(p = 0, 1，"·η且n>10);进一步假设在时间序列tQ、1v"tn Mn上有 一景地表温度数据G与极轨卫星地表温度数据M获取时间上对应，则实现时空一体化融合 包括以下步骤： 步骤一、对地表温度数据M和Gtp(p = 0, 1，…η)进行预处理，包括几何配准、裁剪、 温度单位转换、空值标记、重采样； 步骤二、基于昼夜变化模型的空间修复：根据当前时相待修复像元在其它时相对 应有效像元值的数量分如下两种情况求解所述昼夜变化模型的参数，并进行空值修复； 情况1)若当前时相的待修复像元在其它时相上有足够非空值像元，直接利用这 些有效值与对应时间构成的数据对按需求分白天和晚上进行所述昼夜变化模型参数拟 合； 情况2)若待修复像元的其它时相有效像元及其对应时间构成的数据对不足以 进行温度昼夜变化模型参数拟合时，则以相邻时相1对应的非空值像元为中心像元开一 个窗口，设窗口大小为waX wa，wa根据G的特征而定，在窗口内筛选相似像元，筛选条件 为：邻域像元与中心像元差的绝对值小于或等于该中心像元相关的函数，用公式表达为：，其中G(i，tx)为时间^时静止卫星地表温度数 据G上窗口内邻域像元，G(wa/2, tx)为静止卫星地表温度数据G上窗口内中心像元，（13为 常量，使得筛选的相似像元满足其它时相上有足够多的有效数据对；利用这些有效数据对 进行所述昼夜变化模型参数的拟合数据数据。 步骤三、构建遥感地表温度数据时空一体化融合模型： 在地表温度数据M上逐个以像元为中心像元开一个大小为wbXwb的窗口，窗口大 小视M特征而定，按以下公式在窗口内筛选相似像元： 其中M(i, t)为时间t时M上窗口内邻域像元，M(wb/2, t)为M上窗口内中心像元， db为常量， 计算相似像元地表温度数据在当前时间和预测时间间的温度变化值与中心像元 所在昼夜变化模型的温度变化值的差异TD1如下， TDi= IG (i, t) -G (i, t x) - (Temp (wb/2, t) -Temp (wb/2, tx)) 其中，G (i，t)为当前时间t的静止卫星地表温度数据G上窗口内第i个邻域像元， G (i，tx)为时间tx时静止卫星地表温度数据G上窗口内第i个邻域像元，i为正整数，其最 大值为窗口大小； 结合TD1及其它参量构建以下相似像元的权重函数表达式： W1= (!/(V1 · SD1)/ Σ !/(V1 · SD1)) 其中SD1是窗口内相似像元与中心像元的相似度、V1是由M上窗口内相似像元与 中心像元温度变化差异、欧式距离以及不同分辨率地表温度数据间尺度差异共同构成；假 设地表温度数据M和Gtn在时间t n时对应，作为基础时间对，其它时间序列t xe {t。、心、… tn J作为预测时间，得到时空一体化融合模型： CN 105184076 A 说明书 3/6 页 其中，M(wb/2, tx)为融合的空间无缝且高时空分辨率的地表温度数据，具有静止 卫星G的时间分辨率和极轨卫星M的空间分辨率。 进一步的，在所述步骤二中，由于考虑到观测条件、传感器误差、反演误差的原因 使得现有地表温度数据存在包含一些噪声、误差较大的离群点，采用M估计理论提高模型 参数拟合精度。 所述步骤二中，当待修复像元的其它时相像元都为空值时，选择当前时相邻域有 效像元值，包括原有的和建模求解的；利用空间内插技术完成修复。 本专利技术充分结合多尺度、多时相传感器间的时空互补信息和静止卫星地表温度的 昼夜变化规律，利用静止卫星高时间分辨率特性建立地表温度昼夜变化模型，进而实现空 值像元的修复，同时与较高分辨率极轨卫星温度数据构建时空融合模型，最终生成空间无 缝的高时空分辨率地表温度数据集。而且不需要任何辅助数据，易于实现，实用价值高，可 用于业务化运行。本专利技术与现有技术相比具有空间修复精确、时空权重函数合理、时空一体 化融合等优点。主要体现在：首先，根据待修复空值像元在时间维上的有效像元数分情况求 解地表温度昼夜变化模型参数，并利用M估计等稳健回归方法使模型参数最优化，从而提 高修复精度；其次，通过引入相似像元本身的相似度差异因子、基于温度昼夜变化模型的时本文档来自技高网...

【技术保护点】
一种遥感地表温度数据的时空一体化融合方法，其特征在于：设有一个空间无缝的极轨卫星地表温度数据M和同一地区的多时相静止卫星地表温度数据Gtp(p＝0,1,…n且n>10)；进一步假设在时间序列t0、t1…tn‑1、tn上有一景地表温度数据G与极轨卫星地表温度数据M获取时间上对应，则实现时空一体化融合包括以下步骤:步骤一、对地表温度数据M和Gtp(p＝0,1,…n)进行预处理，包括几何配准、裁剪、温度单位转换、空值标记、重采样；步骤二、基于昼夜变化模型的空间修复：根据当前时相待修复像元在其它时相对应有效像元值的数量分如下两种情况求解所述昼夜变化模型的参数，并进行空值修复；情况1)若当前时相的待修复像元在其它时相上有足够非空值像元，直接利用这些有效值与对应时间构成的数据对按需求分白天和晚上进行所述昼夜变化模型参数拟合；情况2)若待修复像元的其它时相有效像元及其对应时间构成的数据对不足以进行温度昼夜变化模型参数拟合时，则以相邻时相tx对应的非空值像元为中心像元开一个窗口，设窗口大小为wa×wa，wa根据G的特征而定，在窗口内筛选相似像元，筛选条件为：邻域像元与中心像元差的绝对值小于或等于该中心像元相关的函数，用公式表达为：其中G(i,tx)为时间tx时静止卫星地表温度数据G上窗口内邻域像元，G(wa/2,tx)为静止卫星地表温度数据G上窗口内中心像元，da为常量，使得筛选的相似像元满足其它时相上有足够多的有效数据对；利用这些有效数据对进行所述昼夜变化模型参数的拟合数据数据。步骤三、构建遥感地表温度数据时空一体化融合模型：在地表温度数据M上逐个以像元为中心像元开一个大小为wb×wb的窗口，窗口大小视M特征而定，按以下公式在窗口内筛选相似像元:|M(i,t)‑M(wb/2,t)|≤db×2M(wb/2,t)其中M(i,t)为时间t时M上窗口内邻域像元，M(wb/2,t)为M上窗口内中心像元，db为常量，计算相似像元地表温度数据在当前时间和预测时间间的温度变化值与中心像元所在昼夜变化模型的温度变化值的差异TDi如下，TDi＝|G(i,t)‑G(i,tx)‑(Temp(wb/2,t)‑Temp(wb/2,tx))|其中，G(i,t)为当前时间t的静止卫星地表温度数据G上窗口内第i个邻域像元，G(i,tx)为时间tx时静止卫星地表温度数据G上窗口内第i个邻域像元，i为正整数，其最大值为窗口大小；结合TDi及其它参量构建以下相似像元的权重函数表达式：Wi＝(1/(Vi·SDi)/∑1/(Vi·SDi))其中SDi是窗口内相似像元与中心像元的相似度、Vi是由M上窗口内相似像元与中心像元温度变化差异、欧式距离以及不同分辨率地表温度数据间尺度差异共同构成；假设地表温度数据M和Gtn在时间tn时对应，作为基础时间对，其它时间序列tx∈{t0、t1、…tn‑1}作为预测时间，得到时空一体化融合模型：M(wb/2,tx)=Σi=1wWi×(M(i,tn)+G(i,tn)-G(i,tx))]]>其中，M(wb/2,tx)为融合的空间无缝且高时空分辨率的地表温度数据，具有静止卫星G的时间分辨率和极轨卫星M的空间分辨率。...

【技术特征摘要】

【专利技术属性】
技术研发人员：吴鹏海，沈焕锋，刘紫涵，曾超，邹照旭，
申请(专利权)人：安徽大学，
类型：发明
国别省市：安徽;34