Suitable sampling technique in contextual fuzzy c-means classification of remotely sensed data for land cover mapping

In the past, researchers tried hard classification techniques with contextual information to improve classification results. While modelling the spatial contextual information for hard classifiers using Markov Random Field it has been found that the Metropolis algorithm is easier to program and it performs better when compared with the Gibbs sampler. In this study, it has been found that in the case of soft contextual classification, the Metropolis algorithm fails to sample from a random field efficiently and the Gibbs sampler performs better than the Metropolis algorithm, due to the high dimensionality of the soft classification outputs.     

Bayesian methods for outliers detection in GNSS time series

This article is concerned with the problem of detecting outliers in GNSS time series based on Bayesian statistical theory. Firstly, a new model is proposed to simultaneously detect different types of outliers based on the conception of introducing different types of classification variables corresponding to the different types of outliers; the problem of outlier detection is converted into the computation of the corresponding posterior probabilities, and the algorithm for computing the posterior probabilities based on standard Gibbs sampler is designed. Secondly, we analyze the reasons of masking and swamping about detecting patches of additive outliers intensively; an unmasking Bayesian method for detecting additive outlier patches is proposed based on an adaptive Gibbs sampler. Thirdly, the correctness of the theories and methods proposed above is illustrated by simulated data and then by analyzing real GNSS observations, such as cycle slips detection in carrier phase data. Examples illustrate that the Bayesian methods for outliers detection in GNSS time series proposed by this paper are not only capable of detecting isolated outliers but also capable of detecting additive outlier patches. Furthermore, it can be successfully used to process cycle slips in phase data, which solves the problem of small cycle slips.     

处理高杠杆异常值的抗隐差型Bayes方法

给出了一种剔除初始子集中高杠杆异常值的方法。首先根据高杠杆异常值在总观测值集中所占的比例选出若干组观测值,使得至少有一组不含高杠杆异常值的概率很高;然后根据残差最小准则从中选出不含高杠杆异常值的那组作为初始子集;最后用这种初始子集确定方法结合Gibbs抽样给出了相应的Bayes多粗差定位算法。     

Supervised SAR Image MPM Segmentation Based on Region-Based Hierarchical Model

This letter presents a novel method of supervised multiresolution segmentation for synthetic aperture radar images. The method uses a region-based half-tree hierarchical Markov random field model for multiresolution segmentation. To form the region-based multilayer model, the watershed algorithm is employed at each resolution level independently. The nodes of a quadtree in the proposed model are defined as regions instead of pixels. The relationship over scale is studied, and the region-based upward and downward maximization of posterior marginal estimations are deduced. The experimental results for the segmentation of homogeneous areas prove the region-based model much better in terms of robustness to speckle and preservation of edges compared to the pixel-based hierarchical model and the Gibbs sampler with the single-resolution model     

Some issues in contextual fuzzy <Emphasis Type="Italic">c</Emphasis>-means classification of remotely sensed data for land cover mapping

Earlier for the hard classification techniques contextual information was used to improve classification accuracy. While modelling the spatial contextual information for hard classifiers using Markov Random Field it has been found that Metropolis algorithm is easier to program and it performs better in comparison to the Gibbs sampler. In the present study it has been found that incase of soft contextual classification Metropolis algorithm fails to sample from a random field efficiently and from the analysis it was found that Metropolis algorithm is not suitable for soft contextual classification due to the high dimensionality of the soft outputs.     

基于MRF随机场和广义混合模型的遥感图像分级聚类

有限混合模型FM的分级聚类已广泛应用于不同领域,然而,它的计算复杂度与观测数据的平方成正比,因此,在海量数据方面的应用就受到了限制。另一方面,多光谱图像数据中同时包含有空间和光谱两类信息,但大多数基于像素的多光谱图像聚类方法,仅使用了其频谱信息而忽视了空间信息。本文提出了一种新的基于广义有限混合模型GFM的分级聚类方法,该算法把MRF随机场和GFM模型结合在一起,分类数可以通过PLIC准则自动确定。算法在执行过程中,采用K均值聚类方式获得过分类图像,分级聚类从过分类图像开始,代替原来从单点类开始的方式,这样可以方便获取GFM模型成分密度的初始参数。最后,采用由Gibbs采样器生成的仿真测试图对算法的精度进行了定量评价,通过与K均值聚类和FM聚类的比较说明了本文算法的优越性,同时用荷兰Flevoland农业地区的极化SAR图像验证了本文算法的有效性。     

Gibbs sampler for computing and propagating large covariance matrices

