GB/T 41281-2022 光合有效辐射遥感产品真实性检验

　　ICS 07 . 040 CCS A 77

　　中 华 人 民 共 和 国 国 家 标 准

　　GB/T 41281—2022

　　光合有效辐射遥感产品真实性检验

　　validationofphotosyntheticallyactiveradiationremotesensingproducts

　　2022-03-09 发布 2022-10-01 实施

　　国家市场监督管理总局国家标准化管理委员会

　　

　　发

　　

　　布

　　GB/T 41281—2022

　　目 次

　　前言 Ⅲ

　　引言 Ⅳ

　　1 范围 1

　　2 规范性引用文件 1

　　3 术语和定义 1

　　4 参考对象 2

　　5 基本要求 2

　　6 检验方法 3

　　7 检验报告 6

　　附录 A(资料性) 全球光合有效辐射及总辐射地面观测数据列表 8

　　附录 B(资料性) 空间代表性评价方法 9

　　附录 C(规范性) 观测样点布设及光合有效辐射地面观测 10

　　附录 D(资料性) 时间一致性转换方法 12

　　附录 E(资料性) 光合有效辐射遥感产品真实性检验信息简表样例 14

　　参考文献 16

　　Ⅰ

　　GB/T 41281—2022

　　前 言

　　本文件按照 GB/T 1 . 1—2020《标准化工作导则 第 1 部分：标准化文件的结构和起草规则》的规定起草。

　　请注意本文件的某些内容可能涉及专利。 本文件的发布机构不承担识别专利的责任。

　　本文件由中国科学院提出。

　　本文件由全国遥感技术标准化技术委员会(SAC/TC 327)归口 。

　　本文件起草单位：中国科学院空天信息创新研究院、中国科学院地理科学与资源研究所、中国气象科学研究院、中国资源卫星应用中心、生态环境部卫星环境应用中心。

　　本文件主要起草人：李丽、辛晓洲、张海龙、闻建光、吴骅、王培娟、高彦华、姚延娟、仲波、余珊珊、柏军华、杜永明、柳钦火、刘照言、于江丰、龚力峰。

　　Ⅲ

　　GB/T 41281—2022

　　Ⅳ

　　引

　　

　　言

　　光合有效辐射是太阳辐射中能被绿色植物用来进行光合作用的能量。 光合有效辐射作为生物量形成的基本能源，直接影响着植物的生长、发育、产量和质量，光合有效辐射遥感产品能够提供全球及区域尺度的光合有效辐射，在光合潜力、作物生长模拟、碳循环、全球变化等研究中具有重要意义。

　　本文件旨在统一不同类型、不同时空分辨率的光合有效辐射遥感产品的真实性检验参考对象、检验流程、检验方法，并形成光合有效辐射遥感产品质量评价的指标与方法，保证光合有效辐射遥感产品在科学研究和行业应用中的一致性和产品质量。

　　GB/T 41281—2022

　　光合有效辐射遥感产品真实性检验

　　1 范围

　　本文件规定了光合有效辐射遥感产品真实性检验的参考对象、基本要求、检验方法和检验报告。

　　本文件适用于光合有效辐射遥感产品的真实性检验。

　　2 规范性引用文件

　　下列文件中的内容通过文中的规范性引用而构成本文件必不可少的条款。 其中，注 日期的引用文件，仅该日期对应的版本适用于本文件;不注 日期的引用文件，其最新版本(包括所有的修改单)适用于本文件。

　　GB/T

　　36296—2018

　　遥感产品真实性检验导则

　　GB/T

　　39468—2020

　　陆地定量遥感产品真实性检验通用方法

　　3 术语和定义

　　GB/T 36296—2018 界定的以及下列术语和定义适用于本文件。

　　3.1

　　太阳辐射 solarradiation

　　太阳发出的电磁辐射。

　　[来源：GB/T 36299—2018 , 2 . 3]

　　3.2

　　光合有效辐射 photosyntheticactiveradiation

　　太阳辐射光谱中可被绿色植物的质体色素吸收、转化并用于合成有机物质的一定波段的辐射能。注：一般把 400 nm~700 nm 的太阳辐射称为光合有效辐射。

　　[来源：GB/T 31163—2014 , 5 . 8]

　　3.3

　　短波辐射 shortwaveradiation

　　波长介于 0.28 μm~3 μm 的电磁辐射。

　　[来源：GB/T 31163—2014 , 5 . 9]

　　3.4

　　直接辐射 directradiation

　　从日面及其周围一小立体角内发出的辐射。

