普通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(插补效果偏低、计算速度快)进行插补。     

A comprehensive analysis of phenological changes in forest vegetation of the Funiu Mountains,China

Journal of Geographical Sciences - This paper reports the phenological response of forest vegetation to climate change (changes in temperature and precipitation) based on Moderate Resolution...     

北斗三号组网卫星数据质量分析及单系统定轨精度初步评估

选取2018-01-23起10 d内16个iGMAS测站观测数据,对北斗三号卫星的观测数据质量及BDS单系统精密定轨精度进行评估。初步结果表明,老信号B1I、B3I北斗三号卫星的信噪比略强于二号卫星,噪声与多路径基本相当,均在0.1 m量级,新卫星不存在星内多路径偏差。新信号B1C/L1/E1频点GPS信噪比最强,Galileo和BDS卫星相当,B2a/L5/E5a和B2b/E5b各系统基本相当;噪声及多路径方面,B1C/L1/E1频点GPS优于BDS、Galileo卫星0.1 m量级,B2a/L5/E5a和B2b/E5b各系统基本相当,均在0.1 m量级,新信号中北斗三号卫星星内多路径偏差基本消失。单系统精密定轨试验中,分别进行有/无GEO卫星策略、太阳光压模型ECOM 五/九参数策略的比较,并使用卫星激光测距数据进行独立检核。初步结果表明,有GEO卫星、ECOM五参数光压模型的定轨精度最好,C19号卫星7个重叠弧段的平均定轨精度在沿迹向、法向、径向的精度分别为32 cm、16 cm、8 cm,与试验卫星的定轨精度基本相当。     

中国耕地土壤相对湿度时空分异

以全国653个农业气象站1993-2013年耕地的土壤相对湿度数据为基础,运用地统计方法,分析中国耕地土壤相对湿度时空分异特征与规律。结果表明:自1993年以来全国耕地的土壤相对湿度呈现波动上升趋势。全国耕地的土壤相对湿度普遍大于60%,分布区域自4月中旬开始随夏季风推移不断向北向西扩大,自10月下旬开始向东、南方向缩小。耕地土壤相对湿度值随土壤深度的增加而增加。年际尺度上,耕地的土壤相对湿度在夏秋季上升速度最快,变化幅度随土层深度增加而变小。土壤相对湿度与降水量存在较强的正相关关系,与潜在蒸发量、气温普遍存在较强的负相关关系。土壤相对湿度与各气象要素的相关性随深度加深而减弱。春、夏、秋季气象因素对旱地土壤相对湿度影响较大,冬季气象因素对水田影响更大。     

模糊聚类定权法对SLR定轨精度的影响

针对卫星激光测距(satellite laser ranging,SLR)精密定轨过程中存在的测站观测数据合理定权问题,将一种改进的模糊聚类算法引入到SLR观测数据定权中。基于国际激光测距服务(International Laser Ranging Service,ILRS)提供的全球SLR测站性能报告,对测站进行近实时滑动分类定权,改变SLR数据处理中权重的经验或者随意性选取模式。经过LAGEOS1卫星2014年1月至2016年12月3年全球SLR实测数据处理的测试。结果表明,当考虑LAGEOS标准点总数、LAGEOS标准点RMS值以及LAGEOS标准点合格率这3项测站质量控制因素确定的测站权值能最大限度地提高卫星定轨精度和观测数据的使用效率,对参与计算的365个3d弧段数据,91.46%弧段精度得到提高,平均提高约3.7mm,且每个测站的定轨残差RMS也得到了降低。这对于正在迈向毫米级测量精度的SLR技术至关重要。     

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产品外推轨道的分析, 验证了地球反照辐射压和太阳光压模型等动力学模型对轨道的差异.     

一种基于Fréchet距离的断裂等高线内插算法

等高线内插是提高地形图精度的一个常用手段,在地图综合等领域中有很重要的作用。在实际地形图数据库中,等高线除了存在完整的计曲线以外,还有一部分断裂的首曲线。已有的等高线内插研究中未考虑对此类数据的处理,为此,本文提出了一种基于Fréchet距离的断裂等高线内插算法。首先提出了等高线度量关系的计算方法,然后对等高线的节点均匀加密,并利用Fréchet距离进行相似度判断以选择参考等高线,最后根据“最近点”的方式进行插值。通过对江苏某地区实际数据的试验,验证了该算法的合理性,对于鞍部地区的内插,有更好的适应性和准确性。     

地月系L_2点中继星激光测距成功率的估算

为了模拟位于地月系L2点的中继星"鹊桥"与月球的位置关系,进而估算中继星激光测距的成功率,按照轨道周期约为14天的要求对中继星所在的晕轨道进行计算,建立了一个综合考虑望远镜抖动、大气抖动和预报轨道横向偏离的模型。从数值上给出了一条轨道周期为14. 78天,X方向(地月连线方向)振幅为12 493 km,Y方向为34 596 km,Z方向(垂直于地月轨道平面方向)为11 916 km的周期轨道。由于晕轨道的最小振幅远大于月球遮挡的临界振幅4 000 km,因此月球对中继星不存在遮挡问题。基于建立的测距成功率模型,根据昆明站(国际编号:7820)的激光测距系统对运行在该轨道上的中继星进行测距成功率分析,结果表明:测距成功率随着中继星横向轨道标准差的增大呈快速降低的趋势。对于中继星到测站的平均距离而言,当中继星没有横向偏离时,探测器产生的光电子数为0. 151,成功率为14. 07%;横向偏离2 km时,光电子数降为0. 035,成功率降为3. 46%。对比最近距离与最远距离的情况,无横向偏离的情况下,探测器产生的光电子数从0. 174降为0. 139,成功率从16. 01%降为13. 02%。该计算结果可为云南天文台1. 2 m望远镜实现中继星激光测距提供参考。     

非滑动式与滑动式拉格朗日插值法在BDS精密星历内插中的比较分析

在高精度卫星定位数据处理中,对精密星历内插处理是一项基础工作。利用拉格朗日多项式插值法,从三方面对北斗导航系统精密星历进行了内插精度分析,分别是非滑动式拉格朗日插值法在不同插值阶数下的插值精度,滑动式拉格朗日插值法在不同插值阶数下的插值精度,以及非滑动式和滑动式拉格朗日插值法在同一插值阶数下的插值精度。算例结果表明,如果将插值阶数作为横坐标,且向右为正,插值精度为纵坐标,向上为正,那么非滑动式拉格朗日插值法的插值精度与插值阶数呈现一个开口向下的抛物线关系;滑动式拉格朗日插值法的插值阶数高于7后,插值精度处于稳定且高的状态,插值误差最大值为1. 1 mm左右,均方差在0. 4 mm左右;在同一插值阶数下,滑动式拉格朗日插值法的精度高于非滑动式拉格朗日插值法的精度。