面向SWOT派生内陆地表水体产品的分析就绪数据处理框架开发

《Frontiers in Environmental Science》：Developing an analysis-ready data processing framework for SWOT-derived inland surface water products

　　

地表水与海洋地形（Surface Water and Ocean Topography, SWOT）任务通过观测地表水范围与水表面高程（Water Surface Elevation, WSE）为内陆水体监测提供了新机遇。然而，直接将SWOT L2_HR_Raster产品用于地表水制图仍受限于破碎水体模式、随机噪声、混合像元及系统性条带伪影。为提升产品可用性，研究人员开发了一种面向SWOT派生内陆地表水产品（基于地表水范围与水掩膜结构而非独立WSE精度验证）的分析就绪数据（Analysis-Ready Data, ARD）处理框架。该框架集成wse_qual、wse_uncert、water_frac、water_area、sig0、dark_frac与layover_impact进行多参数滤波与空间结构优化，其中WSE相关变量用作产品原生质量指示符。研究选取洞庭湖为研究区，仅以Sentinel-1合成孔径雷达（Synthetic Aperture Radar, SAR）派生水掩膜作为水体范围一致性与空间结构评估的参考数据。结果表明，ARD处理适度提升了SWOT派生水掩膜与Sentinel-1参考间的水体范围一致性：交并比（Intersection over Union, IoU）由0.2926增至0.4081，精确率（Precision）由0.3203增至0.5154，F1分数由0.4527增至0.5796，而召回率（Recall）由0.7718降至0.6621，表明在减少假阳性与保留参考水体像元间存在权衡。空间破碎性大幅降低：连通域由18,187减至557，小斑块面积比由5.63%降至0，边界密度由0.0048降至0.0034。沿轨条带伪影的持续存在构成当前ARD框架的核心约束，表明几何相关系统误差需超越像元级滤波与空间优化的额外校正。总体而言，所提框架提升了SWOT派生内陆地表水产品在地表水制图与洪水监测中的分析就绪性。

研究背景方面，地表水资源是全球水循环关键组分，传统光学遥感（如Landsat、Sentinel系列）受云、气溶胶与光照限制，而SWOT卫星采用宽幅干涉测高可同步获取高分辨率水范围与WSE，但原始L2_HR_Raster产品存在破碎水体、随机噪声、混合像元、暗水体低后向散射不稳定、复杂地形layover与shadow效应及沿轨条带伪影等问题，且单一质量标志筛选不足，缺乏将WSE质量变量（wse_qual、wse_uncert）与water_frac、water_area、sig0、dark_frac、layover_impact联合约束的可复现工作流，导致“可获取但难直接用”。为此研究人员提出面向内陆地表水的SWOT ARD框架，以洞庭湖东洞庭山子区为案例，用2024年7月29日SWOT L2_HR_Raster（100 m，UTM 49R）为主源，Sentinel-1 SAR GRD（10 m，EPSG:32649，VV≤?13 dB阈值提取参考掩膜并聚合至100 m网格）作参考，评估水体范围一致性与空间结构。论文发表于《Frontiers in Environmental Science》。

关键技术方法上，研究人员构建多参数联合约束像元级筛选：以有限值WSE定义初值有效掩膜，wse_qual≤2排除低质、wse_uncert≤1.0 m排除高不确定像元，water_frac≥0.3抑制强混合像元，water_area≥3000 m²压低不稳定小检测，sig0保留场景特异P2–P98分位以去极端后向散射，dark_frac≤0.8抑暗水异常，layover_impact≤5.0减几何污染；各约束取逻辑与得M_ARD，满足者保留WSE否则置NaN。随后做空间结构优化：8邻域连通域分析、移除小于20像元孤立小块。评估在共同有效100 m网格上，以Sentinel-1聚合分数≥0.5判参考水，计算IoU、Precision、Recall、F1及连通域数、小斑块数、小斑块面积比、边界密度。

研究结果部分，4.1节原始SWOT数据问题指出，原始产品中部沿轨条带异常（空值/异常值/破碎像元交错）、东部南部边缘散斑噪声与湖滨过渡带边界破碎，呈“可观测不可直接分析”状态。4.2节ARD产品在洞庭湖流域结果显示，主湖体与连通河道保留，条带区伪影减弱但未消除，湖滨与河道边界连续性提升，噪声区孤立斑块大幅减少，低噪声区稳定；视觉更干净但可能误删真实小水体，故仅证空间相干性提升而非证明被删皆错。4.3节水体范围一致性与空间结构评估表明，相较原始SWOT，ARD的IoU 0.2926→0.4081、Precision 0.3203→0.5154、F1 0.4527→0.5796，Recall 0.7718→0.6621（假阳性大减但参考水像元亦失）；空间上小斑块17,418→0、连通域18,187→557、边界密度0.0048→0.0034、小斑块面积比5.63%→0，破碎性降低部分源于流程自身连通域滤波，但相对Sentinel-1的Precision/IoU提升佐证书删碎片多与参考不一致。

讨论部分，5.1节分析质量问题的合理解释：沿轨条带可能源于观测几何、轨道缝数据不稳、干涉处理不确定与局部缺测，属系统伪影而非随机噪，故局滤难除；局部碎片与暗水、混合像元、滨线几何效应相关。5.2节析ARD效限：框架减像元不稳与碎片但保主水文结构，对条带仅能减弱不能复原缺失区，故研究避免跨条带插值而保留空值。结论翻译为：研究人员开发了面向SWOT派生内陆地表水产品的ARD处理框架，强调水范围与水掩膜结构而非独立WSE验证；框架通过多参数滤波与空间优化适度提升与Sentinel-1参考的范围一致性并大幅降破碎性；沿轨条带伪影持续存在，指明几何相关系统误差需超像元级滤波的附加校正；该可复现工作流提升了SWOT内陆水产品在制图与洪水监测中的分析就绪性，但非独立去噪算法，亦未含元数据标准化、不确定传播与WSE独立验证等未来扩展。