基于PRISM和泰森多边形的地形要素日降水量空间插值研究

以黑河流域河西走廊中段地区为例,利用该研究区年、月降水与地形间较强的相关性特点,在PRISM方法的基础上对该地区日降水量进行了空间插值计算。文章提出了以月降水量的PRISM空间插值结果为该月逐日降水空间分布的参考本底,利用泰森多边形方法确定空间日降水的概率,从而实现黑河流域河西走廊中段地区日降水的空间制图方法,并对该方法得到的日降水时空数据集进行了误差分析和评估。分析结果表明,这一方法简单可靠,满足分布式水文模型或相关陆表过程分布式模拟对分布式日降水数据时空精度的要求。     

人口统计数据空间化的一种方法

人口空间分布信息在环境健康风险诊断、自然灾害损失评估和现场抽样调查比较等地理学和相关学科研究中占有重要的地位。目前随着对地观测技术和地理信息科学的飞速发展, 如何精确地进行人口数据空间化成为了研究的难点和热点。针对采用传统方法解决人口空间化问题所遇到的困难和不足, 设计了遗传规划(genetic programming, GP)、遗传算法(genetic algorithms, GA) 和GIS 相结合的方法, 以GIS 确定量化影响因子权重, 以GP 建立模型结构, 以GA 优化模型参数, 成功建立研究区-山西省和顺县的人口数据格网分布表面。实验证明与传统建模方法(如逐步回归分析模型和重力模型)相比, 所提方法建模过程更为智能化与自动化, 模型结构更为灵活多样, 而且数据拟合精度更高。     

抗差估计的多波束测深数据内插方法

针对传统多波束测深数据网格化内插方法抗差性不足的问题,该文将抗差估计理论应用于网格化内插方法中。基于不确定度指标建立了距离不确定度联合加权内插模型,在该模型的权函数中引入了等价权,同时给出了等价权设计方案及迭代初值确定原则,以增强模型的抗差性;介绍了一种节点邻域参考点合理选取方法,给出了改进内插法实现的基本步骤,最后通过实测数据对反距离加权法、联合加权法及改进方法内插水深质量进行比较。结果表明,改进方法内插结果要优于另外两种方法,可应用于多波束测深数据网格化处理中。     

普通Kriging插值算法研究

地统计学是进行空间变异分析的有力工具.文章说明了地统计学的基本原理及其实现算法,并应用该算法对研究区域的草句白浆土有效磷的空间变异特征进行了分析,通过计算得出的地统计学参数说明了土壤有效磷具有中等空间相关性,并生成了该养分的空间分布图.     

Integrating ASTER and GLASS broadband emissivity products using a multi-resolution Kalman filter

In this study, the multi-resolution Kalman filter (MKF) algorithm, which can handle multi-resolution problems with high computational efficiency, was used to blend two emissivity products: the Global LAnd Surface Satellite (GLASS) (BBE) product and the Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) narrowband emissivity (NBE) product. The ASTER NBE product was first converted into a BBE product. A new detrending method was used to transfer the BBEs into a process suitable for the MKF. The new detrending method was superior to the two existing methods. Finally, both the de-trended GLASS and ASTER BBE products were incorporated into the MKF framework to obtain the optimal estimation at each scale. Field measurements collected in North America were used to validate the integrated BBEs. Visually, the fusion map showed good continuity, with the exception of the border areas, and the quality of the fusion map was better than that of the original maps. The validation results indicate that the MKF improved the BBE product accuracy at the coarse scale. In addition, the MKF was capable of recovering missing pixels at a finer scale.     

基于GIS的坡面泥石流危险性评价

综合不同的泥石流评价因子对泥石流危险性进行了研究,利用地理信息系统（ GIS）来计算潜在的泥石流危险区域。选取虹口乡高原村与庙坝村共处坡面为研究区域,在ArcGIS中生成研究区域的坡面立体图。选取适当的评级指标,确定了流域内主要的泥石流影响因子。采用层次分析法得出各主要因子的权重值,将各因子划分为五个不同等级,并分别对各个因子进行数值化和归一化处理。利用ArcGIS形成各评价因子专题层,再分别给不同评价栅格赋权重值,最终运用叠置分析与插值分析得出研究坡面的不同色带所表示的危险等级区划图,从而实现单坡面泥石流危险性的评价。     

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

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

MATHEMATICAL GEOLOGY

正20141850 Chen Dongyue(School of Earth Sciences and Resources,China University of Geosciences,Beijing 100083,China);Chen Jianping On 3D Ore Prospecting Modeling of Comprehensive Information for Huangshaping Polymetallic Deposit(Journal of Geology,ISSN1674-3636,CN32-1796/P,37(3),2013,p.489-495,12 illus.,12 refs.) Key words:polymetallic ores,data bases,Hunan Province     

When do we need a trend model in kriging?

Under usual estimation practice with local search windows for data and for interpolation situations, universal kriging and ordinary kriging yield the same estimates, using a data set with apparent trend, for both the unknown attribute and its trend component. Modeling the trend matters only in extrapolation situations. Because conditions of the case study presented arise most frequently in practice, the simpler ordinary kriging is the preferred option.     

Foliation fields and 3D cartography in geology: Principles of a method based on potential interpolation

A modeling method that takes into account known points on a geological interface and plane orientation data such as stratification or foliation planes is described and tested. The orientations data do not necessarily belong to one of the interfaces but are assumed to sample the main anisotropy of a geological formation as in current geological situations. The problem is to determine the surfaces which pass through the known points on interfaces and which are compatible with the orientation data. The method is based on the interpolation of a scalar field defined in the space the gradient in which is orthogonal to the orientations, given that some points have the same but unknown scalar value (points of the same interface), and that scalar gradient is known on the other points (foliations). The modeled interfaces are represented as isovalues of the interpolated field. Preliminary two-dimensional tests carried-out with different covariance models demonstrate the validity of the method, which is easily transposable in three dimensions.