---
title: Hierarchical Volume-Preserving Maps
url: https://www.emergentmind.com/topics/hierarchical-volume-preserving-maps
type: topic
---

# Hierarchical Volume-Preserving Maps

A hierarchical volume-preserving map is a smooth, bijective reparameterization of a structured space (such as high-dimensional data or geometric solids) into multiple scales or hierarchical components, with the property that the transformation preserves local and global volume. The defining characteristic is that the determinant of the Jacobian matrix at every point is identically one, ensuring that the original density structure (in data or space) is not locally contracted or expanded. This property enables exact likelihood evaluation and geometrically faithful discretizations in applications ranging from generative modeling to geometric grid construction [2501.06999, 1504.01517].

## 1. Formal Definition and Mathematical Properties

Let $\mathcal{X}$ denote the original space (e.g., $\mathbb{R}^d$ for imaging), and suppose $x\in\mathcal{X}$ is mapped into $S$ hierarchical scales $(z^{(1)},\ldots,z^{(S)})$, where each $z^{(s)}$ belongs to $\mathcal{Z}^{(s)}$. A hierarchical volume-preserving map is a $C^1$-diffeomorphism $f: \mathcal{X} \to \mathcal{Z}^{(1)} \times \cdots \times \mathcal{Z}^{(S)}$ satisfying
$$
\det\big[\nabla_x f(x)\big] = 1\quad\forall\,x\in\mathcal{X} \ .
$$
Invertibility is required so that reconstruction of $x$ from $(z^{(1)},...,z^{(S)})$ is always possible. The volume preservation guarantees that, in probabilistic models, the change-of-variables formula simplifies:
$$
p(x) = p(z^{(1)},...,z^{(S)})\cdot|\det J_x|^{-1} = p(z^{(1)},...,z^{(S)}) \ .
$$
This transformation introduces no local density distortion, facilitating tractable likelihoods and robust geometric discretizations [2501.06999].

## 2. Concrete Realizations: Laplacian Pyramid and Haar Wavelet

Two canonical examples illustrate the concept in data-analytic settings:

**Laplacian Pyramid:**  
Let $D: \mathbb{R}^{j\times j}\to \mathbb{R}^{j/2\times j/2}$ be a norm-preserving downsampler (e.g., bilinear with scaling), and $U$ the corresponding upsampler. Recursively decompose $y^{(s)}=D(y^{(s+1)})$ and $z^{(s+1)}=y^{(s+1)}-U(D(y^{(s+1)}))$ with base $z^{(1)}=D(y^{(2)})$ and inverse $y^{(s+1)}=U(y^{(s)})+z^{(s+1)}$ up to $x=y^{(S)}$. Because $D$ and $U$ are norm-preserving, the overall transform is a tight frame: it satisfies Parseval’s identity and all singular values of its Jacobian equal one. Thus, $\det J_f=1$ [2501.06999].

**Haar Wavelet Transform:**  
Decompose the input using orthonormal filters $L=[1/\sqrt{2},1/\sqrt{2}], H=[-1/\sqrt{2},1/\sqrt{2}]$ into subbands at each scale via convolutions and striding. The transform matrix $W$ is orthonormal ($W^\top W=I$), implying $|\det J_W|=1$. The multiscale (hierarchical) transform is exactly volume-preserving and distance-preserving [2501.06999].

| Map Type           | Construction Principle               | Volume Preservation Mechanism      |
|--------------------|-------------------------------------|------------------------------------|
| Laplacian Pyramid  | Norm-preserving linear tight frame   | Parseval’s identity $\Rightarrow$ $\det J_f=1$ |
| Haar Wavelet       | Orthonormal matrix decomposition     | $\det J_W = \pm1$ (orientation set to $+1$)   |

## 3. Geometric Instances on Polyhedral Domains

