HJ 1098-2020 水华遥感与地面监测评价技术规范.pdf
《HJ 1098-2020 水华遥感与地面监测评价技术规范.pdf》由会员分享,可在线阅读,更多相关《HJ 1098-2020 水华遥感与地面监测评价技术规范.pdf(15页珍藏版)》请在麦多课文档分享上搜索。
1、 中华人民共和国国家环境保护标准 HJ 1098-2020 水华遥感与地面监测评价技术规范 (试行) Technical specifications for monitoring and evaluating algal bloom based on remote sensing and field monitoring ( 发布 稿 ) 本电子版为发布 稿。请以中国环境出版集团出版的正式标准文本为准。 2020-02-12 发布 2020-04-12 实施 生态 环 境 部 发 布 i 目 次 前 言 . ii 1 适用范围 . 1 2 规范性引用文件 . 1 3 术语和定义 . 1 4
2、水华遥感监测 . 3 5 水华地面监测 . 7 6 水华程度评价 . 9 ii 前 言 为贯彻中华人民共和国环境保护法和中华人民共和国水污染防治法,保护 生态 环境,保障人体健康,规范和指导我国淡水水体藻类水华监测和评价工作,制定本标准。 本标准规定了淡水水 体 藻类水华的遥感监测方法、地面监测方法和水华程度评价方法等 内容。 本标准为首次发布。 本标准由 生态 环境部 生态 环境监测司 、 法规与 标准司组织制订。 本标准起草单位:中国环境监测总站、 生态环境部 卫星环境应用中心 。 本标准 生态 环境部 2020年 2月 12日批准。 本标准自 2020年 4月 12日起实施。 本标准由
3、生态 环境部解释。 1 水华遥感与地面监测评价技术规范 1 适用范围 本标准规定了淡水水体 蓝藻 水华遥感监测方法、地面监测方法和水华程度评价方法等内 容 。 甲藻、硅藻及其他藻类水华监测 与评价 可参考 使用 本标准。 本标准适用于我国淡水 湖库 水华监测、评价与管理 。 其他 淡水 水 体 的水华监测与评价可 参考使用本标准。 2 规范性引用文件 本标准内容引用了下列文件或其中的条款。凡是不注明日期的引用文件,其最新版本适 用于本标准。 GB/T 14950 摄影测量与遥感术语 GB/T 30115 卫星遥感 影像植被指数产品规范 3 术语和定义 GB/T 14950 和 GB/T 301
4、15 界定的以及下列术语和定义适用于本标准 。 3.1 水华 algal bloom 淡水水体中藻类大量繁殖的一种自然生态现象,表观特征为水体表面有藻类聚集或藻类 颗粒悬浮在水体当中。 3.2 藻密度 algal density 单 位体积水样中藻类的细胞个数。单位:个 /L。 3.3 水华面积 algal bloom area 藻类漂浮在水面时覆盖水体的面积。 3.4 水华 面积 比例 algal bloom area proportion 水华面积占监测水体面积的比例,百分比( %)。 3.5 水华程度 algal bloom level 反映水华发生的强弱特征,根据水体中藻密度的高低或水
5、华面积比例的大小进行判定。 2 3.6 标准 假彩色合成图像 standard false color composite image 根据加色法彩色合成原理,将遥感图像的近红、红、绿波段分别通过红、绿、蓝滤光系 统合成产生的彩色图像。 3.7 空间 分辨率 spatial resolution 指 像元代表的地面范围的大小。 3.8 时间分辨率 temporal resolution 传感器能够重复获得同一地区影像的最短时间间隔。 3.9 辐射定标 radiometric calibration 根据遥感器的定标方程和定标系数,将其记录的量化数字灰度值转换成对应现场的表观 辐亮度的过程。 3
6、.10 几何校正 geometric correction 为消除影像的几何畸变而进行 的 投影变换和不同波段影像的套合等校正工作。 3.11 大气校正 atmosphere correction 消除或减弱卫星遥感影像获取时在大气传输中因吸收或散射作用引起的辐射畸变。 3.12 数字灰度值 digital number( DN) 由遥感器各波段获取的反射或辐射能量量化而成的灰度等级。 3.13 表观 辐亮度 apparent radiance 大气层顶辐亮度,卫星遥感器入瞳处的辐射亮度。 3.14 表观 反射率 apparent reflectance 表观辐亮度与无大气水平场景绝对白体假设
7、下遥感器应获得入瞳辐亮度之间的比值 。 3.15 地表 反射率 surface reflectance 地物表面反射能量与到达地物表面的入射能量的比值。 3.16 归一化植被指数 normalized difference vegetation index( NDVI) 近红外波段与可见光红波段数值之差和这两个波段数值之和的比值。 3 3.17 数字灰度 值归一化植被指数 normalized difference vegetation index of digital number( NDVIDN) 基于数字灰度值( DN)计算得到的归一化植被指数。 3.18 大气层 顶归一化植被指数 no
8、rmalized difference vegetation index of top of atmospheric reflectance( NDVITOA) 利用大气层顶表观反射率计算得到的归一化植被指数。 3.19 大气层 底归一化植被指数 normalized difference vegetation index of bottom of atmospheric reflectance( NDVIBOA) 利用经过大气校正的地表反射率计算得到的归一化植被指数。 3.20 常规 监测 routine monitoring 以系统掌握监测 水 体 藻类状况为目的的 例行 监测 行为 。
9、3.21 预警监测 warning monitoring 在水华 敏感 期,以及时掌握监测 水 体 藻类生长和发展状况,判断 水华发生风险 为目的的 监测行为。 3.22 应急 监测 emergency monitoring 在水华暴发期,以跟踪掌握监测 水 体 藻类水华状况为目的的监测 行为 。 4 水华遥感监测 监测目的 4.1 采用 卫星遥感影像数据开展水体的蓝藻水华监测,获得水体中蓝藻水华的空间 分布 位置, 计算蓝藻水华面积及其所占水体面积比例,据此评价水华程度。 监测原理 4.2 采用 卫 星遥感 影像 数据 监测蓝藻水华主要是基于正常水体 光谱 与发生 水华 水体 光谱 的 差异
10、。 蓝藻水华暴发时,藻类 聚 集 在水体表面, 因其 对 红光波段的强吸收导致 产生 的红光波 段反射率 较低 ,在近红外波段具有类似于植被光谱曲线特征的“反射峰平台效应”, 近红外 波段反射率较高。 而正常水体 对 近红外 波段 有强烈的吸收 作用, 导致反射率 较低 。 因此,通 过计算植被指数可以区分水华和正常水 体。 监测内容 4.3 水体面积 、 蓝藻水华面积和 空间 分布位置。 监测频次 4.4 4 根据实际工作需要确定监测频次 ,条件允许时建议每天开展监测 。 监测技术流程 4.5 水华遥感监测技术流程如图 1 所示。 遥 感 影 像 预 处 理 ( 几 何 校 正 , 空 间
11、裁 剪 , 辐 射 定 标 , 表 观 反 射 率 , 大 气 校 正 , N D V I 计 算 ) 卫 星 遥 感 影 像 D N 值 水 体 掩 膜 云 掩 膜 水 草 掩 膜 水 华 二 值 图 , 0 , 1 归 一 化 植 被 指 数 ( N D V I D N 或 N D V I T O A 或 N D V I B O A ) 阈 值 分 割 水 体 区 域 的 N D V I 影 像 水 华 面 积 水 华 面 积 比 例 水 体 面 积 图 1 基于卫星遥感影像数据监测蓝藻水华技术流程 监测方法 4.6 4.6.1 遥感数据选择 选取 具有红光( 630 nm 690 nm)
12、和近红外( 760 nm 900 nm)波段的卫星遥感 影像 数据 。 空间分辨率应取决于监测水体的大小,水体面积越小,空间分辨率要求越高。 遥感影 像数据至少 覆盖 监测水体 90%以上面积;在监测水体内,云覆盖率小于 50%。 水陆掩膜、 云掩膜、水草掩膜及观测几何条件等辅助数据坐标系及投影应与遥感 影像 数据保持一致。 4.6.2 几何校正 基于参考影像或空间几何信息,开展 影像 几何精 度 校正和空间投影转换,精度控制在 1 个像元内。 4.6.3 空间裁剪 依据监测水体区域的左上角和右下角经纬度,对所选取的遥感影像进行空间裁剪,范围 要略大于水体区域。 4.6.4 辐射定标 辐射定标
13、根据遥感器的定标公式及各波段 定标系数进行。辐射定标公式的一般形式为: 5 = ( ) + ( ) ( 1) 式中: L 表观辐亮度,单位为 Wm-2sr-1m-1; DN 数字灰度值; Gain 绝对定标系数的增 益 量,单位为 Wm-2sr-1m-1; Offset 绝对定标系数的偏移量,单位为 Wm-2sr-1m-1,空缺值 0。 遥感器不同波段各有其定标系数。随着遥感器的老化和功能衰退,定标系数会有所改变。 4.6.5 大气层顶表观反射率计算 根据辐射定标获得的各波段表观辐亮度,按照公式 ( 2) 计算 各波段表观反射率: ( 2) 式中: TOA 大气层顶表观反射率(无量纲); 常量
14、(球面度 sr); L 表观辐亮度,由公式( 1)计算获得 , 单位为 Wm-2sr-1m-1; D 实际日地距离与平均日地距离的比值; F0 平均日地距离处的 大气层顶太阳谱照度, 单位为 Wm-2m-1; s 太阳天顶角。 L 和 F0 均与遥感器的具体波段有关,而 D 和 s 可以通过卫星过境的日期和时间以及监 测水体的经纬度计算获得。 4.6.6 大气校正 可见光波段和近红外波段的卫星遥感大气校正主要针对大气分子散射、气溶胶散射及水 汽吸收的影响。可采用基于辐射传输模型的方法(如 6S 模型、 Flaash 大气校正模型等)进 行大气校正,得到各波段地表反射率,也称为大气层底反射率 B
15、OA。 4.6.7 归一化植被指数计算 ( 3) 式中: 当 i 为 DN 时,公式为基于 DN 值的归一化植被指数, DN(NIR) 近红外波段的灰度值; DN(R) 红波段的灰度值; NDVIDN 的取值范围为 -1, 1。 当 i 为 TOA 时, 公式为 大气层顶归一化植被指数 , TOA(NIR) 近红外波段的大气层顶表观反射率; TOA(R) 红波段的大气层顶表观反射率; NDVITOA 的取值范围为 -1, 1。 2 T O A 0 c o s sLDF OffsetDNGainL 6 当 i 为 BOA 时,公式为大气层底归一化植被指数, BOA(NIR) 近红外波段 经过大气
16、校正的地表 反射率; BOA(R) 红波段 经过大气校正的地表 反射率; NDVIBOA 的取值 范围为 -1, 1。 归一化 植被指数 ( NDVI) 的 计算可以 采用 三种定量化 级别 的图像 作为 输入,分别是 未 经定量化处理的 原始 DN 值 、 经过 部分定量化处理的大气层顶表观反射率 、 经过 全面定量化 处理的 地表 反射率。 定量化 程度越高, 对应 的 NDVI 的水华阈值 越 稳定 , 但 处理 难度 以及 对 相关输入参数的要求越高。因此 , 在 相关 参数比较完备 , 处理经验比较丰富的情况下,尽量 使用定量化程度 较 高的 NDVI 数据 ;反之 则使用定量化程度
17、较低的 NDVI 数据。 4.6.8 水陆分离、云识别和水草识别 为避免蓝藻水华的误判,应结合实际情况准确识别被监测的水体区域,剔除云、水草等 的干扰,形成监测水体的水体掩膜、云掩膜和水草掩膜等辅助数据,获得包含蓝藻水华的水 体区域所有像元的 NDVI值,按不同处理方式可分 别获得像元的 NDVIDN、 NDVITOA和 NDVIBOA。 4.6.9 阈值分割和水华二值图 针对水体区域的 NDVIDN 和 NDVITOA 数据,蓝藻水华的 NDVI 值要高于正常水体,由于 未经过大气校正,判别蓝藻水华的 NDVI 阈值会随着不同影像发生变化,需要结合遥感假彩 色合成图像目视识别或 NDVI 灰
18、度直方图统计加以确定, NDVI 值 高于阈值的为蓝藻水华 。 针对经过准确大气校正的水体区域的 NDVIBOA 数据,判别水华区别于正常水体的阈值可以设 为 0, NDVI 值 高于 0 的像元为蓝藻水华。 蓝藻水华像元赋值为 1,其他地物(包括正常水体、云、陆地等)像元赋值为 0,获得 水华二值图。 4.6.10 水华面积和水华面积比例计算 基于水华二值图统计水华像元总数,乘以每个像元对应的实际面积,计算得出水华面积; 水华面积除以水体面积获得水华面积比例(百分比)。 ( 4) 式中: A1 水华面积; A 水体面积 (推荐 采用 以当年无云影像获取的水陆掩膜计算 出 的水体面积); P
- 1.请仔细阅读文档,确保文档完整性,对于不预览、不比对内容而直接下载带来的问题本站不予受理。
- 2.下载的文档,不会出现我们的网址水印。
- 3、该文档所得收入(下载+内容+预览)归上传者、原创作者;如果您是本文档原作者,请点此认领!既往收益都归您。
下载文档到电脑,查找使用更方便
5000 积分 0人已下载
下载 | 加入VIP,交流精品资源 |
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- HJ 1098-2020 水华遥感与地面监测评价技术规范 1098 2020 遥感 地面 监测 评价 技术规范
