Comparison of interpolation models for estimating heavy metals in soils under various spatial characteristics and sampling methods

Heavy metal pollution in soils has become increasingly challenging, especially in developing countries. Estimating the spatial distribution of heavy metals in soils is essential to preventing their build‐up. This article aims to identify the effects of spatial scales, spatial autocorrelation, sampling methods, and proportion on interpolation models in estimating the distribution of heavy metals in soils. Six interpolation models (area‐and‐point kriging, AAPK; inverse distance weighting, IDW; local polynomial interpolation, LP; ordinary kriging, OK; simple kriging, SK; and thin plate spline, TPS), three sampling methods (random, stratified, and systematic sampling), and five sampling proportions (1, 5, 10, 15, and 20%) are considered in this study using sets of simulated data, and the real situation was tested for verification. The results show that, in general, with the increase of spatial autocorrelation or the sampling percentage, the accuracy and stability of different interpolation models gradually increase; however, the various interpolation models have their own specific characteristics and application conditions. The best application conditions of the interpolation models compared with other models under the same situation are summarized and explained in theory. These conclusions have implications for future work.     

遥感图像分类精度的点、群样本检验与评估

遥感专题分类结果在使用前,必须进行客观可靠的精度验证和分析,以保持遥感分类结果的可靠性.本文利用不同分辨率遥感数据获取的同一地区土地利用/覆盖信息,进行了简单随机抽样、系统抽样和分层抽样三种不同抽样组织方式下的点样本和群样本检验分析,评估了不同抽样方式下的点样本和群样本检验效果.研究结果表明：（1）抽样方式对遥感分类精度评价结果的影响是客观存在的,不同抽样方式下的点样本和群样本检验结果都存在一定的随机性,但同一种抽样方式下,点样本检验精度评估结果的波动幅度小于群样本检验,稳定性比群样本检验要好;（2）不同抽样方式下的多次点样本和群样本检验的平均精度检验结果基本上都能够反映分类图像的精度特征,其中,点样本检验中,分层随机抽样点样本检验效果较好;群样本检验中,系统抽样群样本检验和分层随机抽样群样本检验的效果优于简单随机抽样群样本检验.     

分层随机抽样下全球30 m火烧迹地产品验证

火烧迹地是全球及区域碳循环和气候变化等研究所需的重要参数之一,卫星遥感技术为快速获取大区域火烧迹地空间分布信息提供了有效手段。中国科学院基于Landsat系列卫星数据研发了首个30 m空间分辨率全球火烧迹地产品GABAM (Global Annual Burned Area Map)。遥感数据产品的精度验证对产品使用具有重要意义,迄今尚未有研究机构对GABAM产品精度进行独立评价和分析。为系统评价GABAM产品精度,利用2010年全球30 m空间分辨率火烧迹地产品(GABAM2010)开展精度验证研究工作,在全球和几个陆地生物群落中估算了产品精度,并探索了全球遥感专题产品精度验证的技术框架。基于分层随机抽样选择80个非重叠的泰森多边形区域TSA (Thiessen Scene Areas),采用误差矩阵和6个精度指标对GABAM2010产品做全面精度评价和分析,以满足火烧迹地产品用户的使用要求。结果表明:在全球范围内,GABAM2010产品的错分率和漏分率分别为24.32%和31.60%,总体精度为97.85%;由于数据质量(如条带、云)等的影响,火烧迹地的范围会被低估,对于较容易发生火灾的生物群落,如热带亚热带草原区域,GABAM2010产品精度较高;在生物群落内部,高密度火烧迹地区域的精度高于低密度火烧迹地区域。     

An evaluation of sampling and full enumeration strategies for Fisher Jenks classification in big data settings

Large data contexts present a number of challenges to optimal choropleth map classifiers. Application of optimal classifiers to a sample of the attribute space is one proposed solution. The properties of alternative sampling‐based classification methods are examined through a series of Monte Carlo simulations. The impacts of spatial autocorrelation, number of desired classes, and form of sampling are shown to have significant impacts on the accuracy of map classifications. Tradeoffs between improved speed of the sampling approaches and loss of accuracy are also considered. The results suggest the possibility of guiding the choice of classification scheme as a function of the properties of large data sets.     

遥感影像预分类精度对地物面积空间抽样估算的敏感性分析

传统的地物面积测量受精度和效率制约,为此引入了结合遥感影像的空间分层抽样方法.首先以遥感影像的预分类结果作为模拟地物的真实分布,在地物外沿等概率随机添加不同比例的错分像元,从而获得准真实地物区的摸拟预分类结果,并依此设定各层等比例取样的样本人层标志,指导地物样本的选取,然后以抽中样本地物的准真实值之和按比例推算出总量.通过比较分析各水平含量的地物类别、不同预分类精度、层内随机和系统抽样下的多次总量估计精度及其稳定性变化情况,结果表明:该方法不需要背景数据库等先验知识,在预分类达到一定精度之上时,依分类区域设立层标志的分层抽样方法所获得的总量估计精度及标准差均好于无分类支持的随机和系统抽样;当预分类精度达到50%以上时,具有较高的成本效率比,其中在60%时,各类地物在0.5%抽样率、95%的置信度下可以保证估计量精度在92%以上.     

