GB/T 38935-2020.On-orbit radiometric characteristics assessment for optical imaging remote sensor-VIS-SWIR.
1范围
GB/T 38935规定了工作波长在可见光到短波红外范围(380nm~2500nm)的星载成像光学遥感器的在轨成像辐射性能评价对象、评价指标和评价方法。
GB/T 38935适用于搭载在卫星平台上,采用线阵列探测器扫描成像的被动光学遥感器在轨运行过程中进行的辐射性能评价。其他类型星载被动光学遥感器的在轨外场辐射性能评价可参考使用。
2规范性引 用文件
下列文件对于本文件的应用是必不可少的。凡是注日期的引用文件,仅注日期的版本适用于本文件。凡是不注日期的引用文件,其最新版本(包括所有的修改单)适用于本文件。
GB/T 33988-2017城镇地物可见光短波红外光谱反射率测量
GB/T 36297-2018光学遥感载荷性能外场测试评价指标
GB/T 36540-2018水体可见光-短波红外光谱反射率测量
3术语和定义、缩略语
3.1 术语和定义
GB/T 36297-2018界定的以及下列术语和定义适用于本文件。
3.1.1
辐射性能 radiometric characteristics
遥感器成像时获取与保持地物目标场景相对或绝对辐射能量分布的能力。
3.1.2
信噪比 signal-to-noise ratio
遥感器输出的有效信号功率与噪声功率之比。
3.1.3
辐射分辨率 radiometric resolution
遥感器感测(敏感)元件在接收波谱辐射信号时能分辨的最小辐射度差,或指对两个不同的辐射源的辐射量的分辨能力。
[GB/T 14950-2009,定义4.103]
3.1.4
动态范围 dynamic range
光学遥感器输出能够随输人发生变化的最小人射辐射量和最大人射辐射量所界定的区间。
3.1.5
非线性度 non-linearity
在动态范围内,响应的实际值偏离对应拟合直线的最大偏差。
本标准规定了工作波长在可见光到短波红外范围（380 nm~2500 nm）的星载成像光学遥感器的在轨成像辐射性能评价对象、评价指标和评价方法。 本标准适用于搭载在卫星平台上，采用线阵列探测器扫描成像的被动光学遥感器在轨运行过程中进行的辐射性能评价。其他类型星载被动光学遥感器的在轨外场辐射性能评价可参考使用。

