A THRESHOLD-FREE FILTERING ALGORITHM FOR AIRBORNE LIDAR POINT CLOUDS BASED ON EXPECTATION-MAXIMIZATION

The expectation maximisation (EM) algorithm is an iterative maximum likelihood procedure often used for estimating the parameters of a mixture model. Theoretically, increases in the likelihood function are guaranteed as the algorithm iteratively improves upon previously derived parameter estimates. The algorithm is considered to converge when all parameter estimates become stable and no further improvements can be made to the likelihood value. However, to reduce computational time, it is often common practice for the algorithm to be stopped before complete convergence using heuristic approaches. In this paper, we consider various stopping criteria and evaluate their effect on fitting Gaussian mixture models (GMMs) to patient length of stay (LOS) data. Although the GMM can be successfully fitted to positively skewed data such as LOS, the fitting procedure often requires many iterations of the EM algorithm. To our knowledge, no previous study has evaluated the effect of different stopping criteria on fitting GMMs to skewed distributions. Hence, the aim of this paper is to evaluate the effect of various stopping criteria in order to select and justify their use within a patient spell classification methodology. Results illustrate that criteria based on the difference in the likelihood value and on the GMM parameters may not always be a good indicator for stopping the algorithm. In fact we show that the values of the difference in the variance parameters should be used instead, as these parameters are the last to stabilise. In addition, we also specify threshold values for the other stopping criteria.

The accurate modeling of uncertain variables, such as wind power and load uncertainty, is an essential foundation for the planning operation and scheduling decisions in power systems. Currently, Gaussian mixture model (GMM) and Dirichlet process mixture model (DPMM) are the primary techniques used to characterize uncertainties in power systems, which are commonly solved by expectation-maximization (EM) algorithm. However, for the massive data of uncertain variables, the algorithm encounters challenges in accurately obtaining GMM and DPMM with a lower time consumption. To address this issue, we propose a method for GMM uncertainty modeling in power systems considering the mutual assistance of latent variables. Specifically, the GMM of uncertain variables is first constructed, and the conditional probability is employed to characterize the mutual assistance of latent variables. Then, an improved EM algorithm is developed to obtain the optimal GMM parameters, in which the expectation step (E-step) and maximization step (M-step) of the algorithm are revised using the conditional probability. Importantly, the closed-form solutions for GMM parameters are rederived in the revised E-step and M-step. Finally, the proposed uncertainty modeling method is compared with the traditional GMM and DPMM on actual wind power and load data from Australia. The proposed method performs the efficiency and accuracy in characterizing uncertainty.

One effective iterative method for estimating parameters in latent variable probabilistic models is the Expectation-Maximization (EM) algorithm, which finds extensive use in classification issues. A versatile probabilistic framework for modeling data from multiple Gaussian distributions, Gaussian Mixture Models (GMMs) employ EM for parameter estimation, which is explored in this article. Problems with class membership estimation (including means, covariances, and mixing coefficients) are the primary focus of this study. Starting with initial parameter estimations, estimating posterior probabilities in the E-step, and repeatedly updating model parameters in the M-step until convergence, the technique entails applying the EM algorithm on both synthetic and real-world datasets. In order to measure performance, we compare it to K-means and other baseline classifiers using measures like classification accuracy, Adjusted Rand Index (ARI), and log-likelihood progression. When data is distributed according to a multimodal Gaussian distribution, the findings show that EM-GMM provides better parameter estimation and better classification accuracy. Initialization procedures have a substantial effect on performance, as convergence study further shows. In addition to reiterating the EM algorithm's usefulness for probabilistic classification, this paper also describes its advantages and disadvantages and offers suggestions for improving its performance in settings with high-dimensional data and noise.

The traditional Expectation-Maximization (EM) training of Gaussian Mixture Model (GMM) is essentially a batch mode procedure which requires the multiple data samples with the sufficient size to update the model parameters. This severely limits the deployment and adaptation of GMM in many real-time online systems since the newly observed data samples are expected to be incorporated into the system upon available via retraining the model. This paper presents a new online incremental EM training procedure of GMM, which aims to perform the EM training incrementally and so can adapt GMM online sample by sample. The proposed method is extended on two kinds of EM algorithms for GMM, namely, Split-and-Merge EM and the traditional EM. Experiments on both the synthetic data and a speech processing task show the advantages and efficiency of the new method.

