A multiresolution hierarchical classification algorithm for filtering airborne LiDAR data

We presented a multiresolution hierarchical classification (MHC) algorithm for differentiating ground from non-ground LiDAR point cloud based on point residuals from the interpolated raster surface. MHC includes three levels of hierarchy, with the simultaneous increase of cell resolution and residual threshold from the low to the high level of the hierarchy. At each level, the surface is iteratively interpolated towards the ground using thin plate spline (TPS) until no ground points are classified, and the classified ground points are used to update the surface in the next iteration. 15 groups of benchmark dataset, provided by the International Society for Photogrammetry and Remote Sensing (ISPRS) commission, were used to compare the performance of MHC with those of the 17 other publicized filtering methods. Results indicated that MHC with the average total error and average Cohen’s kappa coefficient of 4.11% and 86.27% performs better than all other filtering methods.     

空间结构函数在北京地区气象观测站网设计中的应用

通过统计1978—2000年北京东南低地形区有关台站在春季、夏季、秋季、冬季的逐日平均气温和水汽压的结构函数, 分季节分析了该地区这两个二类气象要素的线段及平面内插精度和台站间距的对应关系, 并根据内插标准误差不超过观测标准误差的原则, 对上述两要素在北京东南低地形区的合理布站方案及间距进行了估算, 可以为2008年北京奥运会气象服务系统建设中气象台站布网建设提供一定的依据。结果表明, 正三角形排列方案为北京东南地区二类气象台站的最佳布站方案, 且布站精度应小于等于16 km。     

引入Himawari-8卫星数据协变量的能见度样条插值方法

在ANUSPLIN薄盘光滑样条插值中,高相关协变量的选取决定了插值结果的精确性。本文选取2017—2019年大雾和霾能见度较差的天气过程,利用183个能见度观测站点对能见度进行插值,引入Himawari-8卫星的通道数据和DEM数据作为协变量对能见度的插值结果进行改进,并对能见度插值结果进行对比分析。研究表明,引入Himawari-8数据和DEM数据作为协变量的能见度插值结果在精度上有显著提高,尤其对雾区和霾区的边界范围和纹理的反演更为准确,基于Himawari-8卫星数据和气象监测站点的观测数据,使用协变量的方法进行能见度插值可以做为能见度监测网格化的一种有效途径。     

基于二维三次卷积插值算法的辛几何射线追踪

射线追踪是地震波走时层析成像的基础,射线空间位置的准确性及射线走时的精度决定了层析成像的可靠性.本文根据哈密尔顿系统可以有效提高程函方程解稳定性的特性,采用辛几何算法（SAM-Symplectic Algorithm Method）及二维三次卷积插值技术进行地震波射线追踪.由于采用了SAM算法,保证了地震波波前精度,提高了射线空间位置的准确性.数值模拟结果表明SAM既能保证哈密尔顿系统的稳定性又具有运算速度快的特点,提高了射线追踪的计算精度.     

Spatio-temporal interpolation of precipitation during monsoon periods in Pakistan

Spatio-temporal estimation of precipitation over a region is essential to the modeling of hydrologic processes for water resources management. The changes of magnitude and space–time heterogeneity of rainfall observations make space–time estimation of precipitation a challenging task. In this paper we propose a Box–Cox transformed hierarchical Bayesian multivariate spatio-temporal interpolation method for the skewed response variable. The proposed method is applied to estimate space–time monthly precipitation in the monsoon periods during 1974–2000, and 27-year monthly average precipitation data are obtained from 51 stations in Pakistan. The results of transformed hierarchical Bayesian multivariate spatio-temporal interpolation are compared to those of non-transformed hierarchical Bayesian interpolation by using cross-validation. The software developed by [11] is used for Bayesian non-stationary multivariate space–time interpolation. It is observed that the transformed hierarchical Bayesian method provides more accuracy than the non-transformed hierarchical Bayesian method.     

Observations of submesoscale eddies using high-frequency radar-derived kinematic and dynamic quantities

The spatio-temporal variability of submesoscale eddies off southern San Diego is investigated with two-year observations of subinertial surface currents [O(1) m depth] derived from shore-based high-frequency radars. The kinematic and dynamic quantities — velocity potential, stream function, divergence, vorticity, and deformation rates — are directly estimated from radial velocity maps using optimal interpolation. For eddy detection, the winding-angle approach based on flow geometry is applied to the calculated stream function. A cluster of nearly enclosed streamlines with persistent vorticity in time is identified as an eddy. About 700 eddies were detected for each rotation (clockwise and counter-clockwise). The two rotations show similar statistics with diameters in the range of 5–25 km and Rossby number of 0.2–2. They persist for 1–7 days with weak seasonality and migrate with a translation speed of 4–15 cm s⁻¹ advected by background currents. The horizontal structure of eddies exhibits nearly symmetric tangential velocity with a maximum at the defined radius of the eddy, non-zero radial velocity due to background flows, and Gaussian vorticity with the highest value at the center. In contrast divergence has no consistent spatial shape. Two episodic events are presented with other in situ data (subsurface current and temperature profiles, and local winds) as an example of frontal-scale secondary circulation associated with drifting submesoscale eddies.     

