---
title: 'DC-Reg: Global Optimal Registration'
url: https://www.emergentmind.com/topics/dc-reg
type: topic
---

# DC-Reg: Global Optimal Registration

Searching arXiv for the target paper and closely related registration work for citations.
DC-Reg is a globally optimal point cloud registration framework for settings in which two point sets may only partially overlap, may be severely misaligned, and may contain many outliers. It formulates registration as the simultaneous estimation of transformation parameters and correspondences, optimizing jointly over a transformation vector \(\boldsymbol{\theta}\) and an assignment matrix \(\mathbf{P}\), rather than treating correspondence as a byproduct of pose estimation. The method’s defining contribution is a branch-and-bound procedure equipped with a “holistic concave underestimator” derived from Difference-of-Convex programming, designed to tighten lower bounds for the coupled transformation-assignment objective and thereby accelerate globally optimal search [2603.25442].

## 1. Problem formulation and scope

DC-Reg addresses the registration problem through the coupled objective
\[
E(\mathbf{P}, \boldsymbol{\theta}) = \sum_{i,j} p_{ij} E_{ij}(\boldsymbol{\theta}),
\]
where \(\boldsymbol{\theta}\in\mathbb{R}^{n_\theta}\) lies in a search box
\[
\mathcal{M} = [\underline{\boldsymbol{\theta}}, \overline{\boldsymbol{\theta}}],
\]
and \(\mathbf{P}=\{p_{ij}\}\) belongs to a feasible assignment set \(\mathcal P\) [2603.25442]. The corresponding marginalized objective is
\[
E(\boldsymbol{\theta}) = \min_{\mathbf{P}\in\mathcal P} \sum_{i,j} p_{ij} E_{ij}(\boldsymbol{\theta}).
\]

In implementation, the feasible correspondence set is
\[
\mathcal{P} = \left\{\mathbf{P} \in \{0,1\}^{N \times M} \;\middle|\; \mathbf{P}\mathbf{1} \le \mathbf{1},\; \mathbf{1}^\top \mathbf{P} \le \mathbf{1}^\top,\; \sum p_{ij} = n_p \right\}.
\]
This makes \(\mathbf{P}\) a binary partial permutation matrix: each source point matches at most one target point, each target matches at most one source point, and exactly \(n_p\) matches are selected [2603.25442]. The mechanism for handling outliers and partial overlap is therefore combinatorial rather than penalty-based: unmatched points are excluded because only \(n_p\) pairs are chosen.

The paper’s central viewpoint is that simultaneous transformation and correspondence estimation is attractive because correspondence flexibility makes the formulation naturally tolerant to partial overlap and some forms of nonrigid disturbance. The method does not directly model nonrigid deformation with deformation parameters; rather, it gains robustness to moderate nonrigid effects because the assignment can select a subset of plausible matches rather than forcing all points to fit a single rigid model [2603.25442]. This suggests that DC-Reg occupies an intermediate position between classical rigid registration and explicitly nonrigid registration: it remains a rigid or similarity optimizer at the transformation level, but with a partial-assignment model that can absorb deviations from exact rigidity.

The paper studies this formulation for a 4-parameter 2D similarity transform and for 3D rigid registration approached through a 3D translation search followed by separate rotation recovery [2603.25442]. A plausible implication is that the framework is most effective when the transformation search space remains low-dimensional, because its lower-bound construction requires vertex evaluation over transformation boxes.

## 2. Transformation models and registration energies

For 2D similarity registration, the transformation is parameterized by \(\boldsymbol{\theta}\in\mathbb{R}^4\), with \([\theta_3,\theta_4]^\top\) the translation and
\[
\theta_1=s\cos\phi,\qquad \theta_2=s\sin\phi
\]
encoding scale \(s\) and rotation angle \(\phi\) [2603.25442]. The registration energy is
\[
E(\mathbf{P}, \boldsymbol{\theta}) = \sum_{i,j} p_{ij} \| \mathbf{J}_i \boldsymbol{\theta} - \mathbf{y}_j \|^2,
\]
with
\[
\mathbf J_i \triangleq \begin{bmatrix} x_i^1 & -x_i^2 & 1 & 0\\ x_i^2 & x_i^1 & 0 & 1 \end{bmatrix}.
\]
Here \(\mathbf J_i\boldsymbol{\theta}\) gives the transformed source point under the 2D similarity parameterization, and the cost is purely spatial least squares [2603.25442].

