北斗卫星精密星历插值精度研究

选用非滑动式拉格朗日(Lagrange)插值法和滑动式Lagrange插值法对北斗卫星导航系统(Beidou Satellite Navigation Syetem,BDS)精密星历进行插值,通过大量实验对北斗卫星导航系统卫星的插值精度与两种方法及插值阶次的关系进行研究.实验结果表明:(1)分别采用非滑动式和滑动式La...     

海底地形数据处理与DEM生成

分析了海底DEM生成经常会遇见的问题,介绍了坐标变换、接边处理、插值加密等数据预处理方法,最终生成海底规则格网DEM。     

普通Kriging插值算法研究

地统计学是进行空间变异分析的有力工具.文章说明了地统计学的基本原理及其实现算法,并应用该算法对研究区域的草句白浆土有效磷的空间变异特征进行了分析,通过计算得出的地统计学参数说明了土壤有效磷具有中等空间相关性,并生成了该养分的空间分布图.     

VIIRS/DNB夜间灯光月度产品插补方法对比——以北京为例

新一代VIIRS/DNB(Visible Infrared Imaging Radiometer Suite Day/Night Band)夜间灯光数据因其具备更精细的时空分辨率、数据不存在饱和现象、不同年份数据可比性强等优点而迅速代替DMSP/OLS(Defense Meteorological Satellite Program's Operational Linescan System)夜间灯光数据成为新的研究热点。但由于杂散光的污染,VIIRS/DNB夏季数据数值缺失严重,致使数据在空间和时间上不连续,因此,插补缺失数据成为后续应用的前提。鉴于此,论文从插补结果异常值、与参照值对比、计算性能三个方面,系统比较了三次样条插值(样条)、三次Hermite插值(Hermite)、灰色预测模型(GM)、三次指数平滑法(指数)4种插补方法的适用性,以期为插补VIIRS/DNB夜间灯光数据提供方法选择的依据。研究结果表明：①异常值比较方面,Hermite法未出现异常值,另外3种算法仅出现少量异常值(0.02%~1.34%);②与参照值的对比方面,Hermite法与参考值接近程度最高,GM接近程度最低,样条法和指数法介于两者中间;③算法性能比较方面,4种方法都具备计算简单、容易编程的特点,但指数法的算法时长是另外3种方法的10倍以上。因此,综合多方表现,当插补月份前后两侧均有足够长的原始数据时,插补效果好、计算速度快、不会出现过冲现象的Hermite法最适宜,样条法次之;当插补月份仅单侧有足够长的数据时,适宜采用指数法(插补效果好、计算速度较慢)或GM(插补效果偏低、计算速度快)进行插补。     

一种基于牛顿插值迭代的SAR影像滤波评价方法

针对SAR影像难以精确计算信噪比的问题,该文以滤波评价指标信噪比为依托,提出了一种基于牛顿插值迭代修正SAR影像近似信噪比的滤波质量评价方法:根据已知值通过牛顿插值得到更多的数据,再利用牛顿迭代法拟合非线性函数修正近似信噪比。以模拟影像和真实影像为数据源,选择目视效果较好的滤波结果作为真实影像的近似无噪声影像,人工对模拟影像和无噪影像进行降质,并使用5种滤波器对影像进行处理,分别利用基于牛顿插值迭代得到的修正近似信噪比、传统的峰值信噪比、等效视数和边缘保持指数评价滤波效果。结果表明:基于牛顿插值迭代修正的近似信噪比可作为衡量SAR影像滤波质量的评价指标,且无需参考原始无噪影像,对真实SAR影像的滤波质量评价有一定的可行性。     

基于梯度依赖客观分析技术的全球Argo网格化数据集:构建及初步应用

采用梯度依赖相关尺度方法构建了1套2004—2017年间,月平均的全球海洋(0~1 500 m)1°×1°的Argo数据集,并在对该数据集进行对比检验的基础上,将其初步应用于中西太平洋黄鳍金枪鱼的渔场分析研究。结果表明,所构建的Argo数据集与WOA13数据集的温、盐偏差在上表层(150 m)稍大,最大值分别约为0.5 ℃和0.1,且偏差均随深度的增加而逐渐减小;其与TAO浮标时间序列的温度偏差,2004—2017年间均小于1 ℃,最大盐度偏差则小于0.5,且大部分海域接近0。中西太平洋海域,黄鳍金枪鱼中心渔场多集中在 28~29 ℃ 等温线范围内,在 22 ℃以下的海域单位捕捞努力量渔获量(catch per unit effort,CPUE)值极小;中心渔场区温跃层上界深度范围在20~120 m之间,且中心渔场在各个深度上形成的频数大体呈正态分布,温跃层上界深度为90 m时,形成中心渔场的可能性达到最大。研究表明所构建的数据集在水文环境分析及资源评估中有一定的应用价值。     

基于含砾Folk三角形的底质范围分析

为解决面向航海应用的底质范围分布问题,提出了一种基于含砾Folk三角形的底质范围生成方法。该方法通过将海图底质注记转化为含砾Folk法底质类型,进而逆向利用Folk(G)三角形,将其拆解为各组分的比例,最后对海区内各组分单独插值并通过地图代数的方法整合各插值结果,生成底质的范围分布。对比实验表明其生成结果较传统的Voronio图法可更好地应用于舰船抛锚等应用场景。     

基于EIGEN-6C4的区域重力异常插值方法

重力异常作为地球重力场的一种基本指标,广泛应用于大地测量、地球物理、地质、地震与海洋领域。为了快速获取高精度的未知点重力异常,本文基于EIGEN-6C4重力场模型解算某110 km×110 km区域内重力异常,构建5′×5′的重力异常格网模型,结合格网节点数据采用6种不同的插值算法拟合未知点的重力异常,最后使用外部检核的方法评价插值算法的精度。实验结果表明,6种插值算法的拟合精度均达到4 mgal以内,自然三次样条法和反距离加权插值法拟合精度最高,得到精度达到2.7 magl的重力异常。     

GNSS精密轨道产品不连续性分析及评估

按照目前的国际规范, 高精度GNSS (Global Navigation Satellite System)轨道产品一般以天为周期进行发布, 提供给用户使用. 连续使用多天的产品存在不同天间的跳变问题. 利用德国地学研究中心(GFZ)、欧洲定轨中心(COD)、欧空局(ESA)、美国喷气试验室(JPL)以及上海天文台(SHA)共5个GNSS分析中心2013---2017年的轨道产品, 分析了轨道跳变的特性. 计算结果表明: GFZ、COD、ESA、SHA和JPL的3维轨道跳变平均分别为7.79cm、1.51cm、7.77cm、11.75cm和2.51cm. 轨道跳变序列的周期特性分析表明: 序列存在90d、120d、340d左右的显著周期项, 对应于海潮对地球自转的影响, 其振幅为数毫米至1cm左右. 表明精密轨道确定需要进一步精化该项模型; GPS的跳变序列还存在与卫星星座相关的175d和352d左右的交点年显著周期项. 此外, 针对COD产品外推轨道的分析, 验证了地球反照辐射压和太阳光压模型等动力学模型对轨道的差异.