---
title: Spherical Voronoi Representation
url: https://www.emergentmind.com/topics/spherical-voronoi-representation
type: topic
---

# Spherical Voronoi Representation

A spherical Voronoi representation decomposes the surface of the sphere into regions based on proximity or optimality conditions with respect to a finite set of generators (“sites”). This construct generalizes planar Voronoi diagrams to the geometry of the $2$-sphere, leveraging geodesic distance as the metric and providing foundational structure for combinatorics, directional data modeling, visualization, and non-Euclidean clustering.

## 1. Mathematical Definitions and Fundamental Properties

Let $S = \{ x \in \mathbb{R}^3 : \|x\| = r \}$ denote the sphere of radius $r > 0$. For a point set $U \subset S$ of $n$ elements in general position (no three colinear on a great circle, no four cocircular), the classical (order-1) spherical Voronoi diagram partitions $S$ into $n$ faces, each defined by
\[
V_i = \{ x \in S : d(x,s_i) \leq d(x,s_j) \;\; \forall j \neq i \}
\]
where $d(\cdot,\cdot)$ is the geodesic (angular) distance: $d(x,y) = \arccos(x \cdot y / r^2)$.

The notion extends to order-$k$ spherical Voronoi diagrams. For $1 \leq k \leq n-1$ and each $k$-subset $P_k \subset U$, the order-$k$ face is
\[
f(P_k) = \left\{ x \in S : \text{the set of $k$ nearest neighbors of $x$ in $U$ is exactly $P_k$} \right\}
\]
The complete order-$k$ diagram $SV_k(U)$ is the union of the boundaries of these faces. The region of site $i$ at order $k$ is
\[
R_k(i) = \left\{ x \in S : i \text{ is among the $k$ nearest neighbors of } x \right\} = \bigcup_{P_k \ni i} f(P_k)
\]
Weighted and centroidal generalizations replace the proximity criteria with additive site weights and optimal centroidal assignment conditions [0912.3974, 1705.02232].

## 2. Geometric Construction and Reduction via Inversion

An important geometric method relates spherical and planar Voronoi diagrams through 3D inversion. Fixing a pole $i \in S \setminus U$, inversion maps $S \setminus \{i\}$ to a plane $S' \cong \mathbb{R}^2$ via
\[
\phi(x) = i + \frac{R^2}{\|x - i\|^2}(x - i)
\]
Under this transformation, sites $U$ are mapped to $U'$ in the plane. Circles on $S$ not passing through $i$ invert to circles; those through $i$ become lines. The upshot is a decomposition
\[
SV_k(U) \cong V_k(U') \cup V_{n-k}(U')
\]
where $V_k(U')$ and $V_{n-k}(U')$ are the planar order-$k$ and order-$(n-k)$ diagrams, respectively, glued along unbounded edges associated to the same bisectors. This enables construction of $SV_k(U)$ in $O(T_{\mathrm{plane}} + n)$ time, where $T_{\mathrm{plane}}$ is the planar diagram construction cost [2207.13788].

## 3. Combinatorial Structure and Analytical Results

Order-$k$ spherical Voronoi diagrams admit precise combinatorial characterization. With $V_k$, $E_k$, $F_k$ denoting the numbers of vertices, edges, and faces:
\[
\begin{align*}
V_k &= 4kn - 4k^2 - 2n \\
E_k &= 6kn - 6k^2 - 3n \\
F_k &= 2kn - 2k^2 - n + 2
\end{align*}
\]
Euler’s relation $V_k - E_k + F_k = 2$ holds. Vertices are split into type I ($2k(n-k-1)$) and type II ($2(k-1)(n-k)$) [2207.13788].

In centroidal settings, a site is iteratively relocated to the centroid of its Voronoi cell using,
\[
C(V) = \frac{ \int_{x \in V} x\, dA(x) }{ \left\| \int_{x \in V} x\, dA(x) \right\| }
\]
where $dA(x)$ is the area element. The Weighted Spherical Centroidal Voronoi Tessellation (WSCVT) algorithm updates both site positions and weights to minimize area error against prescribed targets, achieving exact area allocation within arbitrary tolerance $\epsilon$ [0912.3974].

