Skip to content

(ICML 2020) Message Passing Least Squares Algorithm for Rotation Averaging

Notifications You must be signed in to change notification settings

yunpeng-shi/MPLS

Folders and files

NameName
Last commit message
Last commit date

Latest commit

 
 
 
 
 
 
 
 
 
 
 

Repository files navigation

MPLS (Minneapolis) Framework for Robust Rotation Averaging

News:

The python3 version of CEMP-SO(3) for rotation averaging is released at CEMP_SO3_Python !

Introduction and References

Message Passing Least Squares (MPLS) is a powerful alternative to the popular iteratively Reweighted Least Squares (IRLS) algorithm for solving the general problem of group synchronization (and rotation averaging as a special case).

This repo contains matlab files for implementing the method of the following paper.

[1] Message Passing Least Squares Framework and Its Applications in Rotation Synchronization, Yunpeng Shi and Gilad Lerman, ICML 2020.

If you would like to use our code for your paper, please cite [1].

A closely related work to [1] is

[2] Robust Group Synchronization via Cycle-Edge Message Passing, Gilad Lerman and Yunpeng Shi, Foundations of Computational Mathematics, 2021.

The main difference is that in [2] we study cycle-edge message passing (not MPLS) for general compact groups (so it has broader applications), and it uses all 3-cycles (so it is deterministic, unlike MPLS that randomly samples cycles). We include its code in CEMP repository.

Another MPLS-type algorithm for permutation synchronization appears in Robust Multi-object Matching via Iterative Reweighting of the Graph Connection Laplacian. Yunpeng Shi, Shaohan Li, Gilad Lerman. NeurIPS 2020. Its code is available in IRGCL repository.

Usage

Download matlab files to the same directory. Checkout and run the following demo code.

Examples/Compare_algorithms.m

The following is a sample output (error in degrees) under Nonuniform and Adversarial corruption, which shows great advantage of CEMP and MPLS:

Algorithms     MeanError     MedianError
___________    __________    ___________

"Spectral"         13.071         3.902 
"IRLS-GM"          34.093        3.3797 
"IRLS-L0.5"        5.4121       0.57141 
"CEMP+MST"     1.9138e-06    1.7075e-06 
"CEMP+GCW"       0.014078     0.0043606 
"MPLS"         3.0585e-06    2.4148e-06 

Here Spectral refers to eigenvector method for approximately solving least squares formulation. See Angular Synchronization by Eigenvectors and Semidefinite Programming, Amit Singer, Applied and Computational Harmonic Analysis, 2011 for details.

IRLS-GM refers to iteratively reweighted least squares (IRLS) that uses Geman-McClure loss function. See Efficient and Robust Large-Scale Rotation Averaging, Avishek Chatterjee and Venu Madhav Govindu, ICCV 2013 for details.

IRLS-L0.5 refers to iteratively reweighted least squares (IRLS) that uses L_1/2 loss function. See Robust Relative Rotation Averaging, Avishek Chatterjee and Venu Madhav Govindu, TPAMI, 2018 for details.

CEMP+MST, CEMP+GCW, MPLS refer to our methods. CEMP+GCW only appears in CEMP paper Robust Group Synchronization via Cycle-Edge Message Passing, Gilad Lerman and Yunpeng Shi, Foundations of Computational Mathematics, 2021. It combines CEMP with a weighted spectral method (using graph connection weight matrix).

A Variety of Corruption Models

We provide 5 different corruption models. 3 for nonuniform topology and 2 for uniform toplogy (see Uniform_Topology.m and Nonuniform_Topology.m). Uniform/Nonuniform toplogy refers to whether the corrupted subgraph is Erdos Renyi or not. In other words, the choice of Uniform/Nonuniform toplogy decides how to select edges for corruption. In Uniform_Topology.m, two nodes are connected with probability p. Then edges are independently drawn with probability q for corruption. In Nonuniform_Topology.m, two nodes are connected with probability p. Then with probability p_node_crpt a node is selected so that its neighboring edges are candidates for corruption. Next, for each selected node, with probability p_edge_crpt an edge (among the neighboring edges of the selected node) is corrupted. This is a more malicious scenario where corrupted edges have cluster behavior (so local coruption level can be extremely high).

