---
title: 'SplitGS-Loc: Gaussian Splatting for Localization'
url: https://www.emergentmind.com/topics/splitgs-loc
type: topic
---

# SplitGS-Loc: Gaussian Splatting for Localization

Searching arXiv for the primary paper and closely related Gaussian-splatting localization work mentioned in the provided data.
SplitGS-Loc is a localization-specialized framework for constructing Gaussian Splatting-based Feature Fields (GSFFs) from a photometrically trained 3D Gaussian Splatting (3DGS) scene, with the explicit objective of supporting robust direct 2D–3D matching and PnP-based camera pose estimation without any per-scene feature training or iterative pose refinement. Its central premise is that photometrically optimized GSFFs are structurally misaligned with the requirements of point-based pose solvers: the volumetric extent of Gaussians induces ambiguous many-to-one pixel-to-point correspondences, and photometric optimization introduces superfluous Gaussians that lack multi-view consistency. SplitGS-Loc addresses both issues by exploiting Gaussian attributes—especially covariance-derived scales, opacity, and rasterization composition weights—to produce compact, discriminative, and multi-view-consistent feature fields for visual localization [2605.07351].

## 1. Problem formulation and the mismatch between photometric GSFFs and PnP

In structure-based visual localization, the task is to estimate camera pose $(\mathbf{R}, \mathbf{t})$ for a query image by establishing 2D–3D correspondences $(\mathbf{u}_k \leftrightarrow \mathbf{X}_k)$ and solving PnP with RANSAC. The optimization target is written as
$$
\min_{\mathbf{R},\mathbf{t}} \sum_k \left\| \pi\big(\mathbf{R}\mathbf{X}_k + \mathbf{t}\big) - \mathbf{u}_k\right\|^2,
$$
where $\pi$ denotes camera projection [2605.07351].

GSFFs attach learnable descriptors $\mathbf{z}_n$ to Gaussians in a 3DGS scene, so that query image features can be matched directly to Gaussians, yielding 2D–3D correspondences. SplitGS-Loc is motivated by the observation that GSFFs obtained by photometric training, or by rendering-based feature losses, are inherently ill-suited to this use case. A single 3D Gaussian models a spatial volume rather than an isolated point; once projected into an image, that Gaussian may cover many pixels. When the Gaussian is later reduced to a single 3D representative, typically its center $\mu_n$, many distinct pixels become associated with the same 3D location. The resulting many-to-one correspondence pattern is geometrically inconsistent with the assumptions of PnP.

The ambiguity is not merely representational. PnP with RANSAC depends on a set of mutually consistent point-to-point correspondences. When multiple pixels that originate from different surfaces or depths are mapped to one Gaussian center, reprojection errors increase, the inlier set becomes less coherent, and RANSAC may reject otherwise plausible hypotheses, require substantially more iterations, or converge to suboptimal poses. The paper further identifies a second source of degradation: photometric 3DGS optimization is free to introduce Gaussians that explain view-dependent effects or localized appearance variations, even when these Gaussians contribute only in a few training views and exhibit weak multi-view support. Such Gaussians are useful for rendering fidelity but are poorly aligned with the requirements of stable novel-view localization.

A common misconception is that a photometrically strong GSFF is automatically suitable for direct 2D–3D localization. SplitGS-Loc explicitly argues against that view. Its contribution lies in treating localization as a distinct design target rather than as a by-product of photometric scene representation.

## 2. 3D Gaussian Splatting and feature-field construction

A 3DGS scene is represented as a set of Gaussians
$$
\mathcal{G} = \{ g_n \}_{n=1}^{N_g},
$$
where each Gaussian $g_n$ has a mean $\mu_n \in \mathbb{R}^3$, covariance
$$
\Sigma_n = R_n S_n S_n^\top R_n^\top,
$$
with rotation $R_n$ and diagonal scale matrix $S_n = \mathrm{diag}(s_x, s_y, s_z)$, opacity $\alpha_n$, and spherical-harmonic color coefficients $\mathbf{c}_n$ [2605.07351].