The use of sampling-based Monte Carlo methods for the computation and propagation of large covariance matrices in geodetic applications is investigated. In particular, the so-called Gibbs sampler, and its use in deriving covariance matrices by Monte Carlo integration, and in linear and nonlinear error propagation studies, is discussed. Modifications of this technique are given which improve in efficiency in situations where estimated parameters are highly correlated and normal matrices appear as ill-conditioned. This is a situation frequently encountered in satellite gravity field modelling. A synthetic experiment, where covariance matrices for spherical harmonic coefficients are estimated and propagated to geoid height covariance matrices, is described. In this case, the generated samples correspond to random realizations of errors of a gravity field model. AcknowledgementsThe authors are indebted to Pieter Visser and Pavel Ditmar for providing simulation output that was used in the GOCE error generation experiments. Furthermore, the NASA/NIMA/OSU team is acknowledged for providing public ftp access to the EGM96 error covariance matrix. The two anonymous reviewers are thanked for their valuable comments.     

A Real–time Adaptive Sampling Method for Field Mapping in Patchy, Heterogeneous Environments

Many environmental studies require detailed maps describing the spatial distribution of various environmental characteristics. These distributions tend to be 'patchy'; that is, their structure and their relationships vary from place to place according to the influences of the local setting. We present a simple sampling method that adapts the sample spacing on a point–by–point basis as the data are collected. The resulting sample is denser in areas of higher variability and sparser in more 'well–behaved' areas, and is collected in a single traverse of the transect. It uses a combination of simple fuzzy functions representing the concepts 'too close' and 'too far' that are adaptively parameterized based on current measurements. The adaptive sampler produced better representations for 47% of simulated reference transects than uniform or random samples of the same size under perfect measurement conditions, increasing to best performance for 71% of the transects when measurement error was severe with only minimal increase in sampling density. The rapid calculations can be easily incorporated into real–time data acquisition software, and the method may be extended to achieve some type of compromise when faced with the need to sample multiple simultaneous variables.     

GNSS模糊度降相关性能的条件方差平稳度评价法

GNSS模糊度降相关通过整数变换优化条件方差的排列顺序,提高搜索效率。降相关和条件方差的关系及其评价是关键问题之一。针对这一问题,本文从理论上分析了排序后模糊度降相关与条件方差之间的数值关系,发现降相关性能与条件方差数值序列的平稳性有关,降相关性能越强,条件方差数值序列越平稳。基于这一理论关系,给出了"条件方差平稳度"定义,并将其作为评价降相关性能的指标。通过模拟和实测数据验证,并采用条件方差变化趋势图和搜索时间来定性和定量评价降相关性能,用以判定条件方差平稳度的合理性。试验结果表明,条件方差平稳度可以较精确直观地衡量模糊度的降相关性能。本文定义的指标揭示了模糊度降相关的本质。     

基于分布函数的模糊知识表示与推理模型

本文基于模糊集与模糊推理理论,给出了一种模糊谓词模式的知识表示与推理模型。以分布函数来描述模糊语言变量,以可能性分布限定算法作为推理机制,为广泛应用于各种控制系统、专家系统以及人工智能各个领域的模糊推理系统提供了有效的问题描述和推理方式,拓展了Horn子句的表示,为模糊推理系统和工具的研制开发提供了基础。     

An unsupervised segmentation method based on MPM for SAR images

An unsupervised segmentation method for synthetic aperture radar (SAR) images is proposed. It alternately approximates the maximization of the posterior marginals estimate of the pixel class labels and estimates all model parameters except the number of classes during segmentation. In this method, a multilevel logistic (MLL) model for the pixel class labels and Gamma distribution for the marginal distribution of each class in the observed SAR image are employed. In our implementation, the expectation-maximization algorithm is used to estimate parameters of the Gamma distributions, and the iterative conditional estimation algorithm is used to estimate the MLL model parameters. The segmentation results for synthetic and real SAR images show that the proposed method has a good performance.     

非参数吉布斯模型和多波段遥感影像纹理分割方法研究

提出了一种新型的基于概率的非参数Gibbs模型，并建立了光谱空间-二维影像空间的映射关系，通过多级光谱编码策略实现了多波段遥感影像的纹理分割。     

基于马尔柯夫随机场的多波段遥感影像纹理分割研究

利用马尔柯夫随机场（MRF），对光谱矢量在二维空间的映射模式予以描述，建立了基于马尔柯夫随机场的多波段纹理模型，提出了在多波段遥感影像中吉布斯（Gibbs）分布势函数和吉布斯参数的计算方法，并提出了一种迭代修正的多波段遥感影像纹理分割方法。     

全波形LiDAR数据分解的可变分量高斯混合模型及RJMCMC算法

