Efficient Compression of Structured and Unstructured Volumes via Learned 3D Gaussian Representation
Landon Dyken, Sharmistha Chakrabarti, Nathan Debardeleben, Steve Petruzza, Qi Wu, Will Usher, and Sidharth Kumar
ArXiv Preprint, 2026

W-VEG represents a volume as a compact set of learned 3D Gaussians that can be queried directly in world space. Unlike prior implicit models for unstructured data, W-VEG+S can encode both the scalar field and the shape of its valid domain, removing the need to retain an exterior mesh.

Abstract
Recent work has shown that implicit neural representations can effectively compress structured and unstructured volume data, allowing direct data querying with a reduced memory footprint. Existing INR methods for unstructured volumes, however, do not encode geometry and therefore require partial mesh storage for later sampling, limiting the achievable compression ratio. We introduce an explicit volume compression model based on learned 3D Gaussian primitives. Collections of 3D Gaussians are reinterpreted as a representation of a volume scalar field, with scalar values reconstructed at arbitrary spatial locations through weighted aggregation of intersecting Gaussians. CUDA-accelerated pipelines support structured and unstructured sampling, continuity-aware losses, and sampling-error-based densification. By encoding both scalar values and domain geometry, the representation removes mesh storage for unstructured volumes while maintaining competitive reconstruction quality and substantially faster training.
Why Geometry Matters
A structured volume stores values on a regular grid, so its domain is implicit in the grid bounds. An unstructured volume can occupy an arbitrary region and must normally retain mesh geometry to determine whether a query lies inside the domain. Existing unstructured INRs compress the scalar field but still depend on an exterior surface for arbitrary queries and visualization.
For the datasets in this study, that surface alone limits effective compression to roughly 7-24x, regardless of how small the learned scalar model becomes. W-VEG+S instead uses the support of its learned Gaussians to encode the valid domain, allowing the model to answer both questions at once: is this position defined, and what is its scalar value?

| Representation | Scalar field | Domain geometry | Exterior surface needed |
|---|---|---|---|
| Unstructured INR | Implicit neural function | No | Yes |
| W-VEG | Explicit 3D Gaussians | No | Yes |
| W-VEG+S | Explicit 3D Gaussians | Yes | No |
Gaussian Volume Representation
Each world-space volume encoding Gaussian stores a center $\mu_i$, scale $S_i$, rotation $R_i$, scalar value $v_i$, and learnable weight $w_i$. For the Gaussians that influence a query point $p$, W-VEG estimates the scalar field with normalized Gaussian kernel regression:
\[\Phi(p) = \frac{\sum_{i=0}^{N} w_i v_i G_i(p)} {\sum_{i=0}^{N} w_i G_i(p)}, \qquad G_i(p) = \exp\left( -\frac{1}{2}(p-\mu_i)^\top R_i^\top S_i^{-2}R_i(p-\mu_i) \right).\]This is a world-space reconstruction rule, not an image-space splatting rule. It makes the compressed model directly sampleable at arbitrary 3D positions and therefore usable by downstream volume algorithms without first reconstructing a dense grid.
Training and Sampling
Seed isotropic Gaussians from a subset of volume vertices and their scalar values.
Evaluate a regular-grid CUDA rasterizer or an irregular-position BVH sampler.
Backpropagate reconstruction, continuity, and optional domain-boundary losses.
Prune weak primitives and add Gaussians at high-error samples until the size target is met.
Structured volumes
For regular grids, W-VEG adapts tile-based Gaussian rasterization to 3D voxel blocks. It computes each Gaussian’s block-aligned bounding box, sorts Gaussian-block intersections by block, and cooperatively loads intersecting Gaussians into shared memory. The backward pass is organized per Gaussian-block intersection, reducing atomic gradient updates by approximately the volume block size.

Unstructured volumes
Irregular sample positions cannot use the regular-grid pipeline. W-VEG instead builds a CUDA-accelerated bounding volume hierarchy over the requested sample positions. One warp cooperatively traverses the BVH for each Gaussian, computes contributions only for intersected leaves, and atomically accumulates values before normalization. The same traversal structure is reused in the backward pass.

Error-guided densification
The model has a fixed compression target but a dynamic spatial allocation. Every 100 iterations after iteration 500, the method can prune low-weight Gaussians and use the freed capacity in poorly reconstructed regions. The proposed strategy creates new Gaussians directly at the $k$ samples with highest reconstruction error, which is especially effective when spatial complexity is highly nonuniform.