基于不同抽样方法的遥感面积测量方法研究

众多研究结果表明,遥感和抽样技术相结合可以有效地进行地物面积的测量。目前,随机抽样、系统抽样和分层抽样方式在 遥感抽样调查技术领域应用比较广泛。本文以遥感图像为基础,从不同角度对随机抽样、系统抽样及分层抽样（包括等样本量、等 面积、等丰度抽样）进行了有益探讨,分析发现： 对于同一地物,从平均误差百分比、标准差和极差3个角度分析,随机和系统抽 样反推得到的总量精度都低于分层抽样精度; 对于不同的地物类型,利用3种分层抽样方法反推的结果与地物所占百分比成正相关 ,地物所占的百分比越大,反推的结果越好; 等样本量、等面积、等丰度分层抽样从平均误差百分比、标准差和极差3个角度分析 各有优势,跟地物所占的百分比也有密切关系。     

不同比例尺下DEM的地形特征空间自相关性差异——以川渝地区和长三角区为例

采用局部Moran’s I模型对川渝地区和长三角区1∶250 000和1∶1 000 000两种尺度DEM以及全国1∶1 000 000 DEM地形信息进行了空间自相关的计算分析与比较,研究发现两种尺度下,高程自相关存在一致性,而坡度自相关差异性较大。1∶250 000地形自相关特性多呈随机分布,反映地势较平稳,具有均...     

A Spatial Conditioned Latin Hypercube Sampling Method for Mapping Using Ancillary Data

For obtaining maps of good precision by the spatial inference method, the distribution of sampling sites in geographical and feature space is very important. For a regional variable with trends, the predicting error comes from trend estimation, variogram estimation and spatial interpolation. Based on the cLHS (conditioned Latin hypercube Sampling) method, a sampling method called scLHS (spatial cLHS) considering all these three aspects with the help of ancillary data is proposed in this article. Its advantage lies in simultaneously improving trend estimation, variogram estimation and spatial interpolation. MODIS data and simulated data were used as sampling fields to draw sample sets using scLHS, cLHS, cLHS with x and y coordinates as covariates, simple random and spatial even sampling methods, and the distribution and prediction errors of sample sets from different methods were evaluated. The results showed that scLHS performed well in balancing spreading in geographic and feature space, and can generate points pairs with small distances, and the sample sets drawn by scLHS produced smaller mapping error, especially when there were trends in the target variable.     

Improvement in Precision of Crop Acreage Estimation by Remote Sensing Using Frequency Distribution Based Stratification

A study was conducted to improve precision of crop acreage adopting stratified random sampling approach. Remotely sensed data was used to classify mustard crop for the states of Rajasthan, Madhya Pradesh, Uttar Pradesh, Gujarat and Haryana covering 81% of mustard area of India. A grid of size 5 × 5 km was super-imposed on classified image of study area and proportion of mustard crop within the grid was ascertained. Crop proportion was used to determine strata. Stratification was done based on equal interval of proportion, equal sample number and cumulative square root of frequency method. Cumulative square root of frequency method gave highest precision in all the cases.     

A multilevel stratified spatial sampling approach based on terrain knowledge for the quality assessment of OpenStreetMap dataset in Hong Kong

With the development of Volunteered Geographical Information (VGI) data, the OpenStreetMap has high research value in terms of project activity, social influence, urban development, application scope, and historical richness and the number of buildings or roads is increasing every day. However, how to evaluate the quality of a large amount OpenStreetMaps efficiently and accurately is still not fully understood. This article presents the development of an approach regarding multilevel stratified spatial sampling based on slope knowledge and official 1:1000 thematic maps as the reference dataset for OpenStreetMap data quality inspection of Hong Kong. This multilevel stratified spatial sampling plan is as follows: (1) The terrain characteristics of Hong Kong are fully considered by dividing grids into quality estimate strata based on the slope information; (2) Spatial sampling for the selection of grids or objects is used; (3) A more reliable sampling subset is made, regarding the representation of the entire OpenStreetMap dataset of Hong Kong. This sampling plan displays a 10% higher sampling accuracy, but without increasing the sample size, particularly as regards building completeness inspection compared with simple random sampling and systematic random sampling. This research promotes further applications of the Open-Street-Map dataset, thus enabling us to have a better understanding of the OpenStreetMap data quality in urban areas.     

数字土地信息中属性数据的质量控制