## 4. Spherical Voronoi as a Differentiable Partition for Directional Data

Spherical Voronoi diagrams serve as explicit, parameter-efficient, differentiable partitions of $S^2$ for modeling functions defined on the sphere, such as view-dependent appearance in 3D graphics. Using $K$ learnable sites $\{s_k\}$ with per-site parameters (including temperature $\tau_k$ and coefficient $c_k$), the soft region weights are computed as
\[
w_k(\omega; \tau) = \frac{ e^{\tau_k s_k \cdot \omega} }{ \sum_{j=1}^K e^{\tau_j s_j \cdot \omega} }
\]
The spherical function is then
\[
f_{\mathrm{SV}}(\omega; \{s_k, \tau_k, c_k\}) = \sum_{k=1}^K w_k(\omega; \tau_k) \cdot c_k
\]
This framework allows gradients to flow for site, temperature, and coefficient optimization, interpolating between uniform and hard-partitioning as $\tau_k$ varies. Region boundaries are smooth and adapt to local function complexity, avoiding artifacts such as SH ringing and capturing high-frequency content robustly [2512.14180].

## 5. Variants: Weighted, Centroidal, and Generalized Spherical Voronoi Tessellations

Weighted spherical Voronoi diagrams alter the assignment criterion to
\[
V_i = \{ x \in S^2 : d(x, s_i)^2 - w_i \leq d(x, s_j)^2 - w_j \;\; \forall j \neq i \}
\]
allowing regions of controlled area or density [0912.3974].

Centroidal variants (WSCVT) demand each site lie at its cell’s centroid, iteratively adjusting positions and weights to match area allocations. This achieves lower variance in region sizes and is essential for high-quality spherical data visualization and partitioning in non-uniform distributions [0912.3974].

Generalized Voronoi diagrams, as in the context of Spherical Wards clustering, use arbitrary symmetric dissimilarity measures and optimize partitions for within-cluster scatter adjusted by cardinalities. Assignments are based on cost differentials for cluster membership, and the resulting regions form a Voronoi-like diagram under non-Euclidean metrics [1705.02232].

## 6. Algorithmic Realizations and Computational Aspects

Classic spherical Voronoi diagrams can be constructed via inversion and planar diagram computation. For weighted/centroidal settings:
- Build the weighted Voronoi diagram (by geodesic distances minus weights)
- Compute actual cell area fractions $a_i$
- Move each site to the centroid of its current cell, update weights to reduce area error
- Iterate until convergence within tolerance $\epsilon$

For differentiable partitions (as in SV for 3DGS), forward and backward passes involve dot product calculation, temperature scaling, softmax evaluation, and parameter updates via chain rule with normalization constraints for sites [2512.14180].

Performance of such algorithms depends on the number of sites, mesh resolution for area approximation, and softness parameter $\tau$. For clustering-based Voronoi diagrams, Hartigan-style iterative reassignment with efficient updates of within-cluster scatter is employed, followed by pruning of insufficiently populated clusters [1705.02232].

## 7. Applications and Impact

Spherical Voronoi representations are foundational in:
- Computational geometry for the analysis of distance-driven partitions on $S^2$, higher-order combinatorics, and geometric optimization [2207.13788]
- Appearance modeling for explicit 3D representations, providing a differentiable, interpretable, and highly expressive partition for view-dependent effects, outperforming global (SH) and local (SG) bases in task accuracy and efficiency [2512.14180]
- Data visualization: WSCVT enables spatially balanced spherical layouts for hierarchical data with exact control over region sizes, outperforming uniform subdivision and minimizing area waste [0912.3974]
- Spherical clustering, integrating geometric, probabilistic, and variational approaches for robust, scale-invariant partitioning under non-Euclidean metrics [1705.02232]

By combining analytic geometric structure with optimization and differentiable programming, spherical Voronoi models extend far beyond pure geometry into data-driven clustering, graphics, and information visualization.

Source: https://www.emergentmind.com/topics/spherical-voronoi-representation