一种基于双视角光学遥感图像的海面风速探测方法技术

本发明专利技术提供了一种新的利用双视角光学遥感图像探测海面风速的方法。该方法利用同步或准同步观测获得的双视角光学遥感图像记录的海面粗糙度几乎不变的特征，通过比值法，建立了一种新的双视角太阳耀光辐射强度海面粗糙度之间的定量关系，实现了海面粗糙度和海面风速的逐像元计算。该方法可用于高分辨率的海面风速信息，是遥感信息技术用于海面风速探测的中项创新，具有极大的实用价值。

【技术实现步骤摘要】

本专利技术属于海洋遥感技术应用领域，具体是一种能够利用双视角光学遥感图像中所含的不同灰度信息，如所希望那样测量海面风速的方法。
技术介绍
海面风速是基础的海洋环境要素，海面风速探测是认识海洋环境基础工作，也是现今海洋学研究的迫切需求，对于应对全球气候变化、保障海洋环境安全都有十分重要的意义。面向全球广阔的海洋的海面风速探测需求，离散分布的浮标和气象站已难以满足要求。遥感方法已成为全球海面风场探测的主要方法，其基本原理是：海面风速驱动下会引起海面粗糙度的变化，不同的海面粗糙度可以被遥感器探测并记录为不同的强度信息。相关业务化工作中已经建立了海面粗糙度和海面风速的定量关系，因此海面风速遥感探测的关键是海面粗糙度的定量探测。微波散射计和微波辐射计是当前最主要的两种风速探测遥感器，两者的基本原理都是电磁散射强度与海面粗糙度相关。这两种方法能业务化获得风场数据，其探测误差为2m/s左右，然而其海面风速的探测分辨率为数十千米，难以满足高分辨率风速的需求。高分辨率合成孔径雷达技术的发展，提供了一种高分辨率海面风速探测的手段，但其探测范围较小。海洋光学遥感器应用中，太阳耀光信号作为一种污染或干扰信息，被掩膜处理掉。前人研究表明海洋光学遥感图像中的太阳耀光信号强度与海面粗糙度相关，但受到太阳耀光辐射传输复杂过程的影响，难以直接利用太阳耀光辐射信息号探测海面粗糙度，也就难以完成海面风速探测。仅有的工作是基于拟合迭代方法开展了区域性的海面风速光学遥感探测实验，但仍然无法实现逐像元的高分辨率风速探测。因此，从检索到的公开资料看，尚未有与本专利技术完全一致的实现基于双视角光学遥感图像的海面风速探测方法。本专利技术利用遥感技术方法的大范围覆盖和高频率重访的能力，基于双视角光学遥感图像中所包含的不同辐射信息，通过信息处理计算获得海面风速信息。
技术实现思路
本专利技术的目的是提供了一种新的基于双视角光学遥感图像的海面风速探测方法。本专利技术通过以下技术方案实现：一种基于双视角光学遥感图像的海面风速探测方法，其特征在于其包括以下步骤：(1)选取由同一遥感平台同时或准同时(相隔时间很短，指间隔时间小于5分钟)成像的双视角光学遥感图像，根据观测方向与遥感平台运动方向的关系，分别记为前视遥感图像(N)和后视遥感图像(B)。(2)根据遥感图像的成像时间确定太阳高度角(θ0)和方位角(3)根据遥感图像文件所记录的遥感图像获取时的轨道倾角(S)和侧视角(P)，分别逐像元地计算前视遥感图像的传感器高度角(θN)和方位角后视遥感图像的传感器高度角(θB)和方位角(4)对前视遥感图像和后视遥感图像分别进行大气校正，获得前视太阳耀光辐射强度(LgN)和后视太阳耀光辐射强度(LgB)。(5)利用太阳耀光辐射传输模型和海面波面倾斜概率模型，分别建立前视太阳耀光辐射强度(LgN)和后视太阳耀光辐射强度(LgB)与海面粗糙度的关系。(6)采用比值法，逐像元地计算获得海面粗糙度。(7)利用海面粗糙度和海面风速的定量关系，计算海面风速，获得整幅图像的风速信息。所述的一种基于双视角光学遥感图像的海面风速探测方法，其特征在于步骤(2)中，根据获取的遥感图像的经纬度位置和成像时间，逐像元计算太阳天顶角(θ0)和方位角太阳天顶角(θ0)的计算方法为：θ0＝90°-h0sinh0＝sinγsinδ0+cosγcosδ0cost0公式中，各参数确定如下：①h0为太阳高度角；②γ为当地纬度，即遥感图像各像元的纬度；③t0为当地太阳时角，由遥感图像时的地方时确定，中午12时为0°，每早一小时差-15°，每晚一小时差15°；④δ0为太阳赤纬角，计算方法为：其中，D为某天在儒略历中该年中的排序，如2016年1月1日，D为1。太阳方位角的计算方法为：所述的一种基于双视角光学遥感图像的海面风速探测方法，其特征在于步骤(3)中根据遥感图像文件所记录的遥感图像获取时的轨道倾角(S)和侧视角(P)，分别逐像元地计算前视遥感图像的传感器高度角(θN)和方位角后视遥感图像的传感器高度角(θB)和方位角计算公式如下：θN＝|n*IFOV+P|其中，n为距离中心像元的像元数，在中心点左边为正值，右边为负值；IFOV为遥感器的视场角。其中，h为遥感器的高度；m为遥感图像的空间分辨率；前视方向和后视方向的夹角；其中，G为前视和后视影像的基高比。所述的一种基于双视角光学遥感图像的海面风速探测方法，其特征在于步骤(4)中对前视遥感图像和后视遥感图像分别进行大气校正，获得前视太阳耀光辐射强度(LgN)和后视太阳耀光辐射强度(LgB)。大气校正的过程是根据步骤(2)确定的太阳高度角(θ0)和方位角步骤(3)确定的前视遥感图像的传感器高度角(θN)和方位角后视遥感图像的传感器高度角(θB)和方位角分别计算前视遥感图像的瑞利散射强度(LrN)和气溶胶散射强度(LaN)、后视遥感图像的瑞利散射强度(LrB)和气溶胶散射强度(LaB)，计算方式采用单次散射函数进行计算：公式中，各参数确定如下：①下标x为r和a，分别代表瑞利散射和气溶胶散射；②下标y为N和B分别代表前视遥感图像和后视遥感图像；③ρ为大气分子单次散射率；④τ为大气光学厚度；⑤F0为大气层外太阳垂直辐照度；⑥α±为散射角，由下式计算：⑦P(α±)为散射相函数；⑧R为菲尼尔反射率。所述的一种基于双视角光学遥感图像的海面风速探测方法，其特征在于步骤(4)中，获得LrN、LaN、LrB)和LaB后，采用以下公式计算前视太阳耀光辐射强度(LgN)和后视太阳耀光辐射强度(LgB)：Lgy＝Lty-Lry-Lay其中下标y为N和B分别代表前视遥感图像和后视遥感图像，Lty为遥感图像记录的辐射强度。所述的一种基于双视角光学遥感图像的海面风速探测方法，其特征在于步骤(5)中利用太阳耀光辐射传输模型和海面波面倾斜概率模型，分别建立前视太阳耀光辐射强度(LgN)和后视太阳耀光辐射强度(LgB)与海面粗糙度的关系：公式中，各参数确定如下：①F0为大气层外太阳垂直辐照度；②T0为太阳光大气直射透过率；③ωN和ωB分别是前视遥感图像和后视遥感图像的反射角，计算方法为：其中，下标y为前视(N)和后视(B)。④R为菲尼尔反射率，计算方法如下：其中ω’为海水的折射角，计算方法为：⑤βN和βB分别为前视和后视的反射面倾角，计算方法为：其中，下标y为前视(N)和后视(B)。所述的一种基于双视角光学遥感图像的海面风速探测方法，其特征在于步骤(6)中将步骤(5)获得的前视太阳耀光辐射强度(LgN)和后视太阳耀光辐射强度(LgB)进行比值计算：通过变形可以获得海面粗糙度所述的一种基于双视角光学遥感图像的海面风速探测方法，其特征在于步骤(7)中利用步骤(6)获得的海面粗糙度信息，逐像元地计算获得10米海面风速(U10)。计算方法海面粗糙度和海面风速之间的定量关系：利用该公式逐像元计算10米海面风速，获得海面风速探测结果。本专利技术的有益效果是：本专利技术的显著特点是利用同步或准同步观测的两个视角光学遥感图像，通过遥感成像几何关系，计算了前视遥感图像的高度角和方位角、后视遥感图像的高度角和方位角。然后经过对前视遥感图像和后视遥感图像进行大气校正获得前视太阳耀光辐射强度和后视太阳耀光辐射强度。基于同步和准同步观测的海面粗本文档来自技高网...