ICS07.040
中华人民共和国国家标准
GB/T38935—2020
光学遥感器在轨成像辐射性能
评价方法
可见光-短波红外
On-orbit radiometric characteristics assessment foroptical imaging remote sensor-VIS-SWIR2020-07-21发布
国家市场监督管理总局
国家标准化管理委员会
2021-02-01实施
规范性引用文件
术语和定义、缩略语
术语和定义
缩略语
评价对象和评价指标
评价对象
评价指标
评价方法
信噪比
辐射分辨率
动态范围
非线性度
盲元率
附录A（资料性附录）
附录B（资料性附录）
参考文献
信噪比规定化方法
GB/T38935—2020
高空间分辨率光学遥感器在轨动态范围、非线性度评价地面靶标布设需求12
本标准按照GB/T1.1-2009给出的规则起草GB/T38935—2020
请注意本文件的某些内容可能涉及专利。本文件的发布机构不承担识别这些专利的责任本标准由中国科学院提出。
本标准由全国遥感技术标准化技术委员会（SAC/TC327)归口本标准起草单位：中国科学院光电研究院、中国资源卫星应用中心、国家卫星气象中心、北京空间机电研究所、中国科学院上海技术物理研究所、中国科学院长春光学精密机械与物理研究所。本标准主要起草人：李传荣、李晓辉、王新鸿、高彩霞、唐伶俐、马灵玲、王宁、傅俏燕、方翔、伏瑞敏马艳华、王钢、李伟、刘照言、赵永光、朱博、张静、朱家佳、刘耀开、钱永刚、邱实、周勇胜、朱小华、任璐m
1范围
光学遥感器在轨成像辐射性能
评价方法可见光-短波红外
GB/T38935—2020
本标准规定了工作波长在可见光到短波红外范围（380nm～2500nm）的星载成像光学遥感器的在轨成像辐射性能评价对象、评价指标和评价方法。本标准适用于搭载在卫星平台上，采用线阵列探测器扫描成像的被动光学遥感器在轨运行过程中进行的辐射性能评价。其他类型星载被动光学遥感器的在轨外场辐射性能评价可参考使用规范性引用文件
下列文件对于本文件的应用是必不可少的。凡是注日期的引用文件，仅注期的版本适用于本文件。凡是不注日期的引用文件，其最新版本（包括所有的修改单)适用于本文件GB/T33988一2017城镇地物可见光-短波红外光谱反射率测量GB/T36297—2018光学遥感载荷性能外场测试评价指标GB/T36540—2018水体可见光-短波红外光谱反射率测量3术语和定义、缩略语
3.1术语和定义
GB/T36297一2018界定的以及下列术语和定义适用于本文件3.1.1
辐射性能radiometriccharacteristics遥感器成像时获取与保持地物目标场景相对或绝对辐射能量分布的能力。3.1.2
信噪比signal-to-noiseratio
遥感器输出的有效信号功率与噪声功率之比。3.1.3
Eradiometricresolution
辐射分辨率
遥感器感测（敏感)元件在接收波谱辐射信号时能分辨的最小辐射度差，或指对两个不同的辐射源的辐射量的分辨能力。
[GB/T14950—2009，定义4.103]3.1.4
动态范围dynamicrange
光学遥感器输出能够随输人发生变化的最小入射辐射量和最大人射辐射量所界定的区间。3.1.5
non-linearity
非线性度
在动态范围内，响应的实际值偏离对应拟合直线的最大偏差1
GB/T38935—2020
注1：一般采用百分比表示，
注2：改写GB/T17444—2013，定义2.37。3.1.6
盲元率blindpixelratio
遥感器的盲像元数占总像元数的百分比。3.2
缩略语
下列缩略语适用于本文件。
VIS-SWIR
地面采样距离（groundsamplingdistance）可见光-短波红外（visibletoshort-waveinfrared）4评价对象和评价指标
4.1评价对象
本标准针对搭载于卫星平台上，工作波长在可见光到短波红外范围（380nm～2500nm），采用线阵列探测器扫描成像的被动光学遥感器（如可见光全色相机、多光谱相机、高光谱成像仪等），利用其在轨运行时获取的图像数据进行遥感器辐射性能分析与监测。4.2
评价指标
星载光学遥感器在轨成像辐射性能评价指标包括：信噪比、辐射分辨率、动态范围、非线性度、盲元率。
5评价方法
5.1信噪比
基于仅经过相对辐射校正处理的遥感器观测图像数据进行评价。具体步骤如下：a）选取遥感器在相近时间段内获取的若干景受大气影响小的、具有不同灰度等级（灰度等级≥5）的大面积均匀场景的图像作为样本图像，如深海、湖泊等水体，南极或格陵兰等地区的大范围冰雪覆盖区域，沙漠、戈壁或均匀浓密植被覆盖区域等。图像中的均匀区域应尽可能大（至少大于50像素×50像素），最好的情况是在沿遥感器阵列方向覆盖阵列所有像元，在沿遥感器阵列扫描方向（垂直遥感器阵列方向)大于或等于100像素。b）对于样本图像的某一波段l（l=1，2，，Nburd，Nhnd为遥感器波段数），选取图像中满足信噪比评测要求的某一均匀区域（区域内所有像素的灰度均值为D.，对应某一灰度等级k，k=1，2，，K，K为参与评测的均匀场景灰度等级数)作为均匀区子图像I.。若遥感器采用阵列摆扫方式成像，则利用式（1）对均匀区子图像I进行转置操作，式中：
pij..=pi.il..
........(
第l波段灰度等级k下均匀区子图像I..第i行第j列的灰度值。pii
pi.i——第/波段灰度等级k下均匀区于图像I.第j行第i列的灰度值。c）利用式（2)计算得到均匀区子图像I.的差值图像dij.i.. =pi+)j.f.pi.j..
式中：
GB/T38935—2020
第I波段灰度等级k下均匀区子图像I.的差值图像的第i行第；列的灰度值；第1波段灰度等级k下均匀区子图像I.第i十1行第i列的灰度值；d）利用式（3)计算第波段灰度等级k下均匀区子图像I的每一列的列噪声。(dij.）
式中：
のj——第l波段灰度等级k下均匀区子图像I.的第i列的噪声；Mi——第I波段灰度等级k下均匀区子图像I.的行数。e）利用式（4)计算第l波段灰度等级k下均匀区子图像I..每一列的信噪比。RsNRj
式中：
第1波段灰度等级k下均匀区子图像I.第j列的列信噪比。（3）
f）利用式（5)计算列信噪比的平均值作为第1波段灰度等级k下均匀区子图像I.的信噪比。RR
RsNRI=
式中：
第1波段灰度等级k下均匀区子图像I的信噪比；N.
一一第1波段灰度等级k下均匀区于图像I的列数（5）
g）按照步骤b)～f)计算样本图像中第1波段所有灰度等级下的均匀区子图像I，的信噪比。h）利用步骤g)得到的信噪比RsNR（k=1，2，，K.K为参与评测的均匀场景灰度等级数目），采用信噪比规定化方法（参见附录A），计算得到第1波段在某一参考入瞳辐亮度Lo.水平下的信噪比RsNR.
利用式（6)将步骤h)得到的信噪比RsNRo.转换为信噪比的功率表示。i）
RDBSNR.=20·Ig(RsNRa.）
式中：
RDBSNRO.t
第波段在某一参考人瞳辐亮度水平L。.下的信噪比的功率表示，单位为分贝（dB）。按照步骤b）i)得到遥感器所有波段的信噪比RpB-SNRa.r（l=1，2，，Nband，Nband为遥感器j
波段数）。
k）以至少1年1次的频率，按照步骤a)～j)计算遥感器各波段的信噪比。采用二维折线图形式（横轴为评价时段，纵轴为信噪比），显示遥感器各波段的信噪比随时间的变化特征与趋势（对于多光谱遥感器，可以将不同波段的信噪比变化折线显示在同一张折线图中）。5.2辐射分辨率
采用噪声等效辐亮度或噪声等效反射率表示，基于仅经过相对辐射校正处理的遥感器观测图像数据进行评价。具体步骤如下：
a）按照与信噪比评价相同的方法和要求，选择样本图像，同5.1中步骤a）。b）对于样本图像的某一波段l（l=1，2，，Nband，Nhand为遥感器波段数），按照与信噪比评价相同的方法计算遥感器第1波段在某一参考人瞳辐亮度水平Lo.下的信噪比RsNRo.！，同5.1中步3
GB/T38935—2020
骤b)～h）。
c）利用式(7)～式（9)计算第1波段的辐射分辨率。Lo.
RNEALn=
RsNRo-t
元Lrd2
po..=E..coso
式中：
.（7)
..（8)
...(9)
遥感器第I波段的噪声等效辐亮度，单位为瓦每平方米球面度微米[W/（m·sr·μm)]；RNEApn-
遥感器第波段噪声等效反射率；Lo.t
RsNRo,2
遥感器第波段的参考人瞳辐亮度，单位为瓦每平方米球面度微米［W/（m2·sr·μm)]；
与Lo,对应的遥感器第1波段的表观反射率值；遥感器第1波段在参考入瞳辐亮度L。，水平下的信噪比；日地距离因子（天文单位）；
遥感器第1波段的大气层顶太阳辐照度，单位为瓦每平方米微米［W/（m·μm)]太阳天顶角，单位为度（°）。d）
重复步骤b)～c)，完成遥感器的所有波段的辐射分辨率评测。以至少1年1次的频率，按照步骤a)～d)计算遥感器各波段的辐射分辨率。采用二维折线图e)
形式（横轴为评价时段，纵轴为辐射分辨率），显示遥感器各波段的辐射分辨率随时间的变化特征与趋势（对于多光谱遥感器，可以将不同波段的辐射分辨率变化折线显示在同一张折线图中）。
5.3动态范围
基于仅经过相对辐射校正处理的遥感器观测图像数据进行评价。具体步骤如下：a）根据遥感器的GSD，选择地势平坦、表面均一性好、具有不同反射率的试验区域（通常为特定的辐射定标场、验证场）或地面靶标（适用于高空间分辨率的遥感器，地面靶标的布设需求参见附录B)t,，t2，....，tk（K≥4）；且符合下列要求：1）遥感器响应输出达到饱和的高反射率试验区域或靶标至少有一个；2）遥感器响应输出在动态范围内、未达到饱和的高反射率试验区域或靶标至少有一个、低反射率试验区域或靶标至少有一个；3）试验区域或靶标的尺寸要尽可能大，其在图像中对应的均匀区域至少大于5像素×5像素。
b）选择晴朗无云、大气于洁（一般情况下能见度应大于23km）天气，在遥感器过境的同时开展星地同步试验，对试验区域或地面靶标进行地面参数数据和大气参量数据同步测量：1）按照GB/T33988—2017和GB/T36540—2018的规定，利用光谱测量仪对试验区域或靶标ti，t2，，tk（K≥4)测量地表反射率光谱数据；利用太阳光度计、地基激光雷达等大气测量设备同步获取遥感器成像时刻的大气水汽含2）
量、气溶胶光学特性参数等大气参量数据。c）根据遥感器成像时刻和卫星轨道相关参数，获取遥感器对试验区域或靶标成像时的太阳天顶角、太阳方位角、观测天顶角、观测方位角、卫星平台高度等辅助参数d）基于上述同步地面试验中测量的试验区域或靶标的地表反射率光谱数据，根据式（10)计算遥感器各波段的等效反射率：
式中：
入，入12
p(a)s,(a)d)
S,a)da
GB/T38935—2020
..（10)
试验区域或靶标t（k=1，2，，K，K≥4）在遥感器第I（l=1，2，\，Nhand，Nhand为遥感器波段数)波段的等效反射率；光谱测量仪测量的试验区域或靶标的地表反射率光谱数据；遥感器第1波段的归一化光谱响应函数；遥感器第1波段的波长范围，单位为微米（μm）。e）基于上述同步试验获取和计算的地表反射率光谱数据、大气参量数据以及辅助参数（太阳天顶角、太阳方位角、观测天顶角、观测方位角、卫星平台高度等），使用大气辐射传输模型模拟计算各试验区域或靶标t1t2，，tk在遥感器入瞳处各波段的模拟辐亮度L1，L2，,L.k（l1，2，，Nband），单位为瓦每平方米球面度微米[W/（m2·sr·μm)]。f）对同步试验中遥感器获取的各试验区域或靶标的观测图像数据的某一波段I（I=1，2，Nband，Nhand为遥感器波段数），利用式（11)计算各试验区域或靶标ti，t2，，tk在图像中对应均匀区于图像I..的平均灰度值D.（l=1，2，，Nhand；k=1，2，，K,K≥4)；Mi.t N!
Di.=M..XN..
式中：
.（11)
第1波段图像数据中第k个试验区域或靶标t对应均匀区子图像I的平均灰度值；pi.j.e
第I波段图像数据中第k个试验区域或靶标t对应均匀区子图像I.的第i行第列灰度值；
第1波段图像数据中第k个试验区域或靶标t对应均匀区子图像I.的行数；-第1波段图像数据中第k个试验区域或靶标t对应均匀区子图像I的列数。对波段1，选取位于载荷动态范围内、未达到饱和的试验区域或靶标t（k三1，，」，JK≥4），以横轴为遥感器人瞳处的模拟辐亮度L，纵轴为试验区域或靶标图像的平均灰度值D,,将步骤e)中得到的各试验区域或靶标的第I波段的模拟人瞳辐亮度（L.1，LI.2，,LI.J），与步骤f)中得到的遥感器所获取的各试验区域或靶标图像的第1波段的平均灰度值（D,,1，D.2，,D.)进行最小二乘线性拟合，在（L,D,)平面得到一条拟合直线LINE,：D,=G，·L，十B，，该拟合直线即为遥感器第I波段的线性响应曲线，而参数G，和B，分别为遥感器第1波段的增益和偏置，如图1所示。D, t
LINE;: D,-G,· L+B,
图1遥感器第1波段线性响应曲线拟合示意图5
GB/T38935—2020
h）对于步骤g)中拟合得到的直线LINE,，利用式（12)计算其相关系数的平方可以表征遥感器第1波段的拟合线性度。
Z(LIe-L,)(D.-D,)
式中：
拟合直线LINE的相关系数的平方；/(D.-D)
.（12)
遥感器响应输出在动态范围内未达到饱和状态、参与直线拟合的试验区域或靶标的数目；使用大气辐射传输模型模拟计算的第k（k=1，2，，J)个试验区域或靶标在遥感器人瞳处的第1波段的模拟辐亮度，单位为瓦每平方米球面度微米［W/（m2·sr·μm)]；L的平均值，单位为瓦每平方米球面度微米[W/（m2·sr·μm)]；D.
第1波段图像数据中第k个试验区域或靶标对应均匀区子图像I的平均灰度值；D,
-Di.的平均值。
在（L，，D，)平面，经过位于遥感器饱和区的点且平行于横坐标轴做一条直线，该直线与步骤g）i）
得到的拟合直线LINE，相交于点（Lsat.，DH.），则Lsat,l为遥感器第I波段的饱和辐亮度，亦即动态范围高端辐亮度值。将拟合直线LINE，延长与横坐标轴相交，交点的横坐标Lmin.即为遥感器第I波段的动态范围低端辐亮度值（此时对应的图像灰度值为DL.=0），表示为Lmin./=-B,
。（Lmin.，La）即为遥感器第I波段的动态范围，单位为瓦每平方米球面度微米G
[W/（m2·sr·μm)]，如图2所示。D,
动态范围区
图2动态范围示意图
饱和区bzxZ.net
j）重复步骤f)～i)，完成遥感器的所有波段的动态范围评测。k）以至少1年1次的频率，在相同试验场区开展遥感器动态范围测试外场试验，按照步骤a)j)计算遥感器不同波段下的动态范围高端与低端。采用二维折线图形式（横轴为评价时间，纵轴为动态范围高端Lat.），显示遥感器各波段的动态范围高端Lsat.随时间的变化特征与趋势（对于多光谱遥感器，可以将不同波段的动态范围高端L：显示在同一张图中）；采用二维折线图形式（横轴为评价时间，纵轴为动态范围低端Lmin.），显示遥感器各波段的动态范围低端Lmin.随时间的变化特征与趋势（对于多光谱遥感器，可以将不同波段的动态范围低端Lmin显示在同一张图中）。
小提示：此标准内容仅展示完整标准里的部分截取内容，若需要完整标准请到上方自行免费下载完整标准文档。