For pixel $j$ in image $i$, rendered color is computed as
$$
\hat{I}_{ij} = \sum_n w(g_n, I_{ij})\, \mathbf{c}_n,
$$
where $w(g_n, I_{ij})$ is the composition weight. This weight reflects both the projected Gaussian footprint in the image plane and the accumulated transmittance along the ray; it is analogous to volumetric alpha compositing weights. GSFFs augment each Gaussian with a feature vector $\mathbf{z}_n$, giving a rendered feature
$$
\hat{\mathbf{f}}_{ij} = \sum_n w(g_n, I_{ij})\, \mathbf{z}_n.
$$

In conventional GSFF construction, $\mathbf{z}_n$ is often optimized with a rendering loss of the form
$$
\mathcal{L}_{\text{feat}} = \sum_{i,j} \|\hat{\mathbf{f}}_{ij} - \mathbf{f}_{ij}\|^2,
$$
where $\mathbf{f}_{ij}$ are 2D image features from an encoder. SplitGS-Loc identifies a limitation of this strategy: gradients from a ray are distributed across multiple Gaussians, which leads to feature blurring and weakens per-Gaussian distinctiveness. The method therefore avoids rendering-based feature optimization and instead performs direct feature lifting from multi-view pixel descriptors, so that descriptors attached to Gaussians remain aligned with strong pixel-level evidence rather than with a diffuse reconstruction objective [2605.07351].

This design decision is conceptually important. It reframes the feature field not as an auxiliary rendering layer but as a localization-oriented structure whose descriptors must support correspondence extraction under PnP.

## 3. Mixture-of-Gaussians-based splitting

The defining mechanism of SplitGS-Loc is Mixture-of-Gaussians-based splitting. The objective is to replace ambiguous many-to-one pixel-to-point mappings with more precise one-to-one correspondences by decomposing each Gaussian into several smaller Gaussians. If a single large Gaussian spans an extended volumetric region, approximating that region with multiple finer Gaussians reduces the number of pixels that each component significantly influences and tightens the geometric interpretation of each 2D–3D match [2605.07351].

The method begins from a Gaussian distribution
$$
p(\mathbf{x}) = \mathcal{N}(\mu, \Sigma)
$$
and replaces it by a mixture
$$
p(\mathbf{x}) = \sum_{k=1}^{K} \lambda_k \mathcal{N}(\mathbf{x}\mid \mu_k, \Sigma_k),
$$
which in the 3DGS setting corresponds to child Gaussians with opacities $\lambda_k \alpha$. To make the construction tractable, the derivation factorizes the Gaussian along its principal axes and performs splitting along the major axis, that is, the axis with the largest scale. In canonical coordinates, the scalar coordinate along this axis is modeled as
$$
\tilde{p}(x) \sim \mathcal{N}(0, s^2),
$$
and approximated by a symmetric three-component mixture
$$
\tilde{q}(x) = \lambda q_{-}(x) + \lambda_o q_o(x) + \lambda q_{+}(x),
$$
with
$$
q_{-}(x) \sim \mathcal{N}(-\beta s, \sigma^2),\quad
q_o(x) \sim \mathcal{N}(0, \sigma_o^2),\quad
q_{+}(x) \sim \mathcal{N}(+\beta s, \sigma^2),
$$
subject to $2\lambda + \lambda_o = 1$.

The parameters are determined by moment matching up to the fourth order, using
$$
\mathbb{E}[Z^2] = v + m^2,\qquad \mathbb{E}[Z^4] = 3v^2 + 6vm^2 + m^4
$$
for $Z \sim \mathcal{N}(m,v)$. Under the balanced-splitting constraint $\sigma_o = \sigma$, the derived solution is
$$
\lambda_o = \frac{2}{3},\qquad \lambda = \frac{1}{6},\qquad
\sigma^2 = \sigma_o^2 = \left(1 - \frac{\beta^2}{3}\right)s^2.
$$
Positive semi-definiteness requires
$$
1 - \frac{\beta^2}{3} > 0 \quad \Rightarrow \quad \beta \in (0,\sqrt{3}),
$$
and the implementation chooses $\beta = 1.4$, which the paper reports as empirically best for localization performance.

