Supervised Fitting of Geometric Primitives to 3 D Point Clouds – Supplementary Material

Lingxiao Li,Minhyuk Sung,Anastasia Dubrovina,Li Yi,Leonidas Guibas

2019-01-01

Abstract:In our differentiable model estimator, we are solving two linear algebra problems: homogeneous least square and unconstrained least square. In both problems, numerical stability issues can occur. We solve the homogeneous least square using SVD to find v, the right singular vector corresponding to the smallest singular value. However, when backpropagating the gradient through SVD [1, 2], the gradient value goes to infinity when the singular values of the input matrix are not all distinct. In our case, such an issue happens only when the output segment (decided by membership matrix Ŵ) becomes degenerate. For instance, when fitting a plane to a segment via SVD, non-distinct singular values correspond to the case where the points in the segment with significant weights concentrate on a line or a single point. Hence if we get good segmentation by minimizing the segmentation loss (Section 3.3 in the paper), then such degenerate cases should not happen. Thus, we handle the issue by simply bounding the gradient in the following way. We implemented a custom SVD layer following [1, 2], and when computing Kij = 1 σi−σj in Equation 13 of [1] where σi, σj are singular values, we instead use Kij = 1 sign(σi−σj)max(|σi−σj |, ) for = 10−10. When solving the unconstrained least square using Cholesky factorization, numerical unstability can happen even when the segmentation is correct, but the type used in the estimator does not match with the segment. For instance, when fitting a sphere to a segment that is almost a flat plane, the optimal sphere is the one with center at infinity. To deal with such a singular case (as well as cases when the segments are degenerate), we add a l2-regularizer to the formulation (Equation 7 in the paper) and solve instead