SCI投稿准备 · 2026/10/22
遥感产品更新了处理版本,前后时间序列怎样比较?
遥感长时间序列遇到产品更新时,应核对场景版本、缩放系数、质量掩膜和重叠期差异。
一篇遥感论文把十年前下载的地表反射率和今年新下载的影像排成一条趋势线,曲线上某年恰好出现转折。作者以为地表植被变了,后来才发现两段数据来自不同的产品集,缩放系数、填充值和质量标记的定义也不一样。同一传感器的名字不能替代产品版本;要把版本留在方法和数据清单里。
先回到每景数据的产品标识。记录任务、传感器、处理等级、产品集或集合版本、获取日期、场景编号及处理日期。USGS的Landsat说明把Collection 1与Collection 2分开,并提供再处理事件记录;同一观测在后续处理后可能生成修订产品。若方法只写“使用Landsat影像”,读者无法重建所用数据集。数据表或代码清单应列出场景标识和版本,仓储里保留下载时的元数据快照。
接下来检查文件里的整数怎样换算成反射率。以Landsat地表反射率为例,USGS明确说明不同集合的缩放系数、加性偏移和填充值有差别。若把旧脚本原封不动用于新集合,数字可能仍落在看似合理的范围,但某些年份会整体偏移。方法段应写明应用了哪一版本文档的缩放规则,并对一景样本人工复算。尤其不要把填充值当作零反射率送入年度平均,这会制造人为低谷。
质量掩膜也要重做核查。云、云影、雪、饱和像元及无效值的编码可能随产品变化;USGS为Collection 2单列了QA波段说明。方法中应注明采用哪些质量位、排除顺序和边界处理;若某年可用像元明显较少,图中标出观测数量或有效覆盖比例。一个趋势转折若同时伴随云筛选方法更换,应先检验质控口径造成的差异,再谈生态或气候机制。
如果产品可得且方法适用,优先用同一产品集重建全时段序列。如果旧版数据无法重新获得,就在版本交界处设置重叠期:对同一批场景用新旧流程各算一次,比较差值分布、季节性和地物类型依赖。发现系统偏差后,说明校正方法与剩余不确定度;无法建立可靠桥接时,分段展示,避免把跨版本跳跃写成真实地表变化。不能只用一个整体相关系数证明两套产品完全可替换。
图表还应暴露时间序列的分母。每个年度或季节点由多少景、多少有效像元构成,采用平均数、中位数还是合成规则,是否有缺测插值,都要能复核。空间范围若由旧行政边界换成新边界,也会改变汇总值,这属于另一个口径变化,应与产品版本分开标识。若跨传感器,需要另查波段响应和观测条件,而非以“均为Landsat”代替方法比较。
交稿前挑三个点来拆算:起点、版本交界和终点,从场景号追到元数据、缩放、质量位和聚合结果。再用同一规则重绘一条不包含低质量像元的敏感性曲线。只有在这两个检查后仍稳固的转折,才适合进入环境解释。论文把产品版本写清,读者将来面对再次再处理的数据,也能理解当时的分析选择。