Next Article in Journal
Pipeline Leakage Detection Using Acoustic Emission and Machine Learning Algorithms
Next Article in Special Issue
A Targetless Method for Simultaneously Measuring Three-Degree-of-Freedom Angular Motion Errors with Digital Speckle Pattern Interferometry
Previous Article in Journal
Behavior and Task Classification Using Wearable Sensor Data: A Study across Different Ages
Previous Article in Special Issue
Advances of Research on Dual-Frequency Solid-State Lasers for Synthetic-Wave Absolute-Distance Interferometry
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

An Adaptive Hybrid Sampling Method for Free-Form Surfaces Based on Geodesic Distance

College of Control Science and Engineering, China University of Petroleum (East China), Qingdao 266580, China
*
Author to whom correspondence should be addressed.
Sensors 2023, 23(6), 3224; https://doi.org/10.3390/s23063224
Submission received: 3 February 2023 / Revised: 14 March 2023 / Accepted: 15 March 2023 / Published: 17 March 2023

Abstract

:
High precision geometric measurement of free-form surfaces has become the key to high-performance manufacturing in the manufacturing industry. By designing a reasonable sampling plan, the economic measurement of free-form surfaces can be realized. This paper proposes an adaptive hybrid sampling method for free-form surfaces based on geodesic distance. The free-form surfaces are divided into segments, and the sum of the geodesic distance of each surface segment is taken as the global fluctuation index of free-form surfaces. The number and location of the sampling points for each free-form surface segment are reasonably distributed. Compared with the common methods, this method can significantly reduce the reconstruction error under the same sampling points. This method overcomes the shortcomings of the current commonly used method of taking curvature as the local fluctuation index of free-form surfaces, and provides a new perspective for the adaptive sampling of free-form surfaces.

1. Introduction

