（高清版）GBT 41281-2022 光合有效辐射遥感产品真实性检验

光合有效辐射遥感产品真实性检验2022-03-09发布国家标准化管理委员会I前言 12规范性引用文件 13术语和定义 14参考对象 25基本要求 26检验方法 37检验报告 6附录A(资料性)全球光合有效辐射及总辐射地面观测数据列表 8附录B(资料性)空间代表性评价方法 9附录C(规范性)观测样点布设及光合有效辐射地面观测 附录D(资料性)时间一致性转换方法 附录E(资料性)光合有效辐射遥感产品真实性检验信息简表样例 参考文献 ⅢGB/T41281—2022本文件按照GB/T1.1—2020《标准化工作导则第1部分：标准化文件的结构和起草规则》的规定起草。请注意本文件的某些内容可能涉及专利。本文件的发布机构不承担识别专利的责任。本文件由中国科学院提出。本文件由全国遥感技术标准化技术委员会(SAC/TC327)归口。本文件起草单位：中国科学院空天信息创新研究院、中国科学院地理科学与资源研究所、中国气象科学研究院、中国资源卫星应用中心、生态环境部卫星环境应用中心。本文件主要起草人：李丽、辛晓洲、张海龙、闻建光、吴骅、王培娟、高彦华、姚延娟、仲波、余珊珊、光合有效辐射是太阳辐射中能被绿色植物用来进行光合作用的能量。光合有效辐射作为生物量形成的基本能源，直接影响着植物的生长、发育、产量和质量，光合有效辐射遥感产品能够提供全球及区域尺度的光合有效辐射，在光合潜力、作物生长模拟、碳循环、全球变化等研究中具有重要意义。本文件旨在统一不同类型、不同时空分辨率的光合有效辐射遥感产品的真实性检验参考对象、检验流程、检验方法，并形成光合有效辐射遥感产品质量评价的指标与方法，保证光合有效辐射遥感产品在科学研究和行业应用中的一致性和产品质量。1光合有效辐射遥感产品真实性检验本文件规定了光合有效辐射遥感产品真实性检验的参考对象、基本要求、检验方法和检验报告。本文件适用于光合有效辐射遥感产品的真实性检验。2规范性引用文件下列文件中的内容通过文中的规范性引用而构成本文件必不可少的条款。其中，注日期的引用文件，仅该日期对应的版本适用于本文件；不注日期的引用文件，其最新版本(包括所有的修改单)适用于本文件。GB/T36296—2018遥感产品真实性检验导则GB/T39468—2020陆地定量遥感产品真实性检验通用方法3术语和定义GB/T36296—2018界定的以及下列术语和定义适用于本文件。太阳发出的电磁辐射。太阳辐射光谱中可被绿色植物的质体色素吸收、转化并用于合成有机物质的一定波段的辐射能。注：一般把400nm～700nm的太阳辐射称为光合有效辐射。短波辐射shortwaveradiation从日面及其周围一小立体角内发出的辐射。注：一般情况下，直接辐射是由视场角约为5°的仪器测定的。因此，它包括日面周围的部分散射辐射，即环日辐射，而日面本身的视场角仅约为0.5°。散射辐射diffuseradiation太阳辐射被空气分子、云和空气中的各种微粒分散成无方向性的、但不改变其单色组成的辐射。2水平面从上方2π立体角范围内接收到的直接辐射和散射辐射之和。光量子通量密度photonfluxdensity单位时间投射到单位面积上的光子数量。注：单位为微摩尔每平方米秒[μmol/(m²·s)]某一时刻的光合有效辐射值。注：单位为瓦每平方米(W/m²)。一天时间(24h)内光合有效辐射的累加值。注：单位为兆焦耳每平方米(MJ/m²)。一个月时间内光合有效辐射的累加值。注：单位为兆焦耳每平方米(MJ/m²)。一年时间内光合有效辐射的累加值。注：单位为兆焦耳每平方米(MJ/m²)。下行短波辐射downwardshortwaveradiation大气层顶的太阳短波辐射经吸收、散射等作用后到达地表的电磁辐射。4参考对象光合有效辐射遥感产品真实性检验的参考对象包括。a)地面观测数据：包括已有观测站点或布设观测样点利用光量子测量表或总辐射表获取的地面观测数据。全球已公开发布的光合有效辐射及总辐射地面观测数据列表见附录A。b)已检遥感产品：已检验的光合有效辐射遥感产品或下行短波辐射遥感产品。5基本要求光合有效辐射遥感产品真实性检验在符合GB/T36296—2018第7章要求的同时，还符合以下a)参考对象为已检遥感产品时，已检遥感产品的空间分辨率，宜与待检光合有效辐射遥感产品的空间分辨率相同或相近；已检遥感产品的时间分辨率，应高于待检光合有效辐射遥感产品的时间分辨率或与待检光合有效辐射遥感产品的时间分辨率一致；3b)参考对象宜具有区域代表性，覆盖待检光合有效辐射遥感产品范围内的不同高程和地理条件；c)参考对象宜具有时间代表性，覆盖一天中不同时刻、一年中不同季节；d)参考对象宜具有天气状况代表性，包括晴天、阴天，为全天候天气条件。6检验方法6.1检验方法选择光合有效辐射遥感产品的真实性检验方法包括直接检验方法和间接检验方法。应首选地面观测数据进行光合有效辐射遥感产品的真实性检验，当有地面观测数据时，应采用直接检验法；当有已检遥感产品时，可采用间接检验法。6.2直接检验法光合有效辐射遥感产品直接检验法检验流程应符合GB/T36296—2018中8.1的规定，主要操作流程见图1,具体步骤如下。待检验的光待检验的光合有效辐射遥感产品否已有观测站点数据断是已有观测站点空间代表性评价否空间代表性判断是已有观测站点观测样点数据数据获取采集·光量子通量密度与光合有效辐射转换·总辐射与光合有效辐射转换空间匹否时间·致性判断瞬时光合有效辐射遥感产品是时间一致性转换年累积光合有效辐射遥感产品准确度评价不确定度分析光合有效辐射计算观测样点布设图1光合有效辐射遥感严品直接检验法操作流程4a)待检验光合有效辐射遥感产品：获取产品的类别、空间分辨率、时间分辨率、空间范围等信息。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)时间一致性转换：——若待检验的光合有效辐射遥感产品为瞬时光合有效辐射产品，则宜选择待检验遥感产品时刻30min内的地面观测数据平均值作为像元尺度相对真值； 若待检验的光合有效辐射遥感产品为日累积光合有效辐射产品，应对地面观测数据进行日累积计算作为像元尺度相对真值(见附录D中D.1);——若待检验的光合有效辐射遥感产品为月累积光合有效辐射遥感产品，应对地面观测数据进行月累积计算作为像元尺度相对真值(见D.1);——若待检验的光合有效辐射遥感产品为年累积光合有效辐射遥感产品，应对地面观测数据进行年累积计算作为像元尺度相对真值(见D.1)。1)准确度评价：按GB/T36296—2018中6.1的准确度评价指标定量表达检验的光合有效辐射遥感产品的准确度，准确度评价指标至少应包括均方根误差、相关系数。m)不确定度分析：按GB/T36296—2018中6.2的不确定度评价指标进行光合有效辐射遥感产品的不确定度分析，不确定度评价指标至少应包括标准差。6.3间接检验法间接检验法检验流程应符合GB/T36296—2018中8.2.1的规定，主要操作流程见图2,具体步骤5待检验的光合有效辐射遥感产品已检光合有效辐射遥感产品选择否空间一致性判断是空间尺度转换否时问一致性判断是时间一致性转换样木选择准确度评价不确定度分析图2光合有效辐射遥感产品间接检验操作流程a)待检验光合有效辐射遥感产品：获取产品的类别、空间分辨率、时间分辨率、空间范围等信息；b)已检遥感产品选择：若已检遥感产品为已检光合有效辐射遥感产品，则进行空间一致性判断；若为已检下行短波辐射遥感产品，按照C.3.2将下行短波辐射遥感产品转换为光合有效辐射遥感产品；c)空间一致性判断：根据待检光合有效辐射遥感产品的空间分辨率，判断已检遥感产品是否空间尺度一致，若一致，则进行时间一致性判断；若空间尺度不一致，应对已检遥感产品进行空间尺度转换，获得空间尺度一致的已检遥感产品，再进行时间一致性判断；d)空间尺度转换：当待检光合有效辐射遥感产品与已检遥感产品的空间分辨率不一致时，应使用空间重采样方法对已检遥感产品进行空间尺度转换至与待检光合有效辐射遥感产品相一致的空间尺度；e)时间一致性判断：根据待检验的光合有效辐射遥感产品的时间尺度，判断已检光合有效辐射遥感产品是否时间尺度一致，若时间尺度一致，则采用已检光合有效辐射遥感产品对应的像元值作为像元尺度相对真值；若时间尺度不一致，应对已检遥感产品进行时间一致性转换；f)时间一致性转换：时间一致性转换方法见D.2,获得与待检光合有效辐射遥感产品时间一致的已检遥感产品，作为像元尺度真值；6g)样本选择：按GB/T39468—2020中5.2的检验样本抽样方法进行样本选择；h)准确度评价：按GB/T36296—2018中6.1的准确度评价指标定量表达检验的光合有效辐射遥感产品的准确度，准确度评价指标至少应包括均方根误差、相关系数；i)不确定度分析：按GB/T36296—2018中6.2的不确定度评价指标进行光合有效辐射遥感产品的不确定度分析，不确定度评价指标至少应包括标准差。7检验报告7.1封面信息检验报告封面至少应包括以下信息·a)检验报告编号；b)检验报告名称；c)检验负责人；d)检验核对人；e)检验签发人；f)检验单位；g)送检单位；h)检验时间。7.2正文信息7.2.1待检验的光合有效辐射产品概述对待检验的光合有效辐射遥感产品进行简要描述，至少应包括：b)产品的空间分辨率、空间覆盖范围、获取时间、时间分辨率；c)产品的光谱范围、单位；d)产品所采用的主要算法以及算法的特点。7.2.2参考对象描述对作为参考对象进行描述，至少应包括：b)参考对象的测量方法、测量仪器、数据格式，单位；c)参考对象的精度描述。7.2.3检验方法及流程对采用的检验方法和检验过程进行描述，至少应包括：a)检验方法概述；b)检验流程；c)检验结果的评价指标；d)检验过程的记录要求；e)检验结果的存档。7.2.4真实性检验结论描述对真实性检验的评价结果进行描述，至少应包括：7a)真实性检验结果总体评价：描述所检验光合有效辐射遥感产品的总体准确度；b)分项指标评价：描述所检验光合有效辐射遥感产品在研究区内不同区域的准确度；描述光合有效辐射遥感产品在不同时间的准确度；描述光合有效辐射遥感产品在不同天气状况的准确度；c)对产品的评价：通过对各项指标的评价和分析，给出所检验光合有效辐射遥感产品的评价结7.3检验报告信息简表光合有效辐射遥感产品的真实性检验报告信息简表编制见GB/T36296—2018附录D,样例见附8(资料性)全球光合有效辐射及总辐射地面观测数据列表表A.1给出了已公开发布的全球光合有效辐射及总辐射地面观测数据。表A.1全球光合有效辐射及总辐射地面观测数据列表数据来源观测量时间分辨率年限范围中国气象数据网总辐射日光合有效辐射、总辐射日AmeriFlux光合有效辐射、总辐射AsiaFlux光合有效辐射、总辐射ChinaFLUX光合有效辐射、总辐射30min/日光合有效辐射、总辐射黑河试验-水文气象和通量观测网光合有效辐射、总辐射2007—2017光合有效辐射、总辐射9以辐射相关强度(Correlation(资料性)空间代表性评价方法Strength,CS)作为光合有效辐射地面观测数据空间代表性的评价指n——距站点位置的像元数；d;——方向[下标j∈(1,2,3…16),d,∈(0°,22.5°,45°……337.5°)];CS(n,d;)—-在d,方向上，距离站点n个像元位置处的辐射相关强度；rad(0,d;)——站点所在像元的辐射值；rad(i,d;)——在d;方向上距站点i个像元处的辐射值。辐射相关强度也可写成公式(B.2):CS(h,d;)=CS(n,d;),h=n·RCS(h,d;)——表示在某个特定方向上，将空间中某一点的观测值在该方向上h(或n·R)范围内的空间代表性定义为站点观测的总辐射值与h范围内各距离处的总辐射的平均累积的相对误差；R——像元的空间分辨率。为直观地描述辐射相关强度随距离的变化情况，使用最优辐射相关强度的概念，由式(B.3)表示：CSmim(n,d;)=min[CS(0,d;),CS(1,d,)…CS(n,d;)]………………(B.3)CSmin(n,d;)——为最优辐射相关强度。以CSmn(x,d)=y为例，它表示在方向d,上，在x范围内辐射相关性强度的最小值为y,观测站点/样点的空间代表性即指CSmin。(规范性)观测样点布设及光合有效辐射地面观测安装场地应当是地平面以上无任何明显障碍，特别是在一年中日出和日落的方位里，不应有5°以上的障碍物。在观测仪器以上，应避免附近存在浅颜色的屏障。C.2地面观测C.2.1观测仪器光合有效辐射的量度有两种计量系统，一种是能量系统，a)总辐射表，观测使用的辐射传感器多为热电型，传感器由感应面与热电堆组成。感应面涂上吸收率高、光谱响应无选择的无光黑漆，紧贴在感应面下部是热电堆，它与感应面应保持绝缘。热电堆工作端(热端)位于感应面下面。参考端(冷端)位于隐蔽处。当辐射源照射到辐射表时，感应面增热，使热电堆冷热接点间产生温度差，进而感应出电动势。b)光量子测量表，观测光量子通量密度的光电型辐射表工作原理是光电池，当一定波长的光照射半导体光电器件时，半导体产生电子移动形成电流，光照射强度与其输出的短路电流成正比。C.2.2观测仪器安装总辐射表或光量子测量表安装在专用的支架上，支架离地面或楼顶平台不低于1.50m,要求牢固，保持水平状态。C.2.3数据采集总辐射或光量子通量密度每分钟采样不低于6次(瞬时值),计算有效样本数据的平均值，用1min的平均值作为每分钟的总辐射或光量子通量密度平均值。C.3数据转换C.3.1光量子通量密度与光合有效辐射的转换两者之间的量子转换系数由式(C.1)表示：UpAR=μ×QPAR……………(C.1)UpAR——光合有效量子通量密度，单位为微摩尔每平方米秒[μmol/(m²·s)7;μ——通常取值为4.55;QpAR——光合有效辐射照度，单位为瓦每平方米(W/m²)。C.3.2总辐射与光合有效辐射的转换光合有效辐射是太阳总辐射或短波辐射的一部分，由式(C.2)表示：QpAR=ηg×Q7a——光合有效辐射所占总辐射的比例数即光合有效系数，ηα是一个在天同作用下不断变化的量；不同区域的光合有效系数ηa可由式(C.3)估算：Q'——天文辐射。GB/T41281—2022(资料性)时间一致性转换方法D.1地面观测数据时间一致性转换方法地面观测的光合有效辐射数据根据发布数据集的不同，具有不同的时间分辨率。需根据待检遥感产品的时间分辨率进行时间一致性转换。由瞬时地面观测值计算某一时间段的累积地面观测值，由式(D.1)表示：E——从t₁到tz时间段内的累积光合有效辐射，单位为焦耳每平方米(J/m²);t₁——瞬时地面观测值在某一时间段累积的起始时刻；t₂——瞬时地面观测值在某一时间段累积的结束时刻；I;——i时刻的瞬时光合有效辐射值，单位为瓦每平方米(W/m²)。D.2光合有效辐射遥感产品时间一致性转换方法D.2.1相邻时刻的瞬时值计算小时累积值由相邻时刻的瞬时值计算小时累积值，由式(D.2)表示：…………(D.2)W、——小时累积光合有效辐射值，单位为焦耳每平方米(J/m²);△T——积分时间，单位为秒(s);I₁——本小时的瞬时光合有效辐射值，单位为瓦每平方米(W/m²);I₂——下一小时的瞬时光合有效辐射值，单位为瓦每平方米(W/m²)。D.2.2由一个时刻的瞬时值计算日累积值直接辐射积分公式由式(D.3),式(D.4)表示：…(D.3)散射辐射积分公式由式(D.5),式(D.6)表示：…………(D.5)………(D.6)Edr——太阳天顶角为θ时的瞬时太阳直射辐射，单位为瓦每平方米(W/m²);Es.dir.n——当地时间正午时的瞬时太阳直射辐射，单位为瓦每平方米(W/m²);θn——当地时间正午时的太阳天顶角；D——日长，单位为秒(s);Eay.dir——日累积的太阳直射辐射，单位为焦耳每平方米(J/m²)Eai——太阳天顶角为θ时的瞬时太阳散射辐射，单位为瓦每平方米(W/m²);Eat.n——当地时间正午时的瞬时太阳散射辐射，单位为瓦每平方米(W/m²);Eay.di——日累积的太阳散射辐射，单位为焦耳每平方米(J/m²)。GB/T41281—2022(资料性)光合有效辐射遥感产品真实性检验信息简表样例表E.1给出了光合有效辐射遥感产品真实性检验报告信息简表的参考样例，内容以直接检验法为例，同时也填写了间接检验法特征描述供参考。表E.1光合有效辐射遥感产品真实性检验报告信息简表样例待检遥感产品特征描述产品名称MuSyQ光合有效辐射遥感产品产品类型●数值型○类别型数据来源(传感器)MODIS/MTSAT地理参考(椭球、投影方式)经纬度投影，WGS84空间分辨率0.05°获取时间(YYYY-MM-DD,HH-MM-SS)2009-08-01至2009-08-312010-01-01至2010-06-30空间范围经度(°’")70°0'O"E137°0'o"E纬度(°'”)57°3'o"N15°0'o"N时间分辨率●瞬时●日累积○月累积〇年累积○其他光谱范围400nm～700nm单位W/m²,MJ/m²观测几何详见光合有效遥感产品元数据文件产品生产算法MuSyQ光合有效辐射遥感反演算法：是一种极轨与静止卫星结合的全天候光合有效辐射反演算法(POGDA算法