　　注：一般情况下，直接辐射是由视场角约为 5°的仪器测定的。因此，它包括 日面周围的部分散射辐射，即环 日辐射，而 日面本身的视场角仅约为 0.5 °。

　　[来源：GB/T 31163—2014 , 5 . 11]

　　3.5

　　散射辐射 diffuseradiation

　　太阳辐射被空气分子、云和空气中的各种微粒分散成无方向性的、但不改变其单色组成的辐射。

　　1

　　GB/T 41281—2022

　　[来源：GB/T 31163—2014 , 5 . 14]

　　3.6

　　总辐射 globalradiation

　　水平面从上方 2π立体角范围内接收到的直接辐射和散射辐射之和。

　　[来源：GB/T 31163—2014 , 5 . 15]

　　3.7

　　光量子通量密度 photonfluxdensity

　　单位时间投射到单位面积上的光子数量。

　　注：单位为微摩尔每平方米秒[μmol/(m2 · s)]

　　3.8

　　瞬时光合有效辐射 instantaneousphotosyntheticactiveradiation

　　某一时刻的光合有效辐射值。

　　注：单位为瓦每平方米(W/m2 )。

　　3.9

　　日累积光合有效辐射 dailyphotosyntheticactiveradiation

　　一天时间(24 h)内光合有效辐射的累加值。

　　注：单位为兆焦耳每平方米(MJ/m2 )。

　　3 . 10

　　月累积光合有效辐射 monthlyphotosyntheticactiveradiation

　　一个月时间内光合有效辐射的累加值。

　　注：单位为兆焦耳每平方米(MJ/m2 )。

　　3 . 1 1

　　年累积光合有效辐射 yearlyphotosyntheticactiveradiation

　　一年时间内光合有效辐射的累加值。

　　注：单位为兆焦耳每平方米(MJ/m2 )。

　　3 . 12

　　下行短波辐射 downwardshortwaveradiation

　　大气层顶的太阳短波辐射经吸收、散射等作用后到达地表的电磁辐射。

　　4 参考对象

　　光合有效辐射遥感产品真实性检验的参考对象包括。

　　a) 地面观测数据：包括已有观测站点或布设观测样点利用光量子测量表或总辐射表获取的地面观测数据。 全球已公开发布的光合有效辐射及总辐射地面观测数据列表见附录 A。

　　b) 已检遥感产品：已检验的光合有效辐射遥感产品或下行短波辐射遥感产品。

　　5 基本要求

　　光合有效辐射遥感产品真实性检验在符合 GB/T 36296—2018 第 7 章要求的同时，还符合以下要求：

　　a) 参考对象为已检遥感产品时，已检遥感产品的空间分辨率，宜与待检光合有效辐射遥感产品的

　　空间分辨率相同或相近;已检遥感产品的时间分辨率，应高于待检光合有效辐射遥感产品的时间分辨率或与待检光合有效辐射遥感产品的时间分辨率一致;

　　2

　　GB/T 41281—2022

　　b) 参考对象宜具有区域代表性，覆盖待检光合有效辐射遥感产品范围内的不同高程和地理条件;

　　c) 参考对象宜具有时间代表性，覆盖一天中不同时刻、一年中不同季节;

　　d) 参考对象宜具有天气状况代表性，包括晴天、阴天，为全天候天气条件。

　　6 检验方法

　　6 . 1 检验方法选择

　　光合有效辐射遥感产品的真实性检验方法包括直接检验方法和间接检验方法。 应首选地面观测数据进行光合有效辐射遥感产品的真实性检验，当有地面观测数据时，应采用直接检验法;当有已检遥感产品时，可采用间接检验法。

　　6 . 2 直接检验法

　　光合有效辐射遥感产品直接检验法检验流程应符合 GB/T 36296—2018 中 8 . 1 的规定，主要操作流程见图 1,具体步骤如下。

　　图 1 光合有效辐射遥感产品直接检验法操作流程

　　3

　　GB/T 41281—2022

　　a) 待检验光合有效辐射遥感产品：获取产品的类别、空间分辨率、时间分辨率、空间范围等信息。

　　b) 已有观测站点数据判断：根据待检光合有效辐射遥感产品的空间范围判断在此范围内是否有已有观测站点数据，已有观测站点数据见附录 A。 若有已有观测站点数据，则进行已有观测站点空间代表性评价;若没有已有观测站点数据，则进行观测样点布设。

　　c) 已有观测站点空间代表性评价：见附录 B 中的空间代表性评价方法计算已有地面观测站点的空间代表性评价指标。

　　d) 空间代表性判断：根据计算的空间代表性评价指标，若空间代表性评价指标大于或等于 0.8,则选择已有观测站点数据作为参考对象;若空间代表性评价指标小于 0 . 8,则进行观测样点布设。