Free-form surface parts are widely used in aerospace, automotive and other high-end equipment manufacturing fields [1]. It is one of the hotspots in the measurement field to measure their precise geometric quantities and ensure that they meet the accuracy requirements. This kind of part has a complex structure and large size, which also puts forward higher requirements for the dexterity and space of measuring equipment. The geometric measurement of a free-form surface can adopt a contact or non-contact measurement [2]. In the process of the geometric measurement of a free-form surface, the geometric error can be evaluated from two aspects: measurement accuracy and measurement efficiency. The measurement accuracy can be quantified by measuring the deviation between the actual value and the real value. When the measurement accuracy meets the measurement requirements, the measurement efficiency can be improved by optimizing the number and location of the measurement points and the measurement path [3,4,5,6].
As shown in Figure 1a, the coordinate measuring machine (CMM) equipped with a contact probe is one of the common contact measuring instruments [7]. It can obtain the point coordinates of the parts to be measured by scanning the measurement point by point. The point coordinates of the parts to be measured are reconstructed into free-form surfaces, and the geometric errors are compared with the original surfaces to accurately evaluate the machining errors generated in the manufacturing process of free-form surface parts. As shown in Figure 1b, the articulated industrial robot equipped with an optical probe is one of the common non-contact measuring instruments. The scanning measurement path is generated based on the measurement points of the parts to be measured with significant geometric characteristics, and the industrial robot drives the optical probe to move along the scanning measurement path. The scanning point cloud is reconstructed into a free surface, and the geometric error is compared with the original surface to determine the machining error of the free surface [8].
When measuring the geometric quantity of the parts to be measured, the number of measuring points is positively related to the measurement accuracy. The more the number of measuring points, the higher the measurement accuracy will be. When the number of measuring points reaches a certain value, the measurement accuracy will not change significantly. However, if the number of measuring points is too large, the measurement time will be significantly increased and the measurement efficiency will be reduced. Therefore, when measuring the components of the free-form surface to be measured, the number and position of the sampling points on the free-form surface should be reasonably distributed to ensure that the measurement efficiency can be improved as much as possible on the premise of meeting the measurement accuracy.
By setting a reasonable sampling plan, the corresponding measurement path is generated, and the free surface is measured using a contact measurement or non-contact measurement. Thus, the high-precision measurement of free-form surfaces can be realized, and the measurement efficiency can be significantly improved.
Blind sampling is a sampling method that does not consider the geometric characteristics of the free-form surface [9]. The number and location of sample points determine the geometric error of the reconstructed free-form surface. Uniform sampling is a simple and efficient sampling method, which is widely used in the field of free-form surface sampling. Dunbar et al. [10,11] introduced a random sampling strategy of fast disk under arbitrary dimension, which can realize blind sampling of a free-form surface simply and efficiently. Woo et al. [12] sampled the surface to be measured using the Hammersley distribution method. The research results show that blind sampling has the characteristics of simple and efficient sampling, but it cannot adaptively change the number and position of the sampling points based on the geometric characteristics of the free-form surface to be measured, thus affecting the measurement accuracy and efficiency.
Adaptive sampling generates sampling points according to the geometric characteristics of a free-form surface. In short, more sampling points are generated in the area with a large fluctuation of the free-form surface, and fewer sampling points are generated in the area with a small fluctuation of the free-form surface [13]. Ren et al. [14] proposed a method of using the curvature change matrix of adjacent points as the change index of a free-form surface. According to the proposed index, the optimal position of the newly added bi-directional curve mesh was determined, and the free-form surface was reconstructed based on the Gordon surface fitting principle. Javad et al. [15] used the optimization method and particle swarm optimization algorithm to optimize the position of the sampling points. On the basis of the initial sampling points, the optimal position of the new sampling points was determined by iterative optimization. When the geometric error between the reconstructed surface and the original surface or the number of sampling points reached the set threshold, the optimization was terminated. Gao et al. [16] took the aeroengine blade as the measurement object, and adaptively generated the sampling points based on the bending moment theory for the regions with different curvatures of the engine blade, so that the number and position of the sampling points could be accurately determined according to the fluctuation of the blade surface. In consideration of the influence of Gaussian curvature on machining error, Sang et al. [17] proposed a scanning line distribution strategy based on star pattern, which classifies peak points and anchor points, and connects the error peak points in different regions and anchor points in the same region, so that more scanning lines are generated in the region with small machining error, and fewer scanning lines are generated in the region with small machining error. Jiang et al. [18] studied a calculation method based on curve chord deviation, and adopted a two-step sampling method. First, based on the radius of the CMM probe ball as the threshold, adaptive Isoparametric sampling was carried out on the leading edge curve and trailing edge surface of the blade. Secondly, adaptive sampling points were selected on the Isoparametric based on the proposed curve chord deviation theory, and finally, the adaptive sampling of the blade was realized. Suleiman et al. [19] proposed a patch based free-form surface adaptive sampling method, which sorts patches according to the Gaussian curvature of each patch, determines the number of sampling points of each patch according to the sorting size, and selects the points of maximum curvature, minimum curvature and average curvature in the patch as the sampling points. If the number of sampling points exceeds these three types of points, the maximum curvature, minimum curvature and average curvature are taken as the sampling points, and so on until the number of sampling points reaches the threshold of the number of sampling points. Mansour [20] studied an adaptive sampling method for reducing the number of measurement points and improving the measurement efficiency with the blade as the measurement object. Based on the least square method, the minimum number of points required for the curve polynomial is found, so that the fitting error expressed by the polynomial curve and cubic curve is minimal. He et al. [21] proposed an adaptive sampling method for a free-form surface based on the machining error model. The machining error model was established based on the curvature of the free-form surface, and the relationship between Gaussian curvature and machining error was obtained. The adaptive sampling of a free-form surface was carried out based on the error model and the Hammersley principle. The number of sampling points was large in the places with large errors and small in the places with small errors. Yi et al. [22] discretized the triangular mesh of a free surface, and simplified the triangular mesh by iteratively shrinking the triangle edges. Based on the principle of minimum quadratic error, the optimal objective vertex under discrete curvature constraint is determined. By limiting the side length of the triangular mesh to control the number of sampling points, the adaptive sampling of free-form surfaces is realized. Yu et al. [23] selected the initial sampling point set, reconstructed the initial point set and solved it with the original surface to obtain the global error. The point with the largest global error is added to the initial point set to obtain the updated initial point set. The reconstructed surface and the original surface are solved to obtain the global error. The point with the maximum global error between the reconstructed surface and the original surface is added to the initial point set as a new sampling point. Cycle in turn until the global error reaches the set precision threshold or the number of sampling points reaches the set number threshold. Mian et al. [24] studied the influence of different sampling strategies on surface reconstruction accuracy, and the influence of workpiece size and machining quality on sampling methods. Gohari et al. [25] used principal component analysis to dynamically generate sampling points, thus reducing the number of sampling points, reducing the measurement cost and improving the measurement efficiency.
For the precise measurement of geometric quantities of free-form surfaces, the number of samples and the position of sampling points should be reasonably set to improve the measurement accuracy and efficiency. In this study, geodesic distance is used as the index of global fluctuation of a free-form surface. The free-form surface is divided into blocks, and the number of samples is adaptively determined according to the changes in the geometric characteristics of the free-form surface. Combined with the Isoparametric distribution, Poisson distribution, Hammersley distribution and NRook distribution, the distribution location of the sampling points is determined. Therefore, an adaptive hybrid sampling method for a free-form surface based on geodesic distance is proposed, which can effectively improve the measurement accuracy and efficiency.
The rest of this paper is organized as follows: Section 2 introduces the free-form surface modeling method; Section 3 introduces the definition of geodesic distance and its solution method in detail, including three sub-steps; in Section 4, the number of sampling points and the distribution of the sampling points are given; in Section 5, the reconstruction errors between the reconstructed surface and the original surface are solved and the results are analyzed; Section 6 summarizes the conclusions and future work prospects.

2. Free-Form Surface Modeling

Non uniform rational B-spline (NURBS) is one of the most commonly used parametric mathematical models in free-form surface geometric modeling. NURBS surfaces are widely used in the field of computer aided geometric design (CAGD), and are widely used in the geometric representation of complex components in aerospace, automobile and other fields.
The free-form surface can be represented by the control points and degrees in the  u  and  v  directions. The local adjustment of the free-form surface can be realized by adjusting the control points and weight coefficients. The NURBS surface is usually obtained by using the tensor product of two NURBS curves with two independent parameters,  u  and  v .
The expression of a NURBS surface is as follows [26].
S u ,   v = i = 0 m j = 0 n w i , j N i , p u N j , q v P i , j i = 0 m j = 0 n w i , j N i , p u N j , q v
where  P i , j  is the control point;  m  is the number of control points in the  u  direction;  n  is the number of control points in the  v  direction;  p  is the degree of the parameter coordinate  u q  is the degree of the parameter coordinate  v N i , p u  is a basis function of order  p N i , p v  is a basis function of order  q U  and  V  are defined knot vectors,  U = 0 , , 0 , u p + 1 , , u r p 1 , 1 , , 1 V = 0 , , 0 , v q + 1 , , v s q 1 , 1 , , 1 , which specify the distribution of parameters  u  and  v w i , j  is the same as the  P i , j  relevant weight coefficient. The basis function is recursively defined by the Cox–deBoor algorithm, where  N i , p u :
N i , 1 = 1   if   u i < u < u i + 1 0 otherwise
N i ,   p u = u u i u i + p 1 u i N i , p 1 u + u i + p u u i + p u i + 1 N i + 1 , p 1 u
The NURBS surface is shown in Figure 2. The blue curve is the isoparm in the  u  direction, the red curve is the isoparm in the  v  direction, and the red surface and the blue curve form are the isoparm mesh of the NURBS surface.

