《集中式地表水饮用水水源保护区农业种植污染遥感调查技术规范》（征求意见稿）

本标准面向水源保护区环境执法与监管场景，针对农业种植面源污染“发现难、定位难、定量难”的突出问题，构建了基于多源遥感影像的种植污染调查技术体系。全文规定了从数据准备、复种分区、作物调查到污染负荷估算与成果归档的完整工作流程，以耕地地块为基本空间单元，最终输出总氮、总磷年流失负荷的空间分布与总量结果，为保护区农业污染风险分级和优先治理区域识别提供技术依据。

数据基础与前期处理

标准对输入数据提出分层要求：高空间分辨率影像（优于或等于10 m）用于提取耕地地块边界，高时间分辨率影像（优于或等于30 m、时间分辨率≤10 d）用于构建覆盖完整自然年的NDVI时间序列。两类影像均须经辐射校正、大气校正、几何精校正等预处理。矢量数据方面，优先采信自然资源部门发布的耕地地块数据，保护区边界须满足一级、二级及准保护区之间拓扑逻辑合理且无重叠。对含云、气溶胶污染的原始时序，采用滑动窗口回归、滤波迭代等方法重建高质量NDVI序列，确保所有耕地地块获得完整有效观测值。

复种类型分区与作物构成调查

分区以耕地地块NDVI时间序列为聚类对象，聚合地块内部像素值后，选用k均值、模糊c-medoids或层次聚类算法，以DTW、SBD等形状相似性距离进行归类。聚类数目预设1至5类，种植结构复杂区域可酌情增加。通过轮廓系数、Davies-Bouldin指数或Dunn指数对候选方案综合评价，确定最优复种分区。作物构成调查基于各分区聚类中心的NDVI波形，先以动态阈值法（推荐20%，可本地化调整至30%或50%）识别年内生长季数量与起止时间，再按NY/T 3527执行目视判读、监督分类或光谱匹配确定各生长季主栽作物。

精度验证与污染负荷估算

复种分区与作物构成结果须基于独立验证样本构建混淆矩阵检验，总体精度不低于95%，不达标则重新调查直至合格。污染负荷估算采用“面积×流失系数”法：年流失系数按省级及以上主管部门发布数据、相似区域研究成果本地化修正、规范实验实测（依据GB/T 41222）三级优先顺序确定。负荷计算依据《排放源统计调查产排污核算方法和系数手册》，分别统计总氮与总磷的保护区总量，并以耕地地块为制图单元生成分级设色空间分布图。标准明确该方法存在不确定性，成果定位于危险性评估、风险分级及优先治理区识别，而非精确模拟。

成果归档要求

遥感影像成果以GeoTIFF格式存储，复种类型与作物构成成果以Shapefile矢量归档，属性须含地块编号、保护区级别、调查年份、复种类型、作物构成等字段。污染负荷成果包括调查表（xlsx/csv）、氮磷负荷空间分布专题图和调查报告（docx/wps），报告须明确优先治理区域。所有空间数据采用2000国家大地坐标系、1985国家高程基准及经差3度分带高斯-克吕格投影，保证成果可比与可追溯。