290 research outputs found
Generalized Weiszfeld algorithms for Lq optimization
In many computer vision applications, a desired model of some type is computed by minimizing a cost function based on several measurements. Typically, one may compute the model that minimizes the L₂ cost, that is the sum of squares of measurement errors with respect to the model. However, the Lq solution which minimizes the sum of the qth power of errors usually gives more robust results in the presence of outliers for some values of q, for example, q = 1. The Weiszfeld algorithm is a classic algorithm for finding the geometric L1 mean of a set of points in Euclidean space. It is provably optimal and requires neither differentiation, nor line search. The Weiszfeld algorithm has also been generalized to find the L1 mean of a set of points on a Riemannian manifold of non-negative curvature. This paper shows that the Weiszfeld approach may be extended to a wide variety of problems to find an Lq mean for 1 ≤ q <; 2, while maintaining simplicity and provable convergence. We apply this problem to both single-rotation averaging (under which the algorithm provably finds the global Lq optimum) and multiple rotation averaging (for which no such proof exists). Experimental results of Lq optimization for rotations show the improved reliability and robustness compared to L₂ optimization.This research has been funded by National ICT Australia
A D.C. Algorithm via Convex Analysis Approach for Solving a Location Problem Involving Sets
We study a location problem that involves a weighted sum of distances to
closed convex sets. As several of the weights might be negative, traditional
solution methods of convex optimization are not applicable. After obtaining
some existence theorems, we introduce a simple, but effective, algorithm for
solving the problem. Our method is based on the Pham Dinh - Le Thi algorithm
for d.c. programming and a generalized version of the Weiszfeld algorithm,
which works well for convex location problems
A Look at the Generalized Heron Problem through the Lens of Majorization-Minimization
In a recent issue of this journal, Mordukhovich et al.\ pose and solve an
interesting non-differentiable generalization of the Heron problem in the
framework of modern convex analysis. In the generalized Heron problem one is
given closed convex sets in \Real^d equipped with its Euclidean norm
and asked to find the point in the last set such that the sum of the distances
to the first sets is minimal. In later work the authors generalize the
Heron problem even further, relax its convexity assumptions, study its
theoretical properties, and pursue subgradient algorithms for solving the
convex case. Here, we revisit the original problem solely from the numerical
perspective. By exploiting the majorization-minimization (MM) principle of
computational statistics and rudimentary techniques from differential calculus,
we are able to construct a very fast algorithm for solving the Euclidean
version of the generalized Heron problem.Comment: 21 pages, 3 figure
Robust Rotation Synchronization via Low-rank and Sparse Matrix Decomposition
This paper deals with the rotation synchronization problem, which arises in
global registration of 3D point-sets and in structure from motion. The problem
is formulated in an unprecedented way as a "low-rank and sparse" matrix
decomposition that handles both outliers and missing data. A minimization
strategy, dubbed R-GoDec, is also proposed and evaluated experimentally against
state-of-the-art algorithms on simulated and real data. The results show that
R-GoDec is the fastest among the robust algorithms.Comment: The material contained in this paper is part of a manuscript
submitted to CVI
- …