一种时序MODIS定量遥感产品时空修复及批处理实现方法技术

本发明专利技术涉及定量遥感产品处理技术领域，公开一种时序MODIS定量遥感产品时空修复及其批处理实现方法，通过将时序MODIS定量遥感产品时空修复流程化，实现MODIS定量遥感产品空值或低质量数据修复的批量处理，得到空间上连续的高质量MODIS定量遥感产品；而且批处理计算效率较高，易于实现，实用价值高，可用于MODIS定量遥感产品质量改善的业务化运行。

【技术实现步骤摘要】
一种时序MODIS定量遥感产品时空修复及批处理实现方法
本专利技术涉及定量遥感产品处理
，具体涉及一种时序MODIS定量遥感产品时空修复及批处理实现方法。
技术介绍
MODIS定量遥感产品在全球对地观测中发挥着重要作用，但由于云、气溶胶等不利天气影响，现有发布的MODIS产品常出现空值或低质量等无效数值，极大限制了该产品正常使用与长时序应用。因此，在已有MODIS产品的基础上，利用技术手段对无效数据进行必要的修复成为一个十分重要的课题。常用的遥感产品中无效数据修复方法可分为基于单景数据修复和基于多景数据修复两种。其中，基于单景数据修复方法有插值、矩匹配、滤波和MAP估计等方法；基于多景数据修复方法通常是在不同时相或不同传感器影像之间建立关系，利用其它景数据上未受到影响的信息来进行修复，常用的有直方图匹配、回归分析和谐波分析等方法。单景数据的修复方法在无效数据较多时难以得到较好的修复，必须引入另外数据的信息；但多景数据修复方法面临选择哪些辅助数据参与修复的问题，且不同数据中的同一地物会受到其本身的光谱/辐射特性、尺度差异、观测时间和大气因素等影响，进而给产品修复带来很大的干扰。此外，上述两种方法都没有充分利用MODIS产品的时序特点，不具备面向具体应用的批量自动修复能力。
技术实现思路
(一)解决的技术问题针对现有技术的不足，本专利技术提供了一种时序MODIS定量遥感产品时空修复及批处理实现方法，用于解决现有发布的MODIS产品常出现空值或低质量等无效数值，极大限制了该产品正常使用与长时序应用等问题；本专利技术通过以下技术方案予以实现：(二)技术方案本专利技术提供了一种时序MODIS定量遥感产品时空修复及批处理实现方法，MODIS定量遥感产品由于云污染等原因出现空值或低质量等无效数据，通过在无效数据与时序中筛选的有效数据之间建立联系，实现无效数据的修复，包含①：MODIS数据预处理、②：最佳辅助数据的筛选与③：MODIS产品空值修复三个步骤：①MODIS数据预处理：对多时相MODIS定量遥感产品进行必要的预处理，所述预处理包括所述目标数据提取与无效像元标记；②最佳辅助数据的筛选：时序数据中寻找最佳辅助数据，依次遍历所有待修复的数据，在时序中找到与待修复时刻T0数据具有互补观测信息的所有MODIS数据作为候选辅助数据序列，根据待修复数据空值情况及其在时序中的位置，共有多种筛选所述最佳辅助数据的情况；③MODIS产品空值修复：基于时空互补信息的MODIS定量遥感产品修复，即把一景含无效数据的MODIS产品作为目标数据，根据MODIS空值区域在不同时相有所差别的特点，利用时序具有对应观测信息的辅助数据对所述目标数据中的无效区域进行修复，根据最佳辅助数据的筛选这一步骤中匹配的待修复数据及其所述最佳辅助数据的数量，共有多种执行产品修复的情况。优选的，在最佳辅助数据的筛选这一步骤中，筛选所述最佳辅助数据的情况共有五种，其分别为：情况1：若待修复时刻T0数据为时序中第一个数据或最后一个数据，且能在时序中找到与T0时刻数据有效像元重叠数足够多的数据，则筛选与T0时刻有效像元平均量化值差异最小的一景MODIS数据作为所述最佳辅助数据；情况2：若待修复时刻T0数据不是时序中第一个数据和最后一个数据，在以T0时刻为分界线的两组所述候选辅助数据序列中，能找到与T0时刻数据有效像元重叠数足够多的数据，则选择与T0时刻有效像元平均量化值差异最小的两景MODIS数据作为候选辅助数据；若这两个平均量化值差异相差较大，则用平均量化值差异小的一景候选辅助数据作为所述最佳辅助数据；情况3：若待修复时刻T0数据不是时序中第一个数据和最后一个数据，在以T0时刻为分界线的两组所述候选辅助数据序列中，能找到与T0时刻数据有效像元重叠数足够多的数据，则选择与T0时刻有效像元平均量化值差异最小的两景MODIS数据作为候选辅助数据；若这两个平均量化值差异相差较小，则这两景候选辅助数据都作为所述最佳辅助数据；情况4：若待修复时刻T0数据不是时序中第一个数据和最后一个数据，但在以T0时刻为分界线的两组所述候选辅助数据序列中，只能在一组候选辅助数据中找到与T0时刻数据有效像元重叠数足够多的数据，则筛选与T0时刻有效像元平均量化值差异最小的一景MODIS数据作为所述最佳辅助数据；情况5：若所述候选辅助数据序列中找到的与T0时刻数据有效像元重叠数较少，无法建立拟合关系，则暂不处理，直接进入MODIS产品空值修复这一步骤，选择T0时刻前后修复后的数据作为所述最佳辅助数据。优选的，在最佳辅助数据的筛选这一步骤中，当存在多个所述最佳辅助数据时，在时序中找到与待修复时刻数据具有互补观测信息的所有MODIS数据，并以T0时刻为分界线得到两组数据序列，在这两组所述序列数据中分别进一步筛选与时刻数据有效像元重叠数足够多的数据，则选择平均温度差异最小的两景MODIS数据作为最优辅助数据。优选的，在所述最佳辅助数据的筛选这一步骤中，由于空值区域的G和B不能直接求出，只能通过周围的有效像素点计算，在建立相似像素与目标像素关系时，考虑到并不是每个相似像素与目标像素的关系都是一样的，因此要根据每个相似像素与目标像素之间的紧密程度计算一个权重。优选的，执行产品修复的情况共有三种情况，其分别为：情况1：若是修复前的一景辅助数据修复，以待修复数据T0和辅助数据Tx为例，设T0上任一需修复的像元为DN_0，辅助数据TX上对应的像元为DN_X；以DN_X为中心，开一个大小为Nw×Nw的窗口，统计窗口内的所有有效像元，计算其标准差S_C；以S_C为阈值，在窗口内搜索DN_X的相似像元，当找到的相似像元数目少于规定的最小数目H时扩大窗口，直到满足条件或达到窗口大小为止；如果达到窗口上限时相似像元仍少于H，取消约束条件，将所有窗口内的点纳入计算；则T0时刻待修复数据中的无效像素修复后的值通过TX时刻辅助影像上对应区域的像素值计算得到恢复后的像素值：用公式表达为：DN_0＝G*DN_X+BG为两幅数据间的增益量，B为偏移量，可通过两幅数据间相似像元加权回归求解：其中，DN_0i和DN_Xi分别表示T0时刻待修复影像和TX时刻辅助影像中参与计算的第i个相似像素，和表示T0时刻待修复影像和TX时刻辅助影像中参与计算的N个相似像素的平均值，wi表示第i个相似像素的权重，由TX时刻辅助影像中相似像元(χi,yi)与中心像元(χt,yt)的温度差异与空间距离决定：Di＝|DN_Xi-DN_X|*((χi-χt)2+(yi-yt)2)情况2：若是修复前的两景最优辅助数据修复：DN_X和DN_Y，每景辅助数据可按照本阶段的情况1各得到一个修复结果，对这两个结果进行时相加权得到最终结果，其时相加权函数为：则修复结果公式为：DN_0＝DN_X*Wx+DN_Y*Wy情况3：若是修复后的两景最优辅助数据修复，当本阶段的情况1和情况2执行完成后，利用T0时刻前后修复好的两景数据TP和Tq作为所述最佳辅助数据时，根据反时间间隔加权计算修复值：DN_0＝DN_P*(Tq-T0)/(Tq-Tp)+DN_Q*(T0-Tp)/(Tq-Tp)其中，DN_P和DN_Q为修复好的辅助数据上与待修复影像上DN_0对应的像元。优选的，在所述MODI本文档来自技高网...

