DB 3705T 44-2024 芦苇遥感监测技术规范

ICS65.020ICS65.020CCSA77东 营 市 地 方 标 准DB3705/T44—2024芦苇遥感监测技术规范TechnicalspecificationsforremotesensingmonitoringofPhragmitesaustralis2024-12-17发布 2025-01-17实施东营市场监理局 发布PAGE\*ROMANPAGE\*ROMANIIDB3705/T44—2024目 次前言 II1范引文件 1语定义 1学础 1作程 1据集 2基数据 2卫遥图像 2现数据 2据理 3卫遥图处理 3现数处理 3感息取 3分区地盖遥感译志立 3芦植空分息提与证 3空分信处析 4空分信矢据 4空分信验证 4覆度感算 4物遥估算 5模构建 5生量算 5测果 6数集 6图件 6监报告 6附录A（料） 附表 7DB3705/T44—2024前 言本文按照GB/T1.1—2020《标化作导则 第部分标准文的结和起规则的起草。请注意本文件的某些内容可能涉及专利。本文件的发布机构不承担识别专利的责任。本文件由东营市自然资源和规划局提出、归口并组织实施。11DB3705/T44—2024芦苇遥感监测技术规范范围本文件适用于芦苇空间分布、覆盖度和地上生物量的遥感监测工作。（HY/T0460.4 4HY/T147.7 DB37/T4747 下列术语和定义适用于本文件。芦苇落 Phragmitesaustralis根据其生态环境，可分为潮滩盐沼芦苇、淡水沼泽芦苇和旱生草甸芦苇三类。2000:19853°分带。。22DB3705/T44—2024数据收集与数据处理数据收集与数据处理遥感解译标志建立现场数据基础数据卫星遥感图像遥感解译标志建立现场数据基础数据生物量遥感估算模型芦苇提取算法生物量遥感估算模型否优化F1分数>90%

否优化平均相对误差<30%是芦苇空间分布信息 生物量遥感估算芦苇空间分布信息生物量遥感估算监测成果潮滩盐沼芦苇淡水沼泽芦苇旱生草甸芦苇监测成果潮滩盐沼芦苇淡水沼泽芦苇报告图件数据集报告图件芦苇覆盖度信息芦苇覆盖度信息图1 芦苇感测作程芦苇遥感监测的基础数据来源于自然资源和规划部门，包括：卫星遥感图像符合以下要求：2m610GeoTIFF。HY/T0460.433DB3705/T44—2024HY/T147.7kg/m2HY/T147.7利用现场调查数据，根据芦苇植被和作业区内其他土地覆盖类型在卫星遥感图像中的色调、光谱、纹理、形状和空间分布等特征差异建立遥感解译标志，填写土地覆盖类型遥感解译信息表（附录A中表A.1）。利用不同的方法对黄河三角洲潮滩盐沼芦苇、淡水沼泽芦苇和旱生草甸芦苇等3种芦苇类型进行信息提取。计算潮滩范围遥感图像的土壤调节植被指数（SAVI）（计算公式见附录B中B.1）和归一化差值水体指数(NDWI)（计算公式见附录B中B.2）。根据现场调查数据中的芦苇空间分布信息，结合计算得到的SAVI、NDWI指数，宜按照HY/T147.7中规定的决策树方法，提取盐沼芦苇植被空间分布栅格图。B中B.3）指HY/T147.744DB3705/T44—2024芦苇空间分布矢量数据符合以下要求：ESRIShapefileGeodatabaseA中表A.2。通过现场数据对提取结果的总体精度进行评价，采用F1分数作为分类精度评价指标，芦苇提取精度的F1分数度低于90%时，需更换相关指数阈值、分类训练样本或者更换遥感分类方法重新开展提取。F1分数计算方法见公式（1）。𝐹1=2×𝑃𝑖𝑖/𝑃𝑖𝑖+𝑃𝑗𝑖×𝑃𝑖𝑖/𝑃𝑖𝑖+𝑃𝑖𝑗………………(1)𝑃𝑖𝑖/𝑃𝑖𝑖+𝑃𝑗𝑖+𝑃𝑖𝑖/𝑃𝑖𝑖+𝑃𝑖𝑗式中：𝑖 𝑗 ———（2）。𝐹𝑣𝑐=𝑆𝑁𝐷𝑉𝐼−𝑆𝑁𝐷𝑉𝐼,𝑠𝑜𝑖𝑙𝑆𝑁𝐷𝑉𝐼,𝑣𝑒𝑔−𝑆𝑁𝐷𝑉𝐼,𝑠𝑜𝑖𝑙…………(2)式中：𝑆𝑁𝐷𝑉𝐼 B中B.3；𝑆𝑁𝐷𝑉𝐼,𝑠𝑜𝑖𝑙 55DB3705/T44—2024𝑆𝑁𝐷𝑉𝐼,𝑣𝑒𝑔 注：选择现场调查多处裸土覆盖度100%位置NDVI值的平均值作为全裸土区域的归一化差值植被指数；选择现场调查芦苇覆盖度100%位置NDVI值作为纯芦苇覆盖区域的归一化差值植被指数。70（3）。𝑊=𝑎+𝑏×𝐹𝑣𝑐+𝑐1×𝑆1+𝑐2×𝑆2+⋯+𝑐𝑛×𝑆𝑛…………(3)式中：𝑊 𝑎 𝑏 𝑛 𝑐1,𝑆1,注：模型常数abcc2cn度，通过统计回归算法确定。30%（4）。𝛿𝑀𝑅𝐸=

