Identification of the geometric design parameters of propeller blades from 3D scanning

This paper presents a new method for the identification of geometric design parameters of ISO 484 class propeller blades from scanned point cloud data. The method can be used for tolerance inspection and in-line measurement of manufactured blades, and for wear assessment and reverse engineering of propellers in service. The geometry of the propeller blades is specified by a set of blade sections stacked radially on concentric cylindrical surfaces. For each of these sections, the chord line and mean camber line is determined, and geometric design parameters such as the pitch, skew, chord length, camber, and thickness distributions are identified. The proposed method is a complete procedure for identifying the design parameters of marine propeller blades from point cloud data, and includes a novel method for the precise identification of the mean camber line based on Voronoi diagrams and Delaunay triangulation. The paper includes validation of the proposed method in experiments where reverse engineering is applied to a propeller blade of the KVLCC2 propeller, and in 3D scanning of a large, high-skew thruster blade where the results were compared to CMM measurements.


Introduction
The surface profile of a propeller blade has a vital impact on the efficiency of the propeller. It is, therefore, important that accurate measurement techniques are available in the manufacturing of high-performance marine propellers. Evaluation of the geometric errors is typically carried out through the various stages of the production process and for the final inspection of products. Due to the highly complex surface geometry of the propeller blades, inspection is often time-consuming and involved.
Coordinate measuring machines (CMM) are typically used for propeller inspection where discrete points on the surface are sensed with a probe. Since the measurements are conducted pointwise, the efficiency of CMM inspection is relatively low. Propeller manufacturing is often characterized by small batch production with frequent changes in product dimensions, which means that time-consuming programming and set-up of coordinate measuring machines 1 3 must be carried out frequently. CMM inspection processes can also be challenging to integrate with other automated production processes, such as robotic grinding or polishing. Consequently, it is interesting to consider alternative inspection methods that are potentially less time-consuming, less expensive, fully automatic, and still offer an adequate level of accuracy.
Optical measurement technologies have proved increasingly applicable for the inspection of manufactured parts. For marine propeller blades, optical measurement systems offer several advantages. The planning and programming of an inspection task can be simplified, and by quickly and efficiently collecting a massive number of data points from the surface of the part, the time required to inspect a manufactured propeller can be significantly reduced. Optical systems of small size can be used in combination with an industrial robot for inspection of large products, which is interesting for flexible integration in automated production systems.
A variety of methods for optical measurement of propeller blade geometry has previously been proposed and studied. Almost three decades ago, Golini [1] proposed using a profilometer to measure the shape of a propeller blade. The system was capable of recording ten measurements per second using laser interferometry. Menna and Troisi [2] studied the possibility of using various low-cost approaches based on photogrammetry and laser scanning for generating complete and accurate 3D CAD models of small propellers. Zhang and Lu [3] proposed a solution for measuring the width of marine propellers using parallel axis stereo vision. Focus variation, a technique mainly used for measuring geometrical defects such as surface waviness and roughness, gave 3D depth readings. As proposed by Abdullah et al. [4], this method could also be used to evaluate geometric and dimensional defects of small propeller blades. Cavada and Fadón [5] gave a thorough assessment of various laser-based techniques for measuring marine propellers. Both direct measurements with laser and combinations of laser measurement and surface probing were considered. Beng and Choong [6] used an unspecified type of 3D scanner and computer-aided design software to characterize the geometry of a small outboard marine propeller. It is worth mentioning that other sectors of the shipbuilding industry also see potential applications of optical inspection methods. Examples of this include the as-built measurement of prototype ship design [7] and reverse engineering of hull design for energy efficiency evaluation [8]. Reverse engineering of an old propeller is described by Hand et al. [9], where structured light and coherent laser radar scanning sensors were used to recreate the propeller of a ship constructed during the American Civil War.
A marine propeller blade is, in principle, a rotating airfoil that produces thrust. Airfoils are also found in fans, compressors, and turbines, which makes it relevant to consider studies on optical blade inspection for other sectors. These studies include extraction of the sectional foil curves [10], reconstruction of mean camber line [11], leading edge and trailing edge localization [12], shape matching [13], and position and orientation error evaluations.
There exist some challenges for optical measurement techniques when applied to the inspection of marine propellers. Difficulties with highly reflective surfaces, requirements for unobstructed line of sight, imprecise stitching when joining results from multiple scans, and lack of structure in the resulting point clouds are all examples of these challenges. However, the benefits of the rapid non-contact acquisition of dense spatial measurements have outweighed these difficulties and inspired continued advancement.
To extract the required design parameters of the propeller blades, precise and efficient processing of the massive amount of data arising from optical 3D scanning is necessary. Several efforts have previously addressed and studied the problem of identifying the most decisive design parameters of propellers from scanned data. Patrikalakis and Bardis [14] presented a set of algorithms for the extraction of gross geometrical features of marine propeller blades represented in terms of B-spline surfaces. Their approach for camber line calculation requires an involved integration of a system of differential equations along with a complicated error evaluation scheme. Jinkerson et al. [15] extended the work of [14] with useful methods for unconstrained and constrained localization and feature extraction for propellers. Their method requires a non-uniform rational basis spline (NURBS) representation of the blade geometry in addition to a set of corresponding spatial point measurements of the manufactured blade surface. The authors also state that tedious and manual preprocessing of the measured data points is necessary for the localization methods to succeed. Allen et al. [16] proposed an approach to this problem where the measured point clouds were converted into a surface representation of the blade, and performed a series of NURBS curve interpolations. As for the methods in Ref. [15], the approach requires accurate NURBS representation of the blade geometry and, in addition, the computation of the signed deviation from the design surface for each measured point. More recently, Yeo and Choong [6] characterized the design parameters for a small outboard marine propeller. However, their method is also based on interpolating a CAD model of the scanned data, which is a complicated step that is prone to errors and involves several manual steps.
To correctly identify the propeller design parameters, a set of particular sectional curves need to be precisely extracted. Intermediate surface or mesh reconstruction steps will introduce an interpolated and imprecise surface, which effectively reduces the accuracy of the inspection of a point by taking in point deviations from connecting surface patches. Direct slicing and identification of the foil sections from the point cloud data would therefore increase the accuracy of the inspection methods, by rejecting the influence from measurement deviations in the regions close to each point. Li et al. [17], Kostas et al. [18], and several others have proposed airfoil models that can be used to reconstruct airfoils and hydrofoils from possibly noisy measurements of foil sections. A comprehensive survey and comparison of airfoil shape parameterizations are given by Masters et al. [19]. Consistent reconstruction of surfaces and boundary curves requires that the measured points take the form of organized point clouds, which is not straightforward to achieve for standard optical measurement techniques.
The geometric shape of propeller blades can be approximated by parametric blade models, commonly employed for parametric design studies and performance optimization. Several parametric propeller blade models have been proposed in the literature, and a recent example is the work of Arapakopoulos et al. [20]. An alternative approach to the problem of identifying the design parameters from scanned data could be to set up and solve a minimization problem, where the deviation between a parametric model and a scanned point cloud of a blade is minimized. However, solving such a minimization problem for the large amount of data generated from 3D scanning is computationally intensive.
In this paper, we propose a new method for the identification of the propeller design parameters using point cloud data originating from 3D cameras. The main idea of the proposed method is to extract the sectional contours of propeller blades directly from the point cloud and to determine the mean camber line, the leading and trailing edges, and the chord line for each of the sections. The identification of these parameters allows for the determination of the original design parameters of the propeller with high precision. The success of the approach is subject to the accuracy of the camber line determination step, which we will cover in detail. Emphasis has been placed on complying with the globally recognized standard for marine propeller tolerances, ISO 484.
The performance of the approach is validated in experiments where reverse engineering is performed for a propeller blade of the KVLCC2 design. Further validation is performed where we identify the design parameters of a large, high-skew propeller blade from a 3D scanning of its surface, using the ATOS ScanBox developed by GOM, and compare the results with the actual propeller design data and CMM measurements.
The organization of the paper is as follows: Sect. 2 discusses the geometry and representation method of marine propellers and gives an overview of the ISO 484 standard. The proposed method for extracting the sectional curves of the propeller blade is presented in Sect. 3, where a discussion on mean camber line determination is included. Case studies and results are presented in Sect. 4. Finally, conclusions and suggestions for further work are given in Sect. 5.

The geometry of marine propellers
This section presents a formalism for the description of marine propeller geometry based on [21]. A global reference frame (X, Y, Z) is defined, as shown in Fig. 1, where the origin is in the center of the propeller hub. The X-axis is along the propeller shaft in the forward direction of the vessel. The Y-axis is positive toward the port side, and the positive Z-axis is pointing in the vertically upward direction of the vessel. In addition, the local reference frame (x r , y r , z r ) is defined as a body-bound frame rotating about the X-axis. A typical marine propeller will have three to six blades attached to a boss or hub, which in turn is mounted on the propeller shaft. Each propeller blade is defined relative to a propeller reference line, which is fixed in the blade with direction normal to the propeller shaft. There is one propeller reference line for each blade. In the description of a propeller blade, the z r axis of the local reference frame can be aligned with the propeller reference line of the blade.
The quantities that define the overall propeller geometry are its diameter, pitch, pitch distribution, chord length distribution, number of blades, camber distribution, and thickness distribution. Additional parameters of importance are the skew and rake of the propeller and the diameter of its hub.

Propeller diameter
The diameter D of a propeller is defined as the diameter of the circle it spans while rotating, seen directly from the back of the vessel. This concept is relatively straightforward for Fig. 1 Global reference frame. The X axis is along the propeller shaft in the forward direction propellers with fixed pitch. For some classes of propellers, the controllable pitch propellers (CPP), however, the spanning circle will vary depending on the mechanically varied pitch of the blades. In this case, the concept of a design pitch is practiced. The diameter of a CPP is then determined by orienting the propeller blade, such that the specified blade section pitch is equal to the design pitch. Figure 2 illustrates the propeller diameter and accompanying radius.

Pitch and pitch distribution
The propeller pitch is the theoretical distance P a propeller will move forward in a solid medium for each rotation, not allowing for slip. The distance is illustrated in Fig. 3. The pitch of a blade section can also be defined as an angle. The relation between the pitch distance P and the pitch angle is given by Typically, the propeller blades are twisted to achieve almost constant pitch of the blades from root to tip. To describe this twist, a pitch distribution P(r) is given, where r is the radius considered. The relationship between pitch angle and pitch distribution can then be expressed as where R = D∕2 is the propeller radius and is the nondimensional radius defined by Figure 4 illustrates the geometric relationship between the pitch distance and angle for a blade section near the hub, and a section near the tip of the blade. The pitch angle will typically be decreasing with increasing radius to keep the pitch distance P nearly constant. It is common to refer to the pitch located at 70% of the total propeller radius, which is at = 0.7.

Fig. 2
The diameter D is the most apparent feature of any propeller and is simply the circular extent that the extreme tips of the propeller blades span as the propeller rotates

Cylindrical foil sections
The geometry of a propeller blade is described by selecting a set of N values 1 , 2 , … , N for the non-dimensional radius = r∕R . Then, the geometry of the propeller blade is given by the N cylindrical cross sections of the blade at the non-dimensional radii = i , where the cylindrical cross sections are concentric with the propeller axis. At each considered radius, the propeller blade has a different cross section representing the distribution of pitch, thickness, and other design parameters. Figure 5 illustrates how a cylindrical propeller section resembles the cross section of a foil.
Let the z r -axis of the rotating local reference frame be along the propeller reference line of the blade. To describe the foil parameters of the blade, the cylindrical sections are unwrapped onto a planar surface parallel to the x r y r plane of the local reference frames, as shown in Fig. 6.
The unwrapped blade sections are also termed expanded blade sections. The shape of the unwrapped blade sections is defined by the geometric foil parameters such as thickness, camber, chord length, and the leading and trailing edges. Figure 7 shows the unwrapped foil, with the main terminology used for describing the foil sections. Precise identification of the leading edge (LE) and trailing edge (TE) is necessary for determining the chord line and hence the center point of the blade section. The blade reference line, which is plotted for the expanded blade in Fig. 6, is the line through the midpoint of the camber line of the sections.

Chord length distribution
At a given radius, the leading and trailing edges are connected via a helical line segment, called the chord line. For the unwrapped sections, the chord line becomes a straight line and is also called the nose-to-tail line. The length of the chord line is termed the chord length l, as shown in Fig. 7. The chord length distribution l(r i ) is a set of chord lengths at predefined sections of a propeller blade at r = r i = i R.
Through the chord length distribution, the unfolded blade area of the propeller can be expressed as where Z is the number of blades, r 0 is the radius of the root section, and R is the propeller radius. The blade area is a parameter that affects the propeller's efficiency and thrust to a large extent.
The midpoint of the chord line is called the mid-chord, and is commonly used for defining the skew and rake of each propeller blade section.

Camber and thickness distribution
The chord line divides the blade section into an upper and lower surface. If we trace the midpoints between the upper and lower surfaces, we get a curve called the mean camber line. The distance between the chord line and the mean camber line is called the camber of the section. In Fig. 7, the concepts of the camber and mean section camber lines are illustrated. A camber distribution f(l) can be developed by measuring the camber at several points along the chord line of a section. The camber is then measured orthogonal to the chord line. The camber ratio is expressed as where f 0 is the maximum camber, and l is the chord length of the foil section.
The thickness distribution t(l) specifies the distance between the pressure and suction sides of a blade at points along the chord line, as illustrated in Fig. 7. The maximum thickness t 0 is usually between 1/3 and 1/2 chord lengths from the leading edge.

Skew and rake
The skew angle s (r) of a particular section is the tangential component of the angle formed on the propeller blade between the propeller reference line (z-axis) and a radial line from the propeller origin, passing through the mid-chord of the considered section at radius r. The purpose of a skewed shape is to avoid cavitation while maintaining propeller efficiency [21]. The centerline of each blade is then swept curvilinearly backward of the direction of rotation, and as a result, the contour of the blade is not radially symmetric around its centerline. The asymmetrical contour of a skewed propeller blade is shown in Fig. 8. The propeller skew angle sp is the largest spanning angle between two lines from the propeller origin running through the various mid-chords of the radial blade sections.
For a given foil section at radius r, the skew can be expressed as where s (r) is the skew angle, and is the pitch angle for the considered foil section. The propeller skew can be divided into forward and aft skew. A propeller for which the blades are swept in the direction of rotation is said to have forward skew. An aft skewed propeller has blades swept opposite of the direction of rotation. These concepts are illustrated in Fig. 8. Propeller rake is defined as the displacement of the radial foil section along the shaft or x-axis. Displacements in the aft direction are considered as a positive rake, and the rake measured at the blade tip is frequently used as a rake measure. The rake comes from two components, namely the initial rake i in and the skew-induced rake i s . The total rake i tot is then The skew-induced rake is an axial displacement caused by the skew imposed on a blade section. For a given skew angle (r) and pitch angle , the skew-induced rake can be determined by In Fig. 9, the relation between the pitch, skew, and types of rake impacting a propeller blade section is illustrated. The figure also shows how the chord line and the mid-chord is central in determining these parameters.

ISO 484
Marine propellers are almost exclusively manufactured, repaired, inspected, and certified according to the International Standard ISO 484 from 1981. The standard has four tolerance classes, regulating different levels of manufacturing accuracy. When a propeller is inspected, vital measurements are made on each blade of a propeller at specified radii from the propeller axis. Most often, these measurements are conducted at: some radii near the hub, then at the non-dimensional radii 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 0.95. For each of the radii, the blade surfaces are measured at a set of prescribed points. ISO 484 explicitly specifies tolerances on the propeller design pitch, pitch distribution, extreme radius of the propeller, thickness distribution, blade surface roughness, and static unbalance. Also, the standard describes the required accuracy of the measuring equipment used for inspection.
Cavada and Fadón [5] pointed out that manually operated measurement systems, where various points on the propeller's surface are measured by direct probing, are generally not accepted by the classification societies, as inspection of the exact points defined by the ISO 484 standard cannot be guaranteed.

Overview of the proposed method
In this section, our proposed method is presented. The objective is to develop a fully automated approach for deriving geometric design parameters from a given point cloud. The fundamental idea of this paper is as follows. First, we extract the radial sections of a propeller blade directly from scanned point cloud data and transform them into a corresponding planar surface. Then, the foil geometry parameters are identified for each radial section by interpolating the boundary curve and finding the mean camber line, which is found by generating the Voronoi diagram and Delaunay triangulation of the section points. The identified geometric parameters are compiled, and the overall design parameters Blade center axis Fig. 9 The relationship between the propeller shaft axis, the propeller reference line, the blade reference line, and the blade center axis for each of its blades. Note how skew and rake affect the distance between the mentioned axes 1 3 of the propeller are developed. The success of the approach depends on the accuracy of the identified foil geometry parameters for each section, with particular stress on accurately locating the mean camber line and leading and trailing edges. Because the data points generally have uneven distribution and are noisy, systematically locating said parameters is a challenging task. Another prerequisite for successful identification is that the propeller or propeller blade is oriented correctly according to its datum.

Section extraction and unwrapping
For each of the selected radial sections, we construct a cylindrical surface of radius r equal to the considered radius and select all points coinciding with it. Since only few points are at the exact radial distance, all data points within a userdefined distance of from the sectional cylinder surface are radially projected onto the cylinder surface. In the propeller-fixed local reference frame, the selected points (x r , y r , z r ) have to fulfill the following condition: This results in a set of n i points (x r i , y r i , z r i ) , i = 1, … , n i for each radius r i = i R.
As described in the previous section, it is necessary to unwrap the radial sections onto a planar surface parallel to the xy plane of the local reference plane to identify the various foil parameters. In other terms, we need to perform a transformation that preserves distances between points. A point (x r , y r , z r ) laying in the extracted radial section with coordinates in the local reference frame will have coordinates (x, y, z) in the unwrapped planar section where The transformation between radial and flat sections is illustrated in Fig. 6.

Interpolation of unwrapped boundary curves
Following the extraction and unwrapping of the blade sections, we now have a set of unorganized and potentially noisy points representing each of the selected radii r i . To be able to identify the leading and trailing edges, but also the thickness distribution and other essential parameters, we need to interpolate a boundary for each of the sections. With z =r.
coordinates given by (10)- (12), the sets of unwrapped points are separately located in planes at the given radial distances from the datum, which means that finding the sectional parameters for a given radius r i becomes a two-dimensional problem.
For each radius r i , we downsample the set of points j = [x j , y j ] T , j = 1, … , n by dividing them into an integer grid G. Then, all points that fall into the same cell D of the grid are replaced by one new point in D. The cells D are rectangles with two sides l 1 and l 2 , such that its height is l 1 and its width is l 2 . The output of the downsampling process is a new set k of m points, where m ≤ n.
Next, we sort the sampled points k for each of the sections according to the angle k which is formed between the positive x axis and the rays to each of the points in the plane as where ̄ = [p x ,p y ] T is the centroid given by Given a sampling grid of reasonable resolution, the downsampled and sorted points now form a boundary that represents the foil section at each of the selected radii. The boundary curve can be described parametrically as where s i ∈ {s 0 , … , s m } are called the breakpoints. The breakpoints of the boundary curve are not necessarily uniformly distributed, resulting in an unbalanced effect on the subsequent identification of the camber line. Therefore, a piecewise linear interpolation is used, which gives for s ∈ [s k , s k+1 ] and k = 0, 1, … , m − 1 . This gives a continuous interpolated boundary curve

Voronoi diagrams
All propeller foil sections may be regarded as a thickness distribution symmetrically superimposed upon a mean camber line. Consistent and accurate identification of the camber line is required to obtain the foil parameters robustly. A method for characterization of the line should be robust to noise and various realizations of the foil features.
(15) L(s i ) = (x(s i ), y(s i )), We propose an approach to this where Voronoi diagrams and Delaunay triangulation [22][23][24] are introduced in the characterization of ship propellers based on point cloud data.
First, the basics of Voronoi diagrams are presented. Consider a set S of n sites 1 , … , n ∈ ℝ 2 , which are points in the plane. Then, the set of points will form the perpendicular bisector of the sites p i and p j , which are the points of equal distance to the two points or sites. This bisector divides The line segments forming the boundaries of a Voronoi cell are referred to as Voronoi edges, whereas the endpoints of these edges are called Voronoi vertices. Two Voronoi cells with a common Voronoi edge are said to be Voronoi neighbors. An illustration showing a Voronoi diagram for a blade section is given in Fig. 10.

Mean camber line from Delaunay triangulation
Delaunay triangulation [24] is closely related to Voronoi diagrams. Given a Voronoi diagram Vor(S) , a Delaunay triangulation DT(S) is established by connecting the sites i of neighboring Voronoi cells with straight lines which are called Delaunay edges. These Delaunay edges form a set of faces that are called Delaunay triangles. The aggregate of these Delaunay triangles is the Delaunay triangulation DT(S) , which tessellates the area defined by the set of points S.
Let the Delaunay triangle T ijk be the triangle where the corners are the three sites i , j , k of the Voronoi diagram. Then, the circumcircle of T ijk , which is the circle that passes through the three corners of T ijk , will not contain any other site of Vor(S) . Moreover, the center ijk of the circumcircle of T ijk is the Voronoi vertex within the triangle.
The usefulness of Delaunay triangulation in connection with the identification of foil section parameters is that the mean camber line can be found efficiently and accurately from the measured points L(s i ) on the boundary curve. It is then used that the center ijk of the circumcircle of the triangle T ijk ∈ DT( ) is at the Voronoi vertex of the sites constituting the corners of the triangle. The Voronoi vertex will be at equal distance from the three corners of T ijk . This means that the points ijk will have equal distance d to the pressure and suction side of the boundary curve. Note that this distance is measured along a line that is normal to the boundary curve. The thickness of the foil section can then be approximated by t d = 2 d .
The mean camber line is typically defined such that the distance c to the boundary curve is the same on the suction side and the pressure side, where c is measured along a line that is normal to the mean camber line, and where t = 2 c is the thickness of the foil section. The difference between the estimate of d and c will be d ≈ c cos , where is the angle of a tangent of the boundary curve relative to the mean camber line. This angle will be limited to a few degrees, and accordingly, the cos term will be sufficiently close to unity to allow for the approximation d ≈ c . This means that the Delaunay triangulation can be used to find the mean camber line as an interpolation of the centers ijk of the circumcircle of the Delaunay triangles T ijk . It is noted that each Delaunay triangle T ijk has three corners. In this case, two points will be on the suction side and one on the pressure side, or vice versa. This will not pose any problem as long as the points of the boundary curve are sufficiently close. Figure 11 shows the computed Delaunay triangle circumcircles and the circle centers as an approximation of the mean camber line of the example propeller blade section. The corresponding Voronoi diagram is shown in Fig. 10.

Identification of edges and chord line
Close to the leading and trailing edges, reliable identification of the camber line and section thickness is challenging. The approach based on Voronoi diagrams and Delaunay triangulation forms a fork in the mean camber line at the section edges if the radius of the Delaunay triangle circumcircles is less than the edge radius, or if the trailing edge has sharp corners. As a consequence, the triangulation quality is dependent on the shape of the foil section edges. To overcome this difficulty while robustly identify the intersection of the section boundary with the mean camber line, we perform a least-squares fit of a small straight line segment at each end of the foil section and with slope coefficients a LE and a TE , and intercept coefficients b LE and b TE . The coefficients are computed from subsets S LE , S TE ∈ S for each foil section. The subsets are selected by performing a k-nearest neighbor search at each end of the initial mean camber line, which is identified as described in the previous section. It is expected that the straight line segments sufficiently approximate the camber line at each end for small distances.
The points representing the leading edge p LE and trailing edge p TE can then be found from the intersections of the section boundary K(l a ) with the best-fit straight lines f LE (x) and f TE (x) , respectively.
The chord length l is found by computing the Euclidean distance between the leading and trailing edges LE and TE . The notable sectional midpoint p M originates as the bisection of the chord line, which is a line drawn between the leading and trailing edges. The relationships linking the chord length, chord line, sectional midpoint, and leading and trailing edges are shown in Fig. 9. The local pitch angle is found from the sectional midpoint M = [p Mx , p My ] T which is given by The local pitch angle is then computed from From the parameters identified by the arrangements described in this and the previous sections, it is possible to derive the camber of the observed blade section as the maximal Euclidean distance between the sectional chord and camber lines.

Experimental results and discussion
The proposed approach for identifying geometric design parameters is validated by experiments performed on synthetic and real measured data. We demonstrate our procedure by identifying the design parameters of two marine propeller blades: a propulsion blade made from synthetically generated data, and a high-skew thruster blade of which we have (24) = arctan2(p My , p Mx ).

Fig. 11
The propeller blade section from Fig. 10 with a subset of the Delaunay triangulation circumcircles, in addition to the mean camber line identified via the circumcenters of the Delaunay triangles performed 3D scans and CMM measurements. We primarily consider two aspects when analyzing the performance, which is the geometrical shape error and its sensitivity to noisy measurements.

MOERI tanker KVLCC2
To test our approach on a model where ground-truth design parameters exist, we chose to use a cloud of points sampled from a commonly studied and available propeller model. The The KVLCC2 is appended with a 4-blade fixed-pitch propeller with a diameter of 9.86 m named KP458. This propeller model has been used in the testing of our approach, with a model scale of 1:58. Hence, the diameter of the scaled propeller model is 170 cm. The design parameters of the scaled propeller model are listed in Table 1.
To adequately mimic point cloud data originating from 3D scanning, we randomly distribute a set of points over one   [25]. On the resulting unorganized point cloud, the approach for parameter identification is conducted as described in Sect. 3. Figure 12 shows the surface representation of a KP458-blade, along with a point cloud representation of the blade made in the sampling step. The figure also illustrates the first step of our procedure, where the 0.16R, 0.25R, 0.30R, 0.40R, 0.50R, 0.60R, 0.70R, 0.80R, 0.90R, and 0.95R sections are extracted from the point cloud.
The resulting identified blade parameters for each radial section are listed in Table 2. An exception is the 1.0R-section located at the very tip of the blade. For this section, it is not practicable to extract a sufficient number of points, making meaningful retrieval of the geometric parameters a highly involved task. For the sections that were extracted and evaluated, the absolute geometric errors between the designed and identified parameters are given in Table 3. In Fig. 13, the identified propeller geometry parameters are compared with the design parameters. Additionally, the derived 0.7R-section of the KP458-blade is shown.

Noise tolerance
Point cloud data acquired from 3D sensors contain measurement uncertainty. Hence, the parameter identification procedure needs to be robust in the face of noisy data. To test our approach for various levels of measurement noise, we study the set of points S N of points Ni = i + i , i = 1, … , n , where the additive noise i is used to mimic the uncertainly in real scanned depth data.
A wide range of depth sensors is commercially available at different levels of accuracy and costs. The additive noise depends on the particulars of each sensor. As it is complex to define precise noise models for all sensors, we select a standard Gaussian noise model where is the identity matrix, and the standard deviation of the noise is found from for a given signal-to-noise ratio ( SNR ) specified in dB. This model is commonly known as zero-mean additive white Gaussian noise (AWGN). Figure 14 illustrates how the tip of the KP458 propeller blade is influenced by the added noise at various levels: SNR = {40, 45, 50, 55, 60} . A higher SNR value implies less distortion of the point cloud. In Fig. 15, the design parameters of the propeller geometry are compared with the identified parameters for the various levels of additive noise. It is worth noting that additive noise with SNR of 45 dB and below would for the given amount of points and their extent correspond to unrealistic levels of noise when considering industrial 3D scanning systems. From the plots in Fig. 14, it can be seen that for the skew angle s and thickness distribution t 0 , the added noise has a

3D scanning of a marine thruster blade
A second case study is provided to illustrate the capacities of the proposed approach when given real data originating from 3D scanning. The propeller blade examined in this experiment was designed for a four-blade, controllable pitch propeller cast in NiAl-Bronze. The total diameter of the propeller is 3600 mm, and the blade is right-handed with a propeller skew angle of 29 • . Thus, it is categorized as a highly skewed propeller. It was manufactured following the tolerances in ISO 484-1:2015(E), for accuracy Class 1. The design parameters of the thruster blade are listed in Table 4.
The 3D scanning of the blade was conducted in a GOM ATOS Scanbox MV 700, where the scanning 3D camera uses the structured light principle. The system has a specified maximum point spacing of 0.213 mm within the measuring volume. During the acquisition, the edges and root of the blade were coated with micro-particle spray to improve the scanning performance. A photograph of the 3D scanning setup during execution is shown in Fig. 16a, with a rendering of the resulting scan of the propeller blade in Fig. 16b. A total of 4.9 million points were sampled on the blade surface.
Based on the acquired point cloud, we effectuated our suggested approach for parameter identification, as described in Sect. 3. In Table 5, the resulting identified geometric parameters are listed for each radial section. A prerequisite for using the thruster blade in an analysis of the proposed parameter identification method is that the blade has been made entirely to the exact tolerances. During manufacturing, several factors contribute to dimensional inaccuracies, resulting in product geometry that does not fully match the design parameters. Independent measurements of the blade geometry are hence required to evaluate the capabilities of the proposed method. In this regard, we conducted CMM measurements of the thruster blade surfaces for comparison. The blade pitch and thickness distribution were measured, and the results are shown in Table 6.
In Fig. 17, the identified geometry parameters of the thruster blade are compared with the design parameters and with CMM measurements where available. Additionally, the derived 0.7R-section of the blade is shown. For the sections that were extracted and evaluated, the absolute geometric errors between the designed and identified parameters and between identified parameters and CMM measurements are given in Table 7.
From the results obtained in this case study, it can be seen that the proposed approach is successful in identifying the pitch P, skew angle s , chord length l, and maximum thickness t 0 of the thruster blade sections. The findings are, however, not as robust for the identification of camber distribution f 0 . Although the deviations are within the specified manufacturing tolerances, it is worth noting that the camber parameter is highly dependent on accurate identification of the leading and trailing edges. Small deviations in locating the blade edges make a significant impact on the following determination of the maximum camber of the sections.
Another critical aspect concerning precise parameter identification for individual propeller and thruster blades is to ensure the correct orientation of the blade relative to its reference plane and axis, i.e., the datum of the blade. The thruster blade inspected in this case study was oriented according to its center and guiding holes. As for the leading and trailing edges, small deviations in the initial orientation of the blade make a significant impact on parameter identification. Precise orientation is of importance both for 3D scanning and for CMM measurements.

Summary and conclusions
In this work, we have presented a procedure for identifying the various geometric design parameters of marine propeller blades based on discrete point cloud data originating from 3D scanning. The proposed approach is based on extracting the radial sections from the scanned point cloud, and then extracting it to a planar surface where the foil parameters can be identified. A central contribution is the use of Voronoi diagrams and Delaunay triangulation for determining the mean camber line of each section robustly. Experimental results from various case studies have demonstrated the effectiveness of the presented approach in extracting and identifying the chord line and mean camber line of the radial propeller foil sections. Consequently, the various geometric design parameters such as the pitch, chord length, skew, thickness, and camber distributions are identified. The presented method is shown to be robust to various levels of measurement uncertainty, and assumed a reliable method for extracting and analyzing the defining foil sections of marine propeller blades.
Future work will focus on practical applications of the suggested approach. Inspection, reconstruction, and parameter extraction are common problems in propeller blade manufacturing, and the proposed method can find several specific applications. We believe that the proposed method has the potential for successful application in the final inspection of manufactured marine propeller blades, in-line measurements in robotic grinding and surface finishing processes, and wear assessment or reverse engineering of propellers in service.