A signal processing approach to fair surface design
Gabriel Taubin
Proposes a linear-time signal processing framework that smooths large polyhedral surfaces of arbitrary topology via low-pass filtering while preventing shrinkage and enforcing geometric constraints.
Digital 3D models extracted from sources such as medical scans and range sensors often produce massive meshes containing hundreds of thousands to millions of vertices. Smoothing these large, arbitrary shapes without distorting their overall volume has historically presented a major computational bottleneck. Existing techniques rely on global optimization or solving large equation systems that scale quadratically with mesh size, making interactive editing unaffordable, while standard local smoothing methods severely shrink the geometry.
The article demonstrates a fast, signal-processing approach to smoothing polyhedral surfaces of arbitrary topology without causing shrinkage. It sets out to provide a scalable algorithm that integrates seamlessly into interactive geometric design workflows.
The proposed method adapts classical Fourier analysis to surfaces by treating vertex coordinates as discrete signals and framing smoothing as low-pass filtering. Rather than performing full frequency decompositions or solving massive linear systems, the approach uses a two-step iterative filter. It pairs a standard positive smoothing step with a compensating negative scaling step. The evaluation demonstrates this technique on synthetic geometries, medical models such as spine scans, and meshes integrated with subdivision algorithms.
Key findings show that the algorithm operates with linear time and memory complexity, requiring only basic sparse matrix multiplications rather than quadratic optimizations. Crucially, the alternating two-step filter removes high-frequency noise and unwanted undulations while preserving the model's original volume and major features. The algorithm also supports precise design constraints; users can lock specific vertices, control surface normals, and create smooth deformations by modifying local neighborhood definitions or solving only very small, localized linear systems.
These results demonstrate that smoothing very large meshes is practically feasible for interactive applications, drastically reducing processing costs and execution times. Interactive computer-aided design tools can adopt this filter to eliminate the trade-off between smoothing quality and geometric shrinkage, outperforming prior subdivision pipelines that introduced surface ripples or collapsed object boundaries.
Engineering and graphics teams should adopt this low-pass filtering formulation to replace Gaussian smoothing steps in subdivision frameworks and interactive mesh editors. Future development should explore expanding this signal-processing foundation to replace computationally heavy variational formulations in surface reconstruction and physics-based shape fitting.
Confidence in the approach is high for standard symmetric mesh structures, where stability and spectral bounds are mathematically guaranteed. However, users should exercise caution when introducing widespread non-symmetric neighborhood constraints, as these alter eigenvalue distributions and require validation to prevent unexpected deformations.
- Paper: The emerging field of signal processing on graphs: Extending high-dimensional data analysis to networks and other irregular domains, David I Shuman et al. (2012). Reviewing this signal processing framework on graphs provides essential background on generalizing harmonic analysis to irregular domains before reading about surface fairing.
- Paper: Surface simplification using quadric error metrics, Michael Garland et al. (1997). This paper builds directly upon surface processing concepts by introducing quadric error metrics for high-fidelity mesh simplification.
- Paper: Progressive meshes, Hugues Hoppe (1996). This work extends surface mesh manipulation techniques by introducing progressive representations and multi-resolution level-of-detail management.
