遥感时序平滑要保留真实转折并显式处理观测缺口

遥感时间序列常被云、阴影、轨道周期和质量筛选打断，保留下来的观测还包含噪声与异常。平滑和缺口填补可以帮助识别季节曲线与变化时点，也可能把突发事件抹掉或在无观测区制造确定曲线。处理目标应是估计可解释趋势，同时保留原始证据与不确定性。

先理解时间轴与质量

序列日期可能是瞬时观测日，也可能代表八日、十六日或月度合成窗口。合成值的时间标签未必是实际入选像元日期。分析物候或灾害前，应读取产品时间定义和质量层，不把规则间隔文件名直接当作等间隔现场观测。

云和阴影不是随机缺失。雨季、山地和高纬地区的空缺具有季节性，缺失时段可能正是现象快速变化期。只在晴空样本上拟合会偏向特定天气条件。应统计每个像元和时期的有效观测数、最长缺口与质量分布。

质量筛选太宽会保留污染值，太严会造成长缺口。可按任务设置多个质量等级，用权重而非简单删留表达差异。原始 QA 位、筛选规则和被排除比例都要保存。

孤立尖峰可能来自残云、雪、配准误差或传感器异常，也可能是火灾、洪水和收割。按固定阈值删除所有突变，会把真正事件当噪声。异常判断应结合邻近波段、质量位、空间邻域和其他日期。

稳健拟合降低极端点权重，适合处理少量污染，但仍可能压低真实短期变化。可以比较原始值、质量加权拟合和不含该点的曲线，并将“被降权”与“无效”分开记录。

空间邻域有助识别大范围事件，但小地块和线状目标可能只影响少数像元。邻域一致不是绝对规则，应与对象尺度和分辨率匹配。时间与空间检查都无法替代独立影像或现场证据。

移动平均简单直观，会模糊峰值和转折，窗口越大越明显。Savitzky-Golay 类局部多项式方法能够较好保留局部形状，仍受窗口、阶数和不规则采样影响。直接把缺失日期当作已有等间隔样本会改变权重。

样条可以构建平滑连续曲线，过多自由度可能在缺口内振荡，过强惩罚则压平季节峰。非对称高斯、双逻辑等季节模型适合单峰或特定物候形态，对多季作物、常绿植被和突发扰动未必合适。

TIMESAT 使用平滑与拟合方法提取卫星序列的季节参数，并考虑噪声与权重。具体算法和参数应按数据、植被类型与时间分辨率选择。软件输出的返青、峰值和衰退日期是曲线规则下的估计，不是地面观测日期。

短缺口可以由前后趋势约束，长缺口则可能包含未观测事件。线性插值假设两端之间匀速变化，季节模型假设形状重复，空间填补假设邻近像元相似。每种方法都引入条件，不能把填充值与观测值使用同一质量标识。

多传感器融合可增加有效日期，但要协调波段响应、空间网格、观测角度和定标。传感器切换可能形成阶跃，若平滑器把阶跃当地表变化，会产生错误转折。融合前应在重叠期比较偏差并保留传感器来源。

缺口跨越目标事件窗口时，可以返回区间或无法确定，而不是强制给出单一日期。面向统计汇总时，应报告有足够观测的面积比例，避免只展示成功拟合区域。

序列开头和结尾只有单侧观测，局部窗口不完整，平滑值和导数更不稳定。复制端点或镜像填充会加入人为形状。物候参数靠近年度边界时，需要使用相邻年份数据或明确降低置信度。

跨年曲线不能总按日历年截断。南半球季节、冬小麦和跨年雨季会把一个生长周期切成两段。分析窗口应依据地区和目标现象，日期编码还要处理闰年和不同时区产品。

长期趋势与季节变化可以分解，但趋势模型会受传感器更换、产品版本和土地覆盖突变影响。先检查数据处理断点，再解释缓慢变化为生态趋势。

返青、峰值、季末和生长季长度通常由曲线阈值、导数或拟合参数定义。不同定义会给出不同日期，跨研究比较前必须统一方法。固定绝对阈值对不同植被背景不公平，相对振幅阈值又会受年度幅度变化影响。

季节振幅很小时，微小噪声就可能造成日期大幅跳动。应设置最低振幅、有效观测和拟合质量门槛，并将失败像元标记为无结果。两个生长季需要多峰模型，强行选一个主峰会丢失种植制度信息。

变化检测与平滑存在张力。平滑有助降低噪声，过强则推迟断点、降低变化幅度。应在已知变化日期或模拟序列上测量检测偏移，并保留一条较少处理的序列供复核。

验证可使用地面物候相机、样地记录、高分辨率影像和独立传感器。时间支持范围要匹配：一次现场访问无法完整验证十六日合成曲线。按地表类型、云量、缺口长度和观测数分层报告误差。

交付应同时包含原始有效观测、拟合曲线、填补标记、权重、参数和质量指标。只保存平滑结果会失去判断过拟合和真实事件的依据。参数、软件版本和边界处理应写入元数据。

MODIS 土地覆盖动态产品等标准产品提供算法定义和质量信息，使用其物候日期时仍要遵守产品含义。自建曲线也应达到同样的可追踪程度。

时序平滑的目标不是让曲线更漂亮，而是在观测不完整时给出受证据约束的趋势。把填补与观测分开、让长缺口保持不确定，才能避免从平滑曲线中读出并不存在的精确事件。

Jönsson、Eklundh：TIMESAT，https://doi.org/10.1016/j.rse.2004.05.006
NASA：MODIS Land Cover Dynamics MCD12Q2 Version 6.1，https://doi.org/10.5067/MODIS/MCD12Q2.061