3. Geodesic Distance

On the Riemannian manifold, the geodesic is defined as the shortest path between the points on the model surface. Geodesic distance is the distance value of geodesic. Compared with the curvature, the advantage of geodesic is that it can describe the free-form surface globally, and the curvature can only be described based on the local part of a free-form surface [27]. Therefore, the geodesic distance is used as the index of the surface fluctuation of the free-form surface in this paper, which can accurately describe the global fluctuation of the free-form surface and overcome the limitation that curvature can only describe the local fluctuation.
The geodesic distance solution method used in this paper is the thermodynamic method proposed by Crane [28]. It can be imagined that a hot needle touches a point x on the surface, which is a hot core point. With the passage of time, the heat on this point x diffuses to the rest of the surface. The heat at point y on the surface can be expressed by a thermal kernel function    k t , x y    to describe how the geodesic distance between any point  x ,   y  on the Riemannian manifold φ x ,   y  can be recovered by the point state transformation of the thermal core, as shown in Formula (4). The geodesic distance can be recovered by solving the direction of thermal motion [29].
φ x , y = lim t 0 4 t log k t , x y
Step 1: Solve the thermal kinematics in Equation (5) by describing the propagation state of heat, and establish the temperature scalar field U. The time dispersion of the thermal propagation equation is
u t = u
The thermal kinematics equation is discretized and sorted to obtain
i d t u t = u 0
where  i d  is the identity matrix,  t  is the time interval,   is the discrete Laplacian operator,  u t  is the thermal state at time t, and  u 0  is the thermal state at the initial time.
Step 2: The thermal gradient direction calculated in step 1 is the same as the gradient direction of the geodesic distance. Since the gradient of the geodesic distance is a unit vector, the gradient of the geodesic distance is obtained by normalizing the thermal gradient direction.
  X = u / u
where  X  is the gradient of the geodesic distance.
Step 3: after the gradient of geodesic distance is obtained through step 2, the geodesic distance is solved through Formula (8):
min ϕ M ϕ X 2
According to the variational method, the minimization of Formula (8) is the Poisson equation.
Δ ϕ = X
where  ϕ  is the geodesic distance between the vertex and the hot core point.
Based on the principle of the NURBS surface in part 2 of this paper, the free-form surface to be measured is shown in Figure 3, with a size of 50 mm × 50 mm. Reduce the dimension of the free-form surface from the three-dimensional space to the two-dimensional parameter domain  u , v ; carry out Isoparametric sampling on the two-dimensional parameter domain  u , v ; map the sampling points on the two-dimensional parameter domain to the three-dimensional space; use the collection of sampling points on the three-dimensional space to replace the free-form surface, and set the number of sampling points to (41 × 41) = 1681.
In order to evaluate the global fluctuation of a free-form surface, it is necessary to segment the surface. The principle of the surface segment is shown in Figure 4 and Figure 5. As shown in Figure 4, the free-form surface to be measured is divided into four surface segments  P i  in the two-dimensional parameter domain. Select the (41× 41 + 1)/2 = 841 sampling points to be used as hot core points (the red point in Figure 4 needs to be mapped to the three-dimensional space). Since the sampling points in the same row and column as the hot core point belong to the overlapping part of the two segments, the sampling points in the same row and column as the hot core point will be deleted. Solve the geodesic distance from the sampling point to the thermal core point in each surface segment  P i  after deleting the duplicate sampling points; sum the geodesic distances from the sampling point to the thermal core point in each segment; obtain the sum of the geodesic distances  d i  of each segment sampling point, and take it as the global fluctuation change index of the free-form surface. The sum of the geodesic distances  d i  of each segment  P i  is shown in Table 1.
As shown in Figure 5, each surface segment  P i  is further subdivided into four surface segments  P i j , and the sum of the geodesic distances  d i j  from the sampling point to the hot core point (red point in Figure 5) in each segment  P i j  is calculated. The results are shown in Table 2. Finally, the sum of the geodesic distances  d i j  of each free-form surface segment is obtained, which is used as an indicator of global fluctuation. Based on this indicator, the sampling quantity of each free-form surface segment is determined.

4. Sampling Strategy

In this section, firstly, the sampling quantity of each surface segment is calculated based on the sum of the geodesic distances of each surface segment solved in Section 3, and then the positions of the sampling points are generated based on the Isoparametric distribution, Poisson distribution, Hammersley distribution and NRook distribution. Finally, the geometric error evaluation method of the reconstructed surface and the original surface is determined.

4.1. Determine the Sampling Quantity