The effect of interpolation methods in surface definition: an experimental study

Mappings of the earth surface and their representation in 3D (three‐dimensional) models are commonly used in most recent research. Modeling research, which starts with classical surveying methods, acquires new dimensions matching the modern technologies. 3D models of any object or earth surface can be used in much visual and scientific research. A digital model of the landscape is an important part within creation of geo‐information systems used in the public administration and in the commercial sphere. It is an important tool in applications such as geomorphology, hydrology, geology, cartography, ecology, mining etc. Values of volume in terrains that do not have regular geometric structure can be obtained more accurately by using 3D models of surfaces with respect to developing technology. Basic data of 3D models must indicate 3D coordinates of the surveyed object in the reference frame. Distribution and intensity of points are important factors in modeling earth surfaces. A minimum number of points is desired in defining an object in 3D. Interpolation methods employing different mathematical models are used to obtain 3D models of terrain surfaces. In this study, the effect of interpolation methods in defining a terrain surface is investigated. For this purpose, a uniform surface, hill‐shaped artificial object with a known volume is employed. The 3D surface and volume are calculated by using 12 different interpolation methods. Point distribution, point intensity and accuracy of point measurements are not considered. The same data set was used for all the interpolation methods. The interpolation methods are compared and evaluated based on the results. Copyright © 2007 John Wiley & Sons, Ltd.     

基于多参量的模拟退火全局优化傅里叶有限差分算子

傅里叶有限差分法(FFD)能够处理复杂地质构造中的波传播问题,但对陡倾角成像仍有明显的误差.优化参数的方法能够在保持计算效率的前提下进一步提高陡倾角的成像精度.本文在有理近似的基础上,将FFD算子展开式中的常系数由两个拓展为四个,然后采用模拟退火算法对这四个参数进行全局优化.本方法除了考虑速度对比度以外,还考虑了频率和延拓步长等参量的影响.理论误差分析和脉冲响应测试均表明该方法能极大地提高FFD算子的精确传播角度.二维SEG/EAGE盐丘模型实验表明本文方法对陡倾角以及盐下构造的成像精度明显高于未优化的FFD法.将本文的方法与交替方向加插值的方法结合应用于三维脉冲响应测试更进一步证实了本文方法的有效性.     

Missing precipitation data estimation using optimal proximity metric-based imputation,nearest-neighbour classification and cluster-based interpolation methods

Abstract

New optimal proximity-based imputation, K-nearest neighbour (K-NN) classification and K-means clustering methods are proposed and developed for estimation of missing daily precipitation records. Mathematical programming formulations are developed to optimize the weighting, classification and clustering schemes used in these methods. Ten different binary and real-valued distance metrics are used as proximity measures. Two climatic regions, Kentucky and Florida, (temperate and tropical) in the USA, with different gauge density and network structure, are used as case studies to evaluate the new methods. A comprehensive exercise is undertaken to compare the performances of the new methods with those of several deterministic and stochastic spatial interpolation methods. The results from these comparisons indicate that the proposed methods performed better than existing methods. Use of optimal proximity metrics as weights, spatial clustering of observation sites and classification of precipitation data resulted in improvement of missing data estimates.

Editor D. Koutsoyiannis; Associate editor C. Onof     

基于回声探测监测沉水植物盖度的插值方法研究:以杭州西湖沉水植物盖度监测为例

沉水植物作为水生态系统的重要组成成分,在水生态系统物质循环和能量流动中发挥着重要作用,其覆盖度和生物量是评价湖泊等浅水水体系统稳定性的关键参数随着高效和无损伤监测的回声探测仪在沉水植物盖度监测中的应用,其精确度算法也受到了越来越多的关注本研究以成功恢复沉水植物的浅水湖泊杭州西湖为研究对象,利用BioSonics便携型回声探测仪——MX采集沉水植物回声样本同时结合人工样方设置,采集与回声探测对应位点的沉水植物样本,验证回声探测结果的精确性通过建立回归模型分析回声探测得到的沉水植物体积百分比(PVI)与人工样方获得的对应平均鲜重关系,结果表明二者具有较好的相关性分别采用普通克里金法、反距离权重法、径向基函数法3种插值方法对同一季节的不同湖泊和同一湖泊的不同季节未采集区域沉水植物的盖度数据进行插值分析,并对插值结果进行交叉验证,以确定方法的精确度交叉验证结果表明,插值精确度反距离权重法径向基函数法普通克里金法研究结果为回声探测与插值分析方法结合在大尺度浅水水体中沉水植物监测应用提供了技术支撑.