In geometric contexts, area- and volume-preserving maps have been constructed for families of convex polyhedra $\mathbb{K}_n(r, \varepsilon)$, defined as the union of a right prism of height $2\varepsilon r$ over a regular $n$-gon, plus two congruent pyramidal caps. The surface $\partial\mathbb{K}_n$ can be mapped bijectively and area-preservingly to a sphere $\mathbb{S}^2(r)$ by explicit algebraic and trigonometric formulas, constructed zone-wise with area-preserving Jacobians. For a corresponding family of solids $\overline{\mathbb{K}}_n(r', \varepsilon)$, a volume-preserving homeomorphism onto the solid ball $\overline{\mathbb{S}^2(r)}$ can be constructed if and only if $(n, \varepsilon)$ satisfy explicit algebraic constraints on the Jacobian; exact homeomorphisms exist for specific $n=3,4,5$ [1504.01517].

## 4. Hierarchical Grid Construction and Refinable Partitions

Hierarchical volume-preserving maps naturally permit construction of uniform, refinable, bijective grids for both analytic and geometric applications.

- **Surface Grids:** Uniform partitions of each pyramid face in $\mathbb{K}_n$ yield $n\cdot p$ quadrilaterals of equal area, extendable to HEALPix-type iso-latitude pixelations for $n=4$, $r=1$, and appropriate $\varepsilon=p/(p+1)$.
- **Solid Grids:** Subdivision of $\overline{\mathbb{K}}_n$ into $4n$ tetrahedra, recursively refining each into four, yields tetrahedral partitions with uniform volume at each refinement level: $N_k = 4n\cdot 4^k$, each of volume $V_0/4^k$. Volume-preserving mapping to the ball ensures exact volumetric discretization [1504.01517].

| Partition Type   | Construction                      | Properties                  |
|------------------|-----------------------------------|-----------------------------|
| Surface          | Equiareal quadrilaterals          | Uniform area, hierarchical  |
| Solid            | Recursive tetrahedral refinement   | Uniform volume, refinable   |

## 5. Cascaded Multi-Scale Models and Likelihood Tractability

In probabilistic modeling, especially cascaded (multi-scale) diffusion models, hierarchical volume-preserving maps enable exact, tractable likelihood evaluation. Classical cascaded models suffer from intractable marginalization over extraneous latent scales. Under a hierarchical volume-preserving reparameterization $f(x) = (z^{(1)}, \dots, z^{(S)})$, the likelihood satisfies
$$
p_\theta(x) = p_\theta(z^{(1)},...,z^{(S)})
$$
without the need for Jacobian correction. The model can then be trained using exact diffusion ELBOs on each scale, with the total objective decomposed as a sum of ELBOs for the unconditional base scale and conditional super-resolution scales [2501.06999].

## 6. Theoretical Guarantees and Connections to Optimal Transport

Hierarchical volume-preserving maps provide strict area- or volume-preservation: the local area-distortion factor is identically one in all zones, and similarly the Jacobian for 3-volume is one under the solid domain mappings. In data-analytic cascaded models, the training objective under such maps is provably equivalent (up to a constant) to minimization of a weighted sum of Earth Mover’s Distances (EMD/Wasserstein-p metrics) between the true and learned score fields:
$$
\alpha\,\mathcal{C}(x) \geq \sum_{k=1}^{T-1} w_k\, \mathbb{E}_{q(x_k|x)}\big[ W_p(\nabla_x \log q(x_k|x), s_\theta(x_k, k)) \big],\quad 0<p\leq1,
$$
where hierarchical transforms such as Laplacian pyramids or Haar wavelets admit linear-time upper bounds on the EMD by summing $\ell_2$-norms of the hierarchical coefficients. Thus, volume preservation not only makes the likelihood tractable but also effects a perceptually meaningful score matching in generative modeling [2501.06999].

## 7. Special Cases and Integration with Established Schemes

The mapping framework subsumes established equal-area discretization schemes such as HEALPix by explicit construction: for $n=4$, $r=1$, and $\varepsilon=p/(p+1)$, the quadrilateral partition of the mapped polyhedron reproduces the iso-latitude, equal-area pixelization standard in cosmological data analysis, with the mapping formulas yielding the HEALPix parametrization on each cap [1504.01517]. This construction demonstrates the unified applicability of hierarchical volume-preserving maps in both geometric grid generation and modern generative modeling architectures.

Source: https://www.emergentmind.com/topics/hierarchical-volume-preserving-maps