提出了一种数字土地信息属性数据质量的检验、度量和分析的方法。首先给出了基于简单随机抽样和分层抽样的属性数据缺陷率度量数学模型 ,基于该统计模型 ,以某工业开发区的农村土地利用现状数据为例 ,探讨了土地利用属性数据的质量抽样方案、质量度量和质量分析的具体思路     

A Spatial Filtering Specification for an Auto-negative Binomial Model of Anopheles arabiensis Aquatic Habitats

This research accounts for spatial autocorrelation by including latent map pattern components as predictor variables in a malaria mosquito aquatic habitat model specification. The data used to derive the model was from a digitized grid-based algorithm, generated in an ArcInfo database, using QuickBird visible and near-infrared (NIR) data. The Feature Extraction (FX) Module in ENVI 4.4^® was used to categorize individual pixels of field sampled aquatic habitats into separate spectral classes, convert remotely sensed raster layers to vector coverages, and classify output layers to vector format as ESRI shapefiles. These data were used to construct a geographic weights matrix for evaluation of field and remote sampled covariates of Anopheles arabiensis aquatic habitats, a major vector of malaria in East Africa. The principal finding is that synthetic map pattern variables, which are eigenvectors computed for a geographic weights matrix, furnish an alternative way of capturing spatial dependency effects in the mean response term of a regression model. The spatial autocorrelation components suggest the presence of roughly 11 to 28% redundant information in the aquatic habitat larval count samples. The presence of redundant information in the models suggest that the sampling configuration of the An. arabiensis aquatic habitats, in the study sites, may cause field and remote observations of aquatic habitats to be dependent, rather than independent, moving data analysis away from the classical statistical independence model. A Poisson regression model, with a non-constant, gamma-distributed mean, can decompose field and remote sampled An. arabiensis data into positive and negative spatial autocorrelation eigenvectors, which can assess the precision of a malaria mosquito aquatic habitat map and the significance of all factors associated with larval abundance and distribution in a riceland agroecosystem.     

多目标微观邻域粒子群算法及其在土壤空间优化抽样中的应用

提出一种基于多目标微观邻域粒子群的土壤空间优化抽样方法。方法面向土壤空间调查的多目标特征,构建了基于最小克里金方差（MKV）和极大熵准则（ME）的粒子群多目标适应度函数,设计了最小样本量限制、样点可达性、采样成本限制和最小空间关联性四类粒子微观邻域操作策略,能高效协调土壤空间抽样精度、代表性、成本、样本量与样点布局等多目标冲突。实验结果表明,相比单目标粒子群算法和模拟退火算法,该方法的目标冲突协同能力强、收敛效率高,所设计抽样方案最优,为土壤质量精确调查与高效监测提供了技术支持。     

Mapping the distributions of C3 and C4 grasses in the mixed-grass prairies of southwest Oklahoma using the Random Forest classification algorithm

The objective of this paper is to demonstrate a new method to map the distributions of C₃ and C₄ grasses at 30 m resolution and over a 25-year period of time (1988–2013) by combining the Random Forest (RF) classification algorithm and patch stable areas identified using the spatial pattern analysis software FRAGSTATS. Predictor variables for RF classifications consisted of ten spectral variables, four soil edaphic variables and three topographic variables. We provided a confidence score in terms of obtaining pure land cover at each pixel location by retrieving the classification tree votes. Classification accuracy assessments and predictor variable importance evaluations were conducted based on a repeated stratified sampling approach. Results show that patch stable areas obtained from larger patches are more appropriate to be used as sample data pools to train and validate RF classifiers for historical land cover mapping purposes and it is more reasonable to use patch stable areas as sample pools to map land cover in a year closer to the present rather than years further back in time. The percentage of obtained high confidence prediction pixels across the study area ranges from 71.18% in 1988 to 73.48% in 2013. The repeated stratified sampling approach is necessary in terms of reducing the positive bias in the estimated classification accuracy caused by the possible selections of training and validation pixels from the same patch stable areas. The RF classification algorithm was able to identify the important environmental factors affecting the distributions of C₃ and C₄ grasses in our study area such as elevation, soil pH, soil organic matter and soil texture.     

基于多源信息与深度学习特征提取的人口空间抽样方法

人口抽样调查是通过人口样本估算区域人口总体的一种手段。由于人口分布通常具有空间差异性,传统的抽样调查理论难以满足日益增长的空间抽样需求,合理高效的人口空间抽样调查方法对于人口统计、研究人类活动、解决城市问题等有重要意义。本文提出一种基于多源信息与深度学习特征提取的人口空间抽样方法。在不透水面信息的辅助下,利用四叉树分割进行分层抽样,初步选择出可能存在人口分布的调查样本,并通过深度学习的常用模型——卷积神经网络估算样本建筑物密度,以辅助最终调查样本的选择与调查方案的制定。研究结果证明,该方法能够有效地筛选与人口分布密切相关的抽样区域,排除大量的无用样本,提高了人口调查的效率,节约了大量调查成本。     

DEM误差的空间自相关特征分析

采用空间自相关分析方法,从空间角度对数字高程数据误差的空间分布特征进行了研究。实验表明,利用双线性曲面表示地形表面时,产生的数字高程数据误差的全局Moran’sI指数趋近于0,在整个区域单元上的分布不存在显著的全局空间自相关,邻近区域单元上高程数据误差之间的关系在整体上既不综合表现为趋同,也不综合表现为趋异,高程数据误差的整体空间格局为随机格局;而且数字高程数据误差在空间上的分布与地形坡度和地表粗糙度有一定的联系,一般情况下,平均坡度、地表粗糙度越大,高程数据的全局Moran’sI指数偏离0稍远一些;否则,距离0近一些,但全局空间自相关仍不显著,在整体上表现为随机格局。     

Combining Metric Aerial Photography and Near‐Infrared Videography to Define Within‐Field Soil Sampling Frameworks

Abstract

This paper investigates the combination of metric aerial photography and near‐infrared (NIR) videography data to improve the design of field‐survey sampling frameworks. Spatial data collection can contribute up to 80% of the cost of deploying a Geographic Information System (GIS) based Decision Support System (DSS). The use of remotely sensed information, field survey using differential Global Positioning System (dGPS) and geostatistical interpolation methods maximises data quality for a given rate of sampling.

Medium‐format colour aerial photography and NIR videography were orthorectified to the national map base and mosaiced using ERDAS Imagine. The green and red layers of the aerial photography were combined with the NIR videography to form a false‐colour composite image. Two sampling strategies were tested. The first stratified sampling on a per field basis, creating four points per hectare, randomly located within each field. The second strategy used the remotely sensed information to identify within‐field variability classes for each field, using red‐green difference or normalised difference vegetation index (NDVI) models. These variability classes were used as a sub‐stratification framework with each class sampled at the same rate of 4 per hectare. For both strategies the sample points were generated within ESRI ArcView and were located in the field using dGPS. Maps of stone content were created using geostatistical methods and validated against samples collected on a 100 metre grid. It was concluded that combining the two image sources to create a within‐field stratification framework improved the precision of the results obtained from field‐survey.     

sandwichr: Spatial prediction in R based on spatial stratified heterogeneity

Spatial prediction is commonly used in social and environmental research to estimate values at unobserved locations using sampling data. However, most existing spatial prediction methods and software packages are based on the assumption of spatial autocorrelation (SAC), which may not apply when spatial dependence is weak or non-existent. In this article, we develop a modeling framework for spatial prediction based on spatial stratified heterogeneity (SSH), a common feature of geographical variables, as well as an R package called sandwichr that implements this framework. For populations that can be stratified into homogeneous strata, the proposed framework enables the estimation of values for user-defined reporting units (e.g., administrative units or grid cells) based on the mean of each stratum, even if SAC is weak or absent. The estimated values can be used to create predicted surfaces and mapping. The framework also includes procedures for selecting appropriate stratifications of the populations and assessing prediction uncertainty and model accuracy. The sandwichr package includes functions to implement each step of the framework, allowing users to implement SSH-based spatial prediction effectively and efficiently. Two case studies are provided to illustrate the effectiveness of the proposed framework and the sandwichr package.     

Design, analysis, and inference for studies comparing thematic accuracy of classified remotely sensed data: a special case of map comparison

Assessing thematic map accuracy is a special type of map comparison that is frequently applied to remote sensing classification problems. For map comparisons in the accuracy assessment setting, one map represents the classified output and the other map represents the true or “reference” condition. Several articles in this special issue describe state-of-the-art map comparison analysis tools that could serve to quantify accuracy of a single map. However, accuracy assessment objectives generally extend beyond describing accuracy of a single map to comparing accuracy of several maps. Consequently, interest focuses on comparing map comparison measures when these measures are used to represent accuracy. The virtual workshop emphasizes the analysis component of map comparisons, but it is also important to examine the underlying study designs generating the data input into these analyses. The study designs for accuracy comparisons implemented in remote sensing practice often investigate only a single test site, thus limiting our ability to generalize the results of these accuracy comparisons. Map accuracy comparison studies can be designed to provide stronger generalizations by incorporating experimental design principles such as replication and blocking, and identifying an experimental unit appropriate for the application. It is also important to recognize the role of statistical hypothesis testing and inference for different objectives that motivate map accuracy comparisons. Deciding which of two maps to use for a particular site can be addressed by enumerative inference and does not require hypothesis testing. For the objective of a more general comparison of classification procedures, analytic inference is appropriate and hypothesis testing plays a more prominent role.