One can also optionally add noise to inlier graph and outlier graph through setting sigma_in and sigma_out for Nonuniform_Topology.m. For Uniform_Topology.m we assume inlier and outlier subgraph have the same noise level sigma.

The argument crpt_type in the two functions determines how the corrupted relative rotations are generated for those selected edges. In Uniform_Topology.m, there are 2 options of crpt_type: uniform and self-consistent. In Nonuniform_Topology.m, there are the following 3 options of crpt_type.

uniform: The corrupted relative rotations are i.i.d follows uniform distribution over the space of SO(3).

self-consistent: The corrupted are relative rotations of another set of absolute rotations. Namely where those absolute rotations are different from the ground truth and are i.i.d drawn from the uniform distribution in the space of SO(3). In this way, the corrupted relative rotations are also cycle-consistent.

adv: Extremely malicious corruption that replaces the underlying absolute rotations from ground truth to . Namely for the corrupted neighboring edges (i,j) of node i. Additional high noise must be added to the outlier-subgraph, otherwise the recovery of the ground truth can be ill-posed. This model is not included in MPLS paper. It was first introdued in Robust Multi-object Matching via Iterative Reweighting of the Graph Connection Laplacian, NeurIPS 2020 for permutation synchronization.

Implementation of MPLS

The demo code uses the following function for implementing MPLS:

MPLS(Ind, RijMat, CEMP_parameters, MPLS_parameters)

Each row of Ind matrix is an edge index (i,j). The edge indices (the rows of Ind) MUST be sorted in row-major order. That is, the edge indices are sorted as for example (1,2), (1,3), (1,4),..., (2,3), (2,5), (2,8),..., otherwise the code may crash when some edges are not contained in any 3-cycles. Make sure that i<j. If some edges have indices (3,1), then change it to (1,3) and take a transpose to the corresponding Rij.

MPLS is not sensitive to its parameters. The general rule for reweighting parameters

CEMP_parameters.reweighting
MPLS_parameters.reweighting

is the following:

As the iteration # increases, both of them increases (or at least be nondecreasing) and then stop at a constant between 20-50 (cannot be arbitrarily large due to numerical instability).

The following parameter determines the fraction of cycle-consistency information that MPLS takes account.

MPLS_parameters.cycle_info_ratio

That is, at each iteration

Estimated corruption level of each edge
= (1-MPLS_parameters.cycle_info_ratio) * Residual + MPLS_parameters.cycle_info_ratio * (Cycle inconsistency measure)

For sparse graphs (in many real scenarios) we recommend (rely more on residuals)

MPLS_parameters.reweighting = CEMP_parameters.reweighting(end);
MPLS_parameters.cycle_info_ratio = 1./((1:MPLS_parameters.max_iter)+1); % in the end ignore 3-cycle information

For very dense graphs, we find the following parameters perform better (rely more on 3-cycle information)

MPLS_parameters.reweighting = 0.1*1.5.^((1:15)-1);    % more and more aggressive reweighting
MPLS_parameters.cycle_info_ratio = 1-1./((1:MPLS_parameters.max_iter)+1); % in the end focus only on 3-cycle information

For real data experiments we used the data provided by the authors of https://openaccess.thecvf.com/content_cvpr_2017/papers/Sengupta_A_New_Rank_CVPR_2017_paper.pdf. The random seed was chosen as 2020 (as it was submitted to ICML 2020).

Dependencies

The following files in folder Utils are dependencies for running Lie-Algebraic Averaging method that were written by AVISHEK CHATTERJEE (revised and included in this repo). See also Robust Rotation Averaging and Efficient and Robust Large-Scale Rotation Averaging for details.

Weighted_LAA.m
Build_Amatrix.m
R2Q.m
q2R.m
BoxMedianSO3Graph.m
L12.m
RobustMeanSO3Graph.m

About

(ICML 2020) Message Passing Least Squares Algorithm for Rotation Averaging

Topics

Resources

Stars

Watchers

Forks

Releases

No releases published

Languages