For 3D registration, DC-Reg uses a rotation-invariant feature based on point norms to search only over translation \(\mathbf t\in\mathbb R^3\):
\[
E(\mathbf{P}, \mathbf{t}) = \sum_{i,j} p_{ij} \left( \|\mathbf{x}_{i} + \mathbf{t}\| - \|\mathbf{y}_{j}\| \right)^2.
\]
Because
\[
\|\mathbf{R}(\mathbf{x}_i+\mathbf t)\|=\|\mathbf{x}_i+\mathbf t\|
\]
for any rotation \(\mathbf R\), the translation search becomes independent of the unknown rotation [2603.25442]. After \(\mathbf t^*\) is found globally, the remaining rotation is recovered by an external rotation search method from the authors’ prior work. The paper presents this as an efficient instantiation of the framework rather than a direct 6D simultaneous rigid search [2603.25442].

This decomposition is significant because it defines the precise sense in which the 3D pipeline is globally optimal. The translation-assignment subproblem in the rotation-invariant feature space is solved globally; the full rigid registration is then completed by a separate rotation solver [2603.25442]. It is therefore inaccurate to describe the 3D version as a direct globally optimal 6D rigid-registration solver. A more precise characterization is a globally optimal translation-assignment framework embedded in a two-stage rigid-registration pipeline.

## 3. Difference-of-Convex lower bounding and branch-and-bound

The key technical contribution is the Difference-of-Convex perspective used to derive a tighter lower bound inside branch-and-bound than earlier term-wise envelope constructions. The paper assumes each pairwise cost admits a DC decomposition
\[
E_{ij}(\boldsymbol{\theta}) = E^{ij}_{\text{cvx}}(\boldsymbol{\theta}) + E^{ij}_{\text{cav}}(\boldsymbol{\theta}),
\]
where \(E^{ij}_{\text{cvx}}\) is convex and \(E^{ij}_{\text{cav}}\) is concave [2603.25442]. Given a search box \(M\subseteq \mathcal M\) and a linearization point \(\boldsymbol{\theta}_0\in M\), the convex part is under-estimated by its first-order affine lower bound,
\[
E^{ij}_{\text{cvx}}(\boldsymbol{\theta}) \ge E^{ij}_{\text{cvx}}(\boldsymbol{\theta}_0)+ \xi^\top(\boldsymbol{\theta}-\boldsymbol{\theta}_0),
\]
yielding the concave underestimator
\[
m_{ij}(\boldsymbol{\theta}) := E^{ij}_{\text{cvx}}(\boldsymbol{\theta}_0)+ \xi^\top(\boldsymbol{\theta}-\boldsymbol{\theta}_0) + E^{ij}_{\text{cav}}(\boldsymbol{\theta}).
\]

The associated auxiliary function is
\[
z(\boldsymbol{\theta})= \min_{\mathbf P\in\mathcal P} \sum_{i,j} p_{ij}\, m_{ij}(\boldsymbol{\theta}).
\]
Because each \(m_{ij}\) is concave, the paper states that \(z(\boldsymbol{\theta})\), as the pointwise minimum of concave functions, is also concave. Therefore minimizing \(z\) over a box \(M\) reduces to vertex evaluation:
\[
LB(M):=\min_{\boldsymbol v\in V(M)} z(\boldsymbol v),
\]
where \(V(M)\) is the set of the \(2^{n_\theta}\) vertices of \(M\) [2603.25442].

The “holistic concave underestimator” is defined in opposition to earlier methods such as RPM-HTB or HBSP, which derive lower bounds by relaxing each bilinear or trilinear term independently, often using McCormick-type convex envelopes [2603.25442]. DC-Reg instead first views each entire pairwise cost \(E_{ij}(\boldsymbol{\theta})\) as a DC function, constructs one underestimator \(m_{ij}\) for that whole pairwise term, and then keeps the resulting transformation-dependent costs coupled with \(\mathbf P\) inside
\[
\min_{\mathbf P\in\mathcal P}\sum_{i,j} p_{ij} m_{ij}(\boldsymbol{\theta}).
\]
The paper’s claim is that this better captures the joint structure of pose and correspondence, leading to tighter bounds and faster pruning [2603.25442].

The branch-and-bound procedure partitions only the transformation space, not the correspondence variables. At each iteration it computes a lower bound \(LB(M)\) and candidate solution \(\boldsymbol{\theta}(M)\) for each newly created box; updates the global upper bound by evaluating the original objective at the best candidate among previous and new ones; prunes boxes satisfying
\[
LB(M)\ge E(\boldsymbol{\theta}^k)-\epsilon;
\]
terminates when no boxes remain; otherwise selects the box with smallest lower bound and bisects its longest edge to create two child boxes [2603.25442]. The paper explicitly states that if all remaining boxes are pruned, then the current best \(\boldsymbol{\theta}^k\) is the global \(\epsilon\)-minimum solution [2603.25442].

