Evaluation of Diverse Geometric and Geostatistical Estimation Methods Applied to Annual Precipitation in Asturias (NW Spain)

Several alternative estimation and interpolation methods for making annual precipitation maps of Asturias are analysed. The data series in this study corresponds to the year 2003. There exists an evident relationship between precipitation and altitude, with a high correlation coefficient of 0.70, that reflects the hillside effect; that is, the increase in the amount of precipitation in more mountainous areas. The direct spatial variability of precipitation and of altitude and the cross variability of precipitation–altitude are defined by two exponential variogram models: one with a short-range structure (15–30 km) that reflects the control exerted by the lesser, local mountain ranges over the amount of precipitation; and another with a long-range structure (80 km) that supposes the influence over precipitation of the major mountainous alignments of the inland areas of the Cantabrian Mountain Range (Cordillera Cantábrica) situated between 60 and 90 km from the coastline. These variogram models had to be validated for coregionalization by the Pardo-Igúzquiza and Dowd method so as to be able to make the cokriging map. The geometric estimation methods employed were triangulation and inverse distance. The geostatistical estimation methods developed were simple kriging, ordinary kriging, kriging with a trend model (universal kriging), lognormal kriging, and cokriging. In all of these methods, a 3 × 3 km² grid was selected with a total of 2580 points to estimate, a circular search window of 60 km, and a relatively small number of samples with the aim of highlighting the local features and variations on isohyet maps. The kriging methods were implemented using the WinGslib software, incorporating two specific programs, Prog2 and Fichsurf, so as to be able then to make isohyet maps using the Surfer software. All the methods employed, apart from triangulation, rendered realistic maps with good fits to the values of the original data (precipitation) of the sample maps. The problem with triangulation lies not in the reliability of the estimates but in the fact that it gives rise to contrived maps because of the tendency of isohyets to present abundant triangular facets. The reliability of the methods was based on cross-validation analysis and on evaluation of the different types of errors, both in their values and in their graphical representations. Substantial differences were not found in the values of the errors that might discriminate some methods from others in an evident way. Bearing the aforesaid in mind, should we have to make an evaluation of the different estimation methods in decreasing order of acceptance, this would be: kriging with a trend model, inverse distance, cokriging, lognormal kriging, ordinary kriging, simple kriging, and triangulation. The application of other estimation methods such as colocated cokriging, kriging with an external drift, and kriging of variable local means (residual kriging) is dependent on the availability of a digital model of the terrain with an altitude grid of the region.     

An evaluation of void‐filling interpolation methods for SRTM data

The Digital Elevation Model that has been derived from the February 2000 Shuttle Radar Topography Mission (SRTM) has been one of the most important publicly available new spatial data sets in recent years. However, the ‘finished’ grade version of the data (also referred to as Version 2) still contains data voids (some 836,000 km²)—and other anomalies—that prevent immediate use in many applications. These voids can be filled using a range of interpolation algorithms in conjunction with other sources of elevation data, but there is little guidance on the most appropriate void‐filling method. This paper describes: (i) a method to fill voids using a variety of interpolators, (ii) a method to determine the most appropriate void‐filling algorithms using a classification of the voids based on their size and a typology of their surrounding terrain; and (iii) the classification of the most appropriate algorithm for each of the 3,339,913 voids in the SRTM data. Based on a sample of 1304 artificial but realistic voids across six terrain types and eight void size classes, we found that the choice of void‐filling algorithm is dependent on both the size and terrain type of the void. Contrary to some previous findings, the best methods can be generalised as: kriging or inverse distance weighting interpolation for small and medium size voids in relatively flat low‐lying areas; spline interpolation for small and medium‐sized voids in high‐altitude and dissected terrain; triangular irregular network or inverse distance weighting interpolation for large voids in very flat areas, and an advanced spline method (ANUDEM) for large voids in other terrains.     

Assessing housing growth when census boundaries change