　　e) 已有观测样点数据获取：根据待检光合有效辐射遥感产品的类别、空间分辨率、时间分辨率以及空间范围信息下载对应的已有观测站点数据。

　　f) 观测样点布设：按附录 C 中 C.1 的要求进行观测样点布设。

　　g) 观测样点数据采集：按 C.2 的要求进行观测样点数据采集。

　　h) 光合有效辐射计算：

　　— 若获取的地面观测数据为光量子测量表测量的观测数据，按公式(C. 1)将地面观测的光量子通量密度数据转换为光合有效辐射数据;

　　— 若获取的地面观测数据为总辐射表测量的观测数据，按公式(C. 2)将地面观测的总辐射数据转换为光合有效辐射数据。

　　i) 空间匹配：将地面观测站点/样点与待检验的光合有效辐射遥感产品像元进行空间匹配。

　　j) 时间一致性判断：根据待检验的光合有效辐射遥感产品的时间尺度，判断地面观测数据的时间尺度是否一致，若时间尺度一致，则作为像元尺度相对真值;若时间尺度不一致，则应进行时间一致性转换。

　　k ) 时间一致性转换：

　　— 若待检验的光合有效辐射遥感产品为瞬时光合有效辐射产品，则宜选择待检验遥感产品时刻 30 min 内的地面观测数据平均值作为像元尺度相对真值;

　　— 若待检验的光合有效辐射遥感产品为 日 累积光合有效辐射产品，应对地面观测数据进行日 累积计算作为像元尺度相对真值(见附录 D 中 D. 1) ;

　　— 若待检验的光合有效辐射遥感产品为月累积光合有效辐射遥感产品，应对地面观测数据进行月累积计算作为像元尺度相对真值(见 D. 1) ;

　　— 若待检验的光合有效辐射遥感产品为年累积光合有效辐射遥感产品，应对地面观测数据进行年累积计算作为像元尺度相对真值(见 D. 1) 。

　　l) 准确度评价：按 GB/T 36296—2018 中 6.1 的准确度评价指标定量表达检验的光合有效辐射遥感产品的准确度，准确度评价指标至少应包括均方根误差、相关系数。

　　m) 不确定度分析：按 GB/T 36296—2018 中 6.2 的不确定度评价指标进行光合有效辐射遥感产品的不确定度分析，不确定度评价指标至少应包括标准差。

　　6 . 3 间接检验法

　　间接检验法检验流程应符合 GB/T 36296—2018 中 8 . 2 . 1 的规定，主要操作流程见图 2,具体步骤如下：

　　4

　　GB/T 41281—2022

　　图 2 光合有效辐射遥感产品间接检验操作流程

　　a) 待检验光合有效辐射遥感产品：获取产品的类别、空间分辨率、时间分辨率、空间范围等信息;

　　b) 已检遥感产品选择：若已检遥感产品为已检光合有效辐射遥感产品，则进行空间一致性判断;若为已检下行短波辐射遥感产品，按照 C. 3 . 2 将下行短波辐射遥感产品转换为光合有效辐射遥感产品;

　　c) 空间一致性判断：根据待检光合有效辐射遥感产品的空间分辨率，判断已检遥感产品是否空间尺度一致，若一致，则进行时间一致性判断;若空间尺度不一致，应对已检遥感产品进行空间尺度转换，获得空间尺度一致的已检遥感产品，再进行时间一致性判断;

　　d) 空间尺度转换：当待检光合有效辐射遥感产品与已检遥感产品的空间分辨率不一致时，应使用空间重采样方法对已检遥感产品进行空间尺度转换至与待检光合有效辐射遥感产品相一致的空间尺度;

　　e) 时间一致性判断：根据待检验的光合有效辐射遥感产品的时间尺度，判断已检光合有效辐射遥感产品是否时间尺度一致，若时间尺度一致，则采用已检光合有效辐射遥感产品对应的像元值作为像元尺度相对真值;若时间尺度不一致，应对已检遥感产品进行时间一致性转换;

　　f) 时间一致性转换：时间一致性转换方法见 D.2,获得与待检光合有效辐射遥感产品时间一致的已检遥感产品，作为像元尺度真值;

　　5

　　GB/T 41281—2022

　　g) 样本选择：按 GB/T 39468—2020 中 5.2 的检验样本抽样方法进行样本选择;

　　h) 准确度评价：按 GB/T 36296—2018 中 6.1 的准确度评价指标定量表达检验的光合有效辐射遥感产品的准确度，准确度评价指标至少应包括均方根误差、相关系数;

