一种波浪水面透射光的水下偏振场模拟方法
【专利说明】一种波浪水面透射光的水下偏振场模拟方法 (一） 所属技术领域
[0001] 本发明涉及一种波浪水面透射光的水下偏振场模拟方法，属于光学遥感领域，对 水下目标探测技术研究和水下偏振导航应用方面具有重要意义。 (二）
【背景技术】
[0002] 波浪水面作为一种复杂的自然场景，其形态在空间上时刻发生变化。而水下透射 光携带强烈的偏振信息，并且在Snell窗中具有规律的分布与变化特征。研究发现一些水 生动物可以利用透射光的偏振特性进行导航迀徙。由于光的偏振特性与光的传播方向有密 切关系，因此波浪水面的形态特征对水下透射光的偏振态分布有较大影响。随着光的偏振 信息在光学遥感领域应用的广泛与深入，本发明中波浪水面透射光偏振特性模拟的技术方 法对研究大气海洋親合系统矢量福射传输过程和水下偏振场的分布都有着重要的意义。
[0003] 目前水下透射光的偏振特性研究主要集中于测量手段和数值模拟的方法。而受 到探测设备和实验环境的限制，测量水下透射光的全偏振分布极其困难；现有的数值模拟 手段以蒙特卡洛方法为主，又存在效率低的缺点。现有的技术方法要么没考虑入射光的偏 振特性，将入射光视为自然光；要么没考虑波浪水面的干扰，仅研究了 Snell窗以内的偏振 场。实际上海洋一直处于波动中，Snell窗内外都存在偏振场，而且随着海面的波动，水下偏 振场始终处于动态变化过程中。现有的研究是对真实物理场描述的一种非常粗略的近似。 迄今为止，尚未有同时考虑入射光的偏振态与水面波动的水下透射光偏振场模拟方法的报 道，本发明建立的波浪水面透射光偏振分布模型填补了这一空白。本发明方法是一种解析 方法，没有数值计算误差，模型参数可调，提高了适用性，可为水下目标偏振探测和水生动 物的偏振场感应研究提供有效的技术方法。 (三）
【发明内容】
[0004] 本发明涉及一种波浪水面透射光的水下偏振场模拟方法。技术解决方案是：建立 瑞利散射天空光的偏振模型，计算天空光的偏振态；建立波浪水面面元斜坡的概率分布模 型，统一入射光和透射光的参考坐标系；考虑风成浪的遮挡，得到波浪水面折射的Mueller 矩阵；结合Stokes参量表示法，最终可以模拟水下透射光的偏振态分布。其具体步骤如下：
[0005] 1 -种波浪水面透射光的水下偏振场模拟方法。其特征在于包含以下步骤：
[0006] (1)建立天空光偏振场坐标系，设定太阳方位参数和观测几何参数；
[0007] (2)在半球天空上对入射天空光方向和观测位置进行空间立体角采样；
[0008] (3)根据太阳入射光与散射光方向的几何关系计算散射角；
[0009] (4)计算入射天空光的偏振态；
[0010] (5)建立波浪水面面元坡度概率分布模型；
[0011] (6)建立波浪水面透射过程的矢量辐射传输模型；
[0012] (7)计算波浪水面透射天空光的偏振态；
[0013] (8)对采样点进行加权平均从而实现波浪水面下行透射光的偏振场的模拟。
[0014] 2步骤（1)中所述的"建立天空光偏振场坐标系，设定太阳方位参数和观测几何参 数"，确定入射光和散射光的0°参考方向，要求E矢量参考平面为入射光或散射光所在的 子午面。
[0015] 3步骤⑵中所述的"在半球天空上对入射天空光方向和观测位置进行空间立体 角采样"，采样点要涵盖整个半球天空，确定入射光以及透射光的天顶角和方位角间隔。
[0016] 4步骤（3)中所述的"根据太阳入射光与散射光方向的几何关系计算散射角"，具 体计算过程如下：
[0017] 第一步：根据天空光偏振场坐标系确定入射光和散射光的天顶角和方位角大小；
[0018] 第二步：根据公式cosy =CosGi 〇)8θ+??ηθ、sin0 cosqH十算散射角γ。其中，Θ 3为 太阳天顶角，Θ为观测天顶角，φ为相对方位角。
[0019] 5步骤（4)中所述的"计算入射天空光的偏振态"，具体计算过程如下：
[0020] 第一步：利用半经验瑞利散射模型计算天空光的偏振度，
[0022] 其中δ _为经验值，最大取1 ;
[0023] 第二步：设POS为散射平面，其中0点为坐标原点，其坐标为（0, 0, 0) ;S点为太阳 位置，
[0024] 其坐标为（Xs，Ys，Zs) ;Ρ点为遥感器观测位置，其坐标为（XP，YP，Zp)，则
[0025] OS = (Xs, Ys, Zs), OP = (XP, YP, ΖΡ)
[0026] 由于E矢量垂直于平面P0S，则E矢量垂直于向量OS和0Ρ，设E = (X，Υ, 1)，则有
[0027] X · Xs+Y · Ys+Zs= 0
[0028] X · Xp+Y · Yp+Zp= 0
[0030] 则偏振角余角的余弦为
[0032] ON为P点所在子午面的法向量；
[0033] 第三步：最后得到计算偏振角的表达式：
[0035] 6步骤（5)中所述的"建立波浪水面面元斜坡概率分布模型"，具体计算过程如下：
[0036] 第一步：建立波浪水面透射过程的坐标系；
[0037] 第二步：根据入射方向、折射方向、面元法向量的空间关系得到波浪面元的坡度在 X，y两坐标轴的分量：
[0038] CN 105181145 A 说明书 3/4 页
[0039] 其中s表示入射向量，r表示折射（观测）向量；
[0040] 第三步：基于Cox-Munk模型计算波浪面元坡度分布概率：
[0042] 7步骤（6)中所述的"建立了波浪水面透射过程的矢量辐射传输模型"，具体计算 过程如下：
[0043] 第一步：计算透射方向的辐射通量：
[0044] (1Φ t= t (ω) E s cos ω sec β AP
[0045] 其中A为面元面积，β为面元法线与天顶之间的夹角，ω为入射角大小，Es表示 入射太阳辐照度，t(co)为入射角为ω时水面的透过率；
[0046] 第二步：考虑波浪面元遮挡问题，提出一个修正因子S ;
[0047] 第三步：旋转参考面，旋转矩阵为：
[0049] 其中i为参考面旋转角度；
[0050] 第四步：最后得到波浪面元透射的Mueller矩阵：
[0052] 其中M为平静水面透射的Mue 11 er矩阵。
[0053] 8步骤（7)中所述的"计算波浪水面透射天空光的偏振态"，具体计算过程：根据步 骤⑷计算出的入射天空光的偏振态，乘以步骤（6)得到的波浪水面透射的Mueller矩阵， 就得到透射光的偏振态。
[0054] 9步骤（8)中所述的"对采样点进行加权平均从而实现波浪水面下行透射光的偏 振场的模拟"，具体计算过程：根据步骤（2)中的天空光采样点，在每一个观测方向上逐个 计算入射天空光发生折射时的Stokes参量，再进行加权平均得到该观测方向上的透射光 的偏振态。逐步计算每个观测方向上透射光的偏振态就实现了水下偏振场的模拟。
[0055] 本发明与现有技术相比的优点在于：
[0056] (1)目前学术界缺少对波浪水