基于空间变异修正普通克里金的矿区复杂地形高精度建模

普通克里金法是构建矿区三维地形模型并揭示地表沉陷形变场动态变化规律的有效方法。但普通克里金法存在平滑效应的问题,导致估计值的空间变异程度小于实际,无法真实反映矿区复杂地表形态的空间变化特征。本文结合矿区实际地形采样数据,提出了一种空间变异修正的普通克里金法,并将其应用于矿区复杂地形的高精度建模,并与已有方法在建模精度和空间变异性复现方面进行了对比分析。研究结果表明,空间变异修正的普通克里金法能够很好地处理平滑效应对建模的影响,具有更高的建模精度和空间变异性复现能力,在矿区复杂地形高精度建模应用中具有较强的适用性,可以作为一种可靠的建模工具用于矿区复杂地形沉陷形变场动态变化规律分析。     

On a stronger-than-best property for best prediction

The minimum mean squared error (MMSE) criterion is a popular criterion for devising best predictors. In case of linear predictors, it has the advantage that no further distributional assumptions need to be made, other then about the first- and second-order moments. In the spatial and Earth sciences, it is the best linear unbiased predictor (BLUP) that is used most often. Despite the fact that in this case only the first- and second-order moments need to be known, one often still makes statements about the complete distribution, in particular when statistical testing is involved. For such cases, one can do better than the BLUP, as shown in Teunissen (J Geod. doi: 10.1007/s00190-007-0140-6, 2006), and thus devise predictors that have a smaller MMSE than the BLUP. Hence, these predictors are to be preferred over the BLUP, if one really values the MMSE-criterion. In the present contribution, we will show, however, that the BLUP has another optimality property than the MMSE-property, provided that the distribution is Gaussian. It will be shown that in the Gaussian case, the prediction error of the BLUP has the highest possible probability of all linear unbiased predictors of being bounded in the weighted squared norm sense. This is a stronger property than the often advertised MMSE-property of the BLUP.     

New optimal smoothing scheme for improving relative and absolute accuracy of tightly coupled GNSS/SINS integration

For mobile surveying and mapping applications, tightly coupled integration of global navigation satellite system (GNSS) and Strap down Inertial Navigation System is usually recommended for direct georeferencing since it can provide position, velocity, and attitude information at higher accuracy and better reliability in a self-contained manner. A post-mission smoothing method is applied to optimally use observation information of both systems and to overcome the shortcomings of Kalman filter in GNSS degraded environments. We propose the revised Rauch–Tung–Streibel Smoother (RTSS) and Forward–Backward combination (FBC) smoothing algorithms for tightly coupled integration. From the analysis and field test, it is found that RTSS smoothing mainly improves the relative accuracy, while FBC mainly contributes to the absolute accuracy. With the complementary characteristics of both smoothing algorithms, an optimal new smoothing scheme combining RTSS with FBC is built. The performance of these three smoothing algorithms is evaluated through a real vehicular test. Compared with RTSS and FBC smoothing algorithms, the new smoothing scheme improves the mean 3D position RMS and the mean 3D attitude RMS by 65.7 and 70%, respectively. It provides better accuracy and smoothness for the position, velocity, and attitude at the same time.     

载波相位平滑时间常数对GBAS精度增强的影响

陆基增强系统(GBAS)是利用载波相位平滑伪距差分修正实现对导航辅助定位的. 其中,平滑时间常数是影响载波相位平滑伪距精度的关键参数. 本文分析研究了不同平滑时间下电离层时间梯度和空间梯度对Hatch滤波的影响. 在结合电离层时空梯度和多径效应引起的滤波总误差方差的基础上,推导出自适应的最优平滑时间常数. 分别对GBAS静态和动态两种环境下的定位误差进行实验,实验结果表明,采用本文推导出的自适应平滑时间常数降低了GBAS伪距测量误差,从而使定位精度得到增强.      

Implication of adaptive smoothness constraint and Helmert variance component estimation in seismic slip inversion

When ill-posed problems are inverted, the regularization process is equivalent to adding constraint equations or prior information from a Bayesian perspective. The veracity of the constraints (or the regularization matrix R) significantly affects the solution, and a smoothness constraint is usually added in seismic slip inversions. In this paper, an adaptive smoothness constraint (ASC) based on the classic Laplacian smoothness constraint (LSC) is proposed. The ASC not only improves the smoothness constraint, but also helps constrain the slip direction. A series of experiments are conducted in which different magnitudes of noise are imposed and different densities of observation are assumed, and the results indicated that the ASC was superior to the LSC. Using the proposed ASC, the Helmert variance component estimation method is highlighted as the best for selecting the regularization parameter compared with other methods, such as generalized cross-validation or the mean squared error criterion method. The ASC may also benefit other ill-posed problems in which a smoothness constraint is required.     

