遥感论文的云掩膜规则，投稿时怎样写到可复现？

遥感论文写“删除有云影像后得到年度植被指数”，却没说删的是整景影像，还是只屏蔽被云覆盖的像元。两种处理留下的观测数量不同，可能改变年度合成结果。以Landsat Collection 2为例，美国地质调查局的QA_PIXEL质量波段用多个比特记录云、扩张云、云影、卷云等状态；方法里只写“按质量波段去云”仍不足以让别人重复。投稿时应交代产品版本、使用的标志位和排除后实际剩余的数据。

先区分整景筛选与逐像元掩膜

整景筛选通常依据影像的总体云量或研究者设定的阈值决定是否纳入；逐像元掩膜则在已纳入影像内剔除特定像元。两者可以同时使用，但要分别记录规则。比如先限制场景云量，再对每景应用QA_PIXEL，最后把可用像元做时间合成。若只写“云量低于20%”，读者不知道研究区内的云像元是否仍被剔除；若只写“使用QA_PIXEL”，也不知道哪种标志被视为不可用。

质量波段不是一个直接代表“好”或“坏”的单一数值。USGS说明Landsat 8与9的QA_PIXEL把不同条件编码在不同位上，且不同传感器的位定义需查对应产品说明。正文或补充方法可列出排除的类别，例如高置信云、扩张云和云影，并写明是否处理卷云、雪和填充值。若研究覆盖Landsat 4至7的水域，还应留意USGS公布的Collection 2“Clear”位已知问题，不能把Landsat 8的处理语句无条件复制过去。

把掩膜前后的数量留下来

对每个年份或季节，至少能回查原始场景数、经整景筛选留下的场景数、研究区内最终可用像元或观测次数。多时相合成若按中位数、最大值或其他规则计算，也需说明像元在缺少有效观测时怎样处理。云多的雨季可能剩下很少观测；此时地图上的空白或填补区域不应被悄悄当作真实地表值。可把有效观测次数图或覆盖率表作为质量信息，显示结果的空间差异。

设想一个自拟湿地研究，某年有12景影像，整景云量筛选后留8景，逐像元掩膜后各位置有效观测数介于2至8次。方法段应同时写出两道筛选及年度合成规则，结果或补充材料给出有效观测分布。“12景”是场景总数，不能当作每个像元的12次有效测量；湿地范围内的覆盖还要单独核对。数字仅作示范，真实项目应由处理日志复算。

投稿前从代码或平台脚本中核对数据集ID、产品级别、传感器、QA位筛选逻辑和合成函数，再与方法文字及图例对照。若处理流程更新过，保留脚本版本和重新运行后的数量，不用旧图配新方法。写清云掩膜规则的目的，是让植被或水体变化解释建立在可见的观测范围之上。

参考信息：美国地质调查局Landsat Collection 2质量波段说明及已知问题；USGS质量波段示例代码。