In 3D, the parent Gaussian is split along its major axis into three child Gaussians located at offsets $\pm \beta s$ and $0$ in canonical coordinates, with reduced variance along that axis, unchanged covariance components along the other axes, and opacities redistributed as $\lambda\alpha$, $\lambda_o\alpha$, and $\lambda\alpha$. Rotating the construction back to world coordinates via the original rotation $R$ yields three child Gaussians that preserve the overall volumetric distribution while narrowing the individual projected footprints. The paper reports that localization accuracy improves monotonically as $\beta$ increases within the feasible interval, which it interprets as evidence that wider separation of child components better reduces pixel–Gaussian ambiguity [2605.07351].

The significance of the splitting step is not simply densification. The method’s own ablations distinguish Mixture-of-Gaussians splitting from geometry-only densification, showing that nearest-neighbor upsampling remains inferior. This suggests that covariance-aware decomposition, rather than an arbitrary increase in point count, is the relevant mechanism.

## 4. Composition weights, salient Gaussian selection, and feature lifting

SplitGS-Loc builds on PlugGS-Loc, a training-free feature-lifting pipeline that uses composition weights from GS rasterization to select Gaussians that contribute significantly and consistently across views and to aggregate discriminative multi-view features to those Gaussians [2605.07351]. The key observation is that composition weights encode not only projected support but also opacity and transmittance, making them an occlusion-aware indicator of how strongly a Gaussian explains a pixel.

For each training image $I_i$ and pixel $I_{ij}$, the method defines the set of primary Gaussians as
$$
\hat{\mathcal{G}}_{ij} = \{ g_n \in \mathcal{G} \mid w(g_n, I_{ij}) \ge \tau \},
$$
with a fixed threshold $\tau = 0.1$. This eliminates weak pixel–Gaussian associations and retains only contributors with strong explanatory power. Rather than storing all pixel-level weights, the method compresses per-image contribution into a representative image–Gaussian weight
$$
\tilde{w}(g_n, I_i) :=
\begin{cases}
\max_j \big(w(g_n, I_{ij})\big), & \text{if } \cup_j \hat{\mathcal{G}}_{ij} \neq \emptyset,\\[4pt]
0, & \text{otherwise,}
\end{cases}
$$
and for each Gaussian collects the informative nonzero weights as
$$
\mathcal{W}_n = \{\tilde{w}(g_n, I_i) \mid \tilde{w}(g_n, I_i) > 0\}.
$$
Gaussians with $\mathcal{W}_n = \emptyset$ are discarded.

Importance is then scored by
$$
\psi(g_n) = \frac{1}{|\mathcal{W}_n|} \sum_{\tilde{w}\in \mathcal{W}_n} \tilde{w},
$$
which rewards high visibility, high opacity or transmittance, and multi-view consistency. To keep the map compact while preserving spatial coverage, Gaussians are grouped into local regions by anchor sampling and $k$-nearest-neighbor grouping in 3D, and the most important Gaussian in each region is retained. This produces a set of anchor Gaussians that are uniformly distributed and strongly supported across views.

Feature construction proceeds by direct multi-view aggregation. Given image features $\mathbf{f}_{ij}$ from an encoder—SuperPoint in the experiments—the method defines
$$
\mathcal{P}_n = \{(i,j^*) \mid w(g_n, I_{ij^*}) = \tilde{w}(g_n, I_i)\},
$$
namely the set of pixels where Gaussian $g_n$ attains maximal weight in each contributing image. With softmax-normalized weights $\bar{w}_{ij^*}$ over $\mathcal{W}_n$ and channel-normalized pixel descriptors $\bar{\mathbf{f}}_{ij^*}$, the Gaussian feature is
$$
\bar{\mathbf{z}}_n = \sum_{(i,j^*)\in \mathcal{P}_n} \bar{w}_{ij^*}\, \bar{\mathbf{f}}_{ij^*}.
$$

