Kathryn P. Drake,Edward J. Fuselier,Grady B. Wright
Abstract:Surface reconstruction from a set of scattered points, or a point cloud, has many applications ranging from computer graphics to remote sensing. We present a new method for this task that produces an implicit surface (zero-level set) approximation for an oriented point cloud using only information about (approximate) normals to the surface. The technique exploits the fundamental result from vector calculus that the normals to an implicit surface are curl-free. By using a curl-free radial basis function (RBF) interpolation of the normals, we can extract a potential for the vector field whose zero-level surface approximates the point cloud. We use curl-free RBFs based on polyharmonic splines for this task, since they are free of any shape or support parameters. Furthermore, to make this technique efficient and able to better represent local sharp features, we combine it with a partition of unity (PU) method. The result is the curl-free partition of unity (CFPU) method. We show how CFPU can be adapted to enforce exact interpolation of a point cloud and can be regularized to handle noise in both the normal vectors and the point positions. Numerical results are presented that demonstrate how the method converges for a known surface as the sampling density increases, how regularization handles noisy data, and how the method performs on various problems found in the literature.
What problem does this paper attempt to address?
This paper attempts to solve the problem of reconstructing implicit surfaces from a set of unordered point - cloud data. Specifically, the authors propose a new method - the Curl - Free Radial Basis Function Partition of Unity (CFPU) method, for reconstructing the zero - level set (i.e., the implicit surface) from point clouds with orientation information. This method utilizes an important conclusion in vector calculus: the normal of an implicit surface is curl - free, which means that the normal can be represented as the gradient of a certain scalar potential function.
### Main Problem Description
In many application fields, such as computer graphics, remote sensing, medical imaging, etc., reconstructing surfaces from a set of scattered point - cloud data is a common requirement. Traditional methods usually need to handle a large number of unordered data points and their corresponding normal vectors, and require these normal vectors to be "oriented", that is, each point is attached with an approximate normal vector. The CFPU method proposed in this paper aims to simplify this process through the following steps:
1. **Utilize the Curl - Free Property**: Since the normal of an implicit surface is curl - free, these normals can be approximated by the Curl - Free Radial Basis Function (RBF) interpolation method.
2. **Extract the Potential Function**: The normal field obtained by Curl - Free RBF interpolation can further extract a scalar potential function, and its zero - level set can approximate the surface represented by the original point cloud.
3. **Partition of Unity Method**: In order to improve computational efficiency and better capture local features, combined with the Partition of Unity (PU) method, the global problem is decomposed into multiple local sub - problems for solution.
### Method Characteristics
- **No Need to Extend Normals or Calculate Divergence**: Different from some methods based on Poisson equations, CFPU does not need to extend normals or calculate divergence, but directly constructs local potential functions from the normal data on the surface.
- **High Parallelism**: Since the Partition of Unity method is adopted, the potential function can be independently solved in each local area, which is very suitable for parallel computing.
- **No Need for Shape Parameters**: Use Curl - Free Polyharmonic Splines (PHS) and Curl - Free vector polynomials to construct local approximations, so there is no need to select shape parameters or support radii, thereby improving numerical stability.
### Applications and Advantages
The CFPU method can not only efficiently process large - scale point - cloud data, but also performs excellently when dealing with noisy data. In addition, this method can adapt to different application scenarios, such as accurate interpolation and regularization for noisy data. Experimental results show that as the sampling density increases, this method can converge well to the known surface and has good robustness to noise.
In summary, the main contribution of this paper is to propose a novel and efficient surface reconstruction method, which is especially suitable for reconstructing implicit surfaces from point - cloud data with orientation information.