CN119475197A 一种基于哨兵3号卫星高度计数据的内孤立波检测方法 （自然资源部第一海洋研究所）

一种基于哨兵3号卫星高度计数据的内孤立开了一种基于哨兵3号卫星高度计数据的内孤立次提取Ku波段和C波段各个坐标点的后项散射系段和C波段的后项散射系数的差分均方斜率；对ASSHA和ASWH；分别计算差分均方斜率以及21.一种基于哨兵3号卫星高度计数据的内孤立波检测方法，其特征在于，包括如下步步骤3，依次计算各个坐标点对应的Ku波段和C波段的后项散射系数的差分均方斜率步骤4，分别计算各个坐标点的差分均方斜2.根据权利要求1所述的一种基于哨兵3号项散射系数0"和of的取值范围设定为8~18，将响应特征参数SSHA的数据范围设定为ISSHAI<1，并且剔除水蒸气含量大于60Kg/m2和液态水浓度大于0.1Kg/m2的异常后项散射系数和SSHA值。3.根据权利要求1所述的一种基于哨兵3征在于，步骤3中，计算坐标点对应的Ku波段和C波段的后项散射系数的差分均方斜率计算如下：;"为Ku波段坐标点i的有效反射率系数，o"为Ku波段坐标点i的后项散射系数，pf为C波段坐4.根据权利要求1所述的一种基于哨兵3号卫星高度计5.根据权利要求1所述的一种基于哨兵3;3i代表第i个坐标点，x[i]代表窗口内第i个坐标点的差分均方斜率;。6.根据权利要求5所述的一种基于哨兵3号卫星;;4[0001]本发明涉及海洋内孤立波自动检测领域，特别涉及一种基于哨兵3号卫星高度计[0002]海洋内孤立波传播过程中携带巨大能量并会诱发突发性的强流和强烈的等密度[0004]早期的高度计为脉冲限制型传统高度计，沿轨分辨率通常为2_3km，采样频率为[0005]为解决上述技术问题，本发明提供了一种基于哨兵3号卫星高度计数据的内孤立率；对响应特征参数SSHA和SWH进行两次平滑处理，得到两个新的变量ΔSSHA和;5步骤4，分别计算各个坐标点的差分均方斜率以及两个新的变量ΔSSHA围设定为ISSHA]<1，并且剔除水蒸气含量大于60Kg/m2和液态水浓度大于0.1Kg/m2的;为Ku波段坐标点i的有效反射率系数，o"为Ku波段坐标点i的后项散射系数，pf为C二次平滑处理将平滑长度设置为30km；随后计算第一次平滑结果与第二次平滑结果的差，得到两个新的变量ΔSSHA和ASWH，ΔSSHA为平滑后的海表面高度异常变化量，ASWH为平滑后的有效波高变化量。[0010]上述方案中，步骤4中，计算坐标点的差分均方斜率以及两个新的变量ΔSSHA和ASWH的局部标准差SD的公式如下：;i代表第i个坐标点，x[il代表窗口内第i个坐标点的差分均方;通过计算得到差分均方斜率、变量ΔSSHA和dswH的局部平均值6标准差、SDsssna和SDaswn。;;根据系数K，分别计算局部标准差、SDsssna和SDaswn对应的阈值[0012]通过上述技术方案，本发明提供的一种基于哨兵3号卫星高度计数据的内孤立波2、传统方法通常仅仅依赖于均方斜率或波形峰值功率来实现内孤立波的自动检算第一次平滑结果与第二次平滑结果的差，得到两个新的变量么SSHA和ASWH。该变量消除了中大尺度海洋现象对数据的影响，过滤了数据中不属于内孤立波的高频扰动，7[0014]图1为本发明实施例所公开的一种基于哨兵3号卫星高度计数据的内孤立波检测图2为2020年4月6日苏禄海海域的Sen图3为与图2的OLCI光学遥感图像匹配的苏禄海海域的Sentilne_3SRAL图4为苏禄海2020年4月6日Sentilne_3SRAL数据的内孤立波检测结果a）为差分均方斜率变量的局部标准差变化曲线b）为海表面高度异常变量图7为南海2023年04月25日Sentilne_3[0016]本发明提供了一种基于哨兵3号卫星高度计数据的内孤立波检测方法，如图1所[0017]实施例中使用的卫星高度计为2017年至2023年9个相对轨道共581个pass的Sentinel_3A数据和2019年至2023年7个相对轨道共376个pass的Sentinel_3B数据构建数[0018]获取Sentinel_3卫星搭载的光学海陆色度仪（OceanandLandColor步观测的特性，从数据集中提取出129个经过了OLCI光学遥感图像验证的含有内孤立波的8和SSHA范围取为：;。[0023]排除陆地和特殊现象对数据的影响后，降雨也是检测方;;除水蒸气含量大于60Kg/m2和液态水浓度大于0.1Kg/m2的异常后项散射系数和ASWH。[0026]Sentinel_3卫星的SRAL是双频雷达高度计，因此数据产品中提供Ku波段和C波段[0027]计算坐标点对应的Ku波段和C波段的后项散射系数的差分均方斜率计算9;偏置参数来补偿TOPEX与Sentinel_3A/B雷达高度计之间C波段的00差异。本实施例中，。此第一次平滑处理将平滑长度设置为1.8km以过滤数据中的高频扰动；中大尺度海洋现象长度设置为30km；随后计算第一次平滑结果与第二次平滑结果的差，得到两个新的变量ΔSSHA和ASWH，ΔSSHA为平滑后的海表面高度异常变[0031]步骤4，分别计算各个坐标点的HA和aswn的局部标准差，首先需要确;其中，x[il代表窗口内第i个坐标点的差分均方斜率、变量ΔSSHA或;;i代表第i个坐标点。[0033]通过计算得到差分均方斜率、变量ΔSSHA和swn的局部平均值标准差、SDsssna和SDaswn。;;[0035]根据系数K，分别计算局部标准差、SDsssna和SDaswn对应的阈值所示。[0038]对所公开的实施例的上述说明，使本领域专业技术人员能够实现或使用本发