The US Census provides the primary source of spatially explicit social data, but changing block boundaries complicate analyses of housing growth over time. We compared procedures for reconciling housing density data between 1990 and 2000 census block boundaries in order to assess the sensitivity of analytical methods to estimates of housing growth in Oregon. Estimates of housing growth varied substantially and were sensitive to the method of interpolation. With no processing and areal‐weighted interpolation, more than 35% of the landscape changed; 75–80% of this change was due to decline in housing density. This decline was implausible, however, because housing structures generally persist over time. Based on aggregated boundaries, 11% of the landscape changed, but only 4% experienced a decline in housing density. Nevertheless, the housing density change map was almost twice as coarse spatially as the 2000 housing density data. We also applied a dasymetric approach to redistribute 1990 housing data into 2000 census boundaries under the assumption that the distribution of housing in 2000 reflected the same distribution as in 1990. The dasymetric approach resulted in conservative change estimates at a fine resolution. All methods involved some type of trade‐off (e.g. analytical difficulty, data resolution, magnitude or bias in direction of change). However, our dasymetric procedure is a novel approach for assessing housing growth over changing census boundaries that may be particularly useful because it accounts for the uniquely persistent nature of housing over time.     

用于强震记录仿真地动位移的高精度累积数值积分方法研究

针对传统梯形积分方法的低精度缺陷,提出了一套高精度累积数值积分方法.首先对被积数据进行分段三次Hermite插值,之后采用积分区间四等分的Newton-Cotes积分公式进行积分.误差分析表明,本方法精度远高于梯形积分方法精度.在对构建信号进行方法验证的前提下,采用2013年4月20日雅安芦山地震和2008年5月12日汶川地震的真实强震动加速度记录对本方法进行测试,测试结果表明,本方法适用于将强震动加速度记录积分为位移记录.     

http://www.sciencedirect.com/science/article/pii/S1674987111001411

正Three-dimensional geological modeling(3DGM) assists geologists to quantitatively study in three-dimensional(3D) space structures that define temporal and spatial relationships between geological objects.The 3D property model can also be used to infer or deduce causes of geological objects.3DGM technology provides technical support for extraction of diverse geoscience information,3D modeling,and quantitative calculation of mineral resources.Based on metallogenic concepts and an ore deposit model, 3DGM technology is applied to analyze geological characteristics of the Tongshan Cu deposit in order to define a metallogenic model and develop a virtual borehole technology;a BP neural network and a 3D interpolation technique were combined to integrate multiple geoscience information in a 3D environment. The results indicate:(1) on basis of the concept of magmatic-hydrothermal Cu polymetallic mineralization and a porphyry Cu deposit model,a spatial relational database of multiple geoscience information for mineralization in the study area(geology,geophysics,geochemistry,borehole,and cross-section data) was established,and 3D metallogenic geological objects including mineralization stratum,granodiorite, alteration rock,and magnetic anomaly were constructed;(2) on basis of the 3D ore deposit model,23,800 effective surveys from 94 boreholes and 21 sections were applied to establish 3D orebody models with a kriging interpolation method;(3) combined 23,800 surveys involving 21 sections,using VC++ and OpenGL platform,virtual borehole and virtual section with BP network,and an improved inverse distance interpolation(IDW) method were used to predict and delineate mineralization potential targets (Cu-grade of cell not less than 0.1%);(4) comparison of 3D ore bodies,metallogenic geological objects of mineralization,and potential targets of mineralization models in the study area,delineated the 3D spatial and temporal relationship and causal processes among the ore bodies,alteration rock,metallogenic stratum,intrusive rock,and the Tongshan Fault.This study provides important technical support and a scientific basis for assessment of the Tongshan Cu deposit and surrounding exploration and mineral resources.     

黄土高原降水量空间插值方法研究

为了研究黄土高原降水量的空间分布特征,利用ARCGIS软件分别采用多种空间内插方法对黄土高原108个气象站点1980-2009年30年间的平均年降水量进行空降插值。经过正交检验以及内插检验得到最优的插值方法为普通克里金法,最终对利用该方法获得的黄土高原降水量内插图进行分析,得出近30年来黄土高原区降水量分布呈现出明显的由东南向西北递减的态势。     

Performance of geostatistical interpolation methods for modeling sampled data with non-stationary mean

The measured ozone pollution peak in the atmosphere of Mexico City region was considered in order to study the effect of a non-stationary mean of the sampled data in geostatistics interpolation methods. With this objective the local mean value of the sampled data was estimated through a linear regression analysis of their values on the monitoring station’s coordinates. The residuals obtained by removing the data trend are considered as a set of stationary random variables. Several interpolation methods used in geostatistics, such as inverse distance weighted, kriging, and artificial neural networks techniques were considered. In an effort to optimize and evaluate its performance, we fit interpolated values to sampled data, obtaining optimal values for the parameters defining the used model, that means, the values of the parameters that give the lowest mean RMSE between the interpolated value and measured data at 20 stations at 1500 hours for a set of 21 days of December 2001, which was chosen as the training set. The training set is conformed by all the days in December 2001 excepting the days (3,6,9,12,...,27,30) which were considered as the testing set. Once the optimal model is obtained, it is used to interpolate the values at the stations at 1500 hours for the testing days. The RMSE between interpolated and measured values at monitoring stations was also evaluated for these testing values and is shown as a percentage in Table 2. These values and the defined generalization parameter G, can be used to evaluate the performance and the ability of the models to predict and reproduce the peak of ozone concentrations. Scatter plots for testing data are presented for each interpolation method. An interpretation of the ozone pollution levels obtained at 1500 hours at December 21 was given using the wind field that prevailed in the region 1 h before the same day.     

基于Kriging估计的GPS高程转换

随着GPS测量精度的不断提高，其平面测量已经在多种领域得到应用，成为测量工作的一个全新的手段与工具，但在GPS测高方面需要改善和研究。Kliging方法是地统计学的基本方法，其理论核心是变异函数，它是一种最好的线性无偏估计法，在地理学、生态学、环境科学等领域都得到了广泛的应用，取得了很好效果。本文在详细介绍kriging法的基础上，建立了Kriging法的GPS高程拟合模型，利用该模型对两个实例的GPS高程数据拟合，得到了正常高，结果表明Kriging法在GPS高程拟合中的应用精度达到厘米级，且得出了一些有益的结论，为以后GPS高程的应用提供一些参考。     

四大流域面雨量监控预报系统

基于中国气象局下发的“关于全国七大流域面雨量预报业务实施方案的通知”和河南省气象局下发的“黄河中下游面雨量预报业务实施细则”，开发了“四大流域面雨量监控预报系统”。系统程序设计采用C＋＋Builder4.0和Visual Fortran 5.0，实现了降水实况自动入库、全省站点雨量预报及四大流域面雨量预报的计算输出，达到了对四大流域面雨量进行实时预报、监控的目的。     

Volume-preserving interpolation of a smooth surface from polygon-related data

 The interpolation of continuous surfaces from discrete points is supported by most GIS software packages. Some packages provide additional options for the interpolation from 3D line objects, for example surface-specific lines, or contour lines digitized from topographic maps. Demographic, social and economic data can also be used to construct and display smooth surfaces. The variables are usually published as sums for polygonal units, such as the number of inhabitants in communities or counties. In the case of point and line objects the geometric properties have to be maintained in the interpolated surface. For polygon-based data the geometric properties of the polygon boundary and the volume should be preserved, avoiding redistribution of parts of the volume to neighboring units during interpolation. The pycnophylactic interpolation method computes a continuous surface from polygon-based data and simultaneously enforces volume preservation in the polygons. The original procedure using a regular grid is extended to surface representations based on an irregular triangular network (TIN). Received: 5 October 2000 / Accepted: 20 December 2000