请留言
slide1

全球视野 媒体聚焦

接受央视报道共计15次 党媒、央媒、地方官媒总传播量达0.5亿人次

0101
新闻分类

基于西光壹号01星和西光壹号05星(天仙配号)高光谱卫星的在轨交叉定标

2025-07-18

1 (1).png

原理介绍  

交叉定标是利用定标精度较高的传感器作为参考,对待定标的传感器进行定标。其原理是选取对同一目标成像的同步或近似同步的影像对,在分析两个传感器光谱响应、观测几何、大气参数等匹配的基础上,建立两个传感器图像数字计数值之间的关系,利用参考传感器已知的辐射定标系数求解待标定传感器的定标系数。通常情况下,卫星遥感器的DN值与入瞳辐亮度(表观辐亮度,Top-of-Atmosphere Radiance, TOA)存在线性关系:L= gain·DN+ offset式中,gain和 ofset分别是定标系数的增益和截距,L表示传感器的入瞳辐亮度,DN表示图像的数字计数值。

数据

 

 

 西光壹号

 01星

 西光壹号

 05星(天仙配号)

 类型

 参考数据

 待标定数据

 波长范围(nm)

 441-854

 431-853

 波段数量

 150

 72

 FWHM(nm)

 1.945

 9.16

 GSD(m)

 40

 40

 轨道高度(km)

 484

 492

 幅宽(km)(@500km)

 80

 80

 流程

 1、影像对选择

1 (2).jpg

根据区域范围、成像时间以及观测角进行约束。具体为每景影像中心点距离在50km内,保证有大部分重叠区域。成像时间差异1天内,观测角差异小于2°。

 2、区域选择

1 (3).png

 西光壹号01星样本区域示例

1 (4).png

 西光壹号05星(天仙配号)样本区域示例

1 (5).png

 西光壹号01星样本区域平均辐亮度曲线示例

1 (6).png

西光壹号05星(天仙配号)样本区域平均DN值曲线示例  

选择柴达木盆地东部区域的均匀地表作为校正目标区域。最终筛选出16组数据。每组数据中,分别在西光壹号01星和   西光壹号05星(天仙配号)   数据中选择相同区域作为样本。  

  3、西光壹号01星光谱校正  

对西光壹号01星使用760nm氧气吸收谷进行光谱校正,校正中心波长位置以进行后续光谱重采样和光谱匹配因子(SMF)计算。  

  4、西光壹号01星辐射校正  

使用场地定标系数进行辐射校正,获取样本区域的表观辐亮度数据,每个样本计算区域内的辐亮度平均值。  

  5、   西光壹号05星(天仙配号)   暗电流校正  

 使用西光壹号05星(天仙配号)   夜间拍摄的沙漠区域数据进行暗电流计算,逐波段计算暗电流,以对原始影像进行暗电流校正,得到样本区域的校正后DN值,每个样本计算区域内的DN值平均值。  

  6、西光壹号01星对   西光壹号05星(天仙配号)   中心波长辐亮度拟合  

参考和待校正数据均为高光谱数据,其光谱响应函数(SRF)均为高斯型,但半高全宽(FWHM)有明显差异,会导致辐射偏差,尤其在光谱特征显著的波段(如吸收带)。因此需要进行光谱重采样和光谱匹配因子校正以弥补FWHM的差距。  

  SRF和光谱重采样:对西光壹号01星1.945nm数据使用高斯卷积模拟   西光壹号05星(天仙配号)   9.16 nm的SRF,减少光谱细节丢失,同时使用三次插值,匹配   西光壹号05星(天仙配号)   中心波长。接下来获取两颗卫星的SRF,选择不同地物标准目标光谱,使用辐射传输模型进行TOA模拟,以计算SMF,校正辐射偏差。  

  7、西光壹号01星表观辐亮度匹配   西光壹号05星(天仙配号)   DN值  

1 (7).png

 表观辐亮度统计表

1 (8).png

 DN值统计表

 使用SMF校正后的西光壹号01星TOA数据和   西光壹号05星(天仙配号)   DN值数据作为线性拟合的因变量和自变量,进行线性拟合,逐波段获取增益和偏置,共计120个样本。计算所有样本对于每个波段的增益和偏置平均值,作为该波段的最终增益和偏置。

  定标系数验证

1 (9).png

1 (10).png

1 (11).png

1 (12).png

 使用拟合出的   西光壹号05星(天仙配号)   定标系数进行辐射校正和大气校正,并与西光壹号01星TOA和地表反射率进行对比,以验证定标系数。

1 (13).png

 西光壹号01星和西光壹号05星(天仙配号)辐亮度对比, 单位W/m²/sr/μm

1 (14).png

 西光壹号01星和西光壹号05星(天仙配号)反射率对比(扩大10000倍)

通过与西光壹号01星对比,表观辐亮度和地表反射率趋势一致,说明交叉定标之后的定标系数可行性较高。

 交叉定标结果影像示例

1 (15).jpg

 西光壹号05星(天仙配号)孟加拉吉大港市附近影像 (2025年1月15日)

1 (16).png

 不同地物辐亮度曲线,单位W/m²/sr/μm

1 (17).png

 不同地物反射率曲线

1 (18).jpg

 西光壹号05星(天仙配号)法国科西嘉岛影像 (2025年7月15日)

1 (19).png

 不同深度海洋采样点示意图

1 (20).png

 不同深度海洋采样点示意图

 不同深度海洋辐亮度曲线,单位W/m²/sr/μm,编号增加表示深度增加

1 (21).png

 不同深度海洋反射率曲线

1 (22).png

 不同地物反射率曲线