Gravity field approximation in the area of Greece with emphasis on local characteristics

The accuracy of the gravity field approximation depends on the amount of the available data and their distribution as well as on the variation of the gravity field. The variation of the gravity field in the Greek mainland, which is the test area in this study, is very high (the variance of point free air gravity anomalies is 3191.5mgal ²). Among well known reductions used to smooth the gravity field, the complete isostatic reduction causes the best possible smoothing, however remain strong local anomalies which disturb the homogeneity of the gravity field in this area. The prediction of free air gravity anomalies using least squares collocation and regional covariance function is obtained within a ±4 ... ±19mgal accuracy depending on the local peculiarities of the free air gravity field. By taking into account the topography and its isostatic compensation with the usual remove-restore technique, the accuracy of the prediction mentioned obove was increased by about a factor of 4 and the prediction results become quite insensitive to the covariance function used (local or regional). But when predicting geoidal heights, in spite of using the smoothed field, the prediction results remain still depend on the covariance function used in such a way that differences up to about 50cm/100km result between relative geoidal heights computed with regional or local covariance functions.     

Simultaneous estimation of cross-validation errors in least squares collocation applied for statistical testing and evaluation of the noise variance components

The cross-validation technique is a popular method to assess and improve the quality of prediction by least squares collocation (LSC). We present a formula for direct estimation of the vector of cross-validation errors (CVEs) in LSC which is much faster than element-wise CVE computation. We show that a quadratic form of CVEs follows Chi-squared distribution. Furthermore, a posteriori noise variance factor is derived by the quadratic form of CVEs. In order to detect blunders in the observations, estimated standardized CVE is proposed as the test statistic which can be applied when noise variances are known or unknown. We use LSC together with the methods proposed in this research for interpolation of crustal subsidence in the northern coast of the Gulf of Mexico. The results show that after detection and removing outliers, the root mean square (RMS) of CVEs and estimated noise standard deviation are reduced about 51 and 59%, respectively. In addition, RMS of LSC prediction error at data points and RMS of estimated noise of observations are decreased by 39 and 67%, respectively. However, RMS of LSC prediction error on a regular grid of interpolation points covering the area is only reduced about 4% which is a consequence of sparse distribution of data points for this case study. The influence of gross errors on LSC prediction results is also investigated by lower cutoff CVEs. It is indicated that after elimination of outliers, RMS of this type of errors is also reduced by 19.5% for a 5 km radius of vicinity. We propose a method using standardized CVEs for classification of dataset into three groups with presumed different noise variances. The noise variance components for each of the groups are estimated using restricted maximum-likelihood method via Fisher scoring technique. Finally, LSC assessment measures were computed for the estimated heterogeneous noise variance model and compared with those of the homogeneous model. The advantage of the proposed method is the reduction in estimated noise levels for those groups with the fewer number of noisy data points.     

基于超分辨率重建的多时相MODIS与Landsat反射率融合方法

提出一种基于超分辨率重建的MODIS与Landsat反射率图像融合方法,以STARFM算法与超分辨率重建为基础,使用观测的MODIS和Landsat地表反射率图像预测给定时刻的Landsat合成反射率图像。该方法利用基于稀疏表示的超分辨率重建方法对MODIS图像进行分辨率增强,实验结果表明这一操作能够增加原MODIS图像的空间细节,有助于提高STARFM算法的预测精度;另一方面,考虑输入两个基时刻图像相差较大时原STARFM算法预测的反射率会存在"时间平滑"的问题,限制每次只使用一个基时刻MODIS和Landsat图像对进行STARFM预测,使用逐图像块选择策略,从由两个基时刻图像分别进行预测得到的两组预测图像中选择最优的预测,同样得到了优于STARFM算法的预测结果。     

基于高斯样条函数的局部重力异常场解析重构