　　i) 不确定度分析：按 GB/T 36296—2018 中 6.2 的不确定度评价指标进行光合有效辐射遥感产品的不确定度分析，不确定度评价指标至少应包括标准差。

　　7 检验报告

　　7 . 1 封面信息

　　检验报告封面至少应包括以下信息：

　　a) 检验报告编号;

　　b ) 检验报告名称;

　　c) 检验负责人;

　　d) 检验核对人;

　　e) 检验签发人;

　　f) 检验单位;

　　g) 送检单位;

　　h ) 检验时间。

　　7 . 2 正文信息

　　7 . 2 . 1 待检验的光合有效辐射产品概述

　　对待检验的光合有效辐射遥感产品进行简要描述，至少应包括：

　　a) 产品的名称、类型、数据源、地理参考、投影方式;

　　b) 产品的空间分辨率、空间覆盖范围、获取时间、时间分辨率;

　　c) 产品的光谱范围、单位;

　　d) 产品所采用的主要算法以及算法的特点。

　　7 . 2 . 2 参考对象描述

　　对作为参考对象进行描述，至少应包括：

　　a) 参考对象的常规信息，包括类型、空间特征、时间特征、光谱范围;

　　b) 参考对象的测量方法、测量仪器、数据格式，单位;

　　c) 参考对象的精度描述。

　　7 . 2 . 3 检验方法及流程

　　对采用的检验方法和检验过程进行描述，至少应包括：

　　a) 检验方法概述;

　　b ) 检验流程;

　　c) 检验结果的评价指标;

　　d) 检验过程的记录要求;

　　e) 检验结果的存档。

　　7 . 2 . 4 真实性检验结论描述

　　对真实性检验的评价结果进行描述，至少应包括：

　　6

　　GB/T 41281—2022

　　a) 真实性检验结果总体评价：描述所检验光合有效辐射遥感产品的总体准确度;

　　b) 分项指标评价：描述所检验光合有效辐射遥感产品在研究区内不同区域的准确度;描述光合有效辐射遥感产品在不同时间的准确度;描述光合有效辐射遥感产品在不同天气状况的准确度;

　　c) 对产品的评价：通过对各项指标的评价和分析，给出所检验光合有效辐射遥感产品的评价结果，包括产品准确度、不确定度。

　　7 . 3 检验报告信息简表

　　光合有效辐射遥感产品的真实性检验报告信息简表编制见 GB/T 36296—2018 附录 D,样例见附录 E。

　　7

　　GB/T 41281—2022

　　附 录 A

　　(资料性)

　　全球光合有效辐射及总辐射地面观测数据列表

　　表 A. 1 给出了已公开发布的全球光合有效辐射及总辐射地面观测数据。

　　表 A.1 全球光合有效辐射及总辐射地面观测数据列表

　　数据来源

　　观测量

　　时间分辨率

　　年限范围

　　中国气象数据网

　　总辐射

　　日

　　1957 —

　　CERN

　　光合有效辐射、总辐射

　　日

　　1998 —

　　AmeriFlux

　　光合有效辐射、总辐射

　　30 min/ 1 h

　　1996 —

　　AsiaFlux

　　光合有效辐射、总辐射

　　30 min/ 1 h

　　1999—2018

　　ChinaFLUX

　　光合有效辐射、总辐射

　　30 min/ 日

　　2002 —

　　GAME

　　光合有效辐射、总辐射

　　30 min/ 1 h

　　1998—2002

　　黑河试验-水文气象和通量观测网

　　光合有效辐射、总辐射

　　10 min/ 30 min

　　2007—2017

　　SURFRAD

　　光合有效辐射、总辐射

　　1 min

　　1995 —

　　8

　　GB/T 41281—2022

　　附 录 B

　　(资料性)

　　空间代表性评价方法

　　以辐射相关强度(Correlation Strength, Cs)作为光合有效辐射地面观测数据空间代表性的评价指

　　标，由式(B. 1) 表示：

　　n

　　Csif n= 0 , Cs = 1 ……

　　式中：

　　n — 距站点位置的像元数;

　　dj — 方向[下标 j∈ (1,2,3...16) ,dj ∈ (0。，22.5。，45。 337.5。)];

　　Cs(n，dj)— 在 dj 方向上，距离站点 n个像元位置处的辐射相关强度;

　　rad(0,dj)— 站点所在像元的辐射值;

　　rad(i，dj)— 在 dj 方向上距站点i个像元处的辐射值。

　　辐射相关强度也可写成公式(B. 2) :

　　Cs(h，dj)=Cs(n，dj)，h=n ·R …………………………( B.2 )

　　式中：