𝑛𝑌𝑌𝑖−𝑍/𝑍𝑖 𝑖

𝑛 （𝑖=1,2,3,…,𝑛）………………(4)式中：𝛿𝑀𝑅𝐸 𝑛 ii若某类型芦苇生物量遥感估算的平均相对误差大于30%，则应重新选择回归模型的遥感植被指数，再次更新模型参数，重新开展模型构建。66DB3705/T44—2024利用构建的芦苇生物量遥感估算模型公式（3），估算芦苇地上生物量。通过通用地理信息系统分析软件自动获取芦苇地上生物量统计结果。数据集图件遥感影像专题图的地图质量要求、输出分辨率要求、要素信息要求按照DB37/T4747中的规定执行。要素信息专题图的输出分辨率要求、要素信息要求按照DB37/T4747中的规定执行。DB37/T474777DB3705/T44—2024附录A（资料性）附表表A.1土地覆盖类型遥感解译标志信息表表A.1给出了土地覆盖类型遥感解译标志信息格式。表A.1土地覆盖类型遥感解译标志信息表序号土地覆盖类型经纬度遥感图像信息土地覆盖类型遥感图像特征土地覆盖类型现场特征气象条件描（度等）遥感图像截图图像特征描述（波段组合、色调、形状、分布以及纹理等）现场照片特征描述（生长形态、覆盖度等）1经度：° ′″E纬度：° ′″N图像类型填写所用遥感图像的名称空间分辨率填写所用遥感图像的空间分辨率成像时间YYYYMMDD制表人： 制表时间： 检查人：88DB3705/T44—2024表A.2给出了芦苇空间分布信息矢量数据属性信息表。表A.2芦苇空间分布信息矢量数据属性信息表序号字段名称数据类型单位备注1面积浮点型hm22中心经度浮点型°3中心纬度浮点型°4中央经线浮点型°5坐标系字符型--6投影方式字符型--7斑块编号整数型--8斑块周长浮点型m9传感器名称字符型--10成像时间整数型--11空间分辨率浮点型--12提取时间日期型--13提取单位字符型--14检查单位字符型--15检查人字符型--16检查时间日期型--17审核单位字符型--18审核人字符型--19审核时间日期型--20备注字符型--99DB3705/T44—2024附录B（）=𝑅𝑁𝐼𝑅−𝑅𝑅𝐸𝐷1+𝐿𝑅𝑁𝐼𝑅+𝑅𝑅𝐸𝐷+𝐿 (B.1)式中：𝑅𝑁𝐼𝑅 𝑅𝑅𝐸𝐷 𝐿 [0,1]。绿光波段与近红外波段反射率的波段运算。用来提取影像中的水体信息。计算公式如下：𝑆𝑁𝐷𝑊𝐼=𝑅𝐺𝑅𝑁−𝑅𝑁𝐼𝑅𝑅𝐺𝑅𝑁+𝑅𝑁𝐼𝑅 (B.2)式中：𝑆𝑁𝐷𝑉𝐼 𝑅𝐺𝑅𝑁 𝑅𝑁𝐼𝑅 近红外波段的反射率减去红光波段的反射率再除以两者之和。计算公式如下：𝑆𝑁𝐷𝑉𝐼=𝑅𝑁𝐼𝑅−𝑅𝑅𝐸𝐷𝑅𝑁𝐼𝑅+𝑅𝑅𝐸𝐷 (B.3)式中：𝑆𝑁𝐷𝑉𝐼 𝑅𝑁𝐼𝑅 𝑅𝑅𝐸𝐷 近红外和红光波段反射率的比值。适用于覆盖度大于50%的植物分布区域。计算公式如下：=𝑅𝑁𝐼𝑅𝑅𝑅𝐸𝐷 (B.4)1010DB3705/T4