Learning the unstructured domain
An explicit Gaussian set can fail in two ways. A false negative occurs when a valid volume position has insufficient Gaussian support; a false positive occurs when the model defines a value outside the original domain. Randomized within-cell sampling and a false-negative loss preserve coverage. W-VEG+S adds samples just outside the surface and a false-positive loss that suppresses support beyond the domain boundary.

Evaluation
The evaluation measures volumetric PSNR and end-to-end training time. Structured experiments compare against InstantVNR (IVNR) and AMGSRN. Unstructured experiments compare against UGINR and separately report real compression ratio, which includes the exterior surface whenever a model still requires one.
Structured-volume results
Across all four datasets and three compression ratios, W-VEG matches the reconstruction quality of substantially longer-trained neural baselines. Its average PSNR is 38.08 dB, compared with 37.76 dB for IVNR 200k and 38.04 dB for AMGSRN.
| Method | Average PSNR | Training schedule | W-VEG speedup |
|---|---|---|---|
| W-VEG | 38.08 dB | Convergence in 1k-4k steps | Reference |
| IVNR 20k | 36.09 dB | 20k steps | 1.52x |
| IVNR 200k | 37.76 dB | 200k steps | 14.91x |
| AMGSRN | 38.04 dB | 30k steps | 17.48x |

Unstructured-volume results
When only the scalar field is encoded, W-VEG achieves the best PSNR and training time on every tested dataset. On the three smaller volumes, it improves PSNR over UGINR by 33.81 dB on average. On Valley and Earthquake, it improves PSNR by 3.13 dB while training 34.9x faster on average. UGINR cannot train on the 17.9 GB Impact dataset because of memory limits.
The table below isolates the most aggressive 1024x target. Surface-dependent methods retain only 7.6-23.2x real compression on the datasets where they run. W-VEG+S removes that storage and preserves the full 1024x ratio.
| Dataset | UGINR PSNR | UGINR real CR | W-VEG PSNR | W-VEG real CR | W-VEG+S PSNR | W-VEG+S real CR |
|---|---|---|---|---|---|---|
| RBL | 20.72 | 9.8x | 75.63 | 9.8x | 58.57 | 1024x |
| Mito | 23.15 | 22.0x | 58.14 | 22.0x | 52.58 | 1024x |
| SF1 | 27.05 | 23.2x | 47.49 | 23.2x | 45.45 | 1024x |
| Valley | 34.32 | 18.9x | 35.83 | 18.9x | 34.75 | 1024x |
| Earthquake | 45.97 | 7.6x | 47.13 | 7.6x | 44.11 | 1024x |
| Impact | OOM | OOM | 44.82 | 18.3x | 44.77 | 1024x |
PSNR values are volumetric measurements. OOM indicates that UGINR could not train within available GPU memory.

Densification ablation
At 256x compression and a fixed 4,000-step budget, sample-error-based densification gives the highest PSNR in every structured, unstructured, and domain-encoding configuration tested. Its advantage is largest on the irregular Earthquake datasets.
| Model | 3DGS heuristic | MCMC | Error-based |
|---|---|---|---|
| Miranda | 35.88 | 33.97 | 36.58 |
| Chameleon | 45.74 | 41.53 | 47.86 |
| Earthquake | 43.79 | 46.97 | 52.25 |
| Earthquake+S | 42.23 | 38.84 | 50.85 |
| Impact | 46.53 | 44.61 | 48.74 |
| Impact+S | 46.89 | 45.27 | 48.16 |
Limitations and Future Work
- Reported sizes are in-memory models. The current method does not yet apply quantization, vector coding, or other 3D Gaussian file-compression techniques.
- W-VEG+S learns unstructured boundaries lossily. Domain accuracy depends on geometric complexity and trades off against scalar reconstruction quality at very small model sizes.
- The model supports arbitrary sampling, but task-specific accuracy for operations beyond volume rendering remains to be studied.
- Time-varying and multivariate volumes, signed-distance-based boundary losses, and a renderer designed directly around W-VEG are promising extensions.
Compression-ratio interpretation. Model CR compares only learned model size with the original volume. Real CR also counts any exterior mesh that must remain available. W-VEG+S is the only evaluated method whose model CR and real CR are identical for unstructured data.
BibTeX
@misc{dyken2026efficient,
title = {Efficient Compression of Structured and Unstructured Volumes via Learned 3D Gaussian Representation},
author = {Dyken, Landon and Chakrabarti, Sharmistha and Debardeleben, Nathan and Petruzza, Steve and Wu, Qi and Usher, Will and Kumar, Sidharth},
year = {2026},
eprint = {2607.01164},
archivePrefix = {arXiv},
primaryClass = {cs.LG},
doi = {10.48550/arXiv.2607.01164},
url = {https://arxiv.org/abs/2607.01164}
}