(高清版)GBT 40038-2021 植被指数遥感产品真实性检验

Validationofvegetationindexremotesensingproducts国家市场监督管理总局国家标准化管理委员会I 1 1 1 2 2 2 2 4 66.1封面信息 66.2正文信息 76.3检验报告信息简表 7附录A(规范性附录)地面二向性反射率因子测量 8参考文献 9ⅢGB/T40038—2021本标准按照GB/T1.1—2009给出的规则起草。本标准由中国科学院提出。本标准由全国遥感技术标准化技术委员会(SAC/TC327)归口。京大学、中国科学院地理科学与资源研究所。GB/T40038—2021植被指数是两个或多个遥感传感器波段数据的线性或非线性组合而形成的系列指数，可用于诊断植被生长状态、绿色植被活力以及反演多种植被参数。植被指数有多种具体的定义方式，如红光波段和近红外波段反射率比值组合的比值植被指数(RVI)和归一化植被指数(NDVI)、降低土壤背景的土壤调节植被指数(SAVI)、对浓密植被区敏感的增强型植被指数(EVI)等。具有红光和近红外波段的传感器反射率常用于计算植被指数，以此为基础多个全球植被指数产品业已发布。由于数据和算法的限制，植被指数产品不可避免地存在误差，需要对其进行真实性检验以获得该产品的精度和不确定性。以往由于验证的数据、流程、评价标准不同，得到的产品精度不具可比性，给用户定量化应用带来困难。因此，有必要制定科学、规范的植被指数遥感产品真实性检验标准。1GB/T40038—2021植被指数遥感产品真实性检验本标准规定了植被指数遥感产品真实性检验的基本要求、检验方法和检验报告。本标准适用于大气层底植被指数遥感产品的真实性检验。2规范性引用文件下列文件对于本文件的应用是必不可少的。凡是注日期的引用文件，仅注日期的版本适用于本文件。凡是不注日期的引用文件，其最新版本(包括所有的修改版)适用于本文件。GB/T36296—2018遥感产品真实性检验导则GB/T39468—2020陆地定量遥感产品真实性检验通用方法GB/T36296—2018和GB/T39468—2020界定的以及下列术语和定义适用于本文件。3.1双向反射率分布函数物体表面特定方向反射辐亮度的微增量与特定方向入射辐照度的微增量之比。[GB/T36299—2018,定义3.6]3.2双向反射率因子物体表面向特定方向的反射辐亮度与相同辐照度条件下理想漫反射体在该方向上的反射辐亮度之比。[GB/T36299—2018,定义3.7]3.3波段反射率spectralreflectivity物体对某一波段电磁波的反射辐射通量与其入射辐射通量之比。3.4波谱响应函数spectralresponsefunction描述遥感器对不同波长处入射辐射响应率的函数。3.5植被指数vegetationindex一种利用遥感影像不同谱段数据的线性或非线性组合而形成的能反映绿色植物的生长状况和分布的特征指数。[GB/T30115—2013,定义3.11]2GB/T40038—20214基本要求植被指数遥感产品真实性检验应符合GB/T36296—2018第7章的规定，同时：a)参考对象应与待检植被指数遥感产品属于同一种植被指数，对应于相同或相近时期(植被指数时间变化率)的同一区域，并已知其数据质量。b)参考对象应包含两种不同植被覆盖度的典型区域：1)低植被覆盖度区域，如稀疏植被覆盖区域；2)高植被覆盖度区域，如浓密植被覆盖区域。c)参考对象应首选地面测量数据进行植被指数遥感产品的真实性检验，宜选择与待检植被指数遥感产品对应波段和角度一致的地面传感器采集的反射率和植被指数。在没有地面测量数据时，可选用与待检植被指数遥感产品对应波段和角度一致的已检植被指数遥感产品作为参考对象。5检验方法5.1检验方法选择根据地表异质性、地面测量视场和待检植被指数遥感产品像元尺度，选择真实性检验方法：a)当待检植被指数遥感产品像元尺度不超过地面测量视场1.5倍时，应采用直接检验法。b)当待检植被指数遥感产品像元尺度大于地面测量视场1.5倍时，若待检像元区域为均质地表，应采用直接检验法；若待检像元区域为非均质地表，应考虑尺度差异，采用基于地面测量与高分辨率遥感数据的多尺度逐级检验法。c)当检验区无地面测量数据时，采用基于已知精度植被指数遥感产品的交叉检验法。5.2直接检验法5.2.1参考对象的地面抽样方法参考对象地面抽样方法如下：a)按照GB/T39468—2020中4.2的规定进行样区地表异质性评价。b)样区选择应具有空间代表性，宜不小于3×3个待检植被指数遥感产品像元。c)样区地面测量时间与待检植被指数遥感产品时间差应依据植被指数时间变化率确定，通常宜不超过5天，变化较快的地表宜不超过2天。d)按照GB/T39468—2020中4.3的规定进行抽样：对于均质样区，宜在区域内进行随机抽样或系统抽样；对于非均质样区，宜根据植被长势和植被类型在样区内进行分层抽样，使得层内相对均质，并在层内进行局部随机抽样以获取各层的植被指数。检验流程应符合GB/T36296—2018中8.1的规定，主要操作流程见图1,具体步骤如下：a)样本选择：按5.2.1的规定进行地面抽样方法设计，确定样本。b)地面二向性反射率因子测量：测量方法见附录A。c)空间匹配：将地面采样点与待检植被指数遥感产品像元进行空间匹配。d)波谱转换：依据待检植被指数遥感产品对应波段的波谱响应函数，将地面测量的二向性反射率因子转换为与待检植被指数遥感产品对应的波段二向性反射率因子。e)角度一致性判断：判断地面测量的二向性反射率因子与待检植被指数遥感产品角度是否一致。若一致，可直接用于下一步检验；否则，利用地面样本的BRDF特性进行角度归一化，将地面3GB/T40038—2021测量的二向性反射率因子校正为与待检植被指数遥感产品角度相一致的二向性反射率因子。g)植被指数相对真值计算：采用待检植被指数遥感产品定义的波段反射率组合方法计算植被指数相对真值。在植被指数计算过程中有经验系数的，应根据测量区域的特征精确赋值。h)准确度评价：按GB/T36296—2018中6.1规定的准确度评价指标定量表达待检植被指数遥感产品的准确度。i)不确定度分析：按GB/T36296—2018中6.2规定的不确定度评价指标定量表达待检植被指数遥感产品的不确定度。产品测量是是否一图1植被指数遥感产品直接检验操作流程4GB/T40038—20215.3间接检验法5.3.1基于地面测量与高分辨率遥感数据的多尺度逐级检验法检验流程应符合GB/T36296—2018中8.1和8.2.1的规定，主要操作流程见图2,具体步骤如下：a)高分辨率反射率影像数据准备：高分辨率反射率影像空间尺度应与地面测量的二向性反射率因子相同或相近，时间应与待检植被指数遥感产品相同或相近。b)样本选择：按5.2.1的规定进行地面抽样方法设计，确定样本。c)地面二向性反射率因子测量：测量方法见附录A。d)使用地面测量的二向性反射率因子验证高分辨率影像反射率：按5.2.2中c)、d)和e)的规定，将地面测量的二向性反射率因子转换成与高分辨率影像一致的参考反射率，进行高分辨率反射率准确度和不确定度评价。若高分辨率反射率准确度和不确定度满足用户既定阈值要求，可用于下一步检验；否则，重新准备高分辨率反射率影像数据并验证。e)高分辨率影像反射率波谱转换：依据待检植被指数遥感产品对应波段的波谱响应函数，将满足用户要求的高分辨率影像反射率转换为与待检植被指数遥感产品对应的波段反射率。f)角度一致性判断：判断高分辨率影像反射率与待检植被指数遥感产品角度是否一致。若一致，可直接用于下一步检验；否则，根据样本BRDF特性进行角度归一化，将高分辨率影像反射率校正为与低分辨率待检植被指数遥感产品角度相一致的反射率。g)空间尺度转换：按GB/T39468—2020中4.4.1的规定，将经过验证的高分辨率影像反射率升尺度至低分辨率待检植被指数遥感产品相一致的空间尺度。h)植被指数相对真值计算：采用待检植被指数遥感产品定义的波段反射率组合方法计算低分辨率植被指数的相对真值。在植被指数计算过程中有经验系数的，应根据测量区域的特征精确赋值。i)准确度评价：按GB/T36296—2018中6.1规定的准确度评价指标定量表达待检植被指数遥感产品的准确度。j)不确定度分析：按GB/T36296—2018中6.2规定的不确定度评价指标定量表达待检植被指数遥感产品的不确定度。5GB/T40038—2021据准备地面抽样方法设计样本选择地面二向性反射率因子测量否直接检验法高分辨率反射率验证植被指数遥感产品是波谱转换角度一致性?是否角度归一化空间尺度转换不确定度分析图2植被指数遥感产品多尺度逐级检验操作流程5.3.2基于已知精度植被指数遥感产品的交叉检验法检验流程应符合GB/T36296—2018中8.2.1的规定，主要操作流程见图3,具体步骤如下：a)参考植被指数遥感产品选择：宜选用与待检植被指数遥感产品时间相同或相近的已检植被指数遥感产品作为参考对象。b)空间匹配：将参考植被指数遥感产品与待检植被指数遥感产品进行几何匹配。c)波谱一致性判断：依据待检植被指数遥感产品对应波段的反射率波谱特征，判断参考植被指数遥感产品对应波段反射率波谱是否一致。若一致，可进行下一步判断；若不一致，则应依据待检植被指数遥感产品和参考植被指数遥感产品的波谱响应函数，进行波谱转换。d)角度一致性判断：依据待检植被指数遥感产品的角度特征，判断参考植被指数遥感产品对应波段反射率角度是否一致。若一致，可进行下一步判断；若不一致，则应依据参考植被指数遥感6GB/T40038—2021产品对应波段反射率的BRDF特性进行角度归一化，将参考植被指数遥感产品校正为与待检植被指数遥感产品角度相一致的产品。e)空间一致性判断：根据待检植被指数遥感产品像元面积大小，判断参考植被指数遥感产品空间尺度是否一致。若一致，可直接采用参考植被指数遥感产品作为检验相对真值；否则，按GB/T39468—2020中4.4.1的规定进行空间尺度转换，获取与待检植被指数遥感产品像元尺度相一致的参考植被指数遥感产品作为相对真值。f)准确度评价：按GB/T36296—2018中6.1规定的准确度评价指标定量表达待检植被指数遥感g)不确定度分析：按GB/T36296—2018中6.2规定的不确定度评价指标定量表达待检植被指数遥感产品的不确定度。选择空间匹配波谱一致性?是是空问一致性?是否否否图3植被指数遥感产品交叉检验操作流程6检验报告6.1封面信息检验报告封面应包括以下信息：7GB/T40038—2021b)检验负责人；检验核对人；检验签发人；检验单位；送检单位；h)检验时间。6.2正文信息6.2.1待检植被指数遥感产品概述对待检植被指数遥感产品进行描述，宜包括：b)产品主要算法以及算法的特点等。对作为参考对象的验证数据集进行描述，宜包括：b)验证数据集的质量评价描述；c)验证数据集的适用性描述。对采用的检验方法和检验过程进行描述，宜包括：a)检验方法概述；b)检验流程；c)检验结果的评价指标；d)检验过程的记录要求；e)检验结果的存档。对真实性检验结果进行描述，宜包括：a)真实性检验结果总体评价：描述所检验植被指数遥感产品的总体精度。b)分项指标评价：描述所检验植被指数遥感产品在不同区域的精度；描述所检验植被指数遥感产品在不同植被类型的精度。c)分析检验过程中的不确定度。对植被指数遥感产品真实性检验过程中的非常规问题进行说明与描述。6.3检验报告信息简表植被指数遥感产品的检验报告信息简表编制参见GB/T36296—2018附录D。8GB/T40038—2021(规范性附录)地面二向性反射率因子测量A.1.1具有计算植被指数相关波段的成像或非成像设备可用于地面样本二向性反射率因子测量，波谱范围应为380nm～1050nm。依据植被指数的定义，可计算植被指数。A.1.2多角度观测架(可选):用于地面样本二