《GB∕T 41279-2022 反照率遥感产品真实性检验.pdf》由会员分享，可在线阅读，更多相关《GB∕T 41279-2022 反照率遥感产品真实性检验.pdf（22页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、书 书 书犐 犆犛 犆犆犛犃 中 华 人 民 共 和 国 国 家 标 准犌犅犜 反照率遥感产品真实性检验犞犪 犾 犻 犱 犪 狋 犻 狅 狀狅 犳犪 犾 犫 犲 犱 狅狉 犲犿狅 狋 犲狊 犲 狀 狊 犻 狀 犵狆 狉 狅 犱 狌 犮 狋 狊 发布 实施国 家 市 场 监 督 管 理 总 局国 家 标 准 化 管 理 委 员 会发 布目次前言引言范围规范性引用文件术语和定义基本要求检验方法 检验方法选择与检验对象要求 直接检验法 间接检验法检验报告 封面信息 正文信息 检验报告信息简表附录（资料性）空间代表性评价方法附录（规范性）上行下行辐射测量附录（资料性）地面测量辐射数据处理 附录（规范性
2、）波段转换计算 附录（规范性）真实反照率计算 附录（资料性）天空散射光比例获取方法 附录（资料性）反照率遥感产品真实性检验报告信息简表样例 参考文献 犌犅犜 前言本文件按照 标准化工作导则第部分：标准化文件的结构和起草规则的规定起草。请注意本文件的某些内容可能涉及专利。本文件的发布机构不承担识别专利的责任。本文件由中国科学院提出。本文件由全国遥感技术标准化技术委员会（）归口。本文件起草单位：中国科学院空天信息创新研究院、北京师范大学、兰州大学、浙江师范大学、中国科学院青藏高原研究所、北京大学、西南大学。本文件主要起草人：游冬琴、闻建光、刘强、吴小丹、林兴稳、唐勇、肖青、柳钦火、李新、焦子锑、范
3、闻捷、马明国、刘照言、王新鸿。犌犅犜 引言地表反照率反映了地球表面对太阳辐射的反射能力。作为地球辐射收支的一个重要参数，地表反照率在地球能量平衡、中长期天气预测和全球变化研究中具有重要意义。目前，通过遥感模型和反演方法对遥感数据进行系列处理，国内外已发布多种全球或局地区域的反照率遥感产品。本文件旨在制定科学、规范的反照率遥感产品真实性检验方法和流程，有利于全面客观地评价其准确度和分析其不确定度，提升反照率遥感产品的质量，推进反照率遥感产品的定量应用。犌犅犜 反照率遥感产品真实性检验范围本文件规定了地表对太阳辐射反射波段的反照率遥感产品真实性检验的基本要求、检验方法和检验报告。本文件适用于陆地地
4、表和水体地表反照率遥感产品的真实性检验。规范性引用文件下列文件中的内容通过文中的规范性引用而构成本文件必不可少的条款。其中，注日期的引用文件，仅该日期对应的版本适用于本文件；不注日期的引用文件，其最新版本（包括所有的修改单）适用于本文件。遥感产品真实性检验导则 陆地定量遥感产品真实性检验通用方法术语和定义 界定的以及下列术语和定义适用于本文件。反照率犪 犾 犫 犲 犱 狅物体表面向半球（）空间反射的辐射通量与半球（）空间入射在物体表面上的辐射通量之比。来源：，黑天空反照率犫 犾 犪 犮 犽狊 犽 狔犪 犾 犫 犲 犱 狅仅考虑直射光照射的情况下，物体表面向半球（）空间反射的辐射通量与入射的直射
5、辐射通量之比。来源：，白天空反照率狑犺 犻 狋 犲狊 犽 狔犪 犾 犫 犲 犱 狅仅考虑散射光照射的情况下，物体表面向半球（）空间反射的辐射通量与入射的散射辐射通量之比。来源：，真实反照率犪 犮 狋 狌 犪 犾犪 犾 犫 犲 犱 狅入射辐射同时包括直射光照射和散射光照射的自然界真实情况下，物体表面向半球（）空间反射的辐射通量与入射的辐射通量之比。天空散射光比例犱 犻 犳 犳 狌 狊 犲狊 犽 狔 犾 犻 犵 犺 狋犳 狉 犪 犮 狋 犻 狅 狀；犱 犻 犳 犳 狌 狊 犲狊 犽 狔 犾 犻 犵 犺 狋狉 犪 狋 犻 狅入射的天空散射辐射通量与入射辐射通量之比。犌犅犜 窄波段反照率狀 犪 狉 狉
6、 狅狑犫 犪 狀 犱犪 犾 犫 犲 犱 狅对应于传感器特定单一波段的反照率。宽波段反照率犫 狉 狅 犪 犱 犫 犪 狀 犱犪 犾 犫 犲 犱 狅覆盖短波波段、可见光波段、近红外波段范围的无传感器波段响应特征的反照率。注：短波波段（），下限可为 ，上限可为 ；可见光波段（）；近红外波段（）。通常可由窄波段反照率经窄宽波段转换得到。基本要求反照率遥感产品真实性检验应符合 中第章的规定，同时：）参考对象与待检反照率遥感产品时间应一致，入射辐射条件应相同，范围覆盖同一区域；）参考对象应首选地面测量的反照率，在无地面测量反照率时，应选用已检反照率遥感产品作为参考对象；）参考对象宜包含待检反照率遥感产品不
7、同值域典型地表反照率，即低反照率（如浓密植被覆盖区域）、中反照率（如裸土区域）、高反照率（如冰雪区域）。检验方法 检验方法选择与检验对象要求 检验方法选择根据待检反照率遥感产品的像元尺度和样本的空间代表性，选择相应检验方法：）当地面测量的样本空间代表性可代表待检反照率遥感产品像元尺度，应采用直接检验法。否则，宜采用间接检验法中的多尺度逐级检验法。）当有已检反照率遥感产品时，可采用间接检验法中的交叉检验法。检验对象要求在实施反照率遥感产品的检验过程中，依据检验方法的不同对检验对象要求如下。）直接检验法中，若待检反照率遥感产品中无真实反照率，应由待检反照率遥感产品中的黑天空反照率和白天空反照率计算
8、真实反照率进行检验。真实反照率的准确度和不确定度代表反照率遥感产品的准确度和不确定度。）间接检验方法中，检验过程中宜尽量对黑天空反照率和白天空反照率都进行检验，对反照率遥感产品进行全面的评价；待检反照率遥感产品的黑天空反照率与参考对象的黑天空反照率对应检验，待检反照率遥感产品的白天空反照率与参考对象的白天空反照率对应检验。直接检验法 参考对象的地面抽样方法参考对象地面抽样方法如下：）样区地表异质性按照 中 的规定进行样区异质性评价；）样区空间范围宜不小于倍待检反照率遥感产品像元区域；）对于均质样区，应在样区内进行随机抽样或系统抽样，对于非均质样区，宜根据地表覆盖在样犌犅犜 区内进行分层抽样，抽
9、样方法见 的附录。检验流程检验流程应符合 中 的规定，主要操作流程见图。图反照率遥感产品直接检验流程操作流程如下。）样本选择：）应首选已有辐射测量站点，进行空间代表性评价（评价方法见附录），若空间代表性可代表待检反照率遥感产品像元尺度，可用于下一步检验；）若无辐射测量站点，应按 的规定进行地面抽样方法设计，确定样本。）辐射通量数据获取：若为已有辐射测量站点，直接获取其已有的上行下行辐射测量数据；否则，按附录测量上行下行辐射通量。）辐射通量数据选择：通过预期的下行辐射阈值或者云信息，选择有效测量数据，用于下一步检验，选择方法见附录中。）时间一致性判断：根据待检反照率遥感产品的时间代表范围，判定上
10、行下行辐射测量值的时间代表范围是否一致。若一致，用上行辐射除以下行辐射计算得到地面测量反照率（见），用于下一步检验；若不一致，应首先进行时间一致性处理。）时间一致性处理：见。）波段一致性判断：判断地面测量反照率数据与待检反照率遥感产品波段是否一致。若波段一犌犅犜 :己有辐射测量站点空评价 实验场地面抽样方法设计le 己有辐射测量站点辐射数据辐射通量数据获取L 斗le 地面抽样点上行/下行辐射测量致则作为像元反照率相对真值进行下一步检验；若不一致，应进行波段转换。）波段转换：按附录将地面测量的光谱反照率转换为与待检反照率遥感产品波段一致的反照率。）时空匹配点的真实反照率计算：按附录，将与样本时空
11、匹配点的选定波段的待检反照率遥感产品黑天空反照率和白天空反照率通过天空散射光比例计算得到真实反照率，天空散射光比例获取方法见附录。）准确度评价：按 中 的准确度评价指标定量计算待检反照率遥感产品的准确度，准确度评价指标应包括但不限于均方根误差、相关系数。）不确定度分析：按 中 的不确定度评价指标定量计算待检反照率遥感产品的不确定度，不确定度指标应包括但不限于标准差。间接检验法 基于地面测量和高分辨率反照率影像的多尺度逐级检验法检验流程应符合 中 和 的规定，主要操作流程见图。图反照率遥感产品多尺度逐级检验方法流程操作流程如下。）高分辨率反照率影像准备：空间尺度应与地面测量的反照率相同或相近；时
12、间应与待检反照率遥感产品一致。犌犅犜 高分辨率反照率影像准各高分辨率反照率影像直接检验空间配准二E空间尺度转换否否波段转换）高分辨率反照率影像直接检验：按 的检验流程对高分辨率反照率影像进行直接检验，计算高分辨率反照率影像的准确度和不确定度。）准确度满足需求判断：依据用户对准确度既定阈值要求（均方根误差 以内为可接受的准确度；也可依据用户具体需求进行阈值设定），若高分辨率反照率影像满足阈值，可用于下一步检验；否则，重新准备高分辨率反照率影像并进行直接检验。）波段一致性判断：判断高分辨率反照率影像与待检反照率遥感产品波段是否一致。若波段一致，直接进行下一步检验；否则，应进行波段转换。）波段转换：
13、按附录将高分辨率反照率转换为与待检反照率遥感产品波段一致的反照率。）空间配准：将高分辨率反照率影像与待检反照率遥感产品进行空间配准。）空间尺度转换：按 中 的规定，将高分辨率反照率升尺度至待检反照率遥感产品相一致的空间尺度，作为待检反照率遥感产品的像元尺度相对真值。）准确度评价：按 中 的准确度评价指标定量计算反照率遥感产品的准确度，准确度评价指标应包括但不限于均方根误差、相关系数。）不确定度分析：按 中 的不确定度评价指标定量计算反照率遥感产品的不确定度，不确定度指标应包括但不限于标准差。基于已检反照率遥感产品的交叉检验法检验流程应符合 中 的规定，主要操作流程见图。图反照率遥感产品交叉检验
14、操作流程操作流程如下。）已检反照率遥感产品选择：应选用与待检产品时间一致的已检反照率遥感产品。犌犅犜 否波段转换=r一）空间配准：将已检反照率遥感产品与待检反照率遥感产品进行空间配准。）波段一致性判断：判断已检反照率遥感产品与待检反照率遥感产品波段是否一致。若波段一致，直接进行下一步检验；否则，应进行波段转换。）波段转换：按附录将已检反照率遥感产品转换为与待检反照率遥感产品波段一致的反照率。）空间一致性判断：根据待检反照率遥感产品像元大小，判断已检反照率遥感产品是否与其一致，若一致，可直接将已检反照率遥感产品作为相对真值；否则，应进行空间尺度转换。）空间尺度转换：按 中 的规定对已检反照率遥感
15、产品进行空间尺度转换，获取与待检反照率遥感产品像元尺度相一致的反照率，作为待检反照率遥感产品的像元尺度相对真值。）准确度评价：按 中 的准确度评价指标定量计算待检反照率遥感产品的准确度，准确度评价指标应包括但不限于均方根误差、相关系数。）不确定度分析：按 中 的不确定度评价指标定量计算待检反照率遥感产品的不确定度，不确定度指标应包括但不限于标准差。检验报告 封面信息检验报告封面应包括但不限于以下信息：）检验报告编号；）检验报告名称；）检验负责人；）检验核对人；）检验签发人；）检验单位；）送检单位；）检验时间。正文信息 待检反照率遥感产品概述对待检验的反照率遥感产品进行描述，应包括但不限于以下信
16、息：）产品的数据源、待检时空覆盖范围、时空分辨率、地理参考；）待检反照率遥感产品类型，如波段、观测几何等；）产品主要算法以及算法的特点。参考对象描述对作为参考对象的验证数据集进行描述，宜包括：）检验数据集的常规信息，如参考对象类型、基本信息地面测量仪器，地面测量仪器空间覆盖范围和记录频率，高分辨率反照率影像空间分辨率和日期（或景数），已检遥感产品时空分辨率、检验场概况、样本分布等；）检验数据集反照率类型，如波段、观测几何等。犌犅犜 检验方法及流程对采用的检验方法和检验过程进行描述，宜包括：）检验方法；）检验流程；）不确定性分析；）检验结果的评价指标；）检验结果。真实性检验结论对真实性检验的评价
17、结果进行描述，宜包括：）真实性检验结果评价：描述所检反照率遥感产品的准确度；）分项指标评价：描述所检反照率遥感产品在不同检验方法中的结果；）分析检验过程中的不确定性。附加信息对反照率遥感产品真实性检验过程中的非常规问题进行说明与描述，可包括对产品准确度或不确定度时空特点评价，对产品算法改进建议，对产品应用过程注意事项描述等。检验报告信息简表反照率遥感产品真实性检验报告简表编制见 的附录，样例见附录。犌犅犜 附录犃（资料性）空间代表性评价方法可利用空间分辨率高于待检反照率遥感产品的反照率遥感数据，通过地统计学中的半变异函数确定地面辐射测量站点的反照率空间代表性。提取待验证像元内的所有高分辨率反照
18、率像元犣（狓），犣（狓），犣（狓犻），在满足二阶平稳条件下，半变异函数由式（）计算：狉（犺）犖（犺）犖（犺）犻犣（狓犻）犣（狓犻犺）（）式中：犺 空间间隔距离，又称滞后距离，单位为米（）；狉（犺）距离相隔为犺的半变异函数值；犖（犺）距离相隔为犺的用于计算半变异函数值的像元对数；犣（狓犻）高分辨率像元位置狓犻的反照率；犣（狓犻犺）高分辨率像元位置狓犻犺的反照率。拟合半变异函数最常用的模型是基台值模型，包括球型模型（）、指数模型（）和高斯模型（）等。衡量变异程度的指标有基台值、变程和块金值三个参数。基台值反映地表反照率的变异强度，基台值越小，说明地表反照率均一性越好，反之，空间异质性越强。变程可反
19、映像元内地表反照率分布的最小异质性单元。块金值可反映地表反照率的不连续性。比较地面辐射测量站点的有效测量半径与待检反照率遥感产品像元分辨率。若地面的有效测量半径与像元分辨率相近，认为该站点的空间代表性较好，测量值可作为像元尺度相对真值。若地面的有效测量半径小于像元分辨率，比较半变异函数的变程与有效测量半径，如果有效测量半径与变程相近，则认为该站点空间代表性较好，测量值可作为像元尺度相对真值。犌犅犜 附录犅（规范性）上行下行辐射测量犅 测量设备犅 可用于测量辐射的仪器（如辐射表，或光谱仪加余弦积分器），测量仪器光谱范围宜为 。测量目标物（陆地地表或水体地表）的上行辐射和下行辐射，两者之比作为地面
20、测量的反照率。辐射表测量计算的反照率为宽波段反照率，光谱仪测量计算法的反照率为光谱反照率。犅 用于固定测量仪器的测量塔架。犅 测量要求犅 测量时间应依据待检反照率遥感产品时间确定，数据采集频率应不低于 一次。犅 测量人员应身着深色服装。犅 测量上行地表辐射与下行辐射的仪器应固定在塔架一端，保持水平。犅 宜尽量避免测量设备和测量人员阴影对目标物测量的影响。犅 测量仪器的架设高度宜距离样本表面大于。犌犅犜 附录犆（资料性）地面测量辐射数据处理犆 辐射数据选择方法犆 下行辐射阈值选择法绘制下行辐射值随时间的变化曲线，确定阈值，阈值以下的观测被认为是有云污染的观测，选择阈值以上的下行辐射观测值。犆 云
21、信息判断选择法利用云产品，提取对应时刻对应位置的云覆盖比例，云覆盖超过 的认为是有云。犆 直射光与散射光比选择法适用于当云量信息不可用，且有地面测量的直射辐射和散射辐射数据。计算地面测量的直射辐射与散射辐射的比值。比值小于认为是有云，选择比值大于时的地面测量下行辐射通量。犆 地表太阳入射辐射与理论晴空下行辐射比值选择法利用辐射传输模型模拟对应时刻对应位置处的晴空下行辐射通量，计算测量的太阳下行辐射通量与理论晴空下行辐射通量的比值。若比值大于，则选择该地面测量的下行辐射通量。犆 地面辐射数据与待检反照率遥感产品时间一致性处理地面测量辐射通量计算的真实反照率与反照率遥感产品保持时间一致性。若待检反
22、照率遥感产品为瞬时反照率，两者时间一致，宜取待检反照率遥感产品中真实反照率或黑天空反照率太阳入射角度下对应测量时刻前后半小时的辐射均值再进行后续处理，代表该时刻地面测量的真实反照率。若反照率遥感产品为合成反照率，两者时间不一致，对地面测量的上行下行辐射通量按遥感反照率产品的合成方法进行辐射通量合成，再后续处理代表地面测量的该合成时间范围的反照率。犌犅犜 附录犇（规范性）波段转换计算犇 光谱反照率与窄波段反照率转换当地面测量为光谱反照率时，按照公式（）计算窄波段反照率：（犻）狑犻（）（）狑犻（）（）式中：（犻）目标传感器波段犻的反照率；狑犻（）目标传感器波段犻的光谱响应函数；（）地面测量的连续光
23、谱反照率；目标传感器波段犻的起始波长，单位为纳米（）；目标传感器波段犻的终止波长，单位为纳米（）。犇 光谱反照率与宽波段反照率转换当地面测量为光谱反照率，按照公式（）计算宽波段反照率：（）犉（）（）犉（）（）式中：（）宽波段反照率；犉（）光谱反照率对应波长太阳下行辐射通量，单位为瓦每平方米纳米（）；（）地面测量的连续光谱反照率；宽波段反照率起始波长，单位为纳米（）；宽波段反照率终止波长，单位为纳米（）。犇 传感器之间波段转换具有各自波段响应特征的两个传感器波段之间转换，如将待转换传感器波段转换为目标传感器波段犻，可将待转传感器所有波段或相关波段通过经验系数加权，按照公式（）计算：（犻）犑犮犼（
24、犼）犫（）式中：（犻）目标传感器波段犻的反照率；（犼）待转换传感器波段犼的反照率；犮犼 待转换传感器波段犼的权重；犑 待转换传感器参与转换的波段总数；犫 常数。注：本公式也适用于窄波段反照率转换为宽波段反照率，可将宽波段反照率范围涵盖的窄波段或相关窄波段通过 犌犅犜 窄宽波段转换系数为权重加权窄波段反照率，其中，（犻）为宽波段反照率，犮犼为窄宽波段转换系数。犌犅犜 附录犈（规范性）真实反照率计算本文件中定义的地表反照率（）是一个广义概念，针对入射辐射的不同情况，可分为黑天空反照率（）、白天空反照率（）和真实反照率（）。黑天空反照率和白天空反照率可由遥感数据直接定量估算，反照率遥感产品一般仅提供
25、黑天空反照率和白天空反照率，不反演真实反照率。而真实反照率可由黑天空反照率和白天空反照率以天空散射光比例（）加权按式（）计算。地面测量计算得到反照率为真实反照率，反照率遥感产品进行直接检验时，应按照式（）由反照率遥感产品计算真实反照率进行检验。（狉）狉（）式中：真实反照率；黑天空反照率；白天空反照率；狉 天空散射光比例。犌犅犜 附录犉（资料性）天空散射光比例获取方法犉 直接测量方法利用朝上的辐射表，通过遮光环在太阳入射方向遮挡太阳直射辐射，测量天空散射辐射（可用太阳追踪装置）。用朝上的辐射表，直接测量总下行辐射。天空散射辐射除以总下行辐射得到天空散射光比例。在开阔空间开展测量，即无障碍物遮挡天
26、空，辐射表保持水平。犉 模型模拟方法通过或等大气校正软件，输入实测或大气遥感产品中的大气参数（水汽、气溶胶等）、太阳角度、地表高程等，模拟得到该大气状态下天空漫散射和下行总辐射，两者之比即为天空散射光比例。犉 经验计算方法通过经验公式计算天空散射光比例，见式（）：狉 （）（）式中：狉 天空散射光比例；太阳入射天顶角，单位为弧度（）。犌犅犜 附录犌（资料性）反照率遥感产品真实性检验报告信息简表样例表 给出了检验报告信息简表的参考样例，内容以直接检验法为例，同时也填写了间接检验法中的多尺度逐级检验法特征描述供参考。表犌 反照率遥感产品真实性检验报告信息简表样例 直接检验法待检反照率遥感产品特征描述
27、产品名称 反照率产品产品类型数值型类别型数据来源（传感器）、地理参考（椭球、投影方式等）坐标系，正弦曲线（）投影空间分辨率 时间分辨率天时间范围 年月日 年 月 日空间范围经度 纬度 产品生产算法采用 算法，天累积 、的晴空多角度地表反射率，反演标准传感器波段的核系数，模型在角度上积分得到标准传感器窄波段的黑天空和白天空反照率，窄宽波段转换得到宽波段反照率反照率类型直接检验方法中，由产品的正午时刻短波波段黑天空反照率和白天空反照率加权得到短波波段真实反照率进行检验；间接检验方法中，检验产品正午时刻短波波段黑天空反照率和短波波段白天空反照率参考对象特征描述类型地面测量数据地面测量数据和高分辨率反
28、照率影像已检反照率遥感产品基本信息地面测量仪器 仪器覆盖空间范围 数据记录频率 次高分辨率反照率影像空间分辨率 高分辨率反照率影像日期（或景数）年全年无云数据共计 景已检反照率遥感产品空间分辨率已检反照率遥感产品时间分辨率检验场概况（空间特征）选择检验场为黑河流域，地表覆盖为作物、果园、荒漠、沙漠、草甸等 犌犅犜 表犌 反照率遥感产品真实性检验报告信息简表样例 直接检验法（续）参考对象特征描述样本空间分布反照率类型地面测量数据为正午时刻短波波段真实反照率；高分辨率反照率影像为环境卫星正午时刻短波波段黑天空反照率，短波波段白天空反照率检验方法特征描述直接检验像元尺度相对真值获取方法通过空间代表性
29、评价选择可代表待检像元空间尺度的站点数据作验证不确定性分析仪器测量误差（全天辐射累积误差）间接检验法多尺度逐级检验法高分辨率反照率影像直接检验描述环境卫星正午时刻短波波段黑天空反照率和短波波段白天空反照率，计算短波波段真实反照率，通过地面测量数据进行直接检验高分辨率反照率影像准确度均方根误差为 ，相关系数 （样本量 ）空间尺度转换方法环境卫星反照率影像 聚合为 产品空间尺度波段一致性是否波段转换方法无不确定性分析环境卫星产品本身具有不确定性，升尺度后的像元尺度相对真值也含有这些不确定性，但不确定性的传递尚不明确交叉检验法已检反照率遥感产品准确度（不确定度）空间一致性判是否空间一致性转换方法波段
30、一致性是否波段转换方法不确定性分析样本量 样本选取方法覆盖不同值域：植被（作物、草甸）、荒漠、雪覆盖 犌犅犜 表犌 反照率遥感产品真实性检验报告信息简表样例 直接检验法（续）准确度评价及不确定度分析直接检验结果准确度评价均方根误差 相关系数 偏差 不确定度分析标准差 多尺度逐级检验结果准确度评价均方根误差 相关系数 偏差不确定度分析标准差 交叉检验结果准确度评价均方根误差相关系数偏差不确定度分析标准差附加信息（评价结果、改进建议等）从结果来看，产品总体准确度较高，由于多源数据引入，增加了观测频率，在观测数据充足，产品与地面的相关性表现非常好，特别雪覆盖时期，多数产品难以捕捉地表变化，该产品的相关系数可达 ；但是多源数据本身的不确定性也限制了产品的准确度，在观测数据不足时，会引入更多的误差。算法需要进一步对多源数据质量进行控制，更进一步提高产品的准确度。用户使用产品过程，可通过质量标识选择质量最优的反演结果，反演质量标识为差的结果准确度可能较低操作人员描述检验人单位联系方式时间年月日备注 犌犅犜 参考文献 光学遥感辐射传输基本术语陈述彭遥感大辞典北京：科学出版社，：，闻建光，刘强，柳钦火，肖青，李小文陆表二向反射特性遥感建模及反照率反演北京：科学出版社，：，犌犅犜