Biomedical Signal Processing and Modeling Complexity of Living Systems 2014
View this Special IssueResearch Article  Open Access
Controllable Edge Feature Sharpening for Dental Applications
Abstract
This paper presents a new approach to sharpen blurred edge features in scanned tooth preparation surfaces generated by structuredlight scanners. It aims to efficiently enhance the edge features so that the embedded feature lines can be easily identified in dental CAD systems, and to avoid unnatural oversharpening geometry. We first separate the feature regions using graphcut segmentation, which does not require a userdefined threshold. Then, we filter the face normal vectors to propagate the geometry from the smooth region to the feature region. In order to control the degree of the sharpness, we propose a feature distance measure which is based on normal tensor voting. Finally, the vertex positions are updated according to the modified face normal vectors. We have applied the approach to scanned tooth preparation models. The results show that the blurred edge features are enhanced without unnatural oversharpening geometry.
1. Introduction
Optical scanning and geometric processing are two critical techniques in dental CAD systems which are responsible for acquiring tooth shapes and designing dental appliances, respectively. Various studies have been published on building dedicated scanning systems [1, 2] and automating the procedure of generating the shapes of dental appliances [3–5]. However, there are still limitations, of which the feature blurring is a prominent one. The feature blurring problem has a significant impact on the cervical line extraction which is a necessary step in modeling various dental restorations. As shown in Figure 1(a), the original scanned tooth preparation model contains blurred feature regions, which makes the automated cervical line extraction unreliable. The problem lies in the limitations of the structuredlight principle. For example, algorithms based on phase analysis [6] confine the data density according to the resolution of the projected fringes. It is difficult to be solved by improving the structuredlight algorithms. As a result, geometric postprocessing is essential to further improve the quality of the scanned surfaces. As shown in Figure 1(c), by sharpening the blurred feature regions, highquality cervical lines are obtained robustly.
(a)
(b)
(c)
(d)
Geometric filtering is a versatile tool to alter the properties of scanned surfaces represented by triangle meshes. It can make scanned surfaces more appropriate for specific visualization and shapebased product design tasks. For example, surface noise [7–10], the most common defect, can be reduced by geometric filtering, and geometric filteringbased feature enhancement can be used to exaggerate the microstructure on the artifacts surface in archeology. In order to emphasize the interesting surface attributes, a variety of filtering approaches have been developed to modify derived differential quantities instead of vertex positions. For example, Laplacian coordinate has been employed for mesh denoising and enhancing [11, 12]; curvature has been prescribed to directly control the shape of the surface in [13]. In comparison with algorithms involving secondorder differential attributes, normal based filtering algorithms [14–16] are more appropriate to process anisotropic features. The reason is that the secondorder differential attributes integrate characteristics in all directions so that they are not flexible to constrain anisotropic features in some directions. Although existing geometric filtering algorithms alleviate the feature blurring problem to some extent, none of them considers the degree of the sharpness. The processed edge features usually show unnatural oversharpening geometry.
In this paper, we focus on the problem of enhancing blurred edge features in a controllable manner. Specifically, the degree of the sharpness or the fillet radius is controlled to avoid oversharpening geometry. We propose a feature distance measure based on normal tensor voting to control the normal filtering process. After the filtering, the vertex positions are updated by fitting the new face normal vectors in the least square sense. In addition to geometric filtering, feature region detection is also important for solving the feature blurring problem since engineering users demand highfidelity scanned surfaces. As a result, the featureless regions should be untouched. We consider this problem as a segmentation to avoid involving a userdefined threshold which is common in most prior researches. We adopt a graphcut method to compute the segmentation. The main contributions of the paper contain three aspects as follows.(1)Unlike most existing mesh sharpening methods which produce oversharpening geometry benefitting highquality visualization, the proposed mesh sharpening method, which controls the sharpness or the fillet radius of edge features, is more appropriate for designing shapes of dental appliances. The essential strategy is also applicable to process scanned models used in mechanical and arts industry.(2)We propose a feature distance measure based on normal tensor voting to control the sharpness of edge features.(3)We cast the feature region detection into a segmentation problem and solve it with a graphcut algorithm.
The remainder of this paper is organized as follows. In Section 2, we review the most relevant previous works. Then an overview of our approach is presented in Section 3. The core algorithms of the feature region segmentation and the controllable mesh sharpening are detailed in Sections 4 and 5, respectively. After discussing the results and the applications of our approach in Section 6, we conclude the paper in Section 7.
2. Related Works
2.1. Mesh Detail Editing
Several mesh denoising algorithms adapt twodimensional signal processing theory to filter vertex positions. Taubin [7] proposed the first lowpass filtering algorithm for mesh smoothing. Desbrun et al. [8] improved the efficiency of the filter through an implicit solver. In order to achieve feature preserving, a variety of methods employed bilateral filters [9, 10] and anisotropic diffusion [20, 21] to reduce noises in flat regions while they maintain discontinuities in high contrast regions. In contrast to directly dealing with vertex positions, several researchers [11, 14–16] found that filtering highorder differential quantities brings obvious advantages in terms of flexibility and effectiveness. Shen and Barner [14] applied the fuzzy filter on normal vectors, and Yagou et al. [16] applied the boost filter on normal vectors. Since the edge features are naturally represented as discontinuity or large variance of normal vectors, the normal vectors are appropriate for modeling sharp edge features. Su et al. [11] first filtered the Laplacian coordinates and then reconstructed vertex positions. With similar ideas, Wang et al. [12] detailed versatile effects based on filtering Laplacian coordinates. Recently, algorithms which involve explicit feature detection [22, 23] and classify vertices into feature and featureless regions have been proposed based on the idea that multiple segments with different attributes should not be blended. Different vertex groups in neighborhood structure are filtered separately.
Edge and corner features are important for CAD and sculpture models used in mechanical and arts industry. Unfortunately, the edge and corner features are commonly degenerated depending on how the models are obtained. As a result, mesh sharpening is required to reconstruct the sharp edge and corner features which do not exist in original mesh surfaces. Attene et al. [24] proposed a twostep method to repair sharp edge features for mesh surfaces extracted from volume data. Wang [17] employed an incremental filter to extend the geometry of smooth region into the feature region. Wang [25] took advantage of the bilateral filter [10] to detect and recover sharp features. Chen and Cheng [26] used a sharpness dependent filter to recover sharp structure in surface holefilling. Chen and Cheng [26] presented a normal filteringbased algorithm to form sharp edge features. Actually, the key idea of prior algorithms is based on the assumption that the sharp features are intersections between smooth regions. Different strategies are taken to extend smooth regions to form sharp features. However, these methods inevitably produce oversharpening geometry which is undesirable for scanned mesh surfaces.
In addition to the above local methods, global optimization methods are also developed, which can take advantage of the integral property of mesh models. For example, Ji et al. [19] proposed a global optimization procedure to enhance mesh surfaces. He and Schaefer [27] proposed optimization to improve the mesh quality. Although global methods provide high quality results, they require high computation time and memory footprint in general. Moreover, the local characteristics can hardly be controlled by the global methods.
2.2. Feature Detection
Sharp features especially edge features play an important role in structureaware shape processing tasks. For example, in reverse engineering, mesh surfaces are separated along feature lines and fitted into surface patches. Most existing approaches focus on extracting feature lines. Rössl et al. [28] extracted feature lines using morphological operators. Yoshizawa et al. [29] detected feature lines based on the differential definition of the valleys and ridges and located the feature lines by using local surface fitting. All of the above methods are based on curvature information. In contrast, Kim et al. [30] took advantage of normal tensor voting to classify the features into different categories and grouped feature regions through means clustering in the feature space. Wang et al. [31] extended the normal tensor voting method to extract feature lines by proposing a neighbor support saliency. In this paper, feature regions are detected to reduce the amount of calculation.
3. Overview
The target models of our mesh sharpening algorithm are scanned surfaces produced by optical scanning systems. They commonly have a great number of triangles, which makes global approaches such as [19] unqualified. Moreover, the scanned surfaces produced by structuredlight scanners can achieve accuracy of about 60 m, which makes mesh denoising unnecessary. With these considerations, the method in this paper consists of three main stages: detect feature regions, filter the normal vectors of triangle faces, and update vertex positions according to the filtered normal vectors. Although the method in [18] has taken similar steps to sharpen mesh surfaces, the improvements of our method include two aspects: we avoid userdefined thresholds through graphcut segmentation and in order to avoid oversharpening geometry, a feature distance measure quantifying the distance away from the smooth region, as illustrated in Figure 2(b), is proposed based on normal tensor analysis.
(a)
(b)
(c)
(d)
(e)
(f)
(g)
(h)
Prior feature region detection algorithms for mesh sharpening [17, 18] commonly analyze normal variance in the local neighborhood of a central face and specify a threshold to identify feature regions. This strategy does not consider the spatial coherence of the detected feature regions. In contrast, we adopt a graphcut algorithm which involves spatial constraint as shown in Figure 2(c).
The key ideas of most effective mesh sharpening algorithms [17, 18] are similar which propagate the geometry from smooth regions to feature regions to form edge intersections. Plane fitting and skeletonisation are used in [17]; normal filtering and greedy propagation are adopted in [18]. However, these approaches inevitably produce oversharpening edge features which are unnatural for scanned surfaces. Our algorithm involves a feature distance measure to control the degree of sharpness. Figures 2(d) and 2(e) show the normal color map before and after the filtering process. Figure 2(f) shows the final result in which the edge features are enhanced but do not suffer from oversharpening defects.
4. Feature Region Detection Using Graph Cuts
A given scanned surface can be represented by a triangle mesh , where and are the vertices and triangle faces, respectively. Here denotes the cardinality of a set. Each face has a normal vector which is denoted by .
4.1. Feature Distance Metric
The normal tensor describes the local structure of a vertex of . As suggested by Kim et al. [30], the normal tensor classifies local geometries into three types of features, namely, smooth surface, edge feature, and corner feature. The normal tensor at is defined as where is the face set in the onering neighborhood of and is the weight for the covariance matrix generated by face normal vector . The difference of the definitions of normal tensor is mainly about the definition of which is defined here as where is the area of and is the maximum triangle area among ; is the barycenter of and is the edge length of the bounding box including . The eigendecomposition of uncovers the local structure of : where are three eigenvalues of and , , and are three corresponding eigenvectors. As shown in Figure 3, the relative values of , , and determine the feature type in the neighborhood of .
(a)
(b)
(c)
Based on above normal tensor framework, we define a feature distance measure in feature space constructed by eigenvectors of . First, we find the feature points corresponding to smooth regions of through means clustering in the feature space. As shown in Figure 4(a), after means clustering, the feature points are separated into different compact groups. The final result does not heavily depend on the parameter which is chosen as 3 in this paper. The group with highest component along the largest eigenvector is denoted by the smooth set. Other feature points outside the smooth set form the feature set. In order to quantify how feature points are far away from the smooth set, the feature distance measure is defined as the Mahalanobis distance from the smooth set: where is the coordinate of a testing feature point, is the covariance matrix of feature points in smooth set, and is the mean of feature points in the smooth set. As shown in Figure 4(b) the proposed feature distance faithfully captures the anisotropic feature regions.
(a)
(b)
The feature distance measure has two functions in our algorithm: one is to provide a distribution model in the feature detection step; the other is to control the normal vector filtering process.
4.2. Feature Region Segmentation
Feature region detection is commonly solved by thresholding some attributes of mesh surfaces. For example, approaches in [17, 18] employ the normal variance in the local neighborhood as the attribute. However, this scheme involves multiple userdefined parameters such as the size of the local neighbor, tolerant normal variance, etc. In order to avoid these parameters, we adopt a graphcut algorithm to separate the feature regions from smooth regions.
Let be the dual graph of where is the nodes of the dual graph, is the edge set of the dual graph, each edge connects two neighboring faces, and is the weights defined on edges. To perform a graphcut segmentation, we add two virtual nodes. One is the source node which represents the smooth regions; the other is the sink node which represents the feature regions. Then the energy function of the graphcut segmentation is defined as where is a labeling for triangles of , is the regional penalty for assigning labels, is the boundary penalty for assigning different labels between neighbor triangles, and is the relative importance of the two terms in (5) which is specified as 1.0. The behaviors of the segmentation depend on the definitions of and . To separate feature regions, we define them as follows: where is the maximal distance of feature points. We employ the algorithm in [32] to optimize the energy defined in (5). The computation is efficient and the spatial coherence is guaranteed as shown in Figure 5.
(a)
(b)
5. Normal Filtering in a Controllable Manner
In the previous section, we have confined the following normal filtering in the feature regions so that unnecessary computations are avoided. To reconstruct sharp edge features, the common strategy is to propagate geometry from smooth regions to feature regions; the difference from previous approaches is the way to predict vertex positions. However, these approaches all result in oversharpening geometry since the filtered geometry is the same with the smooth regions where the propagation begins. In contrast, we adopt the feature distance measure defined in (4) to control the normal filtering process: where is the triangles in the onering neighborhood of . The feature distance weights make the triangles at feature region tend to maintain its original normal vectors, which is defined as where , , and are positions of the three vertices of . The parameter controls the sharpness of the edge feature region. A larger value of corresponds to a high degree of the sharpness. The impact of different values of parameter is demonstrated in Figure 6. For processing tooth preparation models, the parameter is experimentally chosen as 0.5 in our tests.
(a)
(b)
(c)
In order to propagate the geometry of the smooth region to feature region, we adopt a greedy process to iteratively filter the face normal vectors using (7). The priority is determined by the feature distance measure. After the desired faces normal vectors have been obtained, we update the vertex positions through least square approximation to the filtered normal vectors. We adopt the energy function used in [33]: where is the vertex positions. We solve (9) using gradient descent method.
6. Results
We have developed an implementation of the proposed mesh sharpening algorithm using C++ language. We present several tests on tooth preparation and mechanical and arts models below. All tests are conducted on a PC with Intel Core i5 CPU, 2 GB main memory, and Windows XP operating system. We compare our method with the most similar approach [18] which also employs normal filtering. Firstly, we present the results on tooth preparation models. As shown in Figure 7, both our approach and the method in [18] successively enhance the blurred edge features. However, our method avoids the oversharpening geometry which makes the scanned surface unnatural. Specifically, the sharpened edge features generated by the method in [18] are singleedge wide which can be identified through dihedral angles. As for modeling dental restorations, the oversharpening geometry may destroy the original morphology of the cervical lines. We further compare the Hausdorff distance between the original scanned tooth preparation and its sharpened versions generated by the method in [18] and ours. As shown in Figure 8, our controllable sharpening algorithm can maintain the shape of the cervical line while enhancing the regions around it.
(a)
(b)
(c)
(a)
(b)
In addition to tooth preparation models, as shown in Figures 9, 10, and 11, our approach is also capable of processing scanned surfaces used in mechanical and arts industry.
(a)
(b)
(c)
(a)
(b)
(c)
(a)
(b)
(c)
Note that prior methods try to directly construct feature lines on mesh surface which can be easily identified through dihedral angles. However, this characteristic is only desired for computergenerated CAD models. For scanned surfaces, the mesh sharpening algorithm should avoid oversharpening geometry. In addition, the scanned models are usually quite large. Therefore, the computational cost is critical for practical applications. The timing statistics of the proposed approach is given in Table 1, from which we can conclude that the time cost is reasonable and approximately linear to the model size.

