一种基于空海协同观测的海洋光学卫星辐射定标方法与流程

1.本发明涉及海洋卫星遥感在轨定标与处理领域，具体涉及一种基于空海协同观测的海洋光学卫星辐射定标方法。


背景技术：

2.随着定量遥感信息技术的高速发展，我国遥感产品质量亟待提高，特别是海洋光学遥感系列卫星。海洋光学遥感主要是通过海水的光谱反射、辐射特征(如反射率、辐射温度、或遥感专题指数)与准同步实测的各种水体参数(叶绿素含量、悬浮泥沙浓度、水深、水温等)的关系，建立一系列相关模型来提取或反演水体参数。海上原位光谱、辐射定标与检验是提高遥感定量评估精度的必要条件，也是提升光学遥感技术在海洋科学研究中应用能力的重要保证。海洋光学卫星的辐射定标方法主要包括实验室定标、星上定标以及场地定标。由于实验室定标与星上定标无法完全确定在轨实际工作过程中的定标参数变化，故场地定标的优势性更加突出。针对海洋光学卫星，其探测的信号比陆地遥感低一个数量级，因此定标场地只能选在海上，且以往的海面场地定标主要通过少量离散走航式站位原位观测，实际定标处理中可利用的匹配点较少，导致辐射定标精度受限。因此，需构建一种基于空海协同观测的海洋光学卫星辐射定标方法，通过点面高精度协同观测得到面分布的参考数据，从而提高海洋光学卫星的辐射定标精度，为后续向国际单位制标准溯源的定量化应用奠定基础。


技术实现要素：

