Reconstruction and representation of 3D objects with radial basis functions
J. C. CarrR. K. BeatsonJ. B. CherrieToby J. MitchellW. R. FrightB. C. McCallumTalfan Evans
Develops a scalable implicit modeling framework using polyharmonic radial basis functions and fast evaluation algorithms to reconstruct, repair, and compress 3D surfaces from noisy or incomplete point clouds containing millions of points.
Industrial design, computer graphics, and digital manufacturing increasingly rely on 3D laser scanners to capture real-world objects. However, raw scan data typically contains millions of unorganized, noisy points with significant holes caused by physical occlusions. Traditional mathematical surface-fitting techniques, such as Radial Basis Functions (RBFs), provide smooth and watertight models, but they historically suffered from severe computational bottlenecks. Standard direct fitting required immense memory and processing power, restricting RBF models to small datasets of only a few thousand points.
The article demonstrates a framework that applies fast evaluation algorithms and data-reduction techniques to fit polyharmonic RBFs to 3D point-cloud datasets containing hundreds of thousands to millions of points. The objective is to evaluate whether these fast methods can reconstruct complex, smooth, watertight geometries, repair incomplete meshes, and filter out noise within reasonable memory and time limits on standard desktop hardware.
The approach models surfaces implicitly as zero-distance fields defined by polyharmonic RBFs, primarily biharmonic splines. Points lying on the surface are paired with projected off-surface normal points to construct a continuous signed-distance function. To overcome computational barriers, the authors adapted the Fast Multipole Method (FMM) to evaluate approximations in hierarchical clusters. They also integrated an iterative greedy algorithm for center reduction, which selects only the most significant data points needed to satisfy a user-specified accuracy threshold. Additionally, an adjustable smoothing parameter was introduced to handle noise, and a modified marching tetrahedra algorithm extracted optimized surface meshes.
The findings show that fast methods reduce memory demands from quadratic to linear scaling and cut computational complexity from cubic to near-linear scaling, allowing datasets of over 500,000 points to fit on a modest personal computer with 512 megabytes of memory. The greedy algorithm achieved dramatic compression, representing complex objects with only 10% to 25% of the original data points without sacrificing detail. In practical tests, the framework smoothly filled large holes in partial scans, correctly extrapolated missing anatomical and structural regions, and successfully smoothed noisy airborne laser (LIDAR) data while preserving separate physical boundaries. It also modeled highly intricate geometries, such as an industrial turbine blade with complex internal channels requiring nearly 600,000 centers.
These results establish that implicit functional modeling is a viable, high-performance alternative to traditional piecewise surface fitting. The continuous mathematical representation inherently defines a solid model, enabling direct calculation of smooth surface normals, seamless integration into solid geometry operations, and flexible mesh extraction at arbitrary resolutions. Consequently, organizations working with 3D scanning can significantly lower data storage costs, streamline mesh-repair workflows, and eliminate tedious manual patching of incomplete scans.
Organizations handling high-resolution 3D scan data should consider implementing fast implicit RBF methods to automate mesh repair, simplify dense meshes, and compress volumetric assets. When processing noisy point clouds, teams should tailor the smoothing stiffness parameter to balance surface fidelity against noise reduction. Looking forward, further research is recommended to parallelize fitting algorithms, accelerate processing for structured scanning grids, and decompose global functional models into localized implicit patches for faster ray-tracing and regional editing.
The primary limitation of this framework is its global mathematical structure, which can make local shape edits computationally expensive compared to traditional localized mesh formats. Fitting times can also range from minutes to several hours depending on the chosen accuracy and center reduction targets. Nevertheless, the underlying algorithms are mathematically well-founded, and the reported performance across diverse physical datasets demonstrates high reliability for engineering, graphics, and medical imaging applications.
- Paper: Surface reconstruction from unorganized points, Hugues Hoppe et al. (1992). Reading this foundational surface reconstruction work provides essential background on extracting smooth manifolds from unorganized point clouds before studying implicit RBF-based alternatives.
- Paper: Fourier Features Let Networks Learn High Frequency Functions in Low Dimensional Domains, Matthew Tancik et al. (2020). This book chapter extends the coordinate-based concepts of radial basis functions by exploring how Fourier feature mappings overcome spectral bias in continuous neural representations.