This article proposes the distribution similarity measure–based Gaussian mixtures model for the contact-state (CS) modelling in force-guided robotic assembly processes of flexible rubber parts. The wrench (Cartesian force and torque) signals of the manipulated object are captured for different states of the given assembly process. The distribution similarity measure–based Gaussian mixtures model CS modelling scheme is employed in modelling the captured wrench signals for different CSs. The proposed distribution similarity measure–based Gaussian mixtures model CS modelling scheme uses the Gaussian mixtures model in modelling the captured signals. The parameters of the Gaussian mixtures models are computed using expectation maximisation. The optimal number of Gaussian mixtures model components for each CS model is determined by considering the classification success rate as an index for the similarity measure between the distribution of the captured signals and the developed models. The optimal number of Gaussian mixtures model components corresponds to the highest classification success rate; hence, object elasticity variation would be accommodated by properly choosing the optimal number of Gaussian mixtures model components. The performance of the proposed distribution similarity measure–based Gaussian mixtures model CS modelling strategy is evaluated by a test stand composed of a KUKA lightweight robot doing peg-in-hole assembly processes for flexible rubber objects. Two rubber objects with different elasticity are considered for two experiments; in the first experiment, an elastic peg of 30 Shore A hardness is considered and that of the second experiment has hardness of 6 Shore A which is even softer than the one used in experiment 1. Employing the proposed distribution similarity measure–based Gaussian mixtures model CS modelling strategy excellent classification success rate was obtained for both experiments. However, more Gaussian mixtures model components are required for the softer one that gives a strong impression of the non-stationarity behaviour increment for softer materials. Comparison is performed with the available CS modelling schemes and the distribution similarity measure–based Gaussian mixtures model is shown to provide the best classification success rate performance with a reduced computational time.

The expectation-maximization (EM) algorithm is a popular approach for parameter estimation of finite mixture model (FMM). A drawback of this approach is that the number of components of the finite mixture model is not known in advance, nevertheless, it is a key issue for EM algorithms. In this paper, a penalized minimum matching distance-guided EM algorithm is discussed. Under the framework of Greedy EM, a fast and accurate algorithm for estimating the number of components of the Gaussian mixture model (GMM) is proposed. The performance of this algorithm is validated via simulative experiments of univariate and bivariate Gaussian mixture models.

In a distribution power network, the load model has no certain pattern or predicted behaviour due to large range of data and changes in energy consumption for end-user consumers. Thus, a powerful analysis based on probabilistic structure is required. For this paper Gaussian Mixture Model (GMM) has been used. GMM is a powerful probability model that allows different types of load distributions to be presented as a combination of several Gaussian distributions. The parameters of GMM is unknown for large random data such as real load data and these parameters can be identified by Expectation Maximization (EM) algorithm. This paper presents a method to evaluate probabilistic load data concerning the time-evolution of any type of distribution load for any duration of time. The proposed method is explained through generated load data of 100 residential houses for duration of one year.

Active curve axis Gaussian mixture models

This paper proposes a Gaussian mixture model of a ground filtering method based on hierarchical curvature constraints. Firstly, the thin plate spline function is iteratively applied to interpolate the reference surface. Secondly, gradually changing grid size and curvature threshold are used to construct hierarchical constraints. Finally, an adaptive height difference classifier based on the Gaussian mixture model is proposed. Using the latent variables obtained by the expectation-maximization algorithm, the posterior probability of each point is computed. As a result, ground and objects can be marked separately according to the calculated possibility. 15 data samples provided by the International Society for Photogrammetry and Remote Sensing are used to verify the proposed method, which is also compared with eight classical filtering algorithms. Experimental results demonstrate that the average total errors and average Cohen's kappa coefficient of the proposed method are 6.91% and 80.9%, respectively. In general, it has better performance in areas with terrain discontinuities and bridges.

Skeleton extraction from 3D plant point cloud data is an essential prior for myriads of phenotyping studies. Although skeleton extraction from 3D shapes have been studied extensively in the computer vision and graphics literature, handling the case of plants is still an open problem. Drawbacks of the existing approaches include the zigzag structure of the skeleton, nonuniform density of skeleton points, lack of points in the areas having complex geometry structure, and most importantly the lack of biological relevance. With the aim to improve existing skeleton structures of state-of-the-art, we propose a stochastic framework which is supported by the biological structure of the original plant (we consider plants without any leaves). Initially we estimate the branching structure of the plant by the notion of β-splines to form a curve tree defined as a finite set of curves joined in a tree topology with certain level of smoothness. In the next phase, we force the discrete points in the curve tree to move toward the original point cloud by treating each point in the curve tree as a center of Gaussian, and points in the input cloud data as observations from the Gaussians. The task is to find the correct locations of the Gaussian centroids by maximizing a likelihood. The optimization technique is iterative and is based on the Expectation Maximization (EM) algorithm. The E-step estimates which Gaussian the observed point cloud was sampled from, and the M-step maximizes the negative log-likelihood that the observed points were sampled from the Gaussian Mixture Model (GMM) with respect to the model parameters. We experiment with several real world and synthetic datasets and demonstrate the robustness of the approach over the state-of-the-art.