传统激光雷达(light detection and ranging,LiDAR)数据处理均采用固定数的波形分解方法,容易遗漏部分重叠的返回波,降低波形拟合精度。为了实现可变数波形分解,本文提出了一种自动确定波形分解数的方法。假定波形数据服从混合高斯分布,并以此建立理想的波形模型;定义用于控制理想模型与实际波形拟合程度的能量函数,用吉布斯分布构建或然率;根据贝叶斯定理构建刻画波形分解的后验概率模型;设计可逆跳转马尔科夫链蒙特卡洛(reversible jump Markov chain Monte Carlo,RJMCMC)算法模拟该后验概率模型,以确定波形分解数并同时完成波形分解。为了验证提出算法的正确性,分别对不同区域的ICESat-GLAS波形数据进行了波形分解试验,定性和定量分析结果验证了本文方法的有效性、可靠性和准确性。     

An iteratively reweighted least-squares approach to adaptive robust adjustment of parameters in linear regression models with autoregressive and <Emphasis Type="Italic">t</Emphasis>-distributed deviations

In this paper, we investigate a linear regression time series model of possibly outlier-afflicted observations and autocorrelated random deviations. This colored noise is represented by a covariance-stationary autoregressive (AR) process, in which the independent error components follow a scaled (Student’s) t-distribution. This error model allows for the stochastic modeling of multiple outliers and for an adaptive robust maximum likelihood (ML) estimation of the unknown regression and AR coefficients, the scale parameter, and the degree of freedom of the t-distribution. This approach is meant to be an extension of known estimators, which tend to focus only on the regression model, or on the AR error model, or on normally distributed errors. For the purpose of ML estimation, we derive an expectation conditional maximization either algorithm, which leads to an easy-to-implement version of iteratively reweighted least squares. The estimation performance of the algorithm is evaluated via Monte Carlo simulations for a Fourier as well as a spline model in connection with AR colored noise models of different orders and with three different sampling distributions generating the white noise components. We apply the algorithm to a vibration dataset recorded by a high-accuracy, single-axis accelerometer, focusing on the evaluation of the estimated AR colored noise model.     

经典平差模型的通用模型

对经典平差模型之间等价性进行了深入研究,并且在附有限制条件的间接平差、附有限制条件的条件平差可作为概括平差模型的基础之上,推导了用附有参数的条件平差的平差原理来解算经典平差模型的计算公式,得出附有参数的条件平差模型可作为经典平差模型的通用模型这一结论。     

利用马尔柯夫随机场图模型的变化像斑类别判定方法

提出了一种利用马尔柯夫随机场(MRF)图模型获取变化像斑变化后类别的方法。用吉布斯(Gibbs)分布描述MRF图模型,通过计算最大后验概率获取变化像斑类别,以此类别为初始值,通过建立像斑类别空间关系概率矩阵,在初始结果基础上进一步判定变化像斑的类别。实验结果证实了该方法的可行性。     

一种闭合条件约束的全局最优多视点云配准方法

针对已有多视点云配准存在的问题,提出了一种严密的闭合条件约束配准方法。该方法首先采用"点-切平面"迭代最近点算法分别求解各对应点云间的坐标转换参数;再以单站点云为配准单元,并将转换参数视为随机观测值构建条件方程,采用条件平差方法对转换参数作改正以达到全局最优。通过地面三维激光扫描仪实测两组点云数据进行试验,验证了该方法有效且可行。试验结果表明,对采样间隔为毫米级和厘米级的点云,增加扫描重叠度能够提高配准精度和可靠性。     

Misspecifications in interaction model distance decay relations: A spatial structure effect

 An exclusively statistical approach is proposed to address the spatial structure effects of general interaction models. It is shown that the spatial heterogeneity in the estimated region-specific distance decay parameters may in part be due to the combination of two factors: (a) a functional mis-specification of the global distance decay relationship; and (b) the heterogeneity in the region-specific conditional distance distributions. A properly specified global distance decay function allows controlling for these spatially induced biases in the local distance decay parameters. However, inherent multicollinearities between the set of region specific distance decay parameters and other estimated model parameters prevent an unambiguous interpretation. A key conclusion is that a proper model specification, in particular, the specification of the global distance decay relationship, is of paramount importance in interaction modeling and for accessibility studies. Received: September 2002 / Accepted: January 2003     

GBGWR模型在降水量空间分布模拟中的应用

近年来,我国大部分地区屡遭洪涝与干旱两种自然灾害侵袭,对重洪涝干旱区域进行空间插值具有重要的意义。针对传统地理加权回归（GWR）模型建模过程中模型识别和参数估计易受观测值异常点影响的问题,本文提出了一种基于吉布斯采样的贝叶斯地理加权回归（GBGWR）方法。运用基于吉布斯采样的马尔可夫链蒙特卡罗贝叶斯方法,估计地理加权回归模型参数,通过平滑函数降低观测值中异常点位数据,最后对湖南省1985-2015年35个观测站点的降水观测数据进行了空间分布模拟。试验结果表明,本文提出的方法相较于GWR模型性能提高了19.8%,相较于BGWR模型性能提高了8.2%,该方法可以有效降低异常值和"弱数据"对回归结果的影响,能够更加真实地模拟湖南省降水量的空间分布。