ICS 07.060 B 18 DB37 山 东 省 地 方 标 准 DB37/T 3791—2019 大蒜长势卫星遥感评估技术规范 Specification for evaluation techniques on satellite remote sensing of garlic growth status 2019 - 12 - 24 发布 山东省市场监督管理局 2020 - 01 - 24 实施 发 布 DB37/T 3791—2019 目 次 前言 ................................................................................ II 1 范围 .............................................................................. 1 2 术语和定义 ........................................................................ 1 3 数据要求 .......................................................................... 2 4 评估时间 .......................................................................... 2 5 评估方法 .......................................................................... 2 6 评估流程 .......................................................................... 4 附录 A（规范性附录） 常用的中高空间分辨率星载传感器及其红光和近红外波段参数 .......... 5 附录 B（规范性附录） 大蒜长势等级地面调查评定指标 .................................... 6 附录 C（资料性附录） 基于 GF-1/WFV 卫星数据的大蒜长势评估示例 ......................... 7 参考文献 ............................................................................. 8 I DB37/T 3791—2019 前 言 本标准按照GB/T 1.1—2009给出的规则起草。 本标准由山东省气象局提出并组织实施。 本标准由山东省气象标准化技术委员会（鲁TC 16）归口。 本标准起草单位：山东省气候中心。 本标准主要起草人：李峰、胡先锋、薛晓萍、张承明、王晓默。 II DB37/T 3791—2019 大蒜长势卫星遥感评估技术规范 1 范围 本标准规定了大蒜长势卫星遥感评估的数据要求、评估时间、评估方法和评估流程。 本标准适用于大蒜长势卫星遥感评估。 2 术语和定义 下列术语和定义适用于本文件。 2.1 红光波段 red band 星载仪器涵盖的0.605 μm～0.700 μm的波长范围。 [QX/T 474—2019,定义2.4] 2.2 近红外波段 near infrared band 星载仪器涵盖的0.76 μm～1.25 μm的波长范围。 注：传感器在近红外波段所接收到的辐射主要是太阳辐射的反射。 [QX/T 474—2019,定义2.5] 2.3 归一化差值植被指数 normalized difference vegetation index；NDVI 用于植被指数计算的近红外波段和可见光波段反射率之差与这两个波段反射率之和的比值。 [QX/T 284—2015,定义2.3] 2.4 幼苗期 seedling stage 从大蒜初生叶展开到鳞芽及花芽开始分化的时期。 2.5 花芽和鳞芽分化期 flower bud and squamous bud differentiation stage 从大蒜花芽、鳞芽开始分化到分化结束的时期。 2.6 花茎伸长期 flower stem elongation stage 从大蒜花芽分化至花茎形成和采收的时期。 1 DB37/T 3791—2019 3 数据要求 3.1 数据源 卫星数据应源自携载红光波段和近红外波段传感器的中高空间分辨率（小于或等于30 m）卫星，常 用的中高空间分辨率星载传感器及其红光和近红外波段参数见附录A。 3.2 数据预处理 在NDVI计算前，应对卫星数据做辐射定标、大气校正、几何精校正、投影转换及云检测处理。 4 评估时间 幼苗期长势评估时间在12月中旬，花芽和鳞芽分化期长势评估时间在3月下旬，花茎伸长期长势评 估时间在4月中旬。 5 评估方法 5.1 单时次归一化差值植被指数（NDVI） 计算单时次卫星数据的NDVI，见公式（1）： I NDVI  RNIR  RRED ................................... (1) RNIR  RRED 式中： INDVI ——某个大蒜像元单时次的NDVI； RNIR ——近红外波段的反射率； RRED ——红光波段的反射率。 5.2 NDVI 最大值合成 计算某个大蒜像元各时次的NDVI，按第4章规定的评估时间，选取长势评估日前20天中的观测最大 值作为该像元多时次合成后的值，见公式（2）： I NDVI (i)  max( I NDVI (i,1), I NDVI (i, 2), , I NDVI (i, t )) ...................... (2) 式中： INDVI(i) i INDVI(i,t) t ——第i个大蒜像元合成后的NDVI； ——评估区域内大蒜像元序号； ——第i个大蒜像元第t时次的NDVI； ——观测总时次，t≥2。 5.3 NDVI 均值 计算评估区域内所有大蒜像元NDVI最大值合成后的算术平均值，见公式（3）：  m I NDVI   I i / m，i  r .................................. (3) i 1 2 DB37/T 3791—2019 式中：  I NDVI ——评估区域NDVI均值； m ——评估区域内大蒜像元总数； i Ii r ——评估区域内大蒜像元序号； ——评估区域内第i个大蒜像元合成后的NDVI； ——评估区域代码。 5.4 NDVI 多年平均值 计算评估区域大蒜某生育期最近5年以上所有年份的NDVI平均值，见公式（4）：  I NDVI (r , s )  式中： I NDVI (r , s) 1 N  I NDVI (r , n, s) .............................. (4) N n1  r s N ——评估区域r在生育期s的NDVI多年平均值； ——评估区域代码； ——大蒜生育期代码； ——统计总年份； n ——年份序号； I NDVI (r, n, s) ——评估区域r在第n年生育期s的NDVI。 5.5 NDVI 距平值 计算评估区域当年大蒜某生育期NDVI与该生育期NDVI多年平均值的差值，见公式（5）：  I NDVI (r, n, s)  I NDVI (r, n, s)  I NDVI (r, s) .......................... (5) 式中： I NDVI (r, n, s) ——评估区域r在第n年生育期s的NDVI距平值； r s ——评估区域代码； n ——年份序号； I NDVI (r, n, s)  I NDVI (r , s) ——评估区域r在第n年生育期s的NDVI； ——大蒜生育期代码； ——评估区域r在生育期s的NDVI多年平均值。 5.6 NDVI 标准差 计算评估区域大蒜某生育期NDVI多年平均值的标准差，见公式（6）：   1 N 1 N  I NDVI (r ,n ,s )  n 2 ......................... (6) 1 式中：  ——NDVI标准差； N n ——统计总年份； I NDVI (r, n, s) ——评估区域r在第n年生育期s的NDVI距平值。 ——年份序号； 3 DB37/T 3791—2019 5.7 长势判定条件 5.7.1 利用大蒜幼苗期、花芽和鳞芽分化期和花茎伸长期的 NDVI 距平值与标准差综合判定大蒜长势， 将其划分为好、中、差 3 个等级，判定条件如下： a) b) 若 I NDVI (r , n, s) ＞  ，则判定为长势好； 若-  ≤ I NDVI (r , n, s) ≤  ，则判定为长势中； 若 I NDVI (r , n, s) ＜-  ，则判定为长势差。 5.7.2 大蒜长势可分为好、中、差，大蒜长势等级地面调查评定指标见附录 B。 c) 6 评估流程 大蒜长势卫星遥感评估应按下列流程操作： a) 读取预处理后的卫星数据； b) 计算单时次 NDVI 和最大值合成值； c) 生成区域 NDVI 均值数据集； d) 计算幼苗期、花芽和鳞芽分化期和花茎伸长期的 NDVI 多年平均值、NDVI 距平值及标准差； e) 按大蒜长势判定条件，评估幼苗期、花芽和鳞芽分化期和花茎伸长期大蒜长势，基于 GF-1/WFV 卫星数据的大蒜长势评估示例参见附录 C； f) 在研究区大蒜种植分布图基础上，按第 4 章规定的评估时间，分别制作幼苗期、花芽和鳞芽分 化期和花茎伸长期大蒜长势专题图。 4 DB37/T 3791—2019 AA 附 录 A （规范性附录） 常用的中高空间分辨率星载传感器及其红光和近红外波段参数 表A.1列出了常用的中高空间分辨率星载传感器及其红光和近红外波段参数。 表A.1 常用的中高空间分辨率星载仪器的红光和近红外波段参数 星载传感器 GF-1/WFV GF-1/PMS GF-6/WFV GF-6/PMS HJ-1/CCD Landsat/TM Landsat/ETM+ Landsat/OLI CBERS/CCD 波段 波段范围 μm 波段描述 星下点分辨率 m 3 0.630～0.690 红光（Red） 16.0 4 0.770～0