叶面积指数遥感产品真实性检验

叶面积指数遥感产品真实性检验

国家标准 GB/T40034一2021 叶面积指数遥感产品真实性检验 Validationofleafareaindexremotesensingproduets 2021-04-30发布 2021-11-01实施 国家市场监督管理总局 发布 国家标涯花管理委员会国家标准
GB/T40034一2021 目 次 前言 范围 2 规范性引用文件 术语和定义 .中 .中 基本要求 检验方法 5.1检验方法选择 5.2直接检验法 5,2.1地面抽样方法 5.2.2 检验流程 5.3间接检验法 5.3.1基于地面测量与高分辨率遥感数据的多尺度逐级检验法 基于已知质量精度叶面积指数遥感产品的交叉检验法 53.2 检验报告 6.1封面信息 6,2正文信息 6.2.1待检叶面积指数遥感产品概述 参考对象描述 6,2.2 6.2.3检验方法及流程 6.2.4真实性检验结论 附加信息 6.2.5 6.3检验报告信息简表 附录A资料性附录)分层抽样方法 附录B(资料性附录叶面积指数地面测量方法 附录c(资料性附录森林参数测量方法 附录D(资料性附录空间代表性评价方法 12 附录E资料性附录准确度评价指标及计算方法 l6 参考文献
GB/T40034一2021 前 言 本标准按照GB/T1.1一2009给出的规则起草 本标准由科学院提出 本标准由全国遥感技术标准化技术委员会(SAC/TC327)归口 本标准起草单位:科学院空天信息创新研究院、福州大学、科学院青藏高原研究所、科 学院地理科学与资源研究所,北京大学、北京师范大学、科学院大学、林业科学研究院资源信息 研究所、农业科学院农业资源与农业区划研究所、西南交通大学、华中农业大学 本标准主要起草人;李静,赵静、邹杰、曾也鲁、柳钦火、李新、方红亮、唐伯惠、范闻捷、屈永华 穆西哈、姜小光、陈尔学吴文斌、董亚冬、王新鸿、刘照言、尹高飞、徐保东
GB/T40034一2021 叶面积指数遥感产品真实性检验 范围 本标准规定了叶面积指数遥感产品真实性检验的基本要求、检验方法和检验报告格式 本标准适用于叶面积指数遥感产品的真实性检验 规范性引用文件 下列文件对于本文件的应用是必不可少的 凡是注日期的引用文件,仅注日期的版本适用于本文 件 凡是不注日期的引用文件,其最新版本(包括所有的修改单)适用于本文件 GB/T362962018遥感产品真实性检验导则 GB/T39468一2020陆地定量遥感产品真实性检验通用方法 术语和定义 GB/T362962018和GB/T39468一2020界定的以及下列术语和定义适用于本文件 3.1 index;LA 叶面积指数leafarea 地面单位投影面积内叶片总面积的一半 3.2 有效叶面积指数efteetieleafareaindex 与观测孔隙率值相等的、假定叶片随机分布条件下的叶面积指数 注:可通过光学仪器法测量或遥感反演间接获得 3.3 聚集指数clumpingindes 描述自然界中叶片在冠层内随机分布、规则分布、聚集分布的空间分布模式 3.4 基本抽样单元elementarysamplingmit;Esu 在地面观测场布设的进行叶面积指数测量的区域单元 3.5 无线传感器网络wirelesssensornetworks;wSN 通过无线通信技术把传感器节点以自由式进行组织与结合进而形成的网络 基本要求 叶面积指数遥感产品真实性检验应符合GB/T36296一2018中第7章的规定,并满足下列要求 检验对象为单天叶面积指数遥感产品时,参考对象与检验对象的获取时间差异宜在2d以内; a 检验对象为叶面积指数多天合成产品时,参考对象获取时间应在叶面积指数遥感产品的合成 时间内
GB/T40034一202 b 参考对象应优先选择均质地表,可包含非均质地表; 参考对象应尽量覆盖待检叶面积指数遥感产品的不同取值范围 5 检验方法 5.1检验方法选择 叶面积指数遥感产品的真实性检验方法包括直接检验法和间接检验法 根据参考对象与待检对象 的尺度差异,以及地表空间异质性特征,选择检验方法 当待检叶面积指数遥感产品与地面测量数据尺度一致,或待检像元为均质地表像元时,应采用 a 直接检验法; 当待检叶面积指数遥感产品像元尺度大于地面测量数据尺度,且待检像元为非均质地表像元 b 时,应基于地面测量与高分辨率遥感数据,采用间接检验法中的多尺度逐级检验法 当检验区无地面测量数据时,应基于已知质量精度叶面积指数遥感产品,采用间接检验法中的 交叉检验法 5.2直接检验法 5.2.1地面抽样方法 5.2.1.1样区的地面抽样方法 在开展叶面积指数遥感产品真实性检验时,样区内应布设多个开展地面测量的基本抽样单元 EsU)样点 样区内样点选取采用的地面抽样方法如下 a 样区选择应具有空间代表性,应覆盖待检叶面积指数遥感产品3×3个像元,即正方形边长为 待检叶面积指数遥感产品空间分辨率3倍 按照GB/T39468一2020中4.2的方法进行样区异质性评价; b 样区均质时,宜采用简单随机抽样模型和系统抽样模型布设EsU样点,并计算EsU样本量， 具体方法参见GB/T39468一2020中B.1和B.2: 样区非均质，且样区内土地利用类型、植被类型、植被长势分布[如归一化差值植被指数 d NDVI)分布图]等先验知识未知时,宜采用克里格模型和非均质表面均值估计模型布设Esu 样点,并计算EsU样本量,具体方法参见GB/T39468一2020中B.4和B.5; 样区非均质,但样区内土地利用类型、植被类型、植被长势分布等先验知识已知时,考虑经济因 素限制样本量,结合样区内非均质特征选择抽样方法,具体方法参见附录A的A.1和A.2,并 计算EsU样本量; fD 样区非均质并采用无线传感器网络(wSN)进行多点连续测量时,应考虑WSN在不同植被类 型、不同长势以及时间序列上的代表性,宜采用针对wsN定点连续测量的抽样方法布设Esu 样点,具体方法参见A.3,并计算EsU样本量 5.2.1.2EsU内部测量位置选取 根据EsU内部植被群落的空间结构非均质性特征,测量点位置确定方法如下 在植被分布均质、密集的地区,宜采用方形抽样模式、十字抽样模式和萎形抽样模式,ESU内 a 部样本量确定方法参见GB/T394682020中B.1.2 方形抽样模式测量位置如图1所示,应 在抽样正方形的4条边及中心布设部分测量样点,在正方形内部随机布设剩余测量样点;十字 抽样模式测量位置如图2所示,沿着图中十字方向均匀布设测量样点;菱形抽样模式测量位置 如图3所示,沿着图中菱形的4条边均匀布设测量样点
GB/T40034一2021 b)在植被分布稀疏不连续的地区,宜采用对角线抽样模式,测量位置如图4所示,其中虚线表示 在对角线上以一定的距离间隔进行样线测量 图1方形抽样模式示意图 图2十字抽样模式示意图 图3萎形抽样模式示意图 + 图4对角线抽样模式示意图 5.2.2检验流程 检验流程应符合GB/T362962018中8.1的规定,具体操作流程见图5,步骤如下 叶面积指数地面测量: a 可采用点观测或EsU方式观测,ESU内测量样点的确定方法见5.2.1.2,取EsU内所有 样点测量值的算术平均作为EsU的相对真值 测量点上不同植被类型的叶面积指数地面测量方法参见附录B,森林参数测量方法参见 附录C; b)空间代表性评价:若地面观测与待检叶面积指数遥感产品尺度不匹配,应使用空间代表性评价
GB/T40034一202 方法(参见附录D)分析地面测量样区在待检产品像元内的空间代表性,选取空间代表性高的 地面测量数据进行真实性检验; 准确度评价:按GB/T36296一2018中6.1规定的平均误差、平均绝对误差、相对误差、平均相 对误差、平均绝对相对误差、,均方根误差、相关系数,以及相对均方根误差和产品精度(参见附 录E)等准确度评价指标定量表达叶面积指数遥感产品的准确度 不确定度分析;按GB/T36296一2018中6.2规定的标准差、,方差、,协方差和标准不确定度等不 d 确定度评价指标进行叶面积指数遥感产品的不确定度分析 叶西积指着地面满路 待检叶面积指 空间代表性课价 数遥感产品 准确度评价 不确定度分析 图5叶面积指数遥感产品直接检验操作流程 5.3间接检验法 5.3.1基于地面测量与高分辨率遥感数据的多尺度逐级检验法 检验流程应符合GB/T38296一2018中8.1和8.2.1的规定,操作流程见图6,具体步骤如下 样区选择:样区长和宽应是待检叶面积指数遥感产品分辨率3倍以上,例如待检产品空间分辨 a 率为1km,则样区至少为3kmx3km. b 叶面积指数地面测量: 宜采用EsU,其大小应与高分辨率遥感数据的空间分辨率一致,EsU位置获取见5.2.1.l; 1 ESU内测量样点的确定方法见5.2.1.2; 3 测量点不同植被类型的叶面积指数地面测量方法参见附录B,森林参数测量方法参见附 录C; 4 取多点测量的算术平均作为ESU的真值; 高空间分辨率叶面积指数相对真值图获取 随机选取2/3的ESU观测,使用高分辨率遥感影像,基于高空间分辨率植被指数[如ND VI、比值植被指数(SR、简化比值植被指数(RSR)等]建立叶面积指数的经验回归模型， 或者基于植被辐射传输模型的叶面积指数反演方法,反演得到高空间分辨率叶面积指数 相对真值分布图; 基于剩余1/3的ESU观测,按照5.2规定的直接检验方法,对反演的高空间分辨率叶面 积指数相对真值分布图进行准确度和不确定度分析,当且仅当高空间分辨率叶面积指数 相对真值结果满足既定阂值要求时,可用于下一步检验; d 待检叶面积指数遥感产品相对真值获取;通过几何配准,选取与待检产品像元对应的高空间分 辨率叶面积指数相对真值图的像元,用逐像元平均的方法获取待检产品像元的相对真值; 准确度评价;按GB/T36296一2018中6.1规定的平均误差、平均绝对误差,相对误差、平均相 对误差、平均绝对相对误差、均方根误差、相关系数,以及相对均方根误差和产品精度(参见附 录E)等准确度评价指标定量表达叶面积指数遥感产品的准确度;
GB/T40034一2021 fD 不确定度分析;按GB/T36296一2018中6.2规定的标准差,方差、协方差和标准不确定度等不 确定度评价指标进行叶面积指数遥感产品的不确定度分析 样区选择 叶面积指数地面测量 高空间分辨率叶面积 指数相对真值获取 待检叶面积指 待检叶面积指数遥感 数遥感产品 产品相对真值获取 准确度评价 不确定度分析 图6基于地面测量与高分辨率遥感数据的多尺度逐级检验法操作流程 5.3.2基于已知质量精度叶面积指数遥感产品的交叉检验法 检验流程应符合GB/T362962018中8.2.1的规定,主要操作流程见图7,具体步骤如下 参考对象选取;宜选取精度及不确定度已知的叶面积指数遥感产品作为参考对象 aa b) 空间尺度匹配:将待检叶面积指数遥感产品与参考对象进行空间尺度配准,采用逐像元平均法 将待检产品与参考对象聚合到相同尺度(两产品原始尺度的最小公倍数); 选取参与评价的像元集合:为全面评价不同空间特征的叶面积指数遥感产品质量,应尽量选取 不同植被类型,不同均质性特征以及不同地形特征的像元,应根据质量描述选取参考对象中反 演质量较高的像元,得到开展质量评价的像元集合 准确度评价;按GB/T36296一2018中6.1规定的平均误差、平均绝对误差,相对误差、平均相 对误差、平均绝对相对误差、均方根误差、相关系数,以及数据缺失百分比、相对均方根误差和 产品精度(参见附录E)等准确度评价指标定量表达叶面积指数遥感产品的准确度 不确定度分析;按GB/T36296一2018中6.2规定的标准差,方差协方差和标准不确定度等不 确定度评价指标进行叶面积指数遥感产品的不确定度分析
GB/T40034一2021 待检叶面积指 参考叶面积指数产" 品选取 数遥感产品 空间尺度匹配 选取参与评价的像元 集合 准确度评价 不确定度分析 图7叶面积指数遥感产品交叉检验操作流程 6 检验报告 6.1封面信息 检验报告封面应包括以下信息: 检验报告编号 a b) 检验报告名称 检验负责人; 山 检验核对人: 检验签发人; D 检验单位 送检单位; 8 h检验时间 6.2正文信息 6.2.1待检叶面积指数遥感产品概述 对待检叶面积指数遥感产品进行描述,宜包括 遥感产品来源,时空分辨率,时间序列跨度,空间覆盖范围投影方式、单位等; a b)产品算法简单描述 6.2.2参考对象描述 对作为参考对象的检验数据集进行描述,宜包括 检验数据集的常规信息,包括名称,植被类型、地形特征(坡度、坡向入,非均质特征、测量方达、 a 测量仪器、数据格式,采集时间,采集负责人等 b检验数据集的质量评价描述; 检验数据集的适用性描述 c 6.2.3检验方法及流程 对采用的检验方法和检验过程进行描述,宜包括
GB/T40034一2021 检验方法和流程概述; a b 检验结果的评价指标; c 检验过程的记录要求; d)检验结果的存档 6.2.4真实性检验结论 对真实性检验结果进行描述,宜包括 真实性检验结果总体评价:描述所检验叶面积指数遥感产品的总体精度; a b 分项指标评价;描述直接检验,间接检验的精度评价结果;描述不同植被类型,不同生长期、不 同非均质地表像元以及不同地形特征的精度评价结果;描述产品空间连续性,如果待检产品为 时间序列产品,评价产品时间连续性; 分析检验过程中的不确定性:包括测量误差、仪器误差、代表性误差等 c d)对产品的评价与建议:通过对各项指标的评价和分析,给出所检验叶面积指数遥感产品的具体 评价结果,包括产品精度、时间连续性、空间连续性等 检验信息;包括检验人、检验单位、检验时间,地点、签字,盖章 6.2.5附加信息 对叶面积指数遥感产品真实性检验过程中的非常规问题进行说明与描述 6.3检验报告信息简表 叶面积指数遥感产品检验报告信息简表的编制参见GB/T362962018附录D.
GB/T40034一2021 附 录 A 资料性附录) 分层抽样方法 A.1基于植被分类图先验知识的抽样方法 该方法适用于检验区域内植被类型较多,特定植被类型内的生长差异不显著,如农田区域 具体步 骤如下 依据植被分类图,将待检区域按植被类型进行分层; a b) 各层样本布设方法及样本量计算参见GB/T39468一2020中B.3分层抽样模型 A.2基于植被指数先验知识的抽样方法 该方法适用于测区内植被类型单一但生长差异明显,如草地、森林地区,且分类图提供的先验知识 有限的区域,具体步骤如下 选取检验区域,获取时相接近的高分辨率遥感影像，经数据预处理(辐射校正、几何校正),生成 a 高分辨率植被指数图及植被分类图; b 选取对叶面积指数值变化较敏感的植被指数,例如低覆盖度区域可选择NDVI,高覆盖度区域 可选择SR,根据高分辨率影像的植被指数频率分布直方图,以各层内部的离散方差加权和取 最小作为约束条件设定阔值,将研究区域分割为若干个不同植被类型与生长水平的类别 根据植被指数分割结果图各层离散方差的大小评估各层内部的变异程度,同时考虑各层的离 散方差与面积比例,采取内曼分配原则确定各层应分配的样点数目; d 根据植被指数图进行布样时,采用最邻近指数(NNI)对样点的空间分布格局进行评估,通过对 照高分辨率植被分类图,使样点在不同的植被类型中都有分布,最终确定样点空间位置 A.3针对WSN定点连续测量的抽样方法 对于wsN,宜采用联合植被分类图与多时相的植被指数图作为先验知识的抽样方法进行布样 该 方法采用多时相的高分辨率植被指数图,目标丽数设定为所选取的样点在不同时相上的植被指数频率 分布与检验区域的差异最小,使样点在检验区域不同时相上具有最高代表性 植被指数图的时相数根 据检验区域异质性在时间序列上的变化特征决定农作物、草地、落叶林地区应至少包含春季、夏季和秋 季3个时相,常绿针叶林地区应至少包含春季和夏季2个时相
GB/T40034一2021 附录 B 资料性附录 叶面积指数地面测量方法 B.1直接测量法 当抽样点个数相对较少,投人实验的人力充足,允许较长的测量周期时,采用直接测量法进行植被 叶面积指数测量,具体包括三种方法 收获测量法;适用于冠层结构较小植株(如农作物、草地等)的叶面积指数测量 通常选取边长 a 为1m或2m的正方形区域,收获区域内所有叶片,采用照相法或叶片面积测量仪器测量收 获的所有叶片的总面积,将叶片总面积除以区域面积计算单位面积的叶面积指数 b 落叶收集法;适用于落叶林 在样区内随机设置一定面积的凋落物收集网,将收集到的凋落物 烘干,分离出叶片称重,得到落叶量 用十字分割法从落叶中取出一定重量的叶片测出总叶面 积,计算出比叶面积,结合落叶收集得到的单位时间单位面积落叶的重量以及单位时间落叶量 所占样地总叶量的百分比,计算得到叶面积指数 异速生长测定法;适用于橡树等叶面积指数与植物器官的几何尺寸相对稳定的特定森林树种 利用不同生长周期植物器官的儿何尺寸,根据其生物量与植物器官的异速生长规则,建立适用 于该植被相应生长时期的异速生长函数,用基本的森林清查数据作为输人,从而计算地面叶面 积指数 B.2采用光学仪器的间接测量法 采用光学仪器的间接测量法按测量原理分为基于辐射测量的方法和基于图像测量的方法 前者通 过测量辐射透过率来计算叶面积指数,测量应在晴朗无云的天气进行;后者通过获取和分析植物冠层的 半球数字图像来计算叶面积指数,测量应在均一的光环境下进行,如黎明、黄昏、阴天等,避免晴天时的 直射光 采用光学仪器进行测量前,应针对测量的植被类型,对仪器进行定标 针对3种植被类型的叶面积指数测量方法如下 森林叶面积指数测量:对多个森林参数进行测量,包括有效叶面积指数、聚集指数、束簇面积比 a 参数,木质总面积比参数,测量方法参见附录c 有效叶面积指数采用基于辐射测量的仪器测 量;测量聚集指数时,应在10"一70"天顶角范围内多次测量并取平均; 农作物叶面积指数测量;采用基于辐射测量的方法测量有效叶面积指数 对于垄行结构明显 b 条对角样线,样点在样线上均质分布 如果垄行中存在明 的行播作物冠层,在两垄行之间拉 显的空隙,应先进行冠层空隙检测 当空隙检测误差小于或等于10%(均质冠层)时,可根据 天气条件选取视角尽可能大的遮盖帽;当空隙检测误差大于10%(非均质冠层)时,应采用45" 遮盖帽,且对角样线的数目应至少为原来的2倍; 草地叶面积指数测量;采用基于辐射测量的方法测量有效叶面积指数 当草的高度较低,传感 器放置在地表会接近或超过草的高度时,采用基于图像测量的方法,以向下测量的方式获取叶 面积指数;当草的高度较高,并夹有枯草时,采用附录C中c.4森林叶面积指数测量估计木质 总面积比参数的方式,剔除树枝、树干等非叶片成分和枯草部分的贡献
GB/T40034一2021 附 录 资料性附录) 森林参数测量方法 C.1有效叶面积指数测量方法 采用基于辐射测量的方法测量叶面积指数,也可采用基于图像测量的方法 测量时按照仪器说明 规定的天气条件,根据5.2.1给出的抽样方案对各样点逐一测量 将传感器置于林下植被的下方时,测 量到的是乔木层和林下植被的叶面积指数之和;将传感器置于林下植被的上方时,测量到的是乔木层的 叶面积指数 C.2冠层基本组分及木质组分聚集指数测量方法 冠层基本组分主要测量方法包括样线测量法,半球摄影法和多光谱冠层成像法 测量时按照仪器 说明规定的天气条件,在10°70"天顶角范围内多次测量并取平均,以消除观测天顶角及方位角的影 响 当需要快速测定冠层基本组分聚集指数时,宜采用样线测量法;当需要高精度冠层基本组分聚集指 数测量时,宜采用多光谱冠层成像法 冠层木质组分聚集指数在天顶角方向呈较大的差异性,采用木质总面积比参数消除木质组分对森 林冠层地面叶面积指数测量结果的影响 C.3束簇面积比参数测量方法 测量针叶林的聚集指数时,需基于束簇面积比对测量的聚集指数进行校正 针叶束表面积测量方 法主要有光学投影法和体积替代法 光学投影法基于针叶束投影面积和形状因子计算针叶簇中所有针 叶束表面积;体积替代法将针叶簇投人含有5%洗洁精的清水中,通过两次测量清水的体积变化计算针 叶束总体积,同时结合针叶束数量和形状因子计算针叶束总表面积 测量方法的选择如下 当样地内分布不同年龄的针叶林树木时,其束簇面积比参数测量应在每个树龄段选择典型针 a 叶树,同时在典型树木各高度及方位区间选择典型针叶簇,数目应不少于12个,以避免测量人 员主观性对最终结果的影响; 在设备及人员条件具备时,优先采用体积替代法开展针叶束表面积测量; b c 当采用光学投影法量时,如针叶簇与标准圆柱体儿何形状差异较大,应增加光学投影在天顶 角方向及方位角方向的抽样次数,以实现对针叶簇全方位抽样并保证测量精度 C.4木质总面积比参数测量方法 木质总面积比参数可采用的测量方法主要有直接测量方法、光学测量方法、多光谱冠层成像法 对 于常绿针叶林、阔叶林,优先采用多光谱冠层成像法 测量方案如图c.1,在天顶角方向0"90°内可每 隔10'测量一次,在方位角方向可每隔60"测量一次,共计60对影像即可完成冠层半球方向抽样 10
GB/T40034一2021 观测方向 00 ,2O 303 40 50" 602 07” 80" 90" 天顶角方向观测 方位角方向观测 图C.1木质总面积比参数影像获取示意图 11
GB/T40034一2021 附 录 D 资料性附录) 空间代表性评价方法 基于高分辨率影像计算低分辨率像元内统计指标的空间代表性评价方法 D.1 评价站点观测在低分辨率像元尺度内空间代表性需要明确低分辨率像元内部空间异质性特征 基 于高分辨率影像能够捕获低分辨率像元内部地物空间分布特征的优势,计算低分辨率像元内部多个高 分辨率像元的统计指标,对低分辨率像元内部地表空间异质性开展定量评价 指标定义如下 相对平均绝对偏差,基于多个高分辨率像元评价低分辨率像元内地表空间异质性,由公式 a D.1)表示: 一" RMAD兰 -×100 D.1) 式中 RMAD 相对平均绝对偏差; 低分辨率像元内高分辨率像元个数 低分辨率像元内第i个高分辨率像元值 ， 低分辨率像元内所有高分辨率像元值得算术平均值 相对空间抽样误差,用于评价站点观测值和多个高分辨率影像像元均值的差异,用于表征站 b 点观测值是否落人异常值内由公式(D.2)表示 .r，一工 (D.2 RSSE ×100 式中: RsSE 相对空间抽样误差; -站点所处位置对应的高分辨率像元值; 低分辨率像元内所有高分辨率像元值得算术平均值 RMAD值越小,代表地表空间异质性程度越低;RSSE值越小,代表站点所处位置落人低分辨率像 元内极端异常值的概率越小 在实际应用中,需要同时保持RMAD和RSSE处于较小值才能保证站点 观测的空间代表性越高 D.2基于高分辨率影像和变异函数的空间代表性评价方法 基于地理学第一定律,考虑到地表空间相关性,可采用变异函数结合低分辨率像元内多个高分辨率 像元评价地表的空间异质性 变异函数由公式(D,3)表示 [c(.r.一c(r,] (D.3 r(h= -Nw,.习 式中 r( 变异函数; 低分辨率像元内高分辨率像元对(r.))的距离; Nh 低分辨率像元内高分辨率像元对(r.，rp)的个数; 12
GB/T40034一2021 低分辨率像元内在r 位置处高分辨率像元值; x(.xa" 低分辨率像元内在工位置处高分辨率像元值 c(.r 基于变异丽数拟合模型(球面模型、指数模型等)拟合低分辨率像元内不同高分辨率像元对在不同 距离的变异丽数值 模型的选择取决于对变异函数值的精度,以常用的球面拟合模型为例,采用拟合模 型后的参数定量评价低分辨率像元内多个高分辨率像元间和高分辨率像元内的空间异质性 球面拟合 模型由公式(D.,4)表示 l+会-a(门<" 4 D.4 式中: 球面拟合模型的变异函数; rh(h 低分辨率像元内高分辨率像元对(r ,r)的距离 拟合模型的块金值,衡量高分辨率像元内部的空间异质性; C0 -拟合模型的基台值,衡量低分辨率像元内部所有高分辨率像元间的空间变异性; 拟合模型的变程,超过变程距离的像元之间没有相关性 通过基台值和块金值来确定低分辨率像元内空间均质程度,块金值越小,说明高分辨率像元内部异 质性较小,高分辨率像元空间尺度能够表征低分辨率像元内最小分辨单元;基台值越小,说明低分辨率 像元尺度内空间均质性越强 D.3基于高分辨率影像和地统计学指标的空间代表性评价方法 借助高空间分辨率影像,利用相对变异系数指标衡量站点在像元尺度内的整体可变性,由公式 D.5)表示 CV.5 R (D.5 式中 相对变异系数 Re CV 高分辨率参考图在像元尺度1.5、内的变异系数; 1.5 CV 高分辨率参考图在像元尺度、内的变异系数; 代表像元尺度大小 利用尺度需求指数指标衡量站点测量可代表的空间范围,由公式(D.6)表示 -、鹰 D.6 RE=exp 式中: -尺度需求指数 RsE -观测数据的空间尺度; g 在像元尺度内高分辨率参考图的变程 a 在像元尺度1.5s内高分辨率参考图的变程 a1.s 利用空间相关性的相对强度指标衡量像元尺度内的空间相关强度,由公式(D.7、公式(D.8)表示 ST， ST 1. D.7 Rsr ST Ye(a Co ST一 D.8 yE(a)一E(O 13
GB/T40034一2021 式中 R -空间相关性的相对强度; ST -空间相关性的强度; YE(a) -对应像元尺度内高分辨率参考图的距离相隔为变程a的半变异函数值; YE(0) -对应像元尺度内高分辨率参考图的距离相隔为0的半变异函数值 -拟合模型的块金值,衡量高分辨率像元内部的空间异质性 c 利用结构变化的相对比例指标衡量像元尺度内的空间结构,由公式(D.9)、公式(D.10)表示 SV SV 5 Rs D.9 SV [YE(h SV dh D.10 式中 R、 -结构变化的相对比例 SV -结构变化的比例 低分辨率像元内高分辨率像元对(r ,工)的距离; YE(h -对应像元尺度内高分辨率参考图的距离相隔为h的半变异函数值; -对应像元尺度的高分辨率参考图的基台值 c 当球状模型可以对半方差函数进行拟合时,可基于Rcv,RsE,Rs和Rs四个指标共同构建的标准 分数指标对站点的空间代表性进行评估;当球状模型不适合对半方差函数进行拟合时,采用基于Rev构 建的一阶分数指标对站点的空间代表性进行评估,ST和RAw分别由公式(D.ll)和公式(D.12) 表示 Rel+Rl|Rs ST D.11 十RsE score SV RAW |2Rc- D.12) score 式中 ST -标准分数; score RAwe -阶分数; ST和RAw的分数越高,说明站点的空间代表性越强 D.4基于高分辨率影像和地表分类图的空间代表性评价方法 借助于高空间分辨率影像,利用相对绝对误差(RAE)指标衡量站点测量对应的高分辨率像元与高 分辨率影像在产品像元尺度内聚合后的值的差异,由公式(D.13)表示 |<(r,)一[r(s)] RAE(s)三 (D.13 ×100 E[r(s] 式中 RAE(s) -像元尺度为s的相对绝对误差; 这(.r， -站点位置对应的高分辨率参考图的值; E[r(s门 -高分辨率参考图在站点位置周边像元尺度、内的算术平均值 借助地表覆盖分类图,计算像元尺度内主要植被类型比例(DVTP)指标定量描述像元尺度的空间 异质性,由公式(D.l4)表示: As DVTP(s) -×100 D.14) ( 习n,( 14
GB/T40034一2021 式中: DVTP(s 像元尺度为、的主要植被类型比例 A(s 站点测量的植被类型在产品像元尺度的面积 像元尺度为、的第i类植被类型的面积 n,(s N(s -像元尺度为、的总地物类别个数 利用基台系数(CS)指标判断像元尺度内的整体空间异质性程度,由公式(D.15)表示 Gac CS(s= ×100 D.15 4(s 式中: CS(s 像元尺度为、的基台系数; e,()+c() 在像元尺度s内高分叶面积指数参考图的的基台值,基台值通过计算区域变异 函数也叫半方差函数)时采用球面模型拟合得到 -像元尺度、内高分叶面积指数参考图的均值 (s) 根据DVTP,RAE和CS三个指标将站点叶面积指数测量在产品像元尺度内空间代表性分为五 级,具体分级规则见表D.1 表D.1站点叶面积指数测量空间代表性质量分级规则 空间代表性指标 代表性 质量等级说明 等级 VIP>60%RAE<阔值cs<阔值 质量最好,站点叶面积指数测量能够代表像元尺度,代表性程 度最高 站点叶面积指数测量与产品像元尺度内均值较一致,但像元内 变异程度较大,代表性稳定程度不高 产品像元尺度内较均一,但站点叶面积指数测量处在像元内异 常点处,代表性程度较低 产品像元内地表类型单一,但植被长势差异较大,站点叶面积 指数测量的代表性程度更低 产品像元内地表类型多样,站点叶面积指数测量代表性程度 最低 注“、/”表示满足指标条件;“义”表示不满足指标条件;“一”表示不考虑 值大小依据代表性等级确定,其取值范围0级为00.l(含),1级为0.l0.2(含),2级为0.20.3(含),3级 为大于0.3. DVTP的阔值设定根据中等分辨率成像光谱仪MODIS)地表分类产品MCD12Q1像元分类规 则,以像元内某一类地物比例达到60%进行像元类别确定 当DvTP大于60%,即空间代表性为 0级3级,则站点测量的植被类型对产品像元尺度的植被类型具有空间代表性;当DVTP小于或等于 0%.,即空间代表性为4级,别站点渊量的体一植被类狸不能代衣整个广品像元. 15
GB/T40034一2021 录 附 资料性附录) 准确度评价指标及计算方法 E.1数据缺失百分比 由公式(E.1)表示 PMD E.1 ×100% N 式中 PMD 数据缺失百分比 N 反演缺失像元个数; N -总反演像元个数 E.2相对均方根误差 由公式(E.2)表示: 习[LAI,(i)一LAI,(i)]/N RRMS E.2 LAI,(i)/N 式中 RRMs 相对均方根误差; N -参与检验的像元数目: LAI,(i) 待检产品第i个检验像元的叶面积指数值; LAI,(i 基准产品第i个检验像元的叶面积指数值 E.3产品精度 由公式(E.3)表示: [LAI,()一LAI(i)]/N PP=[1 -]×100% 习LAI(/N E.3 .( 式中 PP 产品精度 参与检验的像元数目; LAI,(i) 待检产品第i个检验像元的叶面积指数值; LAI,(i 基准产品第i个检验像元的叶面积指数值 16
GB/T40034一2021 参 考文献 [1]GB/T3358.2一2009统计学词汇及符号第2部分;应用统计 [[2]GB/T149502009摄影测量与遥感术语

和叶面积指数遥感产品真实性检验类似的标准

GB/T40034-2021

叶面积指数遥感产品真实性检验

2022/7/16 2:22:50 现行