本发明专利技术公开了一种融合近海多颗卫星海表高度异常数据的方法，其特征在于，包括以下步骤：1）获取多颗近海卫星沿轨海表高度异常观测值数据；2）利用经验模分解方法进行卫星沿轨海表高度异常观测值的空间虑波；3）对空间虑波后的数据，利用最优插值方法将不规则空间、时间采样点上的海表高度异常观测值插值到预设空间分辨率、时间分辨率的规则时间、空间网格上。本方法产生的数据更接近沿轨资料及测潮站资料，降低了误差。
【技术实现步骤摘要】
本专利技术公开，属于海洋卫星资料 处理
技术介绍
海表高度异常（Sea Level Anomaly)是表征海洋状况的一个重要指标。从 1992 年 Topex/Poseidon 卫星发射成功后（Fu et al.，1994, TOPEX/POSEIDON mission overview, J. Geophys. Res.，99, 24369-24381)，卫星海表高度计的观测技术越来越成熟，目 前可以同时得到多颗卫星海表高度异常的观测。图1显示了 2012年1月1日至1月11日 11天内三颗卫星观测的中国近海（10° S-50° N，90° E-160。E)海表高度异常资料的空 间分布，这三颗卫星是Jason-1、Jason-2和CryoSat-2。从图1可以注意到卫星海表高度 异常观测在空间上采样不均匀，沿卫星轨道（沿轨）通常可以达到6km，但轨道间的距离为 200-300km。卫星海表高度资料在时间采样上也有其自身特点，比如Jason-1及Jason-2卫 星轨道重复采样的间隔为近10天，CryoSat-2时间采样上约为35天重复。 在科学研宄及社会应用中，时间及空间上采样均匀的海表高度异常资料 会对应用有极大的益处并且方便使用。DUACS (Developing Use of Altimetry for Climate Studies)及 AVISO (Archiving, Validation and Interpretation of Satellite Oceanographic data)是法国的机构，专门从事海表高度异常资料在全球 范围内的处理，产生规则网格点（0.25° X0.25° )、时间分辨率为1天的全球海表 高度异常资料产品 ?ibarboure et al. ，2011，Jason_2in DUACS:Updated system description, first tandem results and impacts on processing and products, Marine Geodesy, 34, 214-241 ；Le Traon et al.,1998,An improved mapping method of multisatellite altimeter data, J. of Atmos, and Oceanic Technology, 15, 522-534)〇 但AVISO资料在产生全球产品时，只用了部分（1/3)的沿轨资料，并且在近海因为使用的 Lanzcos空间滤波方法会造成大的误差。这样便降低了海表高度资料在中国近海的精度，给 资料应用带来一些不便。
技术实现思路
本专利技术针对AVISO全球海表高度异常资料处理中近海资料精度降低的问题，提供 一种新颖、有效、实用的融合中国近海多颗卫星海表高度异常数据的方法。 为解决上述技术问题，本专利技术提供一种融合近海多颗卫星海表高度异常数据的方 法，其特征在于，包括以下步骤： 1)获取多颗近海卫星沿轨海表高度异常观测值数据； 2)利用经验模分解（Empirical Mode Decomposition)方法进行沿轨海表高度异 常观测值的空间虑波； 3)对空间虑波后的数据，利用最优插值方法将不规则空间、时间采样点上海表高 度异常观测值插值到预设空间分辨率、时间分辨率的规则时间、空间网格上。 本专利技术所达到的有益效果：本专利技术的方法针对AVISO的方法进行了改进，提出了 一种融合中国近海多颗卫星海表高度异常数据的方法，利用三颗卫星的沿轨海表高度异常 观测产生规则空间、时间网格上中国近海海表高度异常产品，与沿轨资料的均方根误差为 3. 49cm，而AVISO产品在中国近海海域与沿轨资料的均方根误差为3. 89cm，。因而本方法产 生的近海海表高度异常数据更接近沿轨资料，将误差降低了 10. 3%。与49个独立的测潮站 相比，本方法产生的近海海表高度异常数据与测潮站资料的均方根误差为10.01cm，AVISO 产品的均方根误差为11. 〇〇cm，因而与测潮站资料相比，本方法的误差降低了 9%。【附图说明】 图1为2012年1月1日至1月11日三颗卫星在中国近海海表高度异常的观测采 样点示意图； 图2为经验模分解方法的示意图； 图3a为模拟的海表高度异常真实值加上一个正态分布随机误差后的示意图； 图3b为用Lanczos滤波方法及经验模分解滤波方法进行空间滤波后的结果示意 图； 图4为2012年AVISO海表高度异常资料产品与沿轨资料的均方根误差及本方法 产生的规则网格点上的海表高度异常资料产品与沿轨资料的均方根误差对比图； 图5为评价海表高度异常资料产品所用的49个测潮站的地理位置示意图； 图6为比较AVISO及本方法生成的海表高度异常资料产品与49个测潮站的均方 根误差示意图。【具体实施方式】 下面结合附图对本专利技术做进一步的说明。 本专利技术的，包括以下步骤： 1.获取三颗近海卫星沿轨海表高度异常观测值数据。 2.利用经验模分解方法进行沿轨资料的空间虑波。 在使用沿轨海表高度异常资料进行空间、时间最优插值前，要对沿轨资料进行空 间上的滤波以去掉一部分观测误差及小尺度噪音（在中炜度为60公里以下）。AVISO中使 用了传统的Lanzcos滤波方法，在中炜度滤掉空间尺度小于40-60km的扰动。这种方法将 沿轨资料进行FFT变换，在谱空间进行滤波，然后将滤波后的谱进行叠加以产生滤波后的 资料。因为这种方法的基函数是正弦及余弦函数，对一般性的信号并不能有效地表达，在沿 轨资料的两端易出现大的误差，这是AVISO资料处理中已被注意到的问题。 为解决这一技术问题，本方法采用了经验模分解方法（Empirical Mode Decomposition,也称为 Hilbert-Huang Transform, Huang and ffu,2008, A review on Hilbert-Huang Transform:Method and its applications to geophysical studies, Rev. Geophys.，46, RG2006, doi : 10. 1029/2007RG000228)对沿轨资料进行空间滤 波。 本方法将一段连续的沿轨海表高度异常观测值依照数据特点分解为多个经验模 态，然后去掉空间尺度最小的第一模态（在经验模分解方法中，第一模态就是空间尺度最 小的模态），将其它空间尺度大的模态叠加起来求和得到滤波后的沿轨海表高度异常观测 值数据。 求取经验模的方法为筛选法，以下仅描述与沿轨海表高度异常观测处理有关的部 分，细节及数学原理可参考Huang和Wu (2008)。图2为经验模分解方法的示意图，x-轴为 空间距离，y_轴为信号的大小。a)显示任意一个随空间变化的信号；b)找出其极大值（菱 型）及极小值（圆圈）；c)用三次样条插值将极大值及极小值连起来，得到极大值包络线及 极小值包络线，求出极大值包络线及极小值包络线的平均，图中以虚线表示；d)求出原始 信号与极大值和极小值包络线平均的差；e)对d)中的差重复步骤b)至d) ;f)当以上过程 结束本文档来自技高网...
【技术保护点】
一种融合近海多颗卫星海表高度异常数据的方法，其特征在于，包括以下步骤：1)获取多颗近海卫星沿轨海表高度异常观测值数据；2)利用经验模分解方法进行卫星沿轨海表高度异常观测值的空间滤波；3)对空间虑波后的数据，利用最优插值方法将不规则空间、时间采样点上的海表高度异常观测值插值到预设空间分辨率、时间分辨率的规则时间、空间网格上。
【技术特征摘要】
【专利技术属性】
技术研发人员：王晓春，
申请(专利权)人：南京信息工程大学，
类型：发明
国别省市：江苏;32
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。