This aggregation strategy is central to SplitGS-Loc’s localization behavior. It uses strong pixel–Gaussian associations as soft assignments, avoids rendering-loss optimization entirely, and thereby prevents the feature blurring associated with rendering-based GSFF learning. The paper’s Gaussian-attribute ablation supports this interpretation: a simple Projection+Average strategy, which ignores composition weights, yields markedly worse localization than the weight-aware construction, indicating that opacity- and transmittance-aware selection is a primary source of robustness [2605.07351].

## 5. End-to-end construction pipeline and PnP-based localization

The full SplitGS-Loc pipeline begins from a standard RGB 3DGS model trained for photometric reconstruction only. That base training takes approximately $6$ minutes per scene. The additional localization-specific construction comprises Mixture-of-Gaussians splitting, pre-filtering, weight aggregation, Gaussian sampling, and feature registration, with total runtime of approximately $64$ seconds per scene; the breakdown is Gaussian Splitting $0.4$ s, Pre-filtering $10.6$ s, Weight Aggregation $9.8$ s, Gaussian Sampling $10.2$ s, and Feature Registration $33.4$ s [2605.07351].

After splitting, sampling must remain aware of parent–child structure. Scoring each child independently and then collapsing to a single representative would negate the benefit of decomposition. The method therefore computes per-child scores and aggregates them into a parent-level quantity
$$
\bar{\psi}(g_n) = \frac{1}{3}\bigl(\psi(g_{n_-}) + \psi(g_{n_o}) + \psi(g_{n_+})\bigr),
$$
then performs spatial grouping and sampling on parent centers $\mu_n$ using $\bar{\psi}(g_n)$, while retaining all three child Gaussians for each selected parent. This preserves the local multi-Gaussian structure that splitting introduced.

At query time, SuperPoint extracts 2D keypoints and descriptors from the input image. Each retained Gaussian is treated as a 3D point at its center $\mu_n$ with descriptor $\bar{\mathbf{z}}_n$, and query descriptors are matched to Gaussian descriptors, following STDLoc’s matching and PoseLib settings. The resulting correspondences have the form
$$
\mathbf{u}_k \leftrightarrow \mathbf{X}_k,\qquad \mathbf{X}_k = \mu_{n(k)}.
$$
PoseLib is then used for PnP with RANSAC, conceptually minimizing
$$
\min_{\mathbf{R},\mathbf{t}} \sum_k \rho\big(\|\pi(\mathbf{R}\mathbf{X}_k + \mathbf{t}) - \mathbf{u}_k\|\big),
$$
where $\rho$ is the inlier-robust loss induced by RANSAC.

The empirical claim of the method is that splitting materially alters the correspondence distribution seen by PnP. On Cambridge, the number of many-to-one correspondences per PnP is reduced from approximately $313$ to $199$. Under different RANSAC settings, the number of valid PnP inliers increases from $362$ to $447$ or from $367$ to $450$ when switching from STDLoc to SplitGS-Loc. Median translation error improves from $11.2$ cm to $8.7$ cm, and the method remains stable under more efficient RANSAC limits: STDLoc degrades from $11.2$ cm to $12.8$ cm, whereas SplitGS-Loc remains at approximately $8.7$–$8.9$ cm [2605.07351].

Query-time performance is correspondingly efficient. With a moderate RANSAC setting of $10^2$–$10^3$ iterations, SplitGS-Loc reaches $8.9$ cm / $0.15^\circ$ in about $55$ ms per query on Cambridge. With a higher-accuracy setting of $10^3$–$10^5$ iterations, it reaches $8.7$ cm / $0.15^\circ$ in about $215$ ms. The method’s pose estimation pipeline therefore remains direct: no dense refinement is required to achieve the reported state of the art, although an optional GS-CPR refinement produces the SplitGS-Loc$^+$ variant on 7Scenes [2605.07351].

## 6. Empirical performance, ablations, and limitations