## 4. Linear assignment reduction and concrete instantiations

Once \(\boldsymbol{\theta}\) is fixed to a box vertex \(\boldsymbol v\), the lower-bound evaluation becomes
\[
z(\boldsymbol v) = \min_{\mathbf P\in\mathcal P} \sum_{i,j} p_{ij}\, m_{ij}(\boldsymbol v),
\]
which is an \(n_p\)-cardinality Linear Assignment Problem with cost matrix
\[
C_{ij}(\boldsymbol v)=m_{ij}(\boldsymbol v)
\]
[2603.25442]. In implementation, the paper uses the Jonker–Volgenant algorithm for these LAPs. Because the search dimension is low in the studied cases, each branch-and-bound node requires only \(2^4=16\) LAP evaluations for 2D similarity and \(2^3=8\) for 3D translation [2603.25442].

For the 2D similarity case, each pairwise cost
\[
E_{ij}(\boldsymbol{\theta})=\|\mathbf J_i\boldsymbol{\theta}-\mathbf y_j\|^2
\]
is convex in \(\boldsymbol{\theta}\), so the DC decomposition is trivial:
\[
E^{ij}_{\text{cav}}=0,\qquad E^{ij}_{\text{cvx}}=E_{ij}.
\]
The underestimator becomes the first-order Taylor lower bound
\[
m_{ij}(\boldsymbol{\theta}) = E_{ij}(\boldsymbol{\theta}_0)+ \nabla E_{ij}(\boldsymbol{\theta}_0)^\top(\boldsymbol{\theta}-\boldsymbol{\theta}_0),
\]
with
\[
\nabla E_{ij}(\boldsymbol{\theta}_0) = 2\mathbf J_i^\top(\mathbf J_i\boldsymbol{\theta}_0-\mathbf y_j).
\]
This yields an especially clean implementation in which the bound reduces to evaluating affine costs at box vertices and solving cardinality-constrained assignments [2603.25442].

For the 3D translation search, the pairwise cost
\[
E_{ij}(\mathbf t)=\left(\|\mathbf x_i+\mathbf t\|-\|\mathbf y_j\|\right)^2
\]
admits the explicit DC decomposition
\[
E^{ij}_{\text{cvx}}(\mathbf t)=\|\mathbf x_i+\mathbf t\|^2,
\]
\[
E^{ij}_{\text{cav}}(\mathbf t)=-2\|\mathbf y_j\|\,\|\mathbf x_i+\mathbf t\|+\|\mathbf y_j\|^2.
\]
Linearizing the convex part at \(\mathbf t_0\in M\) gives
\[
E^{ij}_{\text{cvx}}(\mathbf t) \ge E^{ij}_{\text{cvx}}(\mathbf t_0) + \nabla E^{ij}_{\text{cvx}}(\mathbf t_0)^\top(\mathbf t-\mathbf t_0),
\]
with
\[
\nabla E^{ij}_{\text{cvx}}(\mathbf t_0)=2(\mathbf x_i+\mathbf t_0),
\]
and thus
\[
m_{ij}(\mathbf t) := E^{ij}_{\text{cvx}}(\mathbf t_0)+ 2(\mathbf x_i+\mathbf t_0)^\top(\mathbf t-\mathbf t_0) + E^{ij}_{\text{cav}}(\mathbf t).
\]
Since translation is 3D, only eight vertices are checked per box, which the paper identifies as the main source of efficiency in the 3D variant [2603.25442].

The dominant cost per branch-and-bound node is therefore
\[
(\text{number of BnB nodes}) \times 2^{n_\theta} \times (\text{cost of one LAP}),
\]
with practical behavior determined jointly by search dimension and bound tightness [2603.25442]. The authors further note that a direct 3D affine version would be impractical because \(2^{12}=4096\) vertex evaluations per box would be required [2603.25442].

## 5. Robustness model, guarantees, and limitations

DC-Reg handles missing correspondences and outliers through the assignment constraints and the fixed cardinality \(n_p\). Because only \(n_p\) matches are selected, surplus points in either set are treated as outliers; the row and column inequality constraints impose one-to-one matching among selected pairs; and the cardinality controls how much overlap is presumed [2603.25442]. There are no dummy nodes, slack variables, truncation penalties, or explicit outlier costs in the formulation.