In this paper, a relative proportion method based on the sum of the geodesic distances of each segment is used to determine the P of each surface segment. The specific solution steps are as follows.
Step 1: According to the results of Section 3, calculate the average value of the sum of the geodesic distances of  P 1 P 2 P 3  and  P 4 ; obtain the ratio  λ i  of the sum of the geodesic distances of  P 1 P 2 P 3  and  P 4  and  d i . The reference sampling number  N i ¯  of each surface segment  P i  is the ratio of the total number of samples  N  and the number of segments  φ . The sampling number of each surface segment  P i  is the product of the number of reference samples  N i ¯  and  λ i . The formula is as follows:
λ i = d i d ¯
N i ¯ = N φ
N i = N i ¯ λ i
Step 2: Each surface segment  P i  continue to subdivide into surface segment  P i j . According to the principle of step 1, calculate the average value  d i ¯  of the sum of the geodesic distances of  P i j ; obtain the ratio  λ i j  of the sum of the geodesic distances of  P i 1 P i 2 P i 3  and  P i 4  and  d i j . The reference sampling number  N i j ¯  of each surface segment  P i  is the ratio of the total number of samples  N i  and the number of segments  φ ; The sampling number of each surface segment  P i j  is the product of the number of reference samples  N i j ¯  and  λ i j . The sampling number of each surface segment  P i j  is the calculation result of step 1. The formula is as follows:
λ i j = d i j d i ¯
N i j ¯ = N i φ
N i j = N i j ¯ λ i j
According to the requirements of measurement accuracy, the total number of sampling points of the free-form surface is determined to be 1600,  φ  = 4. Calculate the number of sampling points of each surface segment as shown in Table 3.

4.2. Determine the Sampling Position

After determining the number of sampling points in each segment, a specific point distribution algorithm is used to distribute the sampling points. In this paper, the sampling points generated based on the Isoparametric distribution, Poisson distribution, Hammersley distribution and NRook distribution are studied.
(a)
Isoparametric distribution
Isoparametric distribution is to map the free-form surface from the three-dimensional space to the two-dimensional parameter domain, sample the two-dimensional parameter domain according to a certain step size to obtain the sampling points in the two-dimensional parameter domain, remap the sampling points in the two-dimensional parameter domain to the three-dimensional space, and obtain the set of sampling points in the three-dimensional space. It is calculated by Formulas (16) and (17):
u = u m i n + i 1 u m a x u m i n N u 1 ; i = 1 , 2 , , N u
v = v m i n + i 1 v m a x v m i n N v 1 ; i = 1 , 2 , , N v
where  u m i n  is the minimum value of the u direction parameter;  u m a x  is the maximum value of the u direction parameter;  N u  is the number of sampling points in the u direction;  v m i n  is the minimum value of the v direction parameter;  v m a x  is the maximum value of the v direction parameter;  N v  is the number of sampling points in the v direction.
(b)
Poisson distribution
Poisson distribution is a more uniform distribution mode compared to Isoparametric distribution. This distribution method can generate a random point set. The distribution method adopted is that every two points are at least on a specified minimum distance. The algorithm takes the range of the  R n  sample domain, the minimum distance  r  between samples and the constant  k  as the inputs. The steps are as follows [30].
Step 1: Initialize the n dimensional background grid, which is used to store variables and spatial search data. The search cell size is  r / n , and each grid cell contains, at most, one sample. Therefore, the grid is an n dimensional integer array. The default index (−1) indicates that there is no sample, and the non-negative integer indicates the index of the sample in the single grid.
Step 2: The initial sample  x 0  is randomly selected from the sample field, inserted into the background grid, and this uses index (0) to initialize the activity list (sample index list).
Step 3: When the active list is not empty, a random index is selected from the list, and evenly selects  k  points from the spherical ring between the surrounding radius  r    and the radius  2 r . For each point, check whether it is within the radius r of the existing sample (use the background grid to test only the nearby samples). If a point is far enough from the existing sample, it is taken as the next sample and added to the existing index. If no such point is found after  k  attempts, the index is removed from the active list.
(c)
Hammersley distribution
The Hammersley distribution is one of the most prominent uniform distribution sampling algorithms at present. This method is based on the computer binary number representation method, which converts a given decimal number into a binary number, inverts the order, multiplies each number on each bit of the binary by a power series with 1/2 as the base number and the corresponding number of digits as the exponent, and cumulatively sums them. The calculation result is placed after the decimal point to form the sampling value. The sampling formula is as follows:
u i = i / N
v i = j = 0 k 1 b i j 2 j 1
k = log 2 N
where  N  is the number of sampling points;  i  is the ith sampling point, and its range is  0 ,   N 1 b i  is the binary representation of the index;  b i j  is the jth bit representing  b i , and the range of  j  is  0 ,   k 1 k  is b number of digits of  b i .
(d)
NRook distribution
NRook distribution is a uniform distribution algorithm based on the principle of a chess board. If a square is divided into  n n  small squares and  n  sampling points are placed inside, there is only one sampling point in each row and column. The algorithm steps are as follows.
Step 1: Initialize the sampling set  Q  and the sampling step  λ  according to the number of samples  N . The calculation formula is as follows:
Q = 0 ,   1 , , N 1
  λ = 1 / N
Step 2: Selecting random samples  n i  from the sample set  Q ; select the random number    p i q i  from [0, 1],  u i  is the sum of the result of multiplying the sampling step  λ  by  n i  and a random number  p i v i  is the sum of the result of multiplying the sampling step  λ  by the sampling times  i  and a random number  q i . The calculation formula is as follows:
u i = λ n i + p i   v i = λ i + q i
According to the basic principles of the four distribution modes of Isoparametric distribution, Poisson distribution, Hammersley distribution and NRook distribution, the distribution sequence of parameter  u 0 , 1  and parameter  v 0 ,   1  is calculated and generated. The number of samples is set to 700. The results are shown in Figure 6.
The sampling points generated by the four distribution methods, i.e., Isoparametric distribution, Poisson distribution, Hammersley distribution and NRook distribution, need to be remapped to the corresponding sections according to the surface segment interval, and the sampling points are reconstructed. This is described in detail in Section 5.