Evaluation is reported on Cambridge Landmarks, comprising five outdoor scenes up to approximately $7600\ \mathrm{m}^2$, and on 7Scenes, comprising seven indoor scenes spanning $1$–$18\ \mathrm{m}^2$. Metrics are median translation error in centimeters and rotation error in degrees per scene, using SfM pseudo-ground truth from Brachmann et al. 2021 rather than the original D-SLAM ground truth [2605.07351].

| Dataset / variant | Map size | Average error |
|---|---:|---:|
| Cambridge Landmarks, SplitGS-Loc | 97 MB | **8.7 cm / 0.15°** |
| 7Scenes, SplitGS-Loc | 88 MB | **1.0 cm / 0.26°** |
| 7Scenes, SplitGS-Loc$^+$ | ~175 MB | **0.6 cm / 0.19°** |

On Cambridge Landmarks, SplitGS-Loc uses a mapping time of $7$ minutes, consisting of $6$ minutes for 3DGS training and approximately $1$ minute for GSFF construction. The per-scene results are GreatCourt $9.3$ cm / $0.06^\circ$, King’s College $14.1$ cm / $0.17^\circ$, Old Hospital $13.8$ cm / $0.29^\circ$, Shop Facade $2.5$ cm / $0.12^\circ$, and St. Mary’s Church $3.8$ cm / $0.13^\circ$. The paper compares these results against classical feature matching, scene coordinate regression, NeRF-based localization, and GS-based methods including GSplatLoc, GSFFs-PR, STDLoc, GSVisLoc, RAP, GS-RelocNet, and PlugGS-Loc, and reports that SplitGS-Loc is competitive with or better than both training-free and per-scene-trained alternatives while avoiding hours of mapping time and dense pose refinement in its base form [2605.07351].

On 7Scenes, the base SplitGS-Loc variant also uses a $7$ minute mapping time and achieves $1.0$ cm / $0.26^\circ$ with an $88$ MB map. The refined SplitGS-Loc$^+$ variant reaches $0.6$ cm / $0.19^\circ$ with a map size of approximately $175$ MB, matching or surpassing the best per-scene-trained GS and NeRF methods listed in the comparison. A plausible implication is that the correspondence quality produced by the base feature field is already sufficiently strong that refinement acts as an optional enhancer rather than as a structural necessity.

The ablation study clarifies the source of the gains. Varying the splitting factor $\beta$ shows monotonically improving accuracy as $\beta$ increases within the feasible range, supporting the interpretation that wider separation of child Gaussians reduces correspondence ambiguity. Varying the primary-Gaussian threshold $\tau$ shows that values that are too low admit too many weak associations and worsen feature blurring, whereas values that are too high reduce coverage; the method fixes $\tau = 0.1$ and reports robustness in the range $0.05$–$0.1$. Component ablations report $13.1$ cm / $0.22^\circ$ for “STDLoc GSFFs + cosine similarity sampling,” $9.6$ cm / $0.17^\circ$ when replacing the GSFFs with PlugGS-Loc’s feature lifting under the same sampling, and $8.7$ cm / $0.15^\circ$ after adding SplitGS-Loc’s weight-based sampling. A separate attribute ablation reports $14.1$ cm / $0.26^\circ$ for Projection+Average, $10.6$ cm / $0.20^\circ$ for PlugGS-Loc, $10.0$ cm / $0.18^\circ$ for nearest-neighbor upsampling, and $8.7$ cm / $0.15^\circ$ for SplitGS-Loc. These results are presented as evidence that both composition weights and covariance-guided splitting are essential, and that increasing Gaussian count alone is insufficient [2605.07351].

The principal limitation identified by the paper is dependence on the quality of the underlying photometric 3DGS. Challenging scenes with severe illumination changes, dynamic content, or very sparse coverage in large-scale environments remain difficult when the base Gaussian scene fails to reconstruct reliably. SplitGS-Loc does not solve that upstream failure mode; it inherits the reconstruction limits of the photometric representation on which it is built. This suggests that future progress is likely to depend on advances in neural scene representations for difficult capture conditions as much as on further improvements in correspondence disambiguation itself [2605.07351].

Source: https://www.emergentmind.com/topics/splitgs-loc