Propagation of vertical and horizontal source data errors into a TIN with linear interpolation

Digital elevation models (DEMs) have been widely used for a range of applications and form the basis of many GIS-related tasks. An essential aspect of a DEM is its accuracy, which depends on a variety of factors, such as source data quality, interpolation methods, data sampling density and the surface topographical characteristics. In recent years, point measurements acquired directly from land surveying such as differential global positioning system and light detection and ranging have become increasingly popular. These topographical data points can be used as the source data for the creation of DEMs at a local or regional scale. The errors in point measurements can be estimated in some cases. The focus of this article is on how the errors in the source data propagate into DEMs. The interpolation method considered is a triangulated irregular network (TIN) with linear interpolation. Both horizontal and vertical errors in source data points are considered in this study. An analytical method is derived for the error propagation into any particular point of interest within a TIN model. The solution is validated using Monte Carlo simulations and survey data obtained from a terrestrial laser scanner.     

斜压罗斯贝变形半径优化的误差相关尺度及其对最优插值效果的改进

利用最优插值方法进行三维温盐场重构时,如何正确估算背景场误差协方差至关重要。针对背景场误差协方差主要取决于误差相关尺度的问题,本文提出了用中国东部海域及其附近海域的最新的高分辨率气候态数据进行斜压罗斯贝变形半径优化其误差相关尺度计算的研究方案和技术途径;对比分析了均一化相关尺度方案和法国ISAS系统尺度方案,讨论了变形半径对最优插值的影响。结果表明:均一化相关尺度方案的均方根误差小于ISAS方案,但温度场过于平滑,难以刻画一些重要的物理现象;相比而言,本文提出的基于变形半径的相关尺度方案在取2倍变形半径时,不仅均方根误差在各水平层较小,且温度场能够更好地刻画四国海盆海域涡旋及黑潮影响的温度场三维结构。由于实际海洋中各层物理过程的尺度存在差异,实际应用时各层的最优尺度设置也应有所不同。     

基于ROMS模式利用集合最优插值同化沿轨海表面高度异常数据的研究

The ensemble optimal interpolation （EnOI） is applied to the regional ocean modeling system （ROMS） with the ability to assimilate the along-track sea level anomaly （TSLA）. This system is tested with an eddy-resolving system of the South China Sea （SCS）. Background errors are derived from a running seasonal ensemble to account for the seasonal variability within the SCS. A fifth-order localization function with a 250 km localization radius is chosen to reduce the negative effects of sampling errors. The data assimilation system is tested from January 2004 to December 2006. The results show that the root mean square deviation （RMSD） of the sea level anomaly decreased from 10.57 to 6.70 cm, which represents a 36.6% reduction of error. The data assimilation reduces error for temperature within the upper 800 m and for salinity within the upper 200 m, although error degrades slightly at deeper depths. Surface currents are in better agreement with trajectories of surface drifters after data assimilation. The variance of sea level improves significantly in terms of both the amplitude and position of the strong and weak variance regions after assimilating TSLA. Results with AGE error （AGE） perform better than no AGE error （NoAGE） when considering the improvements of the temperature and the salinity. Furthermore, reasons for the extremely strong variability in the northern SCS in high resolution models are investigated. The results demonstrate that the strong variability of sea level in the high resolution model is caused by an extremely strong Kuroshio intrusion. Therefore, it is demonstrated that it is necessary to assimilate the TSLA in order to better simulate the SCS with high resolution models.     

Improvement of the semi-Lagrangian advection scheme in the GRAPES model: Theoretical analysis and idealized tests

ABSTRACT The Global/Regional Assimilation and PrEdiction System （GRAPES） is the newgeneration numerical weather predic- tion （NWP） system developed by the China Meteorological Administration. It is a fully compressible non-hydrostatical global/regional unified model that uses a traditional semi-Lagrangian advection scheme with cubic Lagrangian interpola tion （referred to as the SL_CL scheme）. The SL_CL scheme has been used in many operational NWP models, but there are still some deficiencies, such as the damping effects due to the interpolation and the relatively low accuracy. Based on Reich＇s semi-Lagrangian advection scheme （referred to as the R2007 scheme）, the Re_R2007 scheme that uses the low- and high-order B-spline function for interpolation at the departure point, is developed in this paper. One- and two-dimensional idealized tests in the rectangular coordinate system with uniform grid cells were conducted to compare the Re..R2007 scheme and the SL_CL scheme. The numerical results showed that： （1） the damping effects were remarkably reduced with the Re_R2007 scheme; and （2） the normalized errors of the Re_R2007 scheme were about 7.5 and 3 times smaller than those of the SL_CL scheme in one- and two-dimensional tests, respectively, indicating the higher accuracy of the Re..R2007 scheme. Furthermore, two solid-body rotation tests were conducted in the latitude-longitude spherical coordinate system with non uniform grid cells, which also verified the Re_R2007 scheme＇s advantages. Finally, in comparison with other global advection schemes, the Re_R2007 scheme was competitive in terms of accuracy and flow independence. An encouraging possibility for the application of the Re_R2007 scheme to the GRAPES model is provided.     