　　Cs(h，dj)— 表示在某个特定方向上，将空间中某一点的观测值在该方向上 h(或 n · R)范围内的空间代表性定义为站点观测的总辐射值与 h范围内各距离处的总辐射的平均累

　　积的相对误差;

　　R — 像元的空间分辨率。

　　为直观地描述辐射相关强度随距离的变化情况，使用最优辐射相关强度的概念，由式(B. 3) 表示：

　　Csmin(n，dj)=min[Cs(0,dj)，Cs(1,dj)...Cs(n，dj)]………………( B.3 )

　　式中：

　　Csmin(n，dj)— 为最优辐射相关强度。

　　以 Csmin(x，dj)=y为例，它表示在方向 dj 上，在 x范围内辐射相关性强度的最小值为y，观测站

　　点/样点的空间代表性即指 Csmin 。

　　9

　　GB/T 41281—2022

　　附 录 C

　　(规范性)

　　观测样点布设及光合有效辐射地面观测

　　C.1 观测样点布设

　　安装场地应当是地平面以上无任何明显障碍，特别是在一年中 日出和 日落的方位里，不应有 5°以

　　上的障碍物。 在观测仪器以上，应避免附近存在浅颜色的屏障。

　　C.2 地面观测

　　C.2 . 1 观测仪器

　　光合有效辐射的量度有两种计量系统，一种是能量系统，用总辐射表观测得到，其单位是瓦每平方米

　　(W/m2 );另一种是量子系统，用光量子测量表观测得到，其单位是微摩尔每平方米秒[μmol/(m2 · s)]。

　　a) 总辐射表，观测使用的辐射传感器多为热电型，传感器由感应面与热电堆组成。感应面涂上吸收率高、光谱响应无选择的无光黑漆，紧贴在感应面下部是热电堆，它与感应面应保持绝缘。热电堆工作端(热端)位于感应面下面。参考端(冷端)位于隐蔽处。 当辐射源照射到辐射表时，感应面增热，使热电堆冷热接点间产生温度差，进而感应出电动势。

　　b) 光量子测量表，观测光量子通量密度的光电型辐射表工作原理是光电池，当一定波长的光照射半导体光电器件时，半导体产生电子移动形成电流，光照射强度与其输出的短路电流成正比。

　　C.2 .2 观测仪器安装

　　总辐射表或光量子测量表安装在专用的支架上，支架离地面或楼顶平台不低于 1.50 m,要求牢固，

　　保持水平状态。

　　C.2 . 3 数据采集

　　总辐射或光量子通量密度每分钟采样不低于 6 次(瞬时值)，计算有效样本数据的平均值，用 1 min

　　的平均值作为每分钟的总辐射或光量子通量密度平均值。

　　C.3 数据转换

　　C.3 . 1 光量子通量密度与光合有效辐射的转换

　　两者之间的量子转换系数由式(C. 1)表示：

　　犝PAR = μ × 犙PAR …………………………( C.1 )

　　式中：

　　犝PAR — 光合有效量子通量密度，单位为微摩尔每平方米秒[μmol/(m2 · s)] ;

　　μ — 通常取值为 4.55 ;

　　犙PAR — 光合有效辐射照度，单位为瓦每平方米(W/m2 )。

　　C.3 .2 总辐射与光合有效辐射的转换

　　光合有效辐射是太阳总辐射或短波辐射的一部分，由式(C. 2)表示：

　　犙PAR =ηQ × 犙 …………………………( C.2 )

　　10

　　GB/T 41281—2022

　　式中：

　　ηQ — 光合有效辐射所占总辐射的比例数即光合有效系数，ηQ 是一个在天文因子、气候因子的共

　　同作用下不断变化的量;

　　Q — 太阳总辐射，单位为瓦每平方米(W/m2 )。

　　不同区域的光合有效系数 ηQ 可由式(C.3)估算：

　　ηQ =a+b …………………………( C.3 )

　　式中：

　　a，b — 待定系数;

　　Q’— 天文辐射。

　　11

　　GB/T 41281—2022

　　附 录 D

　　(资料性)

　　时间一致性转换方法

　　D.1 地面观测数据时间一致性转换方法

　　地面观测的光合有效辐射数据根据发布数据集的不同，具有不同的时间分辨率。 需根据待检遥感产品的时间分辨率进行时间一致性转换。

　　由瞬时地面观测值计算某一时间段的累积地面观测值，由式(D. 1)表示：

　　E Iidt …………………………( D.1 )

　　式中：

　　E — 从 t1 到 t2 时间段内的累积光合有效辐射，单位为焦耳每平方米(J/m2 ) ;