一种多角度偏振水色遥感器卫星大气校正方法与流程

1.本发明属于遥感图像处理领域，具体涉及一种多角度偏振水色遥感器卫星大气校正方法。


背景技术：

2.大气是影响水色遥感定量分析与应用的重要因素，在大气顶所接收到的辐射强度信息中，大气信号占据90%左右，因此消除大气效应是决定水色遥感分析精度的重要前提。偏振图像的大气校正，对于获取水体真实的偏振信息有着重要的意义。
3.h.xianqiang等人在the atmospheric correction algorithm for hy-1a/cocts, spie remote sensing (spie, 2005), vol. 5977一文叙述了我国 hy-1a/cocts海洋遥感卫星所采用的大气校正方法，该方法的主要的缺陷在于仅考虑了离水辐射强度分量的大气校正过程，没有考虑离水辐射偏振分量在大气中传输时所对应的大气漫射透过率，然而离水辐射偏振分量对于充分利用离水辐射信息探究水体组分以及海况情况非常重要，如果针对其的大气校正过程不够准确，会极大的影响后续的相关水色组分探究。


技术实现要素：

4.针对现有技术中存在的问题，本发明提供了一种多角度偏振水色遥感器卫星大气校正方法，针对离水辐射线偏振分量进行了准确的大气校正。具体技术方案如下：一种多角度偏振水色遥感器卫星大气校正方法，主要包括以下步骤：s1：利用辐射传输模型，模拟计算海水的光学特性为全吸收性水体时的海面之上的线偏振分量s
1,0+
以及大气顶处的线偏振分量s
1,toa
；s2：利用辐射传输模型，模拟计算不同的影响因子取值条件下的水面之上的线偏振分量s
2,0+
，及大气顶处的线偏振分量s
2,toa
；s3：根据步骤s1、s2的模拟计算结果，按照公式（1）计算得到不同的影响因子取值状态下离水辐射线偏振分量的大气漫射透过率ts，从而构建由不同的影响因子取值状态以及相应的离水辐射线偏振分量的大气漫射透过率ts构建形成的索引查找表；
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
（1）s4：应用上述索引查找表对多角度偏振水色遥感器采集的遥感数据进行大气校正。
5.其中，辐射传输模型选自osoaa辐射传输模型、pcoart辐射传输模型或momo辐射传输模型，也可以使用其他的辐射传输模型，使用的辐射传输模型应当考虑了偏振情况，并且考虑了海气耦合。索引查找表主要围绕中低浑浊度水体进行构建，遍历中低浑浊度水体条件下所选取影响因子的不同取值状态，从而使所述索引查找表具有实际操作价值。
6.离水辐射线偏振分量的大气漫射透过率ts可能由多种影响因子决定。但发明人通
过研究发现，所述影响因子主要包括气溶胶模式sfmodel、气溶胶光学厚度τa、太阳-遥感器观测几何（包括遥感器观测天顶角θv，太阳天顶角θ0，以及太阳-传感器观测相对方位角φ）和模拟波长λ。因此，优选以上述6项参数作为影响因子，并计算得到各种取值情况下相应的的ts，从而构建所述索引查找表。
7.具体的，所述步骤s4中应用上述索引查找表对多角度偏振水色遥感器采集的遥感数据进行大气校正，至少包括以下步骤：s41：导入多角度偏振水色遥感器采集的多个波段的stokes矢量辐射强度，扣除水面白沫反射辐亮度，并进行臭氧双层吸收修正，得到经修正的stokes矢量辐射强度l
′
t
(λ)，其中λ表示模拟波长；s42：利用辐射传输模型，得到rayleigh散射辐亮度lr(λ)；s43：计算大气程辐射辐亮度l
path
(λ)；s44：确定所述影响因子的取值，根据所述索引查找表确定对应的离水辐射线偏振分量的大气漫射透过率ts，获得水面之上的归一化离水辐射强度的线偏振分量。
8.在校正得到归一化离水辐射强度的线偏振分量后，即可根据根据经验公式，计算并发布水体悬浮颗粒物浓度以及气溶胶光学厚度、浓度以及粒径谱分布产品。
9.通过本发明提供的大气校正方法，在校正过程中充分考虑了离水辐射线偏振分量的大气漫射透过率，对离水辐射线偏振分量进行了准确的大气校正，从而使卫星遥感探测的图像质量更高，探测结果更加准确。
附图说明
10.通过结合附图对本公开示例性实施例进行更详细的描述，本公开的上述以及其它目的、特征和优势将变得更加明显，其中，在本公开示例性实施例中，相同的参考标号通常代表相同部件。
11.图1示出了不同的水体无机颗粒物浓度下线偏振分量大气漫射透过率的误差变化情况；图2示出了实施例中构建的索引查找表的插值精度示意图。
具体实施方式
12.下面将参照附图更详细地描述本公开的实施例。虽然附图中显示了本公开的实施例，然而应该理解，可以以各种形式实现本公开而不应被这里阐述的实施例所限制。相反，提供这些例举性的实施例是为了使本公开更加透彻和完整，并且能够将本公开的范围完整地传达给本领域的技术人员。
13.在本文中使用的术语“包括”及其变形表示开放性包括，即“包括但不限于”。除非特别申明，术语“或”表示“和/或”。术语“基于”表示“至少部分地基于”。术语“连接”、“连通”表示直接或通过其他组件间接的连接或连通。术语“第一”、“第二”等可以指代不同的或相同的对象，但并不直接表示先后顺序或重要程度的不同。下文还可能包括其他明确的和隐含的定义。
14.本发明提供了一种多角度偏振水色遥感器卫星大气校正方法，该校正方法的第一个方面主要在于构建由不同的影响因子取值状态以及相应的离水辐射线偏振分量的大气
漫射透过率ts构建形成的索引查找表。以下用osoaa辐射传输模型为例进行说明：1. 导入osoaa辐射传输计算模型;2. 设置sea.depth=0.05，sea.bottype=1，sea.botalb=0，phyto.chl=0，sed.csed=0，det.abs440=1e5，ys.abs440=1e5，计算全吸收性水体(black ocean)下的水面之上以及大气顶的stokes矢量辐射强度，其中各参数的含义如下：sea.depth：水深 (m)；sea.bottype：底质的反照率类型，取值为1时表示由用户自己设定具体的底质反照率值；sea.botalb：底质反照率具体的值，0表示全吸收，不发生底质的反射；phyto.chl：海表叶绿素浓度 (mg/m3)；sed.csed：海表无机颗粒物浓度 (mg/l)；det.abs440：有机碎屑440 nm吸收系数 (m-1
)；ys.abs440：黄色物质 440 nm吸收系数(m-1
)。
15.3. 设置sea.depth=100，sea.bottype=1，sea.botalb=0，phyto.chl=0.05，sed.csed=0，det.abs440=0，ys.abs440=0，计算大洋清洁水体下的水面之上以及大气顶的stokes矢量辐射强度；4. 设置sea.depth=100，sea.bottype=1，sea.botalb=0，phyto.chl=2.0，sed.csed=0，det.abs440=0.5，ys.abs440=0.5，计算富营养化水体下的水面之上以及大气顶的stokes矢量辐射强度；5. 设置sea.depth=100，sea.bottype=1，sea.botalb=0，phyto.chl=0，sed.csed=500，det.abs440=0，ys.abs440=0，计算浑浊水体下的水面之上以及大气顶的stokes矢量辐射强度；6. 利用以上步骤计算得到的stokes矢量辐射强度以及步骤2得到的stokes矢量辐射强度，根据ts的计算公式，计算得到清洁、富营养化、浑浊水体下的离水辐射偏振分量大气漫射透过率分布；7. 改变大气气溶胶参数条件，重复步骤2-6，计算不同的大气条件下，清洁、富营养化、浑浊水体下的离水辐射偏振分量大气漫射透过率分布情况，并由此发现离水辐射偏振分量大气漫射透过率的主要影响因子是气溶胶模式(sfmodel)、气溶胶光学厚度(τa)、太阳-传感器观测几何(θ0，φ，θv)、模拟波长(λ)；8. 遍历中低浑浊度水体条件下的离水辐射偏振分量大气漫射透过率的主要影响因子范围，利用合适的影响因子步长，构建针对中低浑浊度水体的离水辐射线偏振分量大气漫射透过率的索引查找表。
16.图1示出了在不同的水体无机颗粒物浓度(c
sed
=1.00, 10.00, 100.00, 300.00 mg l-1
) 下，stokes矢量辐射强度的q和u线偏振分量大气漫射透过率相对于c
sed = 10 mg l-1
时的相对误差变化情况，其中，re
q|q1
表示q线偏振分量的相对误差，re
u|u1
表示u线偏振分量的相对误差，c50气溶胶模式表示相对湿度50%的大陆型气溶胶；m99气溶胶模式表示相对湿度99%的海洋型气溶胶模式。可见，当无机颗粒物浓度发生变化时，stokes矢量辐射强度的q和u线偏振分量大气漫射透过率发生明显的变化，无机颗粒物浓度是离水辐射线偏振分量的主要影响因子之一。
17.选择正确的影响因子，并利用合适的影响因子步长，可以构建形成具有实际操作价值的查找表。如图2所示，通过以上具体实施方式构建形成的索引查找表具有很高的插值模拟精度。其中，横轴表示查找表插值结果，纵轴表示osoaa辐射传输模型反演计算结果。图
2给出了在441nm（图2中的a、c）和549nm（图2中的b、d）下的反演结果。可以发现，本算法所构建的针对中低浑浊度水下的离水辐射线偏振分量大气漫射透过率查找表精度较高，1000次monte-carlo模拟结果rmse《0.05，同时r2》0.95，查找表插值结果和osoaa辐射传输模型反演计算结果具有很高的一致性。
18.本发明提供的校正方法的第二个方面在于应用上述索引查找表对多角度偏振水色遥感器采集的遥感数据进行大气校正，主要包括以下步骤：s41：导入多角度偏振水色遥感器采集的多个波段的stokes矢量辐射强度，扣除水面白沫反射辐亮度，并进行臭氧双层吸收修正，得到经修正的stokes矢量辐射强度l
′
t
(λ)，其中λ表示模拟波长；s42：利用辐射传输模型，得到rayleigh散射辐亮度lr(λ)；s43：计算大气程辐射辐亮度l
path
(λ)；s44：确定所述影响因子的取值，根据所述索引查找表确定对应的离水辐射线偏振分量的大气漫射透过率ts，获得水面之上的归一化离水辐射强度的线偏振分量。
19.本具体实施方式以osoaa模型为例进行示例性的说明：步骤s41中的水面白沫反射辐亮度具体通过公式（2）确定：
ꢀꢀꢀ
（2）其中，t
iv
(λ)l
wc
(λ)为水面白沫反射辐亮度，l
wc
(λ)为水面白沫反射在水面之上所对应的辐亮度，c(λ)为水体吸收修正系数，r
wc
(λ)为白沫反射率，f0(λ)为大气顶太阳辐照度，θ0为太阳天顶角，θv为遥感器观测天顶角，t
iv
(λ)为stokes矢量辐射强度的i分量从水面到遥感器的大气漫射透过率，t
is
(λ)为i分量从大气顶到水面的大气漫射透过率；其中，水体吸收修正系数c(λ)的取值如表1所示：[表1] 白沫反射水体吸收修正系数波段波段12345678c值1.01.01.01.01.00.8892250.7600460.644950白沫反射率r
wc
(λ)由公式（3）计算得到：
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
（3）式中，w为水面10 m高处的风速，单位为m/s；t
iv
(λ)和t
is
(λ)由公式（4）计算得到：
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
（4）式中，τ
oz
(λ)为表示臭氧吸收光学厚度，τr(λ)为rayleigh散射光学厚度。
[0020]
具体的，步骤s41中的臭氧双层吸收修正通过公式（5）进行修正：
ꢀꢀ
（5）式中，l
t
(λ)表示多角度偏振水色遥感器采集的多个波段的stokes矢量辐射强度，l
′
t
(λ)表示经修正的stokes矢量辐射强度。
[0021]
具体的，所述步骤s42中通过以下方法得到rayleigh散射辐亮度lr(λ)：利用预先生成的rayleigh散射查找表，根据当前像元的海面风速w、太阳天顶角θ0、遥感器观测天顶角θv和太阳-传感器观测相对方位角φ，插值并进行气压修正，得到相应的rayleigh散射辐亮度lr(λ)。
[0022]
具体的，所述步骤s43中通过以下方法计算得到大气程辐射辐亮度l
path
(λ)：预设近红外两个波段λ
nir1
和λ
nir2
的离水辐亮度为0，其中λ
nir1
小于λ
nir2
，则此两波段的大气程辐射辐亮度和rayleigh散射辐亮度的比值为：
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(6)利用预先生成的气溶胶散射查找表，根据当前像元的太阳天顶角θ0、遥感器观测天顶角θv和太阳-传感器观测相对方位角φ，插值获得20种气溶胶模式对应的l
path
(λ)/lr(λ)展开为气溶胶光学厚度τa的二次多项式系数a(λ)、b(λ)和c(λ)：
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(7)根据r(λ
nir2
)，计算获得20种气溶胶模式对应的λ
nir2
波段气溶胶光学厚度τa(λ
nir2
)：
ꢀꢀꢀꢀꢀꢀ
(8)由τa(λ
nir2
)计算20种气溶胶模式对应的第λ
nir1
波段气溶胶光学厚度τa(λ
nir1
)：
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(9)式中，c
ext
为气溶胶衰减截面；由τa(λ
nir1
)和系数a(λ
nir1
)、b(λ
nir1
)、c(λ
nir1
)，根据公式(7)计算20种气溶胶模式对应的第λ
nir1
波段大气程辐射和大气rayleigh散射辐亮度比值r
model
(λ
nir1
)；根据r(λ
nir1
)值，从20种气溶胶模式中找到r
model
(λ
nir1
)最接近r(λ
nir1
)值的两种气溶胶模式mod1和mod2，并由公式(10)计算第一种气溶胶模式的权重：
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(10)根据公式(9)，由τa(λ
nir2
)分别计算气溶胶模式mod1和mod2对应的可见光波段气溶胶光学厚度τ
amod1
(λ)、τ
amod2
(λ)，再根据公式(7)计算获得对应的r
mod1
(λ)、r
mod2
(λ)值；利用两种气溶胶模式的权重，插值获得可见光波段的气溶胶光学厚度τa(λ)以及大气程辐射与rayleigh散射辐亮度的比值r(λ)：
ꢀꢀꢀꢀꢀꢀꢀ
(11)
由r(λ)及lr(λ)，计算得到可见光波段大气程辐射辐亮度l
path
(λ)。
[0023]
具体的，步骤s44中确定所述影响因子的取值，根据所述索引查找表确定对应的离水辐射线偏振分量的大气漫射透过率ts，是指：根据遥感器观测天顶角及两种气溶胶模式mod1、mod2对应的气溶胶光学厚度τ
amod1
(λ)、τ
amod2
(λ)，查找所述索引查找表并插值获得气溶胶模式mod1和mod2的上行大气漫射透过率t
vmod1
(λ)、t
vmod2
(λ)，并利用两种气溶胶模式的权重插值获得上行大气漫射透过率tv(λ)：
ꢀꢀꢀꢀꢀꢀꢀꢀ
(12)根据太阳天顶角θ0及两种气溶胶模式mod1、mod2对应的气溶胶光学厚度τ
amod1
(λ)、τ
amod2
(λ)，查找所述索引查找表并插值获得气溶胶模式mod1和mod2的下行大气漫射透过率t
smod1
(λ)、t
smod2
(λ)，并利用两种气溶胶模式的权重插值获得下行大气漫射透过率ts(λ)：
ꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(13)将所述上行大气漫射透过率tv(λ)和下行大气漫射透过率ts(λ)作为离水辐射线偏振分量的大气漫射透过率ts。
[0024]
具体的，步骤s44中通过以下公式获得水面之上的归一化离水辐射强度的线偏振分量：
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(14)式中，l
wn
(λ)表示水面之上的归一化离水辐射强度的线偏振分量，r是瞬时日地距离，r是平均日地距离，lw(λ)通过公式（15）确定：
ꢀꢀꢀꢀꢀꢀꢀ
（15）。
[0025]
在得到归一化离水辐射强度的线偏振分量后，即可根据经验公式，按照本领域已有的方法计算并发布水体悬浮颗粒物浓度以及气溶胶光学厚度、浓度以及粒径谱分布产品，本发明在此不作赘述。
[0026]
以上已经描述了本公开的各实施例，本领域技术人员应当理解上述说明仅仅是示例性的，并非穷尽性的，并且也不限于所披露的各实施例。在不偏离所说明的各实施例的范围和精神的情况下，对于本技术领域的普通技术人员来说许多修改和变更都是显而易见的。本文中所用术语的选择，旨在最好地解释各实施例的原理、实际应用或对现有技术的改进，或者使本技术领域的其它普通技术人员能理解本文披露的各实施例。