重力匹配辅助导航理论大都建立在离散场的基础上的,为了研究基于连续场重力匹配算法以克服传统匹配算法的局限,必须建立精度高且具有良好解析性质的局部重力异常场解析模型。利用斐波那契数列寻优方法对一维高斯样条函数插值进行最优化,在此基础上提出了基于斐波那契数列寻优的二维高斯样条函数逼近局部重力异常场方法。为了提高寻优算法运算速度,将二维准则函数解耦为X方向和Y方向两个独立的一维准则函数,分别采用斐波那契数列寻优方法对这两个准则函数进行寻优以获取X方向和Y方向最优参数,最终得到高精度逼近局部离散格网数据的局部重力异常场连续解析模型。仿真实验中采用五组不同的参数对变化范围为-51.185mGal~86.1819mGal的重力异常场进行逼近。从最后的仿真实验结果可以看出采用最优参数时逼近绝对误差均值达到0.00069,相对误差均值更达到10-6级,能较好的满足了匹配导航要求,其逼近精度较采用其它非最优参数时均有较大提高,由此验证了文中提出的重构算法有效性。     

顾及卫星钟随机特性的抗差最小二乘配置钟差预报算法

为了更好地反映钟差特性并提高其预报精度,采用抗差最小二乘配置方法建立一种能够同时考虑星载原子钟物理特性、钟差周期性变化与随机性变化特点的钟差预报模型。首先使用附有周期项的二次多项式模型进行拟合提取卫星钟差的趋势项与周期项,然后针对剩余的随机项及其可能存在的粗差,采用抗差最小二乘配置的原理进行建模,其中最小二乘配置的协方差函数通过对比协方差拟合的方法并结合试验进行确定。使用IGS精密钟差数据进行预报试验,将本文方法与二次多项式模型、灰色模型进行对比,预报精度分别提高了0.457 ns和0.948 ns,而预报稳定性则分别提高了0.445 ns和1.233 ns,证明了本文方法能够更好地预报卫星钟差,同时说明本文的协方差函数确定方法的有效性。     

Non-stationary covariance function modelling in 2D least-squares collocation

Standard least-squares collocation (LSC) assumes 2D stationarity and 3D isotropy, and relies on a covariance function to account for spatial dependence in the observed data. However, the assumption that the spatial dependence is constant throughout the region of interest may sometimes be violated. Assuming a stationary covariance structure can result in over-smoothing of, e.g., the gravity field in mountains and under-smoothing in great plains. We introduce the kernel convolution method from spatial statistics for non-stationary covariance structures, and demonstrate its advantage for dealing with non-stationarity in geodetic data. We then compared stationary and non- stationary covariance functions in 2D LSC to the empirical example of gravity anomaly interpolation near the Darling Fault, Western Australia, where the field is anisotropic and non-stationary. The results with non-stationary covariance functions are better than standard LSC in terms of formal errors and cross-validation against data not used in the interpolation, demonstrating that the use of non-stationary covariance functions can improve upon standard (stationary) LSC.     

Noisy Data Smoothing in DEM Construction Using Least Squares Support Vector Machines

Since spatial datasets are subject to sampling errors, a smoothing interpolation method should be employed to remove noise during DEM construction. Although least squares support vector machines (LSSVM) have been widely accepted as a classifier, their effect on smoothing noisy data is almost unknown. In this article, the smoothness of LSSVM was explored, and its effect on smoothing noisy data in DEM construction was tested. In order to improve the ability to deal with large datasets, a local method of LSSVM has been developed, where only the neighboring sampling points around the one to be estimated are used for computation. A numerical test indicated that LSSVM is more accurate than the classical smoothing methods including TPS and kriging, and its error surfaces are more evenly distributed. The real‐world example of smoothing noise inherent in lidar‐derived DEMs also showed that LSSVM has a positive smoothing effect, which is approximately as accurate as TPS. In short, LSSVM with a high efficiency can be considered as an alternative smoothing method for smoothing noisy data in DEM construction.     

残差自适应回归MGM(1,n)

鉴于MGM(1,n)在预测过程中,不同变量的拟合及预测残差差距较大,本文将自适应回归模型引入到不同变量的残差估计中,对整体的残差起到平滑的作用,从而抑制了残差的上扬趋势,经过实例分析,残差自适应回归MGM(1,n)的预测精度得到了明显的提高。     

GIDS空间插值法估算云下地表温度

选用陆面区域温度最佳空间插值法—梯度距离平方反比法(GIDS),为近似估算云下地表温度提供了可能。实验选取暖季南京江宁地区ETM+影像和ASTERGDEMV1高程数据,探索分析GIDS估算云下地表温度的可行性和可信性。对14种空间大小云覆盖区实验研究表明:利用GIDS插值估算云下地表温度具有可行性,且估算误差随着云覆盖区范围增大而增加,其最大MAE<0.9℃,最大RMSE<1.2℃,并在云覆盖区小于100×100像元时,最大MAE<0.8℃、RMSE<1℃;插值精度与最近邻无云像元典型代表性、区域内空间复杂度和地表覆盖类型均有关,存在不稳定性和动态性;云下NDVI均方差与MAE、RMSE有着一致变化趋势,借助NDVI均方差指示云下地表空间异质性及NDVI–LST负相关性,可对插值结果进行可信性评判,以避免插值结果盲目应用,推进和提升地表温度产品应用价值。     