This design yields a combinatorially clean partial-assignment model compatible with LAP solvers, but it also introduces a modeling dependence on \(n_p\). The paper explicitly studies sensitivity to this choice by using \(n_p=1/1\) and \(n_p=1/2\) of the ground-truth inlier count in experiments [2603.25442]. A plausible implication is that the method is most attractive when a reasonable inlier-count prior is available, even if that prior is somewhat underestimated.

The theoretical statements are property-based rather than theorem-heavy. The central validity arguments are that \(m_{ij}\) is a valid underestimator because a convex function dominates its tangent affine approximation; \(z(\boldsymbol{\theta})\) is concave as a pointwise minimum of concave functions; and minimizing a concave function over a box is achieved at a vertex, making \(LB(M)\) valid [2603.25442]. Combined with a standard branch-and-bound scheme over a compact initial box and exact lower bounds for each box, the method inherits the usual global \(\epsilon\)-optimality guarantee of deterministic branch-and-bound [2603.25442].

The authors also argue qualitatively that the bound is tighter than previous term-wise relaxations because it preserves variable interactions that McCormick-style envelopes ignore, but they do not provide a formal dominance theorem over prior bounds in the text shown [2603.25442]. This is an important interpretive boundary: the empirical acceleration claim is strong, but the comparative tightness claim remains qualitative in the paper’s formal development.

The principal limitations identified in the paper are twofold. First, computational overhead still grows with branching and transformation dimension [2603.25442]. Second, the framework does not directly model richer nonrigid deformations [2603.25442]. The authors suggest future work on extending the DC paradigm to nonrigid models and accelerating branch-and-bound with GPU-parallel pruning [2603.25442].

## 6. Empirical evaluation and relation to prior global registration methods

DC-Reg is evaluated on synthetic and real 2D tests, synthetic 3D tests, and the 3DMatch benchmark [2603.25442]. In 2D synthetic experiments on “fish” and “character” shapes under deformation, positional noise, mixed outliers, separate outliers, and occlusion plus outliers, baselines are RPM-HTB, RPM-PA, and RPM-CAV [2603.25442]. The evaluation metric is RMSE between transformed source inliers and corresponding target inliers. The paper reports that DC-Reg is substantially more stable under deformation and positional noise, and especially strong under mixed and separated outliers, where it maintains near-zero error while baselines degrade; even when \(n_p\) is underestimated to half the true inlier count, it remains competitive and often superior [2603.25442].

On 2D real-image experiments using Canny edge point sets from Caltech-256 and VOC2007 images, with source sets rotated by \(180^\circ\), the paper reports near-perfect overlap on examples including bicycles, motorbikes, guns, and the Eiffel Tower, whereas baseline methods often fail under heavy clutter or occlusion [2603.25442]. The authors interpret this as evidence that the bound better captures global spatial structure.

For 3D synthetic evaluation, the baselines include RPM-HTB, GO-ICP, FRS, GORE, and TEASER++ [2603.25442]. DC-Reg with \(n_p=1/1\) is reported to have the lowest registration error in deformation and noise tests, and to remain strongest in mixed-outlier settings up to outlier ratio \(0.5\) [2603.25442]. The paper is more nuanced for separated outliers and heavy occlusion: for separated outliers, DC-Reg is only moderate and can be outperformed by GO-ICP or GORE; for occlusion plus outliers, it remains stable but is not the best, with TEASER++ and GORE showing stronger resilience to large missing-data regimes [2603.25442]. The authors attribute some of this limitation to the rotation-invariant feature based on norms, which is more vulnerable when clutter is spatially structured.

On the 3DMatch benchmark, using standard D3Feat scan pairs across five RGB-D reconstruction datasets, the paper reports average registration error versus overlap ratio [2603.25442]. The main qualitative conclusion is that DC-Reg maintains low error across datasets, remains particularly stable in low-overlap regimes, significantly outperforms RPM-HTB throughout, and achieves parity with strong robust solvers such as TEASER++ and GORE on real-world data [2603.25442]. No detailed numeric tables are provided in the supplied text, so the most faithful summary remains qualitative.

In the broader registration literature, DC-Reg is best understood as a globally optimal partial-assignment method positioned against deterministic global rigid solvers such as GO-ICP [1601.03702], robust correspondence-driven estimators such as TEASER++ [2001.07715], and outlier-rejection methods such as GORE [1711.02588]. Its distinguishing feature is not merely global search, nor merely correspondence flexibility, but the specific lower-bounding strategy that preserves transformation-assignment coupling within the assignment minimization itself [2603.25442]. This suggests that DC-Reg is particularly relevant where global optimality, poor initialization, partial overlap, and highly ambiguous correspondences are simultaneous requirements.

Source: https://www.emergentmind.com/topics/dc-reg