4.3. Error Comparison Method between the Reconstructed Surface and Original Surface

The reconstructed surface and the original surface after surface reconstruction need to be evaluated for the geometric error of reconstruction. The method adopted in this paper is shown in Figure 7. The Isoparametric method is used to take a series of sampling points on the original surface, and make a straight line perpendicular to the z axis of the oversampling point. The straight line has an intersection with the reconstructed surface. The x and y coordinate values of the intersection and the sampling point are the same. The difference between the z coordinate values of the sampling point and the intersection is calculated, which is the reconstructed geometric accuracy of the reconstructed surface and the original surface.
By solving the z coordinate value difference between each sampling point on the original surface and the intersection point on the reconstructed surface, the root mean square error (RMSE) and the global maximum error (ME) are used as the reconstruction accuracy indexes of the reconstructed surface and the original surface. The formula is as follows:
MSE = 1 N i = 1 N z i z i 2
RMSE = 1 N i = 1 N z i z i 2  
ME = max z i z i  
where  N  is the number of sampling points;  z i  is the z coordinate value of the original surface sampling point;  z i  is the coordinate value of the reconstructed surface sampling point.

5. Experiment and Discussion

The free-form surface is reconstructed according to the different sampling strategies. The reconstructed surface is analyzed based on the reconstruction accuracy indexes, RMSE and ME, of the reconstructed surface and the original surface to evaluate the impact of the different sampling strategies on the reconstruction’s accuracy.
The parameter interval obtained by the four distribution methods of Isoparametric distribution, Poisson distribution, Hammersley distribution and NRook distribution is  0 ,   1 , and the length of the parameter u and v interval of the surface segment is 0.25. Therefore, it is necessary to map the boundary range of the divided surface to the corresponding interval. The interval range of each surface segment is shown in Figure 8.
According to the sampling number of each surface segment in Figure 8, four distribution methods, i.e., Isoparametric distribution, Poisson distribution, Hammersley distribution and NRook distribution, are used to calculate and generate sampling points, and map them to the corresponding interval range of each surface segment in Figure 8. The original free-form surface sampling point set obtained is used for the free-form surface reconstruction. The reconstructed surface is shown in Figure 9, where the green point represents the sample point, and the brown surface represents the reconstructed surface. Different sampling methods produce a different distribution of the sample points, and the reconstructed surface is also different.
The reconstructed surface is compared with the original surface in terms of the reconstructed geometric error, and the ME of each point  x , y , z  on the reconstructed surface is converted into the corresponding color. The obtained error distribution results are shown in Figure 10. Figure 10a is the reconstructed geometric error without surface blocking and other parameter distribution; Figure 10b is the reconstructed geometric error of Poisson distribution; Figure 10c is the reconstructed geometric error of Hammersley distribution; and Figure 10d is the reconstructed geometric error of NRook distribution.
Table 4 and Figure 11 shows the error comparison results between the reconstructed surface obtained by surface reconstruction and the original surface after the sampling points of the free-form surface are obtained by the different distribution methods. Comparing Isoparametric distribution from the aspect of reconstruction ME, the reconstruction error of the hybrid adaptive sampling method, based on geodesic distance combined with Poisson distribution, Hammersley distribution and NRook distribution, is reduced by 47.12%, 41.56% and 20.22%, respectively, compared with that of the Isoparametric distribution sampling method. Compared with the adaptive sampling method proposed in Reference [24], the reconstruction error of the hybrid adaptive sampling method based on the geodesic distance combined with Poisson distribution, Hammersley distribution and NRook distribution, is reduced by 69.83%, 66.67% and 54.49%, respectively. The reconstruction error of the adaptive hybrid sampling method, based on the geodesic distance combined with Poisson distribution, Hammersley distribution, NRook distribution and the adaptive sampling method used in Reference [14], is reduced by 62.46%, 62.15% and 52.61%, respectively, compared with Poisson distribution, Hammersley distribution and NRook distribution.
The adaptive hybrid sampling method based on geodesic distance significantly reduced the reconstruction ME compared with the general Isoparametric distribution sampling method and the adaptive sampling method in Reference [24], and the reconstruction accuracy has been greatly improved. The reconstructed ME based on the adaptive hybrid sampling method of geodesic distance is significantly improved compared with the sampling methods of Poisson distribution, Hammersley distribution and NRook distribution that are not based on the adaptive hybrid sampling method of geodesic distance, which verifies the effectiveness of the adaptive hybrid sampling method based on geodesic distance in reducing the reconstructed ME.
From the aspect of reconstruction RMSE, the adaptive hybrid sampling method based on geodesic distance combined with Poisson distribution, Hammersley distribution and NRook distribution, reduces the reconstruction RMSE by 24.06%, 28.95% and 20.00%, respectively, compared with the sampling method of Isoparametric distribution. Compared with the adaptive sampling method proposed in Reference [24], the reconstructed RMSE of the hybrid adaptive sampling method based on geodesic distance combined with Poisson distribution, Hamersley distribution and NRook distribution is reduced by 58.44%, 61.11% and 56.17%, respectively. The reconstructed RMSE of the adaptive hybrid sampling method based on geodesic distance combined with Poisson distribution, Hammersley distribution and NRook distribution is reduced by 56.75%, 73.31% and 82.63%, respectively, compared with the sampling methods of Poisson distribution, Hammersley distribution and NRook distribution.
The RMSE reconstruction of the adaptive hybrid sampling method based on the geodesic distance combined with Poisson distribution, Hammersley distribution and NRook distribution is greatly reduced compared with the sampling method of Isoparametric distribution and the adaptive sampling method in Reference [24], and the reconstruction accuracy is significantly improved. The RMSE of the adaptive hybrid sampling method based on geodesic distance combined with Poisson distribution, Hammersley distribution and NRook distribution is significantly reduced compared with the sampling method based on Poisson distribution, Hammersley distribution and NRook distribution without geodesic distance, and the reconstruction accuracy is greatly improved.
In order to further verify the effectiveness of the proposed method, this paper uses an ABB IRB1200 robot and Creaform MetraSCAN-R BLACK for experimental verification. The measurement accuracy is  25   µ m , and the measurement depth is  250   mm , as shown in Figure 12.
Set the number of sample points to 100, and use the self-developed path algorithm to obtain the robot scanning measurement path, as shown in Figure 13. The free-form surface is scanned and measured. The measured results are shown in Figure 14.
In order to further verify the effectiveness of the sampling algorithm, repeated measurement experiments are carried out. Taking the first measurement data as the error judgment standard, analyze the error of the two measurement experiments. The error results are shown in Figure 15. The average deviation is 0.001 mm and the standard error is 0.013 mm, which verifies the effectiveness and robustness of the sampling method proposed in this paper.

