中高纬度山区气温空间化的方法比较研究——以大兴安岭北麓为例

为比较和探讨中高纬度山区多种气温空间插值方法的精度及适用性,本文利用大兴安岭山脉北段及其周边区域气象站点实测气温数据,以平均绝对误差(MAE)和均方根误差(RMSE)作为评价指标对六种气温空间插值方法进行精度比较.研究结果表明:(1)反距离权重插值法(IDW)、普通克里金插值法(OK)、样条函数插值法(Spline)三...     

基于RS和GIS的人口估计方法研究综述

随着RS和GIS技术的迅速发展, RS和GIS成为进行人口估计的主要手段。根据人口估计的目标和使用的数据源不同,可以把基于RS和GIS的人口估计方法分成两大类：面插值方法和统计模型方法。面插值方法根据插值过程中是否使用辅助数据可以进一步分成无辅助数据的面插值和有辅助数据的面插值两种。统计模型方法根据模型中自变量的不同可以分成建成区面积估计法、土地利用密度法、居住单元估计法、图像像元特征估计法和自然和社会经济特征综合估计法五种。本文按照上述分类标准综述了基于RS和GIS的各种人口估计方法,分析了各种方法的应用条件、优缺点和研究实例。最后提出了在基于RS和GIS进行人口估计方面需要进一步研究的问题。     

地理数据尺度转换方法研究进展

尺度问题在地理学、生态学和水文科学等众多领域都具有非常重要的地位。近年来,随着对地观测技术和地理信息科学的飞速发展,解决地理数据的尺度转换问题成为目前地理信息科学及相关研究中的热点和难点问题之一。在地理信息科学相关领域中,地图学和遥感科学研究人员从尺度概念的理解到尺度转换的理论和方法都做了大量的研究,对解决地理数据空间特征的尺度转换做出了重要的贡献。在地理数据属性特征的尺度转换研究领域,地理信息科学研究者提出的面域插值方法是解决此问题的主要方法之一。同时,在社会经济科学领域,“小区域统计学”也发展了一套相关的理论和方法,试图解决统计单元的变更问题。文章在全面回顾和比较不同研究领域解决“尺度转换”问题方法的基础上,重点介绍面域插值方法和小区域统计学的基本原理及其典型应用。     

城市土地级别划分的方法研究——以北京城市土地定级为例

城市土地定级是制订基准地价的基础，也是有效进行土地资源合理配置和土地管理的重要手段，文章主要探讨空间插值方法在北京城市土地定级中的应用，认为空间插值算法是一种通过已知点的数据推求同一区域其它未知点数据的计算方法，在土地市场发达，土地交易案例较多的城镇，可以适用于运用市场交易地价划分土地级别。     

空间软数据及其插值方法研究进展

由于对地观测技术的迅速发展,空间数据的种类和数量增长迅猛,由空间数据反演得到的各种信息日趋膨胀,这些反演结果中的信息不少以软数据的形式出现。在实际应用中,这些软数据往往与空间插值的目标变量具有一定的相关性,甚至成为控制目标变量空间分布特征的重要因素。然而,由于这些数据通常表示为非数值形式,在计算和处理上存在着一定困难,以致被传统的插值方法所忽视,从而造成信息浪费。近来出现的空间软插值方法是一种利用空间软数据作为辅助信息并以改善插值效果的方法,能够较好的处理并利用软数据所隐含的信息,具有较好的应用发展前景。本文根据空间软数据的特点及其分类,系统综述了空间软插值方法及其应用领域。首先分析了空间数据软硬性质的根本区别,论述了软数据的分类和“硬化”方法,然后介绍空间插值模型中对空间软数据的集成方法和原理,最后对空间软插值方法及其应用研究领域进行了展望。     

降雨信息空间插值的不确定性分析

文章以潮白河流域为样区,根据58个雨量站1990年的降雨观测数据,采用反距离权重法、克立格法、样条函数法、趋势面法等插值方法,分析了站点数量变化、时间尺度变化、栅格像元的尺度变化、插值方法的差异对降雨数据空间插值结果的影响,剖析降雨插值中的不确定性。结果表明：（1）插值站点数量越大,区域降雨插值的不确定性越小;（2）像元尺度在50m~1000m间变化对降雨插值的不确定性只有微弱的影响;（3）对应于时间尺度由年到月到日的变化,降雨插值的不确定性随时间尺度的减小而显著增大;（4）不同插值方法影响到降雨空间插值的不确定性水平。为了减少降雨信息空间插值的不确定性,根本途径是要引入第三方相关变量,并将其整合到现有的插值算法中。高相关性变量的选取及其与插值模型的整合方式将成为降雨插值研究的主导方向。     

小区域癌症数据典型空间插值方法比较研究

癌症已成为危害全球居民健康的重大民生问题,选取合适的空间插值方法分析小区域癌症数据的空间特征可对区域性癌症防控工作的有效开展提供依据。本研究以湖南省苏仙区2012和2016年以村为单位的肺癌死亡率数据为研究对象,以平均误差和均方根误差为评价指标,对反距离加权（IDW）、普通克里金（OK）、趋势面分析（TSA）、多元线性回归（MLR）与协同克里金（CK）五种典型空间插值方法进行精度效果对比及参数优选,并结合不同插值方法的优缺点,确定癌症数据的最优插值方法。结果表明：插值精度方面,CK法的均方根误差最小、插值精度最高,OK、IDW（幂值=1）和MLR次之,TSA（阶数=5）最低;插值效果方面,五种插值方法的实测值和预测值均显著相关,除CK外,其它四种方法均对死亡率低估程度较大,CK和OK插值结果的空间分布效果更好。同时考虑空间因素和影响因子的CK方法是小区域苏仙区2012年、2016年肺癌死亡率最优插值方法,应用该方法可对区域性癌症防控工作的有效开展提供最优的技术支撑。本论文的研究思路也可为小区域癌症数据空间插值方法及参数优选提供参考。     

一种基于超曲面样条函数的三维空间插值方法

将二维曲面样条函数插值法进一步引申,提出一种基于超曲面样条函数进行三维空间插值的新方法,介绍超曲面样条函数的构造方法,给出利用该方法进行三维空间插值的数值实现过程。最后以插值形成三维空间规则数据场为例,证明利用超曲面样条函数法进行三维空间插值的可行性和准确性。     

复杂地形下降水的高空间分辨率插值方法研究

复杂地形山区降水格局在多种地形要素的综合影响下呈现出显著异质性特征,弱化了降水站点观测资料的空间代表性,限制了遥感及再分析产品的适用性及传统插值方法的准确性。常用的PRISM降水插值算法通过提取并权重化地形要素,借助加权最小二乘算法对降水站点观测数据进行空间插值,被广泛应用于降水产品制备。本文针对PRISM算法对中小地形地貌刻画能力不足的问题,在解析影响复杂地形降水模式的地形要素的基础上,改进了PRISM的地形要素计算与权重化过程,同时,考虑到实际日降水量的随机性,将改进后的PRISM嵌入到“平均态日降水—比值”插值框架,构建了适用于复杂地形的日降水插值算法MPRISMR。随后,以具有复杂地形特征的元江流域为例,通过交叉验证及与ERA5-Land和TRMM_3B42降水产品的对比分析,发现该算法具有较高的准确性与可靠性。结果表明,在元江流域23个气象站点上,MPRISMR的插值结果与观测值的相关系数、相对偏差的中位数分别为0.72、0.98%,总体上优于ERA5-Land和TRMM_3B42日降水产品。另外,MPRISMR插值结果的精度随时间变化更小,更为稳定。最后,研究基于MPRISMR制备了空间分辨率约3 km的元江流域日尺度降水格网数据。本文可为复杂地形的陆面模式或流域水文水质模型提供高精度降水驱动场数据产品,从而支持流域可持续管理决策。     

不同时间尺度、季节的气温数据空间平稳特征及其对插值结果的影响

不同时间尺度、季节的气温数据表现出不同的空间平稳特征。为探讨分析空间平稳性对气温插值的影响规律,采用趋势线法对气温数据进行空间平稳性探索,并对比分析不同空间平稳性条件下,普通线性回归、普通克里格、回归克里格的气温插值精度及插值结果的空间分布特点。结果显示：冬季日均、月均气温与年均气温呈现空间非平稳,插值精度随时间序列的增长而提高,随着气温数据逐渐趋于稳定,精度提高的幅度逐渐下降;夏季日均、月均气温呈现空间平稳,随时间序列的增长,插值精度的提高并不显著;夏季日均气温各插值方法的插值精度普遍高于冬季日均气温。与普通克里格相比,回归克里格能有效提高空间非平稳数据的插值精度。时间序列的增长削弱了不同插值算法之间的插值精度差异和插值结果空间分布差异。     

A Point-Based Intelligent Approach to Areal Interpolation

Areal interpolation is the data transfer from one zonal system to another. A survey of previous literature on this subject points out that the most effective methods for areal interpolation are the intelligent approaches, which often take two-dimensional (2-D) land use or one-dimensional (1-D) road network information as ancillary data to give insight on the underlying distribution of a variable. However, the 2-D or 1-D ancillary information is not always applicable for the variable of interest in a specific study area. This article introduces a point-based intelligent approach to the areal interpolation problem by using zero-dimensional (0-D) points as ancillary data that are locationally associated with the variable of interest. The connection between zonal variables and point locations can be modeled with a linear or a nonlinear exponential function, which incorporates the distribution of the variables in the transferring of the information from the source zone to the target zone. An experimental study interpolating the population data at a suburbanized area suggests that the proposed method is an attractive alternative to other areal interpolation solutions based on the evaluation of its resulting accuracy and efficiency.     

Fine-resolution population mapping using OpenStreetMap points-of-interest

Data on population at building level is required for various purposes. However, to protect privacy, government population data is aggregated. Population estimates at finer scales can be obtained through areal interpolation, a process where data from a first spatial unit system is transferred to another system. Areal interpolation can be conducted with ancillary data that guide the redistribution of population. For population estimation at the building level, common ancillary data include three-dimensional data on buildings, obtained through costly processes such as LiDAR. Meanwhile, volunteered geographic information (VGI) is emerging as a new category of data and is already used for purposes related to urban management. The objective of this paper is to present an alternative approach for building level areal interpolation that uses VGI as ancillary data. The proposed method integrates existing interpolation techniques, i.e., multi-class dasymetric mapping and interpolation by surface volume integration; data on building footprints and points-of-interest (POIs) extracted from OpenStreetMap (OSM) are used to refine population estimates at building level. A case study was conducted for the city of Hamburg and the results were compared using different types of POIs. The results suggest that VGI can be used to accurately estimate population distribution, but that further research is needed to understand how POIs can reveal population distribution patterns.     

Stochastic model-based methods for handling uncertainty in areal interpolation

Handling of uncertainty in the estimation of values from source areas to target areas poses a challenge in areal interpolation research. Stochastic model-based methods offer a basis for incorporating such uncertainty, but to date they have not been widely adopted by the GIS community. In this article, we propose one use of such methods based in the problem of interpolating count data from a source set of zones (parishes) to a more widely used target zone geography (postcode sectors). The model developed also uses ancillary statistical count data for a third set of areas nested within both source and target zones. The interpolation procedure was implemented within a Bayesian statistical framework using Markov chain Monte Carlo methods, enabling us to take account of all sources of uncertainty included in the model. Distributions of estimated values at the target zone level are presented using both summary statistics and as individual realisations selected to illustrate the degree of uncertainty in the interpolation results. We aim to describe the use of such stochastic approaches in an accessible way and to highlight the need for quantifying estimation uncertainty arising in areal interpolation, especially given the implications arising when interpolated values are used in subsequent analyses of relationships.     

Evaluating geo-located Twitter data as a control layer for areal interpolation of population

Control data are critical for improving areal interpolation results. Remotely sensed imagery, road network, and parcels are the three most commonly used ancillary data for areal interpolation of population. Meanwhile, the open access geographic data generated by social networks is emerging as an alternative control data that can be related to the distribution of population. This study evaluates the effectiveness of geo-located night-time tweets data as ancillary information and its combination with the three commonly used ancillary datasets in intelligent areal interpolation. Due to the skewed Twitter user age, the other purpose of this study is to test the effect of age bias control data on estimation of different age group populations. Results suggest that geo-located tweets as single control data does not perform as well as the three other control layers for total population and all age-specific population groups. However, the noticeable enhancement effect of Twitter data on other control data, especially for age groups with a high percentage of Twitter users, suggests that it helps to better reflect population distribution by increasing variation in densities within a residential area delineated by other control data.     

Enhancing areal interpolation frameworks through dasymetric refinement to create consistent population estimates across censuses

To assess micro-scale population dynamics effectively, demographic variables should be available over temporally consistent small area units. However, fine-resolution census boundaries often change between survey years. This research advances areal interpolation methods with dasymetric refinement to create accurate consistent population estimates in 1990 and 2000 (source zones) within tract boundaries of the 2010 census (target zones) for five demographically distinct counties in the US. Three levels of dasymetric refinement of source and target zones are evaluated. First, residential parcels are used as a binary ancillary variable prior to regular areal interpolation methods. Second, Expectation Maximization (EM) and its data-extended version leverage housing types of residential parcels as a related ancillary variable. Finally, a third refinement strategy to mitigate the overestimation effect of large residential parcels in rural areas uses road buffers and developed land cover classes. Results suggest the effectiveness of all three levels of dasymetric refinement in reducing estimation errors. They provide a first insight into the potential accuracy improvement achievable in varying geographic and demographic settings but also through the combination of different refinement strategies in parts of a study area. Such improved consistent population estimates are the basis for advanced spatio-temporal demographic research.     

人口统计数据的空间转换

在经济和社会研究中,所要研究的区域之上经常没有数据,而这些数据需要由已知区域的数据求得,即统计数据需要空间转换,这就通常涉及到面积内插。本文从GIS的角度研究如何解决人口内插问题,认为面积内插和GIS中的叠加分析是一致的。在传统的面积内插方法的基础上是提出了基于人口真实分布的面积内插方法,并推导出了公式。同时提出了人口密度的递归算法,即把居住区分为人口稀疏地区和人口稠密地区,估计出人口稀疏地区的人口密度,就可以求出人口密集地区的人口密度;再把人口密集区分为新的人口稀疏区和密集区,此过程反复直至求出接近于人口真实分布的人口模型。     

Interpolating U.S. Decennial Census Tract Data from as Early as 1970 to 2010: A Longitudinal Tract Database

Differences in the reporting units of data from diverse sources and changes in units over time are common obstacles to analysis of areal data. We compare common approaches to this problem in the context of changes over time in the boundaries of U.S. census tracts. In every decennial census, many tracts are split, consolidated, or changed in other ways from the previous boundaries to reflect population growth or decline. We examine two interpolation methods to create a bridge between years, one that relies only on areal weighting and another that also introduces population weights. Results demonstrate that these approaches produce substantially different estimates for variables that involve population counts, but they have a high degree of convergence for variables defined as rates or averages. Finally, the article describes the Longitudinal Tract Database (LTDB), through which we are making available public-use tools to implement these methods to create estimates within 2010 tract boundaries for any tract-level data (from the census or other sources) that are available for prior years as early as 1970.     

Disaggregating population density of the European Union with CORINE land cover

This article describes and compares six disaggregation methods used to produce a dasymetric population density grid of the European Union at a 100 m resolution. Population data were initially available at commune level. The main ancillary information source was CORINE land cover, a land cover map distributed by the European Environment Agency. Information from the Eurostat point survey, land use/cover area frame survey, was also integrated in the parameter estimation of some of the approaches tested. Accurate population data for 1 km cell grids were provided by the Statistical Offices of Austria, Denmark, Finland, the Netherlands, Northern Ireland, Estonia and Sweden. These data provided the basic reference to quantify the accuracy of each method. The best results were obtained with a modified version of the limiting variable method (Eicher, C. and Brewer, C., 2001. Dasymetric mapping and areal interpolation: implementation and evaluation. Cartography and Geographic Information Science, 28, 125–138) that could be implemented, thanks to the national reference grids. For other methods the parameters could be estimated without using the reference grids; among them a method based on logit regression gave the best results. Compared with the traditional choropleth maps that represent a homogeneous density in each commune, the accuracy improvement of the disaggregated maps ranged between 20% and 67% (between 46% and 67% for the best method).     

薄盘光滑样条插值中三种协变量方法的降水量插值精度比较

在薄盘光滑样条插值中,高相关协变量的选取决定了插值结果的精确性。以2001-2009年全国728个气象站点日降水为数据源,提取年降水量数据,在分析多年平均降水量与两协变量高程（DEM）和距海岸线距离（DCL）的空间相关性基础上,利用ANUSPLIN软件,比较不同协变量下降水量插值结果精度在全国尺度以及区域尺度上的差异。以DEM、DCL及DEM-DCL分别为协变量对降水量数据进行空间插值发现：①在全国尺度上,DEM法的平均绝对误差（MAE）为47.79,略低于DEM-DCL法（48.90）,但显著低于DCL法（55.54）;且DEM法的平均相对误差和均方根误差也明显低于其它两种方法。②在区域尺度上,除西藏地区外的其他7个区域,3种方法的插值误差与全国尺度上相一致。西藏地区降水插值结果以DCL法的精度最高,而DEM法则较差。研究建议除在西藏地区的降水量插值研究中采用DCL法,在全国其他大部分区域采用DEM法。