We further compare with the mesh enhancing method in [19], which optimizes all vertex positions through moving the vertices in flat regions to highcurvature regions. As shown in Figure 12, all the models have the same number of vertex samples. The result of the method in [19] modifies all the vertex positions, leading to dense sampling in highcurvature regions. In contrast, our result only filters the vertex samples around the edge features. In addition, the time cost of the method in [19] is 196 seconds with our implementation.
(a)
(b)
(c)
7. Conclusions
In this paper, we have proposed a novel mesh sharpening algorithm which enhances edge features of scanned surface models in a controllable manner. The main components of the proposed approach consist of two factors: detecting feature regions and propagating the geometry from the smooth regions to the feature regions. By introducing a feature distance measure based on normal tensor analysis, we obtain naturallyenhanced edge features on scanned surfaces like tooth preparation and mechanical and arts models.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgment
This work was supported by the Science and Technology Plan of Zhejiang Province (Grant no. 2011C13009).
References
 H. Cui, N. Dai, W. Liao, and X. Cheng, “Intraoral 3D optical measurement system for tooth restoration,” Optik—International Journal for Light and Electron Optics, vol. 124, no. 12, pp. 1142–1147, 2013. View at: Publisher Site  Google Scholar
 M. Chang and S. C. Park, “Automated scanning of dental impressions,” Computer Aided Design, vol. 41, no. 6, pp. 404–411, 2009. View at: Publisher Site  Google Scholar
 H. T. Yau, C. Y. Hsu, H. L. Peng, and C. C. Pai, “Computeraided framework design for digital dentistry,” ComputerAided Design and Applications, vol. 5, no. 5, pp. 667–675, 2008. View at: Publisher Site  Google Scholar
 T. Steinbrecher and M. Gerth, “Dental inlay and onlay construction by iterative laplacian surface editing,” in Proceedings of the Symposium on Geometry Processing (SGP '08), pp. 1441–1447, 2008. View at: Google Scholar
 N. Qiu, R. Fan, L. You, and X. Jin, “An efficient and collisionfree holefilling algorithm for orthodontics,” The Visual Computer, vol. 29, no. 6–8, pp. 577–586, 2013. View at: Publisher Site  Google Scholar
 S. S. Gorthi and P. Rastogi, “Fringe projection techniques: whither we are?” Optics and Lasers in Engineering, vol. 48, no. 2, pp. 133–140, 2010. View at: Publisher Site  Google Scholar
 G. Taubin, “Signal processing approach to fair surface design,” in Proceedings of the 22nd Annual ACM Conference on Computer Graphics and Interactive Techniques, pp. 351–358, August 1995. View at: Google Scholar
 M. Desbrun, M. Meyer, P. Schroder, and A. H. Barr, “Implicit fairing of irregular meshes using diffusion and curvature flow,” in Proceedings of the 26th Annual Conference on Computer Graphics and Interactive Techniques (SIGGRAPH '99), pp. 317–324, 1999. View at: Publisher Site  Google Scholar
 S. Fleishman, I. Drori, and D. CohenOr, “Bilateral mesh denoising,” ACM Transactions on Graphics, vol. 22, no. 3, pp. 950–953, 2003. View at: Publisher Site  Google Scholar
 T. R. Jones, F. Durand, and M. Desbrun, “Noniterative, feature preserving mesh smoothing,” ACM Transactions on Graphics, vol. 22, no. 3, pp. 943–949, 2003. View at: Publisher Site  Google Scholar
 Z. Su, H. Wang, and J. Cao, “Mesh denoising based on differential coordinates,” in Proceedings of the IEEE International Conference on Shape Modeling and Applications (SMI '09), pp. 1–6, Beijing, China, June 2009. View at: Publisher Site  Google Scholar
 H. Wang, H. Chen, Z. Su, J. Cao, F. Liu, and X. Shi, “Versatile surface detail editing via Laplacian coordinates,” The Visual Computer, vol. 27, no. 5, pp. 401–411, 2011. View at: Publisher Site  Google Scholar
 M. Eigensatz, R. W. Sumner, and M. Pauly, “Curvaturedomain shape processing,” Computer Graphics Forum, vol. 27, no. 2, pp. 241–250, 2008. View at: Publisher Site  Google Scholar
 Y. Shen and K. E. Barner, “Fuzzy vector medianbased surface smoothing,” IEEE Transactions on Visualization and Computer Graphics, vol. 10, no. 3, pp. 252–265, 2004. View at: Publisher Site  Google Scholar
 H. Yagou, Y. Ohtake, and A. G. Belyaev, “Mesh smoothing via mean and median filtering applied to face normals,” in Proceedings of the Geometric Modeling and Processing, pp. 124–131, 2002. View at: Google Scholar
 H. Yagou, A. Belyaevy, and D. Weiz, “Highboost mesh filtering for 3D shape enhancement,” Journal of Three Dimensional Images, vol. 17, pp. 170–175, 2003. View at: Google Scholar
 C. C. L. Wang, “Incremental reconstruction of sharp edges on mesh surfaces,” Computer Aided Design, vol. 38, no. 6, pp. 689–702, 2006. View at: Publisher Site  Google Scholar
 J. G. Shen and Z. Y. Chen, “Mesh sharpening via normal filtering,” Journal of Zhejiang University Science A, vol. 10, pp. 546–553, 2009. View at: Google Scholar
 Z. Ji, L. Liu, B. Wang, and W. P. Wang, “Feature enhancement by vertex flow for 3D shapes,” ComputerAided Design and Applications, vol. 8, no. 5, pp. 649–664, 2011. View at: Publisher Site  Google Scholar
 M. Desbrun, M. Meyer, P. Schroder, and A. Barr, “Anisotropic featurepreserving denoising of height fields and bivariate data,” in Proceedings of the Graphics Interface, pp. 145–152, 2000. View at: Google Scholar
 C. L. Bajaj and G. Xu, “Anisotropic diffusion of surfaces and functions on surfaces,” ACM Transactions on Graphics, vol. 22, no. 1, pp. 4–32, 2003. View at: Publisher Site  Google Scholar
 H. Fan, Y. Yu, and Q. Peng, “Robust featurepreserving mesh denoising based on consistent subneighborhoods,” IEEE Transactions on Visualization and Computer Graphics, vol. 16, no. 2, pp. 312–324, 2010. View at: Publisher Site  Google Scholar
 J. Wang, X. Zhang, and Z. Yu, “A cascaded approach for featurepreserving surface mesh denoising,” Computer Aided Design, vol. 44, no. 7, pp. 597–610, 2012. View at: Publisher Site  Google Scholar
 M. Attene, B. Falcidieno, M. Spagnuolo, and J. Rossignac, “Sharpen&Bend: recovering curved sharp edges in triangle meshes produced by featureinsensitive sampling,” IEEE Transactions on Visualization and Computer Graphics, vol. 11, no. 2, pp. 181–192, 2005. View at: Publisher Site  Google Scholar
 C. C. L. Wang, “Bilateral recovering of sharp edges on featureinsensitive sampled meshes,” IEEE Transactions on Visualization and Computer Graphics, vol. 12, no. 4, pp. 629–639, 2006. View at: Publisher Site  Google Scholar
 C. Y. Chen and K. Y. Cheng, “A sharpnessdependent filter for recovering sharp features in repaired 3D mesh models,” IEEE Transactions on Visualization and Computer Graphics, vol. 14, no. 1, pp. 200–212, 2008. View at: Publisher Site  Google Scholar
 L. He and S. Schaefer, “Mesh denoising via L_{0} minimization,” ACM Transactions on Graphics, vol. 32, no. 4, article 64, 2013. View at: Publisher Site  Google Scholar
 C. Rössl, L. Kobbelt, and H. P. Seidel, “Extraction of feature lines on triangulated surfaces using morphological operators,” in Proceedings of the AAAI Symposium on Smart Graphics, 2000. View at: Google Scholar
 S. Yoshizawa, A. Belyaev, and H. Seidel, “Fast and robust detection of crest lines on meshes,” in Proceedings of the ACM Symposium on Solid and Physical Modeling (SPM 2005), pp. 227–232, June 2005. View at: Publisher Site  Google Scholar
 H. S. Kim, H. K. Choi, and K. H. Lee, “Feature detection of triangular meshes based on tensor voting theory,” Computer Aided Design, vol. 41, no. 1, pp. 47–58, 2009. View at: Publisher Site  Google Scholar
 X. Wang, J. Cao, X. Liu, B. Li, X. Q. Shi, and Y. Sun, “Feature detection of triangular meshes via neighbor supporting,” Journal of Zhejiang University Science C, vol. 13, no. 6, pp. 440–451, 2012. View at: Publisher Site  Google Scholar
 Y. Boykov and V. Kolmogorov, “An experimental comparison of mincut/maxflow algorithms for energy minimization in vision,” IEEE Transactions on Pattern Analysis and Machine Intelligence, vol. 26, no. 9, pp. 1124–1137, 2004. View at: Publisher Site  Google Scholar
 X. Sun, P. L. Rosin, R. R. Martin, and F. C. Langbein, “Fast and effective featurepreserving mesh denoising,” IEEE Transactions on Visualization and Computer Graphics, vol. 13, no. 5, pp. 925–938, 2007. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2014 Ran Fan and Xiaogang Jin. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.