灰色动态组合模型及其在大坝变形预测中的应用

提出了先对序列进行平滑处理,再建立灰色动态组合模型,分析了实测数据的计算结果,实现了较高的预测精度。     

地形因子与小微地震相关度的空间异质性分析

针对传统最小二乘回归未能顾及数据的空间特性,且无法度量模型自变量与因变量相关性的空间变异特性的问题,本文提出利用地理加权回归方法分析小微地震频次与地形因子相关度的空间异质性。以四川地区的地震监测资料、DEM为实验数据,选取地形复杂度、坡度变率、坡向变率和地面曲率为自变量,地震发生频次为因变量,构建地理加权回归模型,并进行回归系数的空间变异分析。实验分析发现,地震频次与地形因子具有一定的相关性:地形复杂度与地震频次相关性最强;坡度变率、沟壑密度、剖面曲率与地震频次的相关性依次减弱;不同空间位置的地形因子和地震频次的相关性具有较明显的空间异质性。实验结果表明,地理加权回归可以有效地度量分析地震频次与地形因子相关度的空间异质性,研究结果可为地震及次生灾害的分析与预报提供辅助决策参考。     

混沌理论在桥梁变形监测中的应用

针对桥梁的非线性下沉问题,引用了混沌理论,首先求取时间序列的两重构参数时间延迟τ和嵌入维数m进行相空间重构;随后进行混沌特性判别,确定该时间序列存在混沌迹象;最后根据所求参数建立加权零阶局域预计模型和RBF神经网络混沌预计模型对观测数据进行预计分析,并与系数为0.9的指数平滑预测模型进行比较,结果显示混沌预计模型值更接...     

Predicting soil organic carbon stocks under commercial forest plantations in KwaZulu-Natal province,South Africa using remotely sensed data

ABSTRACT

Commercial forest plantations are increasing globally, absorbing a large amount of carbon valuable for climate change mitigation. Whereas most carbon assimilation studies have mainly focused on natural forests, understanding the spatial distribution of carbon in commercial forests is central to determining their role in the global carbon cycle. Forest soils are the largest carbon reservoir; hence soils under commercial forests could store a significant amount of carbon. However, the variability of soil organic carbon (SOC) within forest landscapes is still poorly understood. Due to limitations encountered in traditional systems of SOC determination, especially at large spatial extents, remote sensing approaches have recently emerged as a suitable option in mapping soil characteristics. Therefore, this study aimed at predicting soil organic carbon (SOC) stocks in commercial forests using Landsat 8 data. Eighty-one soil samples were processed for SOC concentration and fifteen Landsat 8 derived variables, including vegetation indices and bands were used as predictors to SOC variability. The random forest (RF) was adopted for variable selection and regression method for SOC prediction. Variable selection was done using RF backward elimination to derive three best subset predictors and improve prediction accuracy. These variables were then used to build the RF final model for SOC prediction. The RF model yielded good accuracies with root mean square error of prediction (RMSE) of 0.704 t/ha (16.50% of measured mean SOC) and 10-fold cross-validation of 0.729 t/ha (17.09% of measured mean SOC). The results demonstrate the effectiveness of Landsat 8 bands and derived vegetation indices and RF algorithm in predicting SOC stocks in commercial forests. This study provides an effective framework for local, national or global carbon accounting as well as helps forest managers constantly evaluate the status of SOC in commercial forest compartments.     

条件随机场模型约束下的遥感影像模糊C-均值聚类算法

遥感影像具有丰富的空间相关信息,而传统的基于像元光谱的聚类算法并不能将空间信息融入聚类,聚类结果往往不好。针对这一问题,本文提出了一种条件随机场模型约束下的模糊C-均值聚类算法,通过邻域像元的分类先验信息对中心像元的类别进行约束从而提取空间相关信息,基于二阶条件随机场将光谱信息和空间相关信息同时融入聚类,并使用环形置信度迭代算法得到像元分类后验概率的全局最优推测。试验证明,本文算法能够有效地保持地物的形状特征,分类精度相比传统算法有所提高。