---
title: Generalized Balanced Power Diagrams
url: https://www.emergentmind.com/topics/generalized-balanced-power-diagrams-gbpds
type: topic
---

# Generalized Balanced Power Diagrams

Generalized Balanced Power Diagrams (GBPDs) constitute a geometric and optimization-based extension of classical tessellation models, unifying Voronoi, Laguerre (power), and anisotropic diagrams under a common framework. GBPDs are widely applied in the modeling of microstructures (e.g., polycrystals), image segmentation, clustering, and spatial resource allocation, where control over partition cell size, anisotropy, and boundary curvature are essential. Key features include flexible metric tensors, additive weights (balancing parameters), nonlinear cell boundaries, and the capacity to enforce or optimize for prescribed cell measures such as volumes, areas, or counts.

## 1. Mathematical Foundations and Definition

Let $\psi = \{ (s_i, M_i, w_i) : s_i \in \mathbb{R}^d, M_i \succ 0, w_i \in \mathbb{R} \}_{i=1}^n$ denote a finite set of generators, where $s_i$ is the site (center), $M_i$ is a symmetric positive-definite metric tensor, and $w_i$ is the balancing weight. For any $x \in \mathbb{R}^d$, the generalized (anisotropic power) distance to generator $i$ is defined as
\[
\operatorname{dist}_G(x; s_i, M_i, w_i) = (x-s_i)^\top M_i (x-s_i) - w_i.
\]
The $i$-th GBPD cell is the region
\[
C_i = \left\{ x \in \mathbb{R}^d : \operatorname{dist}_G(x; s_i, M_i, w_i) \leq \operatorname{dist}_G(x; s_j, M_j, w_j),\ \forall j \right\}.
\]
Special cases:
- $M_i \equiv I,\ w_i \equiv 0$ yields the ordinary Voronoi diagram.
- $M_i \equiv I$ and general $w_i$ gives the Laguerre (power) diagram.
- $w_i$ equal and general $M_i$ yields the anisotropic Voronoi diagram [2601.18593][2012.11772][1411.4535][2303.15275].

The boundaries between GBPD cells correspond to quadric surfaces; in 2D, the bisector between two generators is a conic (ellipse, hyperbola, or degenerate line) and in higher dimensions, a quadric hypersurface [2303.15275][2601.18593]. The defining equations guarantee that any uniform shift in weights or common scaling of $(M_i, w_i)$ leaves the tessellation unchanged.

## 2. Optimization Formulations and Balanced Assignments

A central property of GBPDs is their relation to balanced least-squares assignments. Let $P = \{y_1, \dots, y_m\}$ be a set of points (e.g., voxels, residents, pixels). Given $k$ sites and quota vector $\{\mu_i\}_{i=1}^k$ (cell sizes or measures), the classic optimization is
\[
\min_{f: P \to [k]} \sum_{y \in P} \| y - s_{f(y)} \|_{M_{f(y)}}^2
\quad \text{subject to} \quad |f^{-1}(i)| = \mu_i.
\]
Alternatively, for pixel/voxel assignments: $\xi_{ij}$ are indicator variables, each point assigned to one site, each site receives prescribed measure, and
\[
\min_{\xi \geq 0} \sum_{i=1}^k \sum_{j=1}^m \xi_{ij} \|y_j - s_i\|_{M_i}^2
\]
with constraints $\sum_i \xi_{ij} = 1$ (unique assignment), $\sum_j \xi_{ij} = \mu_i$ (balancing), $\xi_{ij} \geq 0$ [2012.11772][2507.14268][1411.4535][1710.03358]. The Lagrange dual provides the weights $w_i$ enforcing these constraints.

In spatial allocation, such as electoral districting, the GBPD constructs districts whose populations differ by at most one, each defined as the intersection of a polygonal (power cell) region with a geographic domain [1710.03358].

For clustering and outlier detection, the GBPD framework—via LPs with margin slacks—can handle soft assignments and quantify departures from balanced assignments, providing optimality certificates on the number of margin errors and support vectors [1307.3949].

## 3. Analytical and Algorithmic Aspects

### Boundary Structure
Between two generators $(s_i, M_i, w_i)$, $(s_j, M_j, w_j)$, the cell boundary is the locus
\[
(x-s_i)^\top M_i (x-s_i) - w_i = (x-s_j)^\top M_j (x-s_j) - w_j,
\]
which is a conic in 2D (general quadric in higher dimensions) [2303.15275][2601.18593]. In the isotropic case or when $M_i = M_j$, this reduces to planar boundaries; otherwise, faces can be curved.

### Affine Invariance and Sections
GBPDs are closed under invertible affine maps. For $x \mapsto Ax + b$, the map is realized as $s_i \mapsto As_i + b$, $M_i \mapsto A^{-\top} M_i A^{-1}$, $w_i$ unchanged; other special cases (translation, rotation, scaling) are subsumed [2601.18593]. Flat sections (intersections with affine subspaces) of GBPDs yield lower-dimensional GBPDs with modified parameters.

### Algorithmic Construction
Discrete GBPDs (over gridded domains) can be constructed by minimizing the assignment or clustering LP, with strong duality ensuring optimal balancing [2012.11772][1411.4535][1710.03358]. Efficient algorithms include:
- Capacitated Lloyd–style iteration: alternating balanced assignment (min-cost flow/LP) and centroid update, guaranteeing monotonic decrease of the k-means cost and finite termination [1710.03358].
- Improved rendering algorithm: two-stage assignment over large grids achieves $O(N\log n)$ complexity under stochastic generator configurations [2601.18593].
- Analytical construction of 2D diagrams: full vertex-edge incidence and topology via explicit conic intersection routines [2303.15275].

## 4. Fitting GBPDs to Data and Methodological Trade-Offs

GBPDs are systematically fitted to high-dimensional structural data (e.g., segmented 3D images of polycrystals) by optimizing generator parameters to maximize agreement with ground-truth grain partitions. Canonical methods include:
- **Linear Programming:** fixing $(s_i, M_i)$ (e.g., by PCA or first and second moments) and solving for $w_i$ subject to exact or bounded volume constraints, yielding globally optimal assignments and integral solutions when the LP matrix is totally unimodular [2507.14268][1411.4535].
- **Gradient Descent (Softmax Relaxation):** simultaneously optimizing all generator parameters over smooth surrogate losses; efficient with GPU acceleration, sensitive to initialization [2507.14268].
- **Stochastic Optimization (Cross-Entropy):** treating generator parameters as random variables and minimizing interface-based discrepancy; robust to local minima but computationally costly [2507.14268].
- **Derivative-Free Local Optimization:** e.g., Subplex, Praxis methods directly fitting generator parameters by minimizing point-to-boundary discrepancies [2507.14268].
  
The choice of method depends on balancing model and algorithmic complexity, convergence guarantees, fidelity to volume/area/topology, and computational resources.

A consolidated table of main fitting methods [2507.14268]:

| Method                | Variables             | Optimality | Constraints/Focus         |
|-----------------------|----------------------|------------|--------------------------|
| LP (weighted assign.) | $w_i$ (weights)      | Global     | Exact/bounded volumes    |
| GD (softmax)          | $(s_i, M_i, w_i)$    | Local      | Volume or shape loss     |
| Cross-entropy (CE)    | $(s_i, w_i)$ or all  | Heuristic  | Interface placement      |
| Neper/DFO             | $(s_i, M_i, w_i)$    | Local      | Boundary agreement       |

## 5. Analytical Representation and Cell Topology

Recent advances include analytic representations of 2D GBPDs, allowing for exact computation of cell boundaries, vertices, and adjacency relations without pixelization [2303.15275]. The bisector between two sites is a general conic $\mathcal{E}_{ij}(x, y) = (x-c_i)^\top A_i (x-c_i) - w_i - (x-c_j)^\top A_j (x-c_j) + w_j = 0$, classified by the inertia of the associated matrix. Vertices are obtained as solutions to pairs of such quadratic equations; edge arcs are assembled by sorting intersection parameters and validating local cell inclusions. This framework enables area and perimeter evaluation via closed-form line integrals, accurate topology analysis, and robust handling of degenerate cases such as lenses and unique-neighbor cells.

## 6. Applications and Empirical Performance

### Polycrystal Modeling
GBPDs efficiently reconstruct polycrystalline grain maps from minimal statistics (center-of-mass, volume, second moments), achieving high accuracy in voxel labeling and topological neighborhood reconstruction, even under anisotropic conditions. In practical tests, GBPDs consistently outperform heuristic Laguerre tessellations and power diagrams on labeling rates and neighborhood accuracy, with only a handful of parameters per grain needed [1411.4535][2507.14268].

### Superpixel Image Segmentation
Power-SLIC applies GBPDs to superpixel segmentation, generating piecewise-quadratic cell boundaries with exact area balancing, competitive boundary recall, undersegmentation error, and robustness to noise. Continuous diagram representations ensure resolution-independence and efficient compression [2012.11772].

### Redistricting and Spatial Allocation
GBPD-based algorithms ensure perfectly balanced, compact, and contiguous partitioning of geographic domains, with cells given by truncations of convex polygons (or more general quadric cells) to the domain boundary. Iterative capacitated Lloyd–style algorithms provide local minima with provable compactness and side count bounds [1710.03358].

### Clustering and Outlier Detection
GBPD formulations enable efficient linear programming algorithms for clustering with balanced assignments, margin maximization, and explicit control over misclassification and support vector counts, facilitating robust outlier detection [1307.3949].

## 7. Theoretical Properties and Invariance

GBPDs inherit several crucial invariance properties:
- **Affine Covariance:** Closed under invertible affine maps; directions and anisotropies transform appropriately, preserving diagram topology [2601.18593].
- **Sectional Stability:** Slices through GBPDs yield lower-dimensional GBPDs, enabling detailed study of microstructure cross-sections or embedding problems [2601.18593].
- **Parameter Redundancy:** Uniform weight shifts and common scaling of $(M_i, w_i)$ leave the diagram unchanged, implying redundancy in parameter specification [2601.18593].
- **Compactness and Cell Topology:** The set-theoretic and geometrical structure of GBPDs restricts the average number of sides (in planar domains) and maintains centroidality between generator positions and resulting cell mass centers when optimization is used [1710.03358][1411.4535].

Robust theoretical guarantees exist in special cases (e.g., volume-matching LPs, support vector counts in soft-margins), while in general, practical optimization is subject to complexity and local minima considerations [2507.14268][1307.3949].

---

**References:**  
[1710.03358]  
[2012.11772]  
[1411.4535]  
[1307.3949]  
[2601.18593]  
[2507.14268]  
[2303.15275]

Source: https://www.emergentmind.com/topics/generalized-balanced-power-diagrams-gbpds