【技术保护点】
1.一种时序MODIS定量遥感产品时空修复及批处理实现方法，其特征在于：通过在无效数据与时序中筛选的有效数据之间建立联系，实现无效数据的修复，包含①：MODIS数据预处理、②：最佳辅助数据的筛选与③：MODIS产品空值修复三个步骤：①MODIS数据预处理：对多时相MODIS定量遥感产品进行必要的预处理，所述预处理包括所述目标数据提取与无效像元标记；②最佳辅助数据的筛选：时序数据中寻找最佳辅助数据，依次遍历所有待修复的数据，在时序中找到与待修复时刻T0数据具有互补观测信息的所有MODIS数据作为候选辅助数据序列，根据待修复数据空值情况及其在时序中的位置，共有多种筛选所述最佳辅助数据的情况；③MODIS产品空值修复：基于时空互补信息的MODIS定量遥感产品修复，即把一景含无效数据的MODIS产品作为目标数据，根据MODIS空值区域在不同时相有所差别的特点，利用时序具有对应观测信息的辅助数据对所述目标数据中的无效区域进行修复，根据最佳辅助数据的筛选这一步骤中匹配的待修复数据及其所述最佳辅助数据的数量，共有多种执行产品修复的情况。

【技术特征摘要】
1.一种时序MODIS定量遥感产品时空修复及批处理实现方法，其特征在于：通过在无效数据与时序中筛选的有效数据之间建立联系，实现无效数据的修复，包含①：MODIS数据预处理、②：最佳辅助数据的筛选与③：MODIS产品空值修复三个步骤：①MODIS数据预处理：对多时相MODIS定量遥感产品进行必要的预处理，所述预处理包括所述目标数据提取与无效像元标记；②最佳辅助数据的筛选：时序数据中寻找最佳辅助数据，依次遍历所有待修复的数据，在时序中找到与待修复时刻T0数据具有互补观测信息的所有MODIS数据作为候选辅助数据序列，根据待修复数据空值情况及其在时序中的位置，共有多种筛选所述最佳辅助数据的情况；③MODIS产品空值修复：基于时空互补信息的MODIS定量遥感产品修复，即把一景含无效数据的MODIS产品作为目标数据，根据MODIS空值区域在不同时相有所差别的特点，利用时序具有对应观测信息的辅助数据对所述目标数据中的无效区域进行修复，根据最佳辅助数据的筛选这一步骤中匹配的待修复数据及其所述最佳辅助数据的数量，共有多种执行产品修复的情况。2.根据权利要求1所述的一种时序MODIS定量遥感产品时空修复及批处理实现方法，其特征在于：在最佳辅助数据的筛选这一步骤中，筛选所述最佳辅助数据的情况共有五种，其分别为：情况1：若待修复时刻T0数据为时序中第一个数据或最后一个数据，且能在时序中找到与T0时刻数据有效像元重叠数足够多的数据，则筛选与T0时刻有效像元平均量化值差异最小的一景MODIS数据作为所述最佳辅助数据；情况2：若待修复时刻T0数据不是时序中第一个数据和最后一个数据，在以T0时刻为分界线的两组所述候选辅助数据序列中，分别能找到与T0时刻数据有效像元重叠数足够多的数据，则选择与T0时刻有效像元平均量化值差异最小的两景MODIS数据作为候选辅助数据；若这两个平均量化值差异相差较大，则用平均量化值差异小的一景候选辅助数据作为所述最佳辅助数据；情况3：若待修复时刻T0数据不是时序中第一个数据和最后一个数据，在以T0时刻为分界线的两组所述候选辅助数据序列中，分别能找到与T0时刻数据有效像元重叠数足够多的数据，则选择与T0时刻有效像元的平均量化值差异最小的两景MODIS数据作为候选辅助数据；若这两个平均量化值差异相差较小，则这两景候选辅助数据都作为所述最佳辅助数据；情况4：若待修复时刻T0数据不是时序中第一个数据和最后一个数据，但在以T0时刻为分界线的两组所述候选辅助数据序列中，只能在一组候选辅助数据中找到与T0时刻数据有效像元重叠数足够多的数据，则筛选与T0时刻有效像元平均量化值差异最小的一景MODIS数据作为所述最佳辅助数据；情况5：若所述候选辅助数据序列中找到的与T0时刻数据有效像元重叠数较少，无法建立拟合关系，则暂不处理，直接进入MODIS产品空值修复这一步骤，选择T0时刻前后修复后的数据作为所述最佳辅助数据。3.根据权利要求1所述的一种时序MODIS定量遥感产品时空修复及批处理实现方法，其特征在于：在最佳辅助数据的筛选这一步骤中中，当存在多个所述最佳辅助数据时，在时序中找到与待修复时刻数据具有互补观测信息的所有MODIS数据，并以时刻为分界线得到两组数据序列，在这两组所述序列数据中分别进一步筛选与T0时刻数据有效像元重叠数足够多的数据，则选择与T0时刻平均量化值差异最小的两景MODIS数据作为最优...

【专利技术属性】
技术研发人员：吴鹏海，朱娴飞，吕永辉，宋旭颖，陆雨婷，马晓双，吴艳兰，
申请(专利权)人：安徽大学，
类型：发明
国别省市：安徽,34