3.本发明所要解决的问题是，针对以往海洋光学卫星场地定标主要采用船舶走航式原位观测，可利用的站位匹配点较少且导致定标精度不高，提出一种基于空海协同观测的海洋光学卫星辐射定标方法。
4.本发明的技术方案为一种基于空海协同观测的海洋光学卫星辐射定标方法，包括以下步骤：
5.步骤1，根据海洋光学遥感卫星的轨道与载荷特性，设计海洋原位-无人机协同观测模式与走航站位布设；
6.步骤2，根据原位观测数据反演无人机影像绝对离水辐射亮度产品；
7.步骤3，计算海洋光学卫星相对辐射校正产品；
8.步骤4，根据无人机绝对离水辐射亮度产品，计算海洋光学卫星绝对辐射校正系数。
9.所述确定海洋原位-无人机协同观测模式包括如下步骤：
10.(11)根据需要进行辐射定标的海洋光学卫星的轨道特性，计算走航船只过境某海域的轨迹与时刻，并利用便携式地物光谱仪，获取不同站位的海面离水辐亮度；
11.(12)根据海洋光学卫星载荷的波段范围与类型，选择相应的无人机相机载荷，以走航船只为中心，规划无人机飞行航线，获取观测海域一定范围海水的遥感反射率影像产
品，并以船舶为参考进行影像拼接。
12.所述便携式地物光谱仪与无人机采用同步协同观测，且与海洋光学卫星的观测时间差异控制在5h以内。
13.所述根据原位观测数据反演无人机影像绝对离水辐射亮度产品包括如下步骤：
14.(21)将便携式地物光谱仪测得的点位原位观测数据与无人机遥感数据进行空间匹配，以点位原位观测数据对应的无人机遥感产品像素点为中心，选取10
×
10的窗口内像素观测值进行平均拟合得到对应的海水遥感反射率；
15.(22)采用线性回归拟合模型构建原位观测数据与无人机影像产品转换关系，得到无人机影像绝对离水辐射亮度产品。
16.所述线性回归拟合模型为：
17.ρ
λi
(b,l)＝a
λi
l
λi
(b,l) b
λi
，
18.其中，ρ
λi
(b,l)表示波段i对应点位的海水遥感反射率，l
λi
(b,l)表示波段i对应点位的离水辐亮度值；a
λi
、b
λi
表示线性回归模型的拟合系数；
19.所述无人机影像绝对离水辐射亮度产品为无人机观测归一化离水辐亮度其中，ρ
wno
为无人机对应测得的海水遥感反射率。即先通过拟合模型得到a
λi
和b
λi
，然后反推得到l
wno
(λ)。
20.所述计算海洋光学卫星相对辐射校正产品包括如下步骤：通过对海洋光学卫星获取的海面观测数据进行分析，拟合出不同探元的光电响应模式并求取相应的辐射校正系数；采用的相对辐射校正公式表示如下：
[0021][0022]
其中，为载荷相对辐射校正后输出观测值，为原始图像第k号探元的dn值，bk为第k号探元的偏移值，ngk为第k号探元的归一化増益。
[0023]
上述公式中，bk、ngk是辐射校正系数，为是已知的数值，是目标数值。
[0024]
将与海洋光学卫星同步获取的无人机绝对离水辐射亮度产品作为外部真实参考，采用替代定标的方法计算绝对辐射校正系数，具体定标模型表示如下：
[0025][0026]
其中，l
t
(λ)表示卫星载荷入瞳处接收到的辐亮度；lr(λ)表示大气分子瑞利散射引起的光辐射；la(λ)表示气溶胶散射引起的光辐射；t(λ)表示大气漫射透过率；l
ru
(λ)表示瑞利散射和气溶胶散射相互作用引起的光辐射；lw(λ)表示离水辐亮度；a、b表示载荷实验室绝对辐射定标系数；表示载荷相对辐射校正后输出观测值；l
wns
(λ)，l
wno
(λ)分别表示海洋光学卫星以及无人机观测归一化离水辐亮度；δ表示系统定标整体限差，c为绝对辐射校正系数。
[0027]
上述公式中，lr(λ)、la(λ)、t(λ)、l
ru
(λ)、a、b、dn(λ)、l
wno
(λ)是已知的，l
wno
(λ)、分别是从步骤(2)、(3)中得到的。
[0028]
与现有技术相比，本发明的有益效果在于：
[0029]
本发明提供的空海协同观测的海洋光学卫星辐射定标方法，可以充分利用点面协同观测得到高精度面分布绝对辐射亮度参考数据，提高定标参数计算的匹配点数量，进一步提高海洋光学卫星辐射定标精度。
附图说明
[0030]
图1为本发明的流程示意图。
具体实施方式
[0031]
除非另作定义，在本说明书和权利要求书中使用的技术术语或者科学术语应当为本发明所属技术领域内具有一般技能的人士所理解的通常意义。本文中列举的所有的从最低值到最高值之间的数值，是指当最低值和最高值之间相差两个单位以上时，最低值与最高值之间以一个单位为增量得到的所有数值。
[0032]
以下将描述本发明的具体实施方式，需要指出的是，在这些实施方式的具体描述过程中，为了进行简明扼要的描述，本说明书不可能对实际的实施方式的所有特征均作详尽的描述。在不偏离本发明的精神和范围的情况下，本领域技术人员可以对本发明的实施方式进行修改和替换，所得实施方式也在本发明的保护范围之内。
[0033]
实施例
[0034]
下面将对本发明的实施例作详细说明，本实施例在以本发明技术方案为前提下进行实施，给出了详细的实施方式和具体的操作过程，但本发明的保护范围不限于下述的实施例。
[0035]
实施例1
[0036]
一种基于空海协同观测的海洋光学卫星辐射定标方法，其流程如图1所示，包括如下步骤：
[0037]
步骤1，根据海洋光学遥感卫星的轨道与载荷特性，设计海洋原位-无人机协同观测模式与走航站位布设。
[0038]
本发明首先根据需要进行辐射定标的海洋光学卫星的轨道特性计算过境某海域的轨迹与时刻，进一步设计走航原位观测站位，原位观测设备采用便携式地物光谱仪，获取不同站位的海面离水辐亮度；根据海洋光学卫星载荷的波段范围与类型，选择相应的无人机相机载荷，以走航船只为中心，规划无人机飞行航线，获取观测海域一定范围海水的遥感反射率影像产品，并以船舶为参考进行影像拼接；海洋原位与无人机采用同步协同观测，与海洋光学卫星的观测时间差异控制在5h以内。
[0039]
步骤2，根据原位观测数据反演无人机影像绝对离水辐射亮度产品。
[0040]
得到某海域一定范围内的原位海水遥感反射辐亮度点位数据以及无人机遥感反射率影像拼接产品，按如下步骤反演无人机遥感绝对离水辐亮度产品。
[0041]
步骤2.1，将点位原位观测数据与无人机遥感数据进行空间匹配，即空间上以原位数据对应的无人机遥感产品像素点为中心，选取10
×
10的窗口内像素观测值进行平均拟合
得到对应的反射率；
[0042]
步骤2.2，采用线性回归拟合模型构建原位观测数据与无人机影像产品转换关系，进一步得到无人机影像绝对离水辐射亮度产品。本发明采用的线性回归拟合模型表示如下：
[0043]
ρ
λi
(b,l)＝a
λi
l
λi
(b,l) b
λi
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(1)
[0044]
其中，ρ
λi
(b,l)，l
λi
(b,l)分别表示波段i对应点位的海水遥感反射率与离水辐亮度值；a
λi
，b
λi
表示线性回归模型的拟合系数。
[0045]
步骤3，计算海洋光学卫星相对辐射校正产品。
[0046]
相对辐射校正，是指为了校正各个探元响应度差异而对海洋光学载荷观测的原始数字计数值进行“再”量化的一种处理过程。由于载荷探元的光电响应函数之间存在一定差异，使得对同一亮度区域成像时，不同探元会输出不同dn(digital number)值，从而造成影像上出现贯穿推扫方向的条带现象。本发明通过对海洋光学卫星获取的海面观测数据进行分析，拟合出不同探元的光电响应模式并求取相应的辐射校正系数或查找表，用于消除由于探元的响应不一致引起的条带效应，使这些条带的影响降低到最小程度或彻底去除。采用的相对辐射校正公式表示如下：
[0047][0048]
其中，为相对辐射校正后第k号探元的dn值，为原始图像第k号探元的dn值，bk为第k号探元的偏移值，ngk为第k号探元的归一化増益，增益值在1.0左右，保证定标后的dn值范围与原始数值的范围基本一致。
[0049]
步骤4，根据无人机绝对离水辐射亮度产品，计算海洋光学卫星绝对辐射校正系数。
[0050]
基于步骤3得到的相对辐射校正结果，进一步计算海洋光学卫星绝对辐射校正系数。绝对辐射定标主要是建立载荷数字量化输出值与入瞳绝对辐射通量的定量关系。本发明将与海洋光学卫星同步获取的无人机绝对离水辐射亮度产品作为外部真实参考，采用替代定标的方法计算绝对辐射校正系数，具体定标模型表示如下：
[0051]
l
t
(λ)＝lr(λ) la(λ) t(λ)lw(λ) l
ru
(λ)
[0052]
l
t
(λ)＝c[a
·
dn(λ) b]
[0053][0054]
其中，l
t
(λ)表示卫星载荷入瞳处接收到的辐亮度；lr(λ)表示大气分子瑞利散射引起的光辐射；la(λ)表示气溶胶散射引起的光辐射；t(λ)表示大气漫射透过率；l
ru
(λ)表示瑞利散射和气溶胶散射相互作用引起的光辐射；lw(λ)表示离水辐亮度；a，b表示载荷实验室绝对辐射定标系数；dn(λ)表示载荷相对辐射校正后输出观测值；l
wns
(λ)，l
wno
(λ)分别表示海洋光学卫星以及无人机观测归一化离水辐亮度；δ表示系统定标整体限差。通过迭代计算绝对辐射校正系数c，以满足定标模型的整体限差要求。
[0055]
上述对实施例的描述是为了便于本技术领域的普通技术人员能理解和应用本技术。熟悉本领域技术的人员显然可以容易地对这些实施例做出各种修改，并把在此说明的
一般原理应用到其它实施例中而不必付出创造性的劳动。因此，本技术不限于这里的实施例，本领域技术人员根据本技术披露的内容，在不脱离本技术范围和精神的情况下做出的改进和修改都在本技术的范围之内。