To model manipulation tasks, we propose a novel method for learning manipulation skills based on the degree of motion granularity. Even though manipulation tasks usually consist of a mixture of fine-grained and coarse-grained movements, to the best of our knowledge, manipulation skills have so far been modeled without considering their motion granularity. To model such a manipulation skill, Gaussian mixture models (GMMs) have been represented using several well-known techniques such as principal component analysis, k-means, Bayesian information criterion, and expectation-maximization (EM) algorithms. However, in this GMM, there is a problem in that when a mixture of fine-grained and coarse-grained movements is modeled as a GMM, fine-grained movements tend to be poorly represented. To resolve this issue, we measure a continuous degree of motion granularity for every time step of a manipulation task from a GMM. Then, we remodel the GMM by weighting a conventional k-means algorithm with motion granularity. Finally, we also estimate the parameters of the GMM by weighting the conventional EM with motion granularity. To validate our proposed method, we evaluate the GMM estimated using our proposed method by comparing it with those estimated by different GMMs in terms of inference, regression, and generalization using a robot arm that performs two daily tasks, namely decorating a very small area and passing through a narrow tunnel.

Computing multiple aggregation levels and contextual features for road facilities recognition using mobile laser scanning data

Ground points of lidar is the redundant information in the perception process of autonomous driving. Removal of ground points is beneficial to the improvement of accuracy and real-time performance in the perception process. Aiming at the problems of under-segmentation and over-segmentation in current methods, proposes a segmentation method based on point cloud normal clustering. First, the original point cloud is preprocessed to reduce the time cost of the whole process; then, the preprocessed point cloud is quickly pre-segmented by the ray threshold method, and the ground point cloud containing some non-ground points is obtained; finally, according to the normal vector features of ground point and non-ground points, we filter out the non-ground points in the ground point cloud to obtain the final ground point cloud. The experimental results show that the average precision and recall rate of this method in various road conditions can reach 90%, and the average time consumption is 50 ms, which has good real-time performance and adaptability.

Surface registration is an inevitable step for solving the correspondence problem in statistical shape modeling and useful for automatic 3D segmentation of objects of a certain shape. Common techniques for a pairwise surface registration rely either on a point-based or a mesh-based representation of the source (moving) and the target (fixed) surface. Iterative closest point (ICP) algorithms operate on unstructured point clouds of surfaces. In the so-called EM-ICP methods, they assume a Gaussian mixture model (GMM) for the registration in order to efficiently estimate the optimal transform and the correspondence probabilities by an iterative expectation-maximization (EM) algorithm. However, point-based techniques cannot preserve the topology of the source surface. Common mesh-based techniques employ different elastic deformations to estimate a parametric model of the optimal registration transform. Both the point-based and the mesh-based schemes can lead to a deformed mesh with flipped or irregular faces. In this paper, we propose a pairwise nonlinear surface registration algorithm that processes both the point clouds and the triangular meshes of surfaces within a unified framework by combing a modified symmetric robust EM-ICP (SR-EM-ICP) algorithm with a mesh optimization mechanism. It improves the registration accuracy and the quality of the deformed mesh in a processing time comparable to the SR-EM-ICP algorithm. These are achieved by an efficient estimation of the inhomogeneous variances of the GMM and the optimal registration transform, reducing the distance between the surfaces and an explicit recovery of the flipped or irregular faces of the deformed mesh in each iteration.

Model-based clustering approaches concern the paradigm of exploratory data analysis relying on the finite mixture model to automatically find a latent structure governing observed data. They are one of the most popular and successful approaches in cluster analysis. The mixture density estimation is generally performed by maximizing the observed-data log-likelihood by using the expectation-maximization (EM) algorithm. However, it is well-known that the EM algorithm initialization is crucial. In addition, the standard EM algorithm requires the number of clusters to be known a priori. Some solutions have been provided in [31], [12] for model-based clustering with Gaussian mixture models for multivariate data. In this paper we focus on model-based curve clustering approaches, when the data are curves rather than vectorial data, based on regression mixtures. We propose a new robust EM algorithm for clustering curves. We extend the model-based clustering approach presented in [31] for Gaussian mixture models, to the case of curve clustering by regression mixtures, including polynomial regression mixtures as well as spline or B-spline regressions mixtures. Our approach both handles the problem of initialization and the one of choosing the optimal number of clusters as the EM learning proceeds, rather than in a twofold scheme. This is achieved by optimizing a penalized log-likelihood criterion. A simulation study confirms the potential benefit of the proposed algorithm in terms of robustness regarding initialization and funding the actual number of clusters.