降雨资料Kriging与IDW插值对比分析—以漓江流域为例

降水空间化信息在很多领域都具有重要意义,而进行空间插值方法形成降水空间化信息是当代较为常用的方式。面对众多的插值方法其插值精度成为我们是否采用的关键,就kirging插值法与IDW(Inverse Distance Weighting)插值进行研究以探讨其插值效果,为此类科学研究提供依据。研究以漓江流域内各气象站点降水统计资料为基础分别采用上述两种方法进行插值处理,通过与预留实测站点比较评判方法的插值效果的优劣。结果表明,多年月平均降水量作为时间步长时,降水丰沛的月份使用kriging插值法较优于IDW插值法,而枯水月份则使用IDW插值法较优于kriging插值法。     

高压线对CSAMT探测结果的影响及数据处理

以干扰严重的某工作区为例,探讨高压线对CSAMT法所测数据的影响,总结了本工作区高压线对CSAMT数据影响的特征;根据其电磁干扰特征,寻求数据处理的方法,运用插值法恢复接近正常场的磁信号,再利用原始电场信号和恢复的磁信号,计算得出视电阻率。这种方法经对比和钻探验证,效果良好。     

土壤粒径空间预测方法研究综述

土壤粒径（砂粒、粉粒和黏粒）是各种陆表过程和生态系统服务评估等模型的关键参数。作为一种土壤成分数据,土壤粒径的空间预测方法有和为1（或100%）等特殊要求,其空间分布精度受预测方法影响较大。本文针对土壤粒径相较于其他土壤属性的特殊性,提出了土壤粒径空间预测方法框架,综述了土壤粒径数据变换、空间插值和精度验证等系列方法,总结了提升土壤粒径空间预测精度的各种途径,包括通过有效的数据变换改善数据分布、结合数据分布特点选择合适的预测方法、结合辅助变量提升制图精度和分布合理性、使用混合模型提升插值精度、使用多成分联合模拟模型提升预测的系统性等。最后,提出了今后土壤粒径空间预测方法研究的未来方向,包括从考虑数据变换原理和机制角度改善数据分布、发展多成分联合模拟模型和高精度曲面建模方法,以及引入土壤粒径函数曲线并与随机模拟结合等。     

遥感影像大气校正通用查找表的设计与插值算法

大气参数查找是遥感图像大气校正中一个重要的步骤。本文针对当前多数大气参数查找表普适性不强、查找表数据存储空间巨大、具体实施步骤不详细的问题,建立了一个通用化的多维大气参数查找表,其光谱分辨率为1 nm,并适用于多种遥感传感器。本文采用二进制文件方式存储查找表,并研究采用多维拉格朗日和多维反距离加权插值方法对输入参数进行插值。结果表明以二进制文件方式存储查找表,不但可以最大化地减少磁盘存储空间,还可以实现查找表数据快速随机读取;相对于多维反距离加权插值方法,多维拉格朗日插值方法的速度更快,精度更高,与MODTRAN4.0模型计算结果(真值)相比,它的误差仅为±1.0%左右。     

真彩色点云三维模型方法研究

通过DLT变换建立三维几何模型上的点与相应文物相片上各彩色像素之间的关系。首先,根据点的扫描坐标,计算该点在真彩色像片上的位置;其次,内插出红、绿、蓝三原色;最后,将三维几何模型与真彩色纹理进行最佳融合。建立文物真彩色真实感觉的数字模型,为数字文物和文物保护提供了理论和技术参考。     

边界虚拟钻孔在复杂地质体3维建模中的引入与确定

在矿山复杂地质体3维建模中,由于数据缺乏致使模型难以控制地质体空间边界特征,直接影响到模型的准确性和实用价值。借助边界虚拟钻孔方法来增加复杂地质体边界数据点密度,从而有效控制3维形体的边界,在一定程度上能缓解这一问题。以地层受断层破坏为例,分别就地层与断层4种典型空间分布组合特征下边界虚拟钻孔的确定方法进行研究,提出相应解决方案。边界虚拟钻孔和普通钻孔在数据存储方式、数据结构和表达方式等方面具有一致性,可用类似的处理方式,因而可以加快3维地质建模,提高建模精度。