Rasterizer for Guassian Splatting using Taichi and PyTorch - embedded in python library. Currently very usable but in active development, so likely will break with new versions!
This work is originally derived off Taichi 3D Gaussian Splatting, with significant re-organisation and changes.
Key differences are the rendering algorithm is decomposed into separate operations (projection, shading functions, tile mapping and rasterization) which can be combined in different ways in order to facilitate a more flexible use, and gradients can be enabled on "all the things" as required for the application (and not when disabled, to save performance).
Using the Taichi autodiff for a simpler implementation where possible (e.g. for projection, but not for the rasterization).
Examples:
- Projecting features for lifting 2D to 3D
- Colours via. spherical harmonics
- Depth covariance without needing to build it into the renderer and remaining differentiable.
- Fully differentiable camera parameters (and ability to swap in new camera models)
A document describing some performance benchmarks of taichi-splatting here. Through various optimizations, in particular optimizing the summation of gradients in the backward gradient kernel. Taichi-splatting achieves a very large speedup (often an order of magnitude) over the original taichi_3d_gaussian_splatting, and is faster than the reference diff_guassian_rasterization for a complete optimization pass (forward+backward), in particular much faster at higher resolutions.
- taichi-nightly
Some bug fixes which have occurred after the 1.7.0 release:
pip install --upgrade -i https://pypi.taichi.graphics/simple/ taichi-nightly
- torch >= 1.8 (probably works with earlier versions, too)
- Install taichi-nightly as above (note that the external repo cannot be listed in the pyproject.toml)
- Clone down with
git clone
and install withpip install ./taichi-spatting
pip install taichi-splatting
There exists a toy optimizer for fitting a set of randomly initialized gaussians to some 2D images fit_image_gaussians
- useful for testing rasterization without the rest of the dependencies.
Fitting an image (fixed points):
fit_image_gaussians <image file> --show --n 20000
Fitting an image (split and prune to target):
fit_image_gaussians <image file> --show --n 1000 --target 20000
See --help
for other options.
There exist benchmarks to evaluate performance on individual components in isolation under taichi_splatting/benchmarks/
Tests (gradient tests and tests comparing to torch-based reference implementations) can be run with pytest, or individually under
taichi_splatting/tests/
A viewer for reconstructions created with the original gaussian-splatting repository can be found here or installed with pip. Has dependencies on open3d and Qt.
A benchmark for a full rendererer (in the same repository as above) with real reconstructions (rendering the original camera viewpoints). Options exist for tweaking all the renderer parameters, benchmarking backward pass etc.
-
Benchmarks with original + taichi_3dgs rasterizer
-
Simple view culling
-
Projection with autograd
-
Tile mapping (optimized and improved culling)
-
Rasterizer forward pass and optimized backward pass
-
Spherical harmonics with autograd
-
Gradient tests for most parts (float64) - including rasterizer!
-
Fit to image training example/test
-
Depth and depth-covariance rendering
-
Compute point visibility in backward pass (useful for model pruning)
-
Example training on images with split/prune operations
-
Novel heuristics for split and prune operations computed optionally in backward pass
-
Depth covariance example
-
3D training code (likely different repository)
-
Backward projection autograd takes a while to compile and is not cached properly
- Exposed all internal constants as parameters
- Switched to matrices as inputs instead of quaternions
- Tile mapping tighter culling for tile overlaps (~30% less rendered splats!)
- All configuration parameters exposed (e.g. tile_size, saturation threshold etc.)
- Warp reduction based backward pass for rasterizer, a decent boost in performance
Transformations are notated T_x_y
, for example T_camera_world
can be used to transform points in the world to points in the local camera by points_camera = T_camera_world @ points_world