【技术保护点】
一种基于双视角光学遥感图像的海面风速探测方法，其特征在于其包括以下步骤：(1)选取由同一遥感平台，在时间间隔小于5分钟内成像的双视角光学遥感图像，根据观测方向与遥感平台运动方向的关系，分别记为前视遥感图像N和后视遥感图像B；(2)根据遥感图像的成像时间确定太阳高度角θ0和方位角(3)根据遥感图像文件所记录的遥感图像获取时的轨道倾角S和侧视角P，分别逐像元地计算前视遥感图像的传感器高度角θN和方位角后视遥感图像的传感器高度角θB和方位角(4)对前视遥感图像和后视遥感图像分别进行大气校正，获得前视太阳耀光辐射强度LgN和后视太阳耀光辐射强度LgB；(5)利用太阳耀光辐射传输模型和海面波面倾斜概率模型，分别建立前视太阳耀光辐射强度LgN和后视太阳耀光辐射强度LgB与海面粗糙度的关系；(6)采用比值法，逐像元地计算获得海面粗糙度；(7)利用海面粗糙度和海面风速的定量关系，计算海面风速，获得整幅图像的风速信息。

【技术特征摘要】
1.一种基于双视角光学遥感图像的海面风速探测方法，其特征在于其包括以下步骤：(1)选取由同一遥感平台，在时间间隔小于5分钟内成像的双视角光学遥感图像，根据观测方向与遥感平台运动方向的关系，分别记为前视遥感图像N和后视遥感图像B；(2)根据遥感图像的成像时间确定太阳高度角θ0和方位角(3)根据遥感图像文件所记录的遥感图像获取时的轨道倾角S和侧视角P，分别逐像元地计算前视遥感图像的传感器高度角θN和方位角后视遥感图像的传感器高度角θB和方位角(4)对前视遥感图像和后视遥感图像分别进行大气校正，获得前视太阳耀光辐射强度LgN和后视太阳耀光辐射强度LgB；(5)利用太阳耀光辐射传输模型和海面波面倾斜概率模型，分别建立前视太阳耀光辐射强度LgN和后视太阳耀光辐射强度LgB与海面粗糙度的关系；(6)采用比值法，逐像元地计算获得海面粗糙度；(7)利用海面粗糙度和海面风速的定量关系，计算海面风速，获得整幅图像的风速信息。2.如权利要求1所述的一种基于双视角光学遥感图像的海面风速探测方法，其特征在于步骤(2)中，根据获取的遥感图像的经纬度位置和成像时间，逐像元计算太阳天顶角θ0和方位角太阳天顶角θ0的计算方法为：θ0＝90°-h0sinh0＝sinγsinδ0+cosγcosδ0cost0公式中，各参数确定如下：①h0为太阳高度角；②γ为当地纬度，即遥感图像各像元的纬度；③t0为当地太阳时角，由遥感图像时的地方时确定，中午12时为0°，每早一小时差-15°，每晚一小时差15°；④δ0为太阳赤纬角，计算方法为：其中，D为某天在儒略历中该年中的排序；太阳方位角的计算方法为：3.如权利要求1所述的一种基于双视角光学遥感图像的海面风速探测方法，其特征在于步骤(3)中根据遥感图像文件所记录的遥感图像获取时的轨道倾角S和侧视角P，分别逐像元地计算前视遥感图像的传感器高度角θN和方位角后视遥感图像的传感器高度角θB和方位角计算公式如下：θN＝|n*IFOV+P|其中，n为距离中心像元的像元数，在中心点左边为正值，右边为负值；IFOV为遥感器的视场角；其中，h为遥感器的高度；m为遥感图像的空间分辨率；前视方向和后视方向的夹角；其中，G为前视和后视影像的基高比。4.如权利要求1所述的一种基于双视角光学遥感图像的海面风速探测方法，其特征在于步骤(4)中对前视遥感图像和后视遥感图像分别进行大气校正，获得前视太阳耀光辐射强度LgN和后视太阳耀光辐射强度LgB；大气校正的过程是根据步骤(2)确定的太阳高度角θ0和方位角步骤(3)确定的前视遥感图像的传感器高度角θN和方位角后视遥感图像的传感器高度角θB和方位角分别计算前视遥感图像的瑞利散射强度LrN和气溶胶散射强度LaN、后视遥感图像的瑞利散射强度LrB和气溶胶散射强度LaB，计算方式采用单次散射函数进行计算：Lxy=ωxτxF0[Px(αy-)+(R(θ0)+R(θy)Px(αy+))]4πcosθy]]>公式中，各参数确定如下：①下标x为r和a...

【专利技术属性】
技术研发人员：张华国，杨康，楼琇林，王隽，郑罡，任林，史爱琴，厉冬玲，
申请(专利权)人：国家海洋局第二海洋研究所，
类型：发明
国别省市：浙江;33