6. Conclusions

In this paper, an adaptive hybrid sampling method for a free-form surface based on geodesic distance is proposed, which can be used for the precise measurement of geometric parameters of free-form surfaces. This method can effectively improve measurement efficiency.
The geodesic distance is introduced as a measure of the global fluctuation of the free-form surface. The free-form surface is divided into multiple surface segments. The sampling number of each surface segment is determined according to the sum of the geodesic distances of each surface segment within the surface segment. The sampling points of each surface segment are generated based on the Poisson distribution, Hammersley distribution and NRook distribution sampling methods. Finally, the sampling points of the overall free-form surface are obtained.
The reconstruction error of the free-form surface obtained by this method was compared with the common sampling strategies. The results show that the adaptive hybrid sampling method of a free-form surface based on geodesic distance can effectively reduce the reconstruction error and significantly improve the reconstruction accuracy of free-form surfaces.
In future work, the sampling distribution method can be self-optimized by combining intelligent algorithms, to further improve the reconstruction accuracy of free-form surfaces.

Author Contributions

Methodology, C.C. and L.Y.; Software, C.C. and L.Y.; Investigation, All; Writing—original draft, C.C.; Writing—review and editing, C.C., H.J. and L.Y.; Supervision, L.Y.; Project administration, L.Y.; Funding acquisition, L.Y. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Natural Science Foundation of China (Grant Nos. 51927811/51875165/52005147), and the National Key Research and Development Program of China(2019YFE0107400/2022YFF0705700).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data available on request due to restrictions, e.g., privacy or ethical.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. He, G.Y.; Sang, Y.C.; Wang, H.L.; Sun, G.M. A profile error evaluation method for freeform surface measured by sweep scanning on CMM. Precis. Eng.-J. Int. Soc. Precis. Eng. Nanotechnol. 2019, 56, 280–292. [Google Scholar] [CrossRef]
  2. Abdullah, A.B.; Sapuan, S.M.; Samad, Z. Profile Measurement Based on Focus Variation Method for Geometrical Defect Evaluation: A Case Study of Cold Forged Propeller Blade. Adv. Mech. Eng. 2014, 6, 874691. [Google Scholar] [CrossRef]
  3. Shen, Y.; Zhang, W.; Zhang, Y.; Huang, N.; Zhang, Y.; Zhu, L. Distributed Particle Swarm Optimization for the Planning of Time-optimal and Interference-free Five-axis Sweep Scanning Path. IEEE Trans. Ind. Inform. 2022, 18, 8703–8713. [Google Scholar] [CrossRef]
  4. Zhang, Y.; Tang, K. Automatic Sweep Scan Path Planning for Five-Axis Free-Form Surface Inspection Based on Hybrid Swept Area Potential Field. IEEE Trans. Autom. Sci. Eng. 2019, 16, 261–277. [Google Scholar] [CrossRef]
  5. Zhang, Y.; Zhou, Z.; Tang, K. Sweep scan path planning for five-axis inspection of free-form surfaces. Robot. Comput.-Integr. Manuf. 2018, 49, 335–348. [Google Scholar] [CrossRef]
  6. Zhang, Y.; Zhu, L.; Zhao, P.; Hu, P.; Zhao, X. Skeleton Curve-Guided Five-Axis Sweep Scanning for Surface With Multiple Holes. IEEE Trans. Autom. Sci. Eng. 2021, 19, 2471–2486. [Google Scholar] [CrossRef]
  7. Zhang, Y.; Shen, Y.; Zhang, W.; Zhang, L.; Zhu, L.; Zhang, Y. Generation of efficient and interference-free scanning path for inspecting impeller on a cylindrical CMM. Measurement 2022, 198, 111352. [Google Scholar] [CrossRef]
  8. Rishi, K.; Gupta, S.K. Planning algorithms for acquiring high fidelity pointclouds using a robot for accurate and fast 3D reconstruction. Robot. Comput.-Integr. Manuf. 2022, 78, 102372. [Google Scholar] [CrossRef]
  9. Huo, X.; Tran, H.D.; Shilling, K.M.; Kim, H. Determination and Optimization of Spatial Samples for Distributed Measurements; Sandia National Laboratories (SNL): Albuquerque, NM, USA; Livermore, CA, USA, 2010.
  10. Dunbar, D.; Humphreys, G. A spatial data structure for fast Poisson-disk sample generation. ACM Trans. Graph. 2006, 25, 503–508. [Google Scholar] [CrossRef]
  11. Wang, T.; Suda, R.; Assoc Comp, M. Fast Maximal Poisson-Disk Sampling by Randomized Tiling. In Proceedings of the High-Performance Graphics (HPG) Conference, Los Angeles, CA, USA, 28–30 July 2017. [Google Scholar]
  12. Woo, T.C.; Liang, R.; Hsieh, C.C.; Lee, N.K. Efficient sampling for surface measurements. J. Manuf. Syst. 1995, 14, 345–354. [Google Scholar] [CrossRef]
  13. Wang, J.; Jiang, X.; Blunt, L.A.; Leach, R.K.; Scott, P.J. Intelligent sampling for the measurement of structured surfaces. Meas. Sci. Technol. 2012, 23, 085006. [Google Scholar] [CrossRef] [Green Version]
  14. Ren, M.; Kong, L.; Sun, L.; Cheung, C. A Curve Network Sampling Strategy for Measurement of Freeform Surfaces on Coordinate Measuring Machines. IEEE Trans. Instrum. Meas. 2017, 66, 3032–3043. [Google Scholar] [CrossRef]
  15. Zahmati, J.; Amirabadi, H.; Mehrad, V. A hybrid measurement sampling method for accurate inspection of geometric errors on freeform surfaces. Measurement 2018, 122, 155–167. [Google Scholar] [CrossRef]
  16. Gao, F.; Pan, Z.; Zhang, X.; Li, Y.; Duan, J. An adaptive sampling method for accurate measurement of aeroengine blades. Measurement 2021, 173, 108531. [Google Scholar] [CrossRef]
  17. Sang, Y.; Yan, Y.; Yao, C.; He, G. A new scanning lines distribution strategy for the form error evaluation of freeform surface on CMM. Measurement 2021, 188, 109578. [Google Scholar] [CrossRef]
  18. Jiang, R.S.; Wang, W.H.; Zhang, D.H.; Wang, Z.Q. A practical sampling method for profile measurement of complex blades. Measurement 2016, 81, 57–65. [Google Scholar] [CrossRef]
  19. Obeidat, S.M.; Raman, S. An intelligent sampling method for inspecting free-form surfaces. Int. J. Adv. Manuf. Technol. 2009, 40, 1125–1136. [Google Scholar] [CrossRef]
  20. Mansour, G. A developed algorithm for simulation of blades to reduce the measurement points and time on coordinate measuring machine (CMM). Measurement 2014, 54, 51–57. [Google Scholar] [CrossRef]
  21. He, G.Y.; Sang, Y.C.; Pang, K.R.; Sun, G.M. An improved adaptive sampling strategy for freeform surface inspection on CMM. Int. J. Adv. Manuf. Technol. 2018, 96, 1521–1535. [Google Scholar] [CrossRef]
  22. Yi, B.; Qiao, F.; Huang, N.; Wang, X.; Wu, S.; Biermann, D. Adaptive sampling point planning for free-form surface inspection under multi-geometric constraints. Precis. Eng.-J. Int. Soc. Precis. Eng. Nanotechnol. 2021, 72, 95–101. [Google Scholar] [CrossRef]
  23. Yu, M.; Zhang, Y.; Li, Y.; Zhang, D. Adaptive sampling method for inspection planning on CMM for free-form surfaces. Int. J. Adv. Manuf. Technol. 2013, 67, 1967–1975. [Google Scholar] [CrossRef]
  24. Mian, S.H.; Al-Ahmari, A.; Alkhalefah, H. Analysis and Realization of Sampling Strategy in Coordinate Metrology. Math. Probl. Eng. 2019, 2019, 9574153. [Google Scholar] [CrossRef] [Green Version]
  25. Gohari, H.; Barari, A. A quick deviation zone fitting in coordinate metrology of NURBS surfaces using principle component analysis. Measurement 2016, 92, 352–364. [Google Scholar] [CrossRef]
  26. Piegl, L.; Tiller, W. The NURBS Book; Springer: Berlin/Heidelberg, Germany, 1997. [Google Scholar]
  27. Quan, L.; Tang, K. Polynomial local shape descriptor on interest points for 3D part-in-whole matching. Comput.-Aided Des. 2015, 59, 119–139. [Google Scholar] [CrossRef]
  28. Crane, K.; Weischedel, C.; Wardetzky, M. Geodesics in Heat: A New Approach to Computing Distance Based on Heat Flow. ACM Trans. Graph. 2013, 32, 1–11. [Google Scholar] [CrossRef] [Green Version]
  29. He, D.; Li, Y.; Li, Z.; Tang, K. Geodesic Distance Field-Based Process Planning for Five-Axis Machining of Complicated Parts. J. Manuf. Sci. Eng.-Trans. ASME 2021, 143, 061009. [Google Scholar] [CrossRef]
  30. Bridson, R. Fast Poisson disk sampling in arbitrary dimensions. In Proceedings of the SIGGRAPH Sketches, San Diego, CA, USA, 5–9 August 2007. [Google Scholar]
Figure 1. Two measurement methods. (a) CMM; (b) Robot optical scanning system.
Figure 1. Two measurement methods. (a) CMM; (b) Robot optical scanning system.
Sensors 23 03224 g001
Figure 2. NURBS surface and its Isoparametric mesh.
Figure 2. NURBS surface and its Isoparametric mesh.
Sensors 23 03224 g002
Figure 3. Free-form surface.
Figure 3. Free-form surface.
Sensors 23 03224 g003
Figure 4. Schematic diagram of the geodesic distance of the first surface segment.
Figure 4. Schematic diagram of the geodesic distance of the first surface segment.
Sensors 23 03224 g004
Figure 5. Schematic diagram of the segment geodesic distance of the second surface.
Figure 5. Schematic diagram of the segment geodesic distance of the second surface.
Sensors 23 03224 g005
Figure 6. Sampling distribution (a) Isoparametric distribution; (b) Poisson distribution; (c) Hammersley distribution; (d) NRook distribution.
Figure 6. Sampling distribution (a) Isoparametric distribution; (b) Poisson distribution; (c) Hammersley distribution; (d) NRook distribution.
Sensors 23 03224 g006
Figure 7. Schematic diagram of the reconstruction error surface.
Figure 7. Schematic diagram of the reconstruction error surface.
Sensors 23 03224 g007
Figure 8. Interval range and sampling quantity of the surface segment.
Figure 8. Interval range and sampling quantity of the surface segment.
Sensors 23 03224 g008
Figure 9. Sampling point distribution of the free-form surface and surface reconstruction: (a) Isoparametric; (b) Poisson; (c) Hammersley; (d) NRook; (e) Adaptive hybrid Poisson; (f) Adaptive hybrid Hammersley; (g) Adaptive hybrid NRook.
Figure 9. Sampling point distribution of the free-form surface and surface reconstruction: (a) Isoparametric; (b) Poisson; (c) Hammersley; (d) NRook; (e) Adaptive hybrid Poisson; (f) Adaptive hybrid Hammersley; (g) Adaptive hybrid NRook.
Sensors 23 03224 g009
Figure 10. Reconstructed surface ME distribution map: (a) Isoparametric; (b) Poissson; (c) Hammersley; (d) NRook; (e) Adaptive hybrid Poisson; (f) Adaptive hybrid Hammersley; (g) Adaptive hybrid NRook.
Figure 10. Reconstructed surface ME distribution map: (a) Isoparametric; (b) Poissson; (c) Hammersley; (d) NRook; (e) Adaptive hybrid Poisson; (f) Adaptive hybrid Hammersley; (g) Adaptive hybrid NRook.
Sensors 23 03224 g010
Figure 11. Comparison of the reconstruction errors using different sampling methods.
Figure 11. Comparison of the reconstruction errors using different sampling methods.
Sensors 23 03224 g011
Figure 12. Measuring the surface using a MetraSCAN laser scanner.
Figure 12. Measuring the surface using a MetraSCAN laser scanner.
Sensors 23 03224 g012
Figure 13. Robot scanning measurement path.
Figure 13. Robot scanning measurement path.
Sensors 23 03224 g013
Figure 14. Point cloud reconstruction results of the free-form surface.
Figure 14. Point cloud reconstruction results of the free-form surface.
Sensors 23 03224 g014
Figure 15. Error comparison results.
Figure 15. Error comparison results.
Sensors 23 03224 g015
Table 1. Sum of the geodesic distances of the first surface segment.
Table 1. Sum of the geodesic distances of the first surface segment.
SegmentDistance/mm
17750
27816
37578
47171
Table 2. Sum of the geodesic distances of the second quadric surface segment.
Table 2. Sum of the geodesic distances of the second quadric surface segment.
SegmentDistance/mm
114207
124305
133983
143714
214174
224186
233824
243728
313546
323442
334059
344004
413366
423427
434069
444163
Table 3. Number of sampling points for each surface segment.
Table 3. Number of sampling points for each surface segment.
SegmentSample Size
11131
12144
13119
14106
21181
22194
23169
24156
3169
3256
3394
3481
416
4219
4331
4444
Table 4. Comparison of the reconstruction errors using different sampling methods.
Table 4. Comparison of the reconstruction errors using different sampling methods.
Sampling MethodME/mmRMSE/mm
Isoparametric0.18550.0266
Poisson0.26130.0467
Adaptive hybrid Poisson0.09810.0202
Hammersley0.28640.0708
Adaptive hybrid Hammersley0.10840.0189
NRook0.31230.1226
Adaptive hybrid NRook0.14800.0213
Adaptive sampling method [24]0.32520.0486
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Chen, C.; Jia, H.; Lu, Y.; Zhang, X.; Chen, H.; Yu, L. An Adaptive Hybrid Sampling Method for Free-Form Surfaces Based on Geodesic Distance. Sensors 2023, 23, 3224. https://doi.org/10.3390/s23063224

AMA Style

Chen C, Jia H, Lu Y, Zhang X, Chen H, Yu L. An Adaptive Hybrid Sampling Method for Free-Form Surfaces Based on Geodesic Distance. Sensors. 2023; 23(6):3224. https://doi.org/10.3390/s23063224

Chicago/Turabian Style

Chen, Chen, Huakun Jia, Yang Lu, Xiaodong Zhang, Haohan Chen, and Liandong Yu. 2023. "An Adaptive Hybrid Sampling Method for Free-Form Surfaces Based on Geodesic Distance" Sensors 23, no. 6: 3224. https://doi.org/10.3390/s23063224

APA Style

Chen, C., Jia, H., Lu, Y., Zhang, X., Chen, H., & Yu, L. (2023). An Adaptive Hybrid Sampling Method for Free-Form Surfaces Based on Geodesic Distance. Sensors, 23(6), 3224. https://doi.org/10.3390/s23063224

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop