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
Ground-truth and W-VEG reconstructions of the Impact volume at 1024x compression
Volume rendering of the 17.9 GB unstructured Impact dataset and its 17.5 MB W-VEG+S reconstruction. At 1024x compression, the model retains 38.43 dB PSNR and 0.984 SSIM while storing both the scalar field and its domain.

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.

38.08 dB average PSNR across the structured-volume benchmarks
14.91x average training speedup over IVNR trained for 200k steps
17.48x average training speedup over AMGSRN
1024x real unstructured compression without exterior surface storage
W-VEG reconstruction of the unstructured Impact dataset
The Impact reconstruction sampled from a 17.5 MB W-VEG+S model. The original unstructured volume occupies 17.9 GB.

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?

Comparison of sampling on regular structured grids and irregular unstructured volume domains
Structured volumes have an implicit regular domain. Unstructured volumes require an explicit geometry and may be undefined at arbitrary query positions.
RepresentationScalar fieldDomain geometryExterior surface needed
Unstructured INRImplicit neural functionNoYes
W-VEGExplicit 3D GaussiansNoYes
W-VEG+SExplicit 3D GaussiansYesNo

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.

\(\mu_i\)World-space center
\(S_i, R_i\)Anisotropic extent and orientation
\(v_i\)Encoded scalar value
\(w_i\)Learned influence on nearby queries

Training and Sampling

01 Initialize

Seed isotropic Gaussians from a subset of volume vertices and their scalar values.

02 Sample

Evaluate a regular-grid CUDA rasterizer or an irregular-position BVH sampler.

03 Optimize

Backpropagate reconstruction, continuity, and optional domain-boundary losses.

04 Adapt

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.

Forward and backward CUDA sampling pipeline for structured volumes
Structured sampling. The forward pass bins Gaussians into grid blocks and estimates cell values; the backward pass reuses the intersection list to accumulate gradients with less atomic contention.

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.

BVH-based sampling pipeline for arbitrary positions in unstructured volumes
Unstructured sampling. A BVH over arbitrary query positions is traversed with each Gaussian's bounding box, avoiding any regular-grid assumption.

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.

Comparison of 3DGS, MCMC, and sample-error-based Gaussian densification
Three density-control strategies evaluated in the paper. Error-based densification places new primitives where the sampled scalar field is reconstructed least accurately.

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.

False-negative and false-positive domain errors in a Gaussian volume model
False negatives create holes inside the valid domain; false positives extend Gaussian support beyond it. W-VEG+S explicitly penalizes both cases.

Evaluation

4 structured datasets537 MB to 32.21 GB
6 unstructured datasets174 MB to 17.93 GB
3 compression targets64x, 256x, and 1024x
NVIDIA A100In-memory model sizes; no file compression

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.

MethodAverage PSNRTraining scheduleW-VEG speedup
W-VEG38.08 dBConvergence in 1k-4k stepsReference
IVNR 20k36.09 dB20k steps1.52x
IVNR 200k37.76 dB200k steps14.91x
AMGSRN38.04 dB30k steps17.48x
Ground truth, W-VEG reconstruction, error, and Gaussian density for Miranda and Richtmyer volumes
Adaptive memory allocation at 1024x compression. Gaussian density remains relatively uniform for Miranda, but concentrates around Richtmyer's mixing surface where the scalar field is most complex.

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.

DatasetUGINR PSNRUGINR real CRW-VEG PSNRW-VEG real CRW-VEG+S PSNRW-VEG+S real CR
RBL20.729.8x75.639.8x58.571024x
Mito23.1522.0x58.1422.0x52.581024x
SF127.0523.2x47.4923.2x45.451024x
Valley34.3218.9x35.8318.9x34.751024x
Earthquake45.977.6x47.137.6x44.111024x
ImpactOOMOOM44.8218.3x44.771024x

PSNR values are volumetric measurements. OOM indicates that UGINR could not train within available GPU memory.

Rendered comparison of UGINR, W-VEG, and W-VEG+S on Earthquake and Mito
Neural volume renderings from 1024x models. Without exterior-surface intersection tests, W-VEG+S alone preserves the domain well enough to produce meaningful Earthquake and Mito renderings.

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.

Model3DGS heuristicMCMCError-based
Miranda35.8833.9736.58
Chameleon45.7441.5347.86
Earthquake43.7946.9752.25
Earthquake+S42.2338.8450.85
Impact46.5344.6148.74
Impact+S46.8945.2748.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}
}