首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 125 毫秒
1.
Surface modeling of DEMs based on a sequential adjustment method   总被引:1,自引:0,他引:1  
A sequential adjustment (SA) method is employed to decrease the computational cost of high-accuracy surface modeling (HASM), and the SA of HASM (HASM-SA) is being developed. A mathematical surface was used to comparatively analyze the computing speed of SA and the classical iterative solvers provided by MATLAB 7.7 for solving the system of linear equations of HASM. Results indicate that SA is much faster than the classical iterative solvers. The computing time of HASM-SA is determined by not only the total number of grid cells but also the number of sampling points in the computational domain. A real-world example of surface modeling of digital elevation models (DEMs) with various resolutions shows that HASM-SA is averagely more accurate and much faster than the commonly used interpolation methods, such as inverse distance weighting (IDW), kriging, and three versions of spline, namely regularized spline (RSpline), thin-plate spline (TPS), and ANUDEM in terms of root mean square error (RMSE), mean absolute error (MAE), and mean error (ME). In particular, the ME of HASM-SA at different spatial resolutions is averagely smaller than those of IDW, kriging, RSpline, TPS, and ANUDEM by 85%, 83%, 83%, 53%, and 19%, respectively. The high speed and high accuracy of HASM-SA can be due to the absence of matrix inversion computation, combined with the perfect fundamental theorem of HASM.  相似文献   

2.
We developed a method for high-accuracy surface modeling (HASM) in terms of the fundamental theorem of surfaces, which has theoretically found a solution for error problems. However, the computing speed of HASM was too slow to be widely applied in practice. Thus, adjustment computation of HASM (HASM-AC) is developed in this article. For comparatively testing HASM-AC accuracy, a mathematical surface is first selected so that the true value is able to be predetermined to avoid uncontrollable data errors. The numerical test indicates that HASM-AC has the highest accuracy and its accuracy is 20.67, 15.67, and 14.67 times higher than the inverse distance weighting (IDW), kriging, and spline, respectively. Then, a 0.4 km × 0.5 km rectangular area is used to test the effects of different spatial resolutions and sampling intervals on accuracy. This real-world test demonstrates that HASM-AC accuracy increases at a much better and stable pace as the spatial resolution becomes finer and sampling intervals get shorter, compared to the classic methods. Finally, the computing speed is tested in an area with 6000 × 6000 grid cells where Qinghai, Gansu, and Sichuan provinces meet. The computing speed of HASM-AC is 11, 8, and 563 times faster than IDW, spline, and kriging, respectively, which makes HASM-AC able to process data in a huge size and make real-time visualization realizable. In short, HASM-AC performs best in both the numerical and real-world tests.  相似文献   

3.
赵娜  岳天祥  赵明伟 《地理研究》2014,33(7):1297-1305
日照百分率作为研究日照时数及太阳辐射等的重要因素之一,其模拟结果的好坏,直接关系到相关领域的研究应用。而高精度曲面建模方法(HASM)是近几年发展起来的用于生态建模的高精度曲面模拟方法。首先对现有的HASM进行改进,给出建立在完整理论基础之上、精度更高的曲面建模方法,并记为HASM.MOD;以高斯合成曲面为数值试验对象,验证HASM与HASM.MOD的模拟精度;最后,根据全国1951-2010 年752 个气象站点的月平均日照百分率数据,运用HASM.MOD研究近60 年月平均日照百分率的分布状况,同时比较了HASM.MOD、HASM、Kriging 和IDW法的插值精度。数值试验和实例验证结果表明,HASM.MOD的模拟精度最高。用该方法所提供的日照百分率数据可作为基础地理数据供相关研究应用。  相似文献   

4.
High accuracy surface modeling(HASM) is a method which can be applied to soil property interpolation.In this paper,we present a method of HASM combined geographic information for soil property interpolation(HASM-SP) to improve the accuracy.Based on soil types,land use types and parent rocks,HASM-SP was applied to interpolate soil available P,Li,pH,alkali-hydrolyzable N,total K and Cr in a typical red soil hilly region.To evaluate the performance of HASM-SP,we compared its performance with that of ordinary kriging(OK),ordinary kriging combined geographic information(OK-Geo) and stratified kriging(SK).The results showed that the methods combined with geographic information including HASM-SP and OK-Geo obtained a lower estimation bias.HASM-SP also showed less MAEs and RMSEs when it was compared with the other three methods(OK-Geo,OK and SK).Much more details were presented in the HASM-SP maps for soil properties due to the combination of different types of geographic information which gave abrupt boundary for the spatial varia-tion of soil properties.Therefore,HASM-SP can not only reduce prediction errors but also can be accordant with the distribution of geographic information,which make the spatial simula-tion of soil property more reasonable.HASM-SP has not only enriched the theory of high accuracy surface modeling of soil property,but also provided a scientific method for the ap-plication in resource management and environment planning.  相似文献   

5.
基于多重网格求解的土壤属性高精度曲面建模   总被引:4,自引:1,他引:3  
高精度曲面建模(HASM)是近几年发展起来的可用于地理信息系统和生态建模的一种较高精度的曲面建模方法.本研究选择南方红壤丘陵区江西省吉安市辖区、吉安县和泰和县为研究区,采集150个表层土壤(0~20cm)样品,分别随机选取60、90和120个点作为模拟数据集,90、60和30个点作为验证数据集,基于多重网格(MG)作为...  相似文献   

6.
高精度曲面建模方法(HASM)作为新一代的曲面模拟方法,经过20多年的发展,其理论体系不断完善,算法体系不断丰富,应用领域不断拓广。然而,目前HASM方法尚未建立科学的分类体系,仅依据求解算法为标准进行简单区分,忽视了HASM所研究问题本质上的差异,阻碍了HASM方法在相关领域的进一步应用。为此,本文在总结HASM方法发展过程的基础上,按照研究问题的本质特点,以输入数据的类型为依据,将HASM分为空间插值方法和空间数据融合方法。其中,HASM空间插值方法是根据离散采样点得到目标曲面,而HASM空间数据融合方法则是融合多源数据,并综合各个数据源优势而得到新曲面的过程。该分类科学、合理,为今后HASM方法的进一步应用提供了理论指导。最后本文叙述了应用两种HASM方法求解问题时的一般步骤,同时还对两种方法的发展前景进行了展望。  相似文献   

7.
中国表层土壤全氮的空间模拟分析   总被引:11,自引:0,他引:11  
基于第二次全国土壤普查5336个典型土壤剖面数据,分析表土全氮(A层)与环境因素的相关关系,利用多元回归模型和HASM模型结合的方法模拟中国国家尺度上表层土壤全氮的空间分布格局。结果表明:对350个检验点模拟结果的平均绝对误差和平均相对误差为0.67g·kg-1和61.06%,与普通克里格法相比分别降低了0.05g·kg-1和17.53%;对样点分布较少的西北地区的模拟结果也更符合实际情况。多元回归模型和HASM模型结合的方法考虑了环境因素的影响,可作为目前模拟大尺度土壤性质空间分布相对有效的方法。  相似文献   

8.
This paper revisits the computation of product combinations for quantification of resources (e.g., pore volume for hydrocarbon reservoir formations). The explanations have been simplified by considering the exhaustive numerical model for properties that are multiplied first and then added or averaged for evaluation. The analysis starts without any probabilistic considerations. Abbreviated (up-scaled) computations in the sum of multiplications are proposed by substituting the individual values by averages for each rock property in the coarser cell resolution model. A result found is that averaged properties can be utilized for estimation at non-sampled locations, instead of individual values; however, covariances and cumulants must also be included in the abbreviated computations. The smoothing effect of kriging is found to be irrelevant if the kriging variance is also included in the up-scaled abbreviated pore volume computations. Thus, the equivalence between computation of resource volumes from kriging estimates and conditional stochastic simulations is established, with the condition that numerical estimation must incorporate the complete covariance and cumulant information as well. An example shows pore volume prediction from a kriging model matches the unbiased result from stochastic simulations.  相似文献   

9.
利用1951-2010年中国711个气象观测站的月降水资料,对多年季平均降水根据中国农业气候类型进行分区模拟。针对中国降水特点,首先分析了影响各分区降水的地理、地形因素及局部地形因素,利用多项式回归和逐步回归的方法对各分区降水进行了趋势拟合;在此基础上,采用改进的高精度曲面建模(HASM)方法,对各模拟区域去掉趋势后的残差进行迭代修正,并比较验证了模拟效果。同时,为保证HASM在边界附近的模拟精度,根据区域内站点间的距离,对每一分区设置一个缓冲区,将HASM实际插值区域扩展为缓冲区内的部分。模拟结果表明:HASM方法的模拟精度在不同区域不同季节内均比经典的插值方法模拟精度高。利用上述方法分析了同一季节各分区降水的分布特点,并模拟了不同季节内多年平均降水的空间分布状况,模拟结果符合我国降水的实际分布特点。  相似文献   

10.
基于地学信息的土壤属性高精度曲面建模   总被引:5,自引:1,他引:4  
史文娇  刘纪远  杜正平  岳天祥 《地理学报》2011,66(11):1574-1581
目前土壤属性的曲面建模面临两大问题,一是缺乏足够多的采样点进行模拟,二是采样点的土壤属性与环境变量间存在非线性关系。高精度曲面建模(HASM)方法在如何融合地学信息对土壤属性进行空间插值,尚需深入研究。本文提出高精度曲面建模与地学信息相结合的土壤属性空间插值方法(High accuracy surface modeling for soil properties, HASM-SP),分别基于研究区的土壤类型、土地利用类型和母岩类型,模拟了江西省典型红壤丘陵区土壤的速效磷(AP)、锂(Li)、pH、碱解氮(AN)、全钾(K) 和铬(Cr) 6 种土壤属性的空间分布。将HASM-SP与普通克立格(OK)、OK与地学信息相结合的方法(OK-Geo) 和分层克立格(SK) 相对比,结果表明:结合地学信息可以使插值方法(HASM-SP和OK-Geo) 模拟结果的平均误差更趋近于0;相对于其他3 种方法,HASM-SP具有较小的平均绝对误差和均方根误差。从模拟的空间分布图上来看,HASM-SP获得了由于地学要素类型突变导致的土壤属性空间变异的细节信息。因此,基于地学信息的土壤属性高精度曲面建模方法(HASM-SP) 与传统插值方法相比,不仅提高了土壤属性空间分布的模拟精度,还能更好地刻画突变边界处土壤属性的空间变异,使土壤属性分布图更好地符合地学规律和实际情况。HASM-SP的提出对于丰富土壤属性的高精度曲面建模理论具有重要意义,并为土壤管理、精准农业的实施以及区域环境规划等提供科学依据。  相似文献   

11.
The shuttle radar topography mission (SRTM), was flow on the space shuttle Endeavour in February 2000, with the objective of acquiring a digital elevation model of all land between 60° north latitude and 56° south latitude, using interferometric synthetic aperture radar (InSAR) techniques. The SRTM data are distributed at horizontal resolution of 1 arc‐second (~30 m) for areas within the USA and at 3 arc‐second (~90 m) resolution for the rest of the world. A resolution of 90 m can be considered suitable for the small or medium‐scale analysis, but it is too coarse for more detailed purposes. One alternative is to interpolate the SRTM data at a finer resolution; it will not increase the level of detail of the original digital elevation model (DEM), but it will lead to a surface where there is the coherence of angular properties (i.e. slope, aspect) between neighbouring pixels, which is an important characteristic when dealing with terrain analysis. This work intents to show how the proper adjustment of variogram and kriging parameters, namely the nugget effect and the maximum distance within which values are used in interpolation, can be set to achieve quality results on resampling SRTM data from 3” to 1”. We present for a test area in western USA, which includes different adjustment schemes (changes in nugget effect value and in the interpolation radius) and comparisons with the original 1” model of the area, with the national elevation dataset (NED) DEMs, and with other interpolation methods (splines and inverse distance weighted (IDW)). The basic concepts for using kriging to resample terrain data are: (i) working only with the immediate neighbourhood of the predicted point, due to the high spatial correlation of the topographic surface and omnidirectional behaviour of variogram in short distances; (ii) adding a very small random variation to the coordinates of the points prior to interpolation, to avoid punctual artifacts generated by predicted points with the same location than original data points and; (iii) using a small value of nugget effect, to avoid smoothing that can obliterate terrain features. Drainages derived from the surfaces interpolated by kriging and by splines have a good agreement with streams derived from the 1” NED, with correct identification of watersheds, even though a few differences occur in the positions of some rivers in flat areas. Although the 1” surfaces resampled by kriging and splines are very similar, we consider the results produced by kriging as superior, since the spline‐interpolated surface still presented some noise and linear artifacts, which were removed by kriging.  相似文献   

12.
An efficient monitoring network is very important in accessing the marine environmental quality and its protection and management. In an estuary, there are fronts that separate distinctly different water masses and affect material transport, nutrient distribution, pollutant aggregation, and diffusion. This stratified heterogeneous surface neither satisfies the stationary requirements of kriging, nor can be handled adequately by removing a spatially continuous trend. This article presents a stratified optimization method for a multivariate monitoring network. In this method, principal component analysis (PCA) was used to reduce the dimensionality of the correlated targets, and the mean of surface with nonhomogeneity (MSN) method was adopted to produce the best linear unbiased estimator for a spatially stratified heterogeneous surface that failed to satisfy the requirements for a kriging estimate. The existing monitoring network in the Yangtze River estuary and its adjacent sea, which was designed by purposive sampling year ago was optimized as an illustration. The optimization consisted of two steps: reduce the redundant monitoring sites and then optimally add new sites to the remaining sites. After optimization, the inclusion of 51 sites in the monitoring network was found to produce a smaller total estimated error than that of the current network, which has 70 sites; moreover, the use of 55 sites can produce a higher precision of estimation for all three principal components (PCs) than that of the current 70 sites. The results demonstrated that the proposed method is suitable for optimizing environmental monitoring sites that have dominant stratified nonhomogeneity and that involve multiple factors.  相似文献   

13.
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 km2 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.  相似文献   

14.
中国土壤温度的空间插值方法比较   总被引:15,自引:1,他引:14  
利用中国698个气象站点1971~2000年的地面气候资料,采用三种不同方法预测中国0cm、20cm和40cm深度年均土壤温度的空间分布,其中普通克里格和泛克里格法直接以年均土壤温度数据为源数据、回归克里格法以中国年均气温数据和中国DEM数据为源数据进行预测。预测结果的准确性通过平均绝对误差(MAE)和均方根误差(RMSE)值来评价。结果表明回归克里格法预测的MAE值和RMSE值均为最小,说明其预测结果的准确性最好、预测的极端误差也最小;其次为泛克里格法;普通克里格法预测的效果最差。回归克里格法预测结果由于采用了中国DEM数据进行修正,在空间特征表达方面能够更好地表达复杂地形地区的局部变异,其平滑效应明显小于泛克里格法和普通克里格法的预测结果。  相似文献   

15.
Spatial interpolation of marine environment data using P-MSN   总被引:1,自引:0,他引:1  
ABSTRACT

When a marine study area is large, the environmental variables often present spatially stratified non-homogeneity, violating the spatial second-order stationary assumption. The stratified non-homogeneous surface can be divided into several stationary strata with different means or variances, but still with close relationships between neighboring strata. To give the best linear-unbiased estimator for those environmental variables, an interpolated version of the mean of the surface with stratified non-homogeneity (MSN) method called point mean of the surface with stratified non-homogeneity (P-MSN) was derived. P-MSN distinguishes the spatial mean and variogram in different strata and borrows information from neighboring strata to improve the interpolation precision near the strata boundary. This paper also introduces the implementation of this method, and its performance is demonstrated in two case studies, one using ocean color remote sensing data, and the other using marine environment monitoring data. The predictions of P-MSN were compared with ordinary kriging, stratified kriging, kriging with an external drift, and empirical Bayesian kriging, the most frequently used methods that can handle some extent of spatial non-homogeneity. The results illustrated that for spatially stratified non-homogeneous environmental variables, P-MSN outperforms other methods by simultaneously improving interpolation precision and avoiding artificially abrupt changes along the strata boundaries.  相似文献   

16.
This paper is concerned with the problem of predicting the surface elevation of the Braden breccia pipe at the El Teniente mine in Chile. This mine is one of the world’s largest and most complex porphyry-copper ore systems. As the pipe surface constitutes the limit of the deposit and the mining operation, predicting it accurately is important. The problem is tackled by applying a geostatistical approach based on closed-form non-stationary covariance functions with locally varying anisotropy. This approach relies on the mild assumption of local stationarity and involves a kernel-based experimental local variogram a weighted local least-squares method for the inference of local covariance parameters and a kernel smoothing technique for knitting the local covariance parameters together for kriging purpose. According to the results, this non-stationary geostatistical method outperforms the traditional stationary geostatistical method in terms of prediction and prediction uncertainty accuracies.  相似文献   

17.
基于HASM的中国森林植被碳储量空间分布模拟   总被引:2,自引:1,他引:1  
赵明伟  岳天祥  赵娜  孙晓芳 《地理学报》2013,68(9):1212-1224
当前区域尺度上森林碳储量估算主要依据森林资源清查数据,整个过程不仅消耗大量人力、物力,而且十分耗时,严重影响了森林碳储量估算的时效性。针对这一问题,本文提出了基于HASM的森林植被碳储量模拟方法,该方法以全球植被动态模型LPJ-Guess 输出的植被碳储量为驱动场,以森林清查样地数据为精度控制点,模拟生成中国陆地森林碳储量分布情况。研究以第7 次中国森林资源清查数据作为精度控制点数据源,同时作为本文模拟方法的精度验证。结果表明,中国森林碳储量为9.2405 Pg,考虑到森林资源清查是基于一定的郁闭度进行的,因此HASM模拟的结果与根据森林资源清查结果计算得出的7.8115 Pg 相比更符合实际情况,西南山区和东北林区是中国森林最主要的碳库,其碳储量分别占中国森林植被碳储量的39.82%和20.46%。同时与之前(1975-1995 年) 相比具有较大幅度的增长,表明近几十年来中国坚持大规模植树造林的碳汇效果显著。同时也表明基于HASM的森林植被碳储量空间分布模拟方法是有效的,模拟结果合理且精度较高,表明该方法在全球尺度上森林植被碳储量模拟及其它生态系统中碳储量模拟中具有应用潜力。  相似文献   

18.
喀斯特地区春季土壤水分空间插值方法对比   总被引:1,自引:0,他引:1  
以杨眉河小流域为研究区,通过土壤水分采样,选取辅助变量,采用普通克里金、协同克里金、回归克里金3种地统计学方法对土壤水分数据进行空间插值。结果表明:1)回归克里金对研究区土壤水分估算误差最小,其次为协克里金,普通克里金的误差最大;2)普通克里金生成的土壤水分表面最为平滑,而回归克里金最大程度反映了研究区实际的土壤水分空间变化;3)对于协同克里金,以湿度指数(WI)样点数据作为辅助变量的估算误差小于将WI栅格数据作为辅助变量的估算误差。总之,在可获得有效辅助变量的条件下,回归克里金对研究区土壤水分估算的效果优于协同克里金与普通克里金。  相似文献   

19.
It is very important in accurately estimating the forests' carbon stock and spatial distribution in the regional scale because they possess a great rate in the carbon stock of the terrestrial ecosystem. Yet the current estimation of forest carbon stock in the regional scale mainly depends on the forest inventory data, and the whole process consumes too much labor, money and time. And meanwhile it has many negative influences on the forest carbon storage updating. In order to figure out these problems, this paper, based on High Accuracy Surface Modeling (HASM), proposes a forest vegetation carbon storage simulation method. This new method employs the output of LPJ-GUESS model as initial values of HASM and uses the inventory data as sample points of HASM to simulate the distribution of forest carbon storage in China. This study also adopts the seventh forest resources statistics of China as the data source to generate sample points, and it also works as the simulation accuracy test. The HASM simulation shows that the total forest carbon storage of China is 9.2405 Pg, while the calculated value based on forest resources statistics are 7.8115 Pg. The forest resources statistics is taken based on a forest canopy closure, and the result of HASM is much more suitable to the real forest carbon storage. The simulation result also indicates that the southwestern mountain region and the northeastern forests are the important forest carbon reservoirs in China, and they account for 39.82% and 20.46% of the country's total forest vegetation carbon stock respectively. Compared with the former value (1975-1995), it mani- fests that the carbon storage of the two regions do increase clearly. The results of this re- search show that the large-scale reforestation in the last decades in China attains a signifi- cant carbon sink.  相似文献   

20.

Interpolation of point measurements using geostatistical techniques such as kriging can be used to estimate values at non-sampled locations in space. Traditional geostatistics are based on the spatial autocorrelation concept that nearby things are more related than distant things. In this study, additional information was used to modify the traditional Euclidean concept of distance into an adjusted distance metric that incorporates similarity in terms of quantifiable landscape characteristics such as topography or land use. This new approach was tested by interpolating soil moisture content, pH and carbon-to-nitrogen (C:N) ratio measured in both the mineral and the organic soil layers at a field site in central Sweden. Semivariograms were created using both the traditional distance metrics and the proposed adjusted distance metrics to carry out ordinary kriging (OK) interpolations between sampling points. In addition, kriging with external drift (KED) was used to interpolate soil properties to evaluate the ability of the adjusted distance metric to incorporate secondary data into interpolations. The new adjusted distance metric typically lowered the nugget associated with the semivariogram, thereby better representing small-scale variability in the measured data compared to semivariograms based on the traditional distance metric. The pattern of the resulting kriging interpolations using KED and OK based on the adjusted distance metric were similar because they represented secondary data and, thus, enhanced small-scale variability compared to traditional distance OK. This created interpolations that agreed better with what is expected for the real-world spatial variation of the measured properties. Based on cross-validation error, OK interpolations using the adjusted distance metric better fit observed data than either OK interpolations using traditional distance or KED.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号