Specular Reflection Effects Elimination in Terrestrial Laser Scanning Intensity Data Using Phong Model

The intensity value recorded by terrestrial laser scanning (TLS) systems is significantly influenced by the incidence angle. The incidence angle effect is an object property, which is mainly related to target scattering properties, surface structures, and even some instrumental effects. Most existing models focus on diffuse reflections of rough surfaces and ignore specular reflections, despite that both reflections simultaneously exist in all natural surfaces. Due to the coincidence of the emitter and receiver in TLS, specular reflections can be ignored at large incidence angles. On the contrary, at small incidence angles, TLS detectors can receive a portion of specular reflections. The received specular reflections can trigger highlight phenomenon (hot-spot effects) in the intensity data of the scanned targets, particularly those with a relatively smooth or highly-reflective surface. In this study, a new method that takes diffuse and specular reflections, as well as the instrumental effects into consideration, is proposed to eliminate the specular reflection effects in TLS intensity data. Diffuse reflections and instrumental effects are modeled by a polynomial based on Lambertian reference targets, whereas specular reflections are modeled by the Phong model. The proposed method is tested and validated on different targets scanned by the Faro Focus3D 120 terrestrial scanner. Results imply that the coefficient of variation of the intensity data from a homogeneous surface is reduced by approximately 38% when specular reflections are considered. Compared with existing methods, the proposed method exhibits good feasibility and high accuracy in eliminating the specular reflection effects for intensity image interpretation and 3D point cloud representation by intensity.


Introduction
Terrestrial laser scanning (TLS) is a powerful technique for the spatial data acquisition and 3D characterization of different targets.This active remote sensing technique allows the direct, illumination-independent measurement of objects in a rapid, contactless, and accurate manner and has been increasingly adopted in various fields in the past two decades.In addition to the conventional precise 3D coordinates, TLS systems simultaneously measure the power of the backscattered laser signal of each scanned point and record it as an intensity value.The intensity value is the momentary amplitude of the return signal, which can be derived from the analog electrical signal output of the photodetector or digitized waveform [1].Theoretically, surfaces of higher reflectance will reflect a greater portion of the incident laser radiation, thereby increasing the backscattered signal power and further the intensity [1].Therefore, the intensity is a source of information closely associated with the reflectance properties of the scanned surface [2][3][4].Intensity, which is insensitive to ambient light and shadowing [5], is initially used to improve point cloud separability.Apart from visualization purposes, intensity data can be used as a major or complementary data source in various studies, such as vegetation and forest investigation [6,7], road traffic marking identification [8,9], water content extraction [6,[10][11][12][13], metro tunnel inspection [10,14], and lithological differentiation [5,15,16].
However, many confounding variables distort the capability of the original intensity to directly retrieve the target characteristics, of which the instrumental mechanism, atmospheric conditions, target surface properties, and data acquisition geometry plays a significant and dominant role [17,18].During one campaign period, instrumental configurations are kept constant and atmospheric attenuation can be ignored.As such, the intensity data recorded by the same TLS system are predominantly influenced by the target reflectance, distance, and incidence angle [19].In order to utilize the intensity for further applications, all influences not related to the material composing the scanned surface should be corrected [5].A number of studies have been successfully conducted to derive a corrected intensity that is merely related to the scattering properties of the scanned target by correcting the effects of incidence angle and distance.The irregular TLS distance effect is strongly dominated by instrumental factors (e.g., aperture size, automatic gain control, amplifier for low-reflective surfaces, and a brightness reducer for near distances [1]) and differs significantly among different systems.The distance effect does not completely follow the inverse square range function from the radar range equation [18][19][20].Empirically, the distance effect was satisfactorily corrected by a polynomial in [21] or based on reference targets in [19].In this study, the polynomial method was adopted to eliminate the distance effect to obtain the distance-independent intensity value.
In contrast to the distance effect that mainly depends on instrumental properties, the incidence angle effect is related to target scattering properties and surface structure [22,23].Theoretically, the incidence angle effect should be individually corrected in terms of target surface properties.To simplify the correction of the incidence angle effect, most existing methods assume that scanned surfaces are diffuse reflectors (i.e., Lambertian).Lambert's cosine law is extensively used to correct the incidence angle effect, regardless of discrepancies in the surface characteristics of various targets [24][25][26].However, natural surfaces do not behave as perfect Lambertian reflectors and Lambert's cosine law leads to an over-correction of the intensity values at large incidence angles [5].A number of reflection models are, therefore, proposed to optimize the Lambertian reflection model.To improve the accuracy of the correction of incidence angle effect, a revised model called the extended Lambertian reflection model was proposed in [22] to modify the contribution of incidence angle by changing the constant exponential value 1 to a parameter which describes the beam divergence.Considering that non-Lambertian behavior occurs in some natural surfaces, a model, which is a linear combination of the Lommel-Seeliger law pertaining to intrinsically dark surfaces and the Lambertian law related to brighter targets, was applied to correct the incidence angle effect in [17,20].Based on the idea that the surface geometry can be modeled by a relief of multiple micro-facets, the Oren-Nayar model was adopted in [5] to correct the incidence angle effect of different rock outcrops.However, these models mainly take into account the diffuse reflections for matte and rough surfaces.
Theoretically, both diffuse and specular reflections exist in all natural surfaces [22,27].The type of reflections influence the direction and strength of the backscattered light [26].Specular reflections can be ignored for rough surfaces.However, for targets with a relatively smooth or glossy surfaces, highlight phenomena (i.e., hot-spot effects in remote sensing [28][29][30][31]) caused by the specular reflections may appear in the intensity data.Highlights are detrimental for 3D visualization of the point cloud and further applications (e.g., target recognition, classification, and feature extraction) based on intensity.In this study, a novel method that considers diffuse and specular reflections, as well as the instrumental effects, is proposed to eliminate the specular reflection effects (i.e., highlight phenomenon).The effects of diffuse reflections and instrumental properties on the incidence angle are the same as that of Lambertian targets and are modeled by a polynomial.Specular reflections are modeled by the Phong model.Compared with existing methods that ignore specular reflections and are only suitable for rough surfaces, the proposed method can be applied to all surfaces.This study specifically aims to: 1.
Explore the mechanism of the specular reflection effects in TLS intensity; 2.
Use the Phong model to eliminate the highlight phenomenon caused by specular reflections; and 3.
Propose a new method to estimate the parameters of the Phong model.
The rest of the study is organized as follows: The proposed method is illustrated in detail in Section 2. Section 3 outlines the experiments.Sections 4 and 5 present the results and discussion, respectively.The conclusions are provided in Section 6.

Specular Reflections and Highlight Phenomenon
The TLS sensor emits laser pulses, illuminates the surface objects, and records the returned laser pulse signals backscattered from the surface objects [32].TLS backscattered light is significantly influenced by the surface scattering characteristics.If the scanned target is an ideal Lambertian reflector, pure diffuse reflection occurs.The entire surface reflects laser shots by Lambert's cosine law and the incident ray is uniformly reflected to a hemisphere isotropically.The radiation observed from a Lambertian surface is directly proportional to the cosine of the angle between the observer's viewing direction and the surface normal.Given that light paths from the emitter to the detector coincide in TLS, the received radiation is proportional to the cosine of the incidence angle θ (0 ≤ θ ≤ π/2) [26]: where L d is the diffuse radiation received by the observer and L in is the incident radiation.When a ray of light strikes a plane mirror (i.e., an ideal specular reflector), the incident ray reflects off the mirror.According to the law of reflection, the incidence and reflection angles are equal and are symmetrically distributed on two sides of the surface normal.However, for a generally-smooth surface constituted by multiple micro-facets with different orientations, specular reflections are distributed in the surroundings of the direction of reflection angle.When the viewpoint is within the vicinity of the direction of the reflection angle, strong specular reflection light can be received by the observer.On the contrary, specular light will be weak, or even disappear, when the viewpoint deviates significantly from the reflection direction [33].
The spatial distribution of specular reflections is directional.Specular reflection intensity depends on the intensity of the incident light and the surface material; it is also related with the viewing direction.Generally, the power function is used to simulate the spatial distribution of the specular reflections for a generally smooth surface [33], as shown by Equation (2): where L s is the specular radiation received by the observer, φ (0 ≤ φ ≤ π/2) is the angle between the directions of specular reflection and observer, and n is the specular highlight coefficient which is related to the object surface glossiness [27].Due to the coincidence of the positions of the emitter and detector, the angle between the directions of specular reflection and observer is the sum of the incidence and reflection angles, i.e., φ = 2θ.Therefore, Equation (2) can be written as: Actually, diffuse and specular reflections exist simultaneously on all natural surfaces [22,27].There are no ideal specular or diffuse targets.Therefore, considering only diffuse reflections is not appropriate when both diffuse and specular reflections exist.The Phong model is a bidirectional reflectance distribution function (BRDF) model [34,35].It describes the way that a surface reflects light as a combination in both diffuse and specular forms and has found wide application in computer graphics and 3D model rendering.In the study of [27], the Phong model was used to attenuate the effects of strong reflections that typically occur over wet or water-dominated areas scanned by airborne laser scanning.In this study, we adopted the Phong model to eliminate the highlight phenomenon in the TLS intensity data of smooth targets.
According to the Phong model [33], the total reflection radiation L t is a combination of diffuse and specular reflections for natural targets [22,27], as shown by Equation (4): where k s (0 ≤ k s ≤ 1) is the specular reflectance parameter which differs significantly for different targets.If k s = 0, the Phong model is completely equivalent to the Lambertian reflection model.On the contrary, the Phong model is exactly the same as the ideal specular reflection model when k s = 0.For natural surfaces, 0 < k s < 1.
According to Equation (3), at incidence angles larger than 45 • , the laser receiver cannot receive specular reflections though it exists in all natural targets (Figure 1a).Only diffuse reflections reach the receiver.However, at incidence angles smaller than 45 • , the received radiations include both diffuse and specular reflections (Figure 1b).The smaller the incidence angle is, more specular reflections can be received by the detector.If specular reflections are dominant over diffuse reflections, then highlight phenomena may exist in the intensity data.Highlight phenomena are those regions with small incidence angles that are brighter in the intensity image than those with large incidence angles where merely diffuse reflections can be received.Strictly, specular reflections should be considered in the intensity correction for all natural targets.k s is usually small for rough targets.As shown in Figure 2a, the diffuse reflections are nearly the same as the total reflections even at small incidence angles.This explains why Lambertian model can be used to approximately correct the incidence angle effect for rough targets [25].However, specular reflections of smooth surfaces are much larger than diffuse reflections at incidence angles smaller than 45 • and must be considered to eliminate the highlight phenomenon (Figure 2b).
attenuate the effects of strong reflections that typically occur over wet or water-dominated areas scanned by airborne laser scanning.In this study, we adopted the Phong model to eliminate the highlight phenomenon in the TLS intensity data of smooth targets.
According to the Phong model [33], the total reflection radiation L is a combination of diffuse and specular reflections for natural targets [22,27], as shown by Equation ( 4): where k (0 ≤ k ≤ 1) is the specular reflectance parameter which differs significantly for different targets.If k = 0, the Phong model is completely equivalent to the Lambertian reflection model.On the contrary, the Phong model is exactly the same as the ideal specular reflection model when k = 0.For natural surfaces, 0 < k < 1.
According to Equation (3), at incidence angles larger than 45°, the laser receiver cannot receive specular reflections though it exists in all natural targets (Figure 1a).Only diffuse reflections reach the receiver.However, at incidence angles smaller than 45°, the received radiations include both diffuse and specular reflections (Figure 1b).The smaller the incidence angle is, more specular reflections can be received by the detector.If specular reflections are dominant over diffuse reflections, then highlight phenomena may exist in the intensity data.Highlight phenomena are those regions with small incidence angles that are brighter in the intensity image than those with large incidence angles where merely diffuse reflections can be received.Strictly, specular reflections should be considered in the intensity correction for all natural targets.k is usually small for rough targets.As shown in Figure 2a, the diffuse reflections are nearly the same as the total reflections even at small incidence angles.This explains why Lambertian model can be used to approximately correct the incidence angle effect for rough targets [25].However, specular reflections of smooth surfaces are much larger than diffuse reflections at incidence angles smaller than 45° and must be considered to eliminate the highlight phenomenon (Figure 2b).

Physical Background of TLS Intensity
The TLS intensity value is theoretically a function of target reflectance ρ, incidence angle θ, and distance R. According to the Lambert's cosine law and radar range equation, the intensity value I can be expressed as [26]:

Physical Background of TLS Intensity
The TLS intensity value is theoretically a function of target reflectance ρ, incidence angle θ, and distance R. According to the Lambert's cosine law and radar range equation, the intensity value I can be expressed as [26]: However, the scattering behaviors of most natural surfaces do not exhibit Lambertian attributes.Additionally, intensity measurement is generally intended to enhance range determination; the raw values of the intensity signal can be modified to obtain an optimal range value [18,36].The range-dependent inverse square model from the radar range equation is insufficient to estimate the accurate intensity at the entire range scale.Considering the fact that the incidence angle and distance effects are theoretically independent of each other and can be solved individually [24], the intensity is written as [19,21]: where f ρ , f cos θ , and f R are functions of reflectance, cosine of incidence angle, and distance, respectively.Therefore, the corrected intensity I that is merely related to the target reflectance can be expressed as: where θ and R are the reference incidence angle and distance, respectively.θ and R can be arbitrarily defined [19].By dividing Equations ( 6) and ( 7), we obtain [3]: Distance-corrected intensity I and incidence angle-corrected intensity I are expressed as [3]: where C = f R and C = f cos θ are two constants.Empirically, the distance and incidence angle effects were satisfactorily corrected by polynomials in [21] or by using reference targets in [19].
In this study, we aimed to eliminate the highlight phenomena that are merely related to the incidence angle effect.First of all, the original intensity is corrected by a polynomial, i.e., f R = ∑ β • R , where β and N are polynomial coefficients, to eliminate the influence of the distance effect and obtain the distance corrected intensity . I only depends on the reflectance and incidence angle, as shown by the right part of the first equation of Equation (9).

Proposed Method
Previous studies [19,21] have indicated that, even for pure Lambertian targets, Lambert's cosine law cannot be applied to correct the incidence angle effect because some commercial undisclosed However, the scattering behaviors of most natural surfaces do not exhibit Lambertian attributes.Additionally, intensity measurement is generally intended to enhance range determination; the raw values of the intensity signal can be modified to obtain an optimal range value [18,36].The range-dependent inverse square model from the radar range equation is insufficient to estimate the accurate intensity at the entire range scale.Considering the fact that the incidence angle and distance effects are theoretically independent of each other and can be solved individually [24], the intensity is written as [19,21]: where f 1 (ρ), f 2 (cos θ), and f 3 (R) are functions of reflectance, cosine of incidence angle, and distance, respectively.Therefore, the corrected intensity I s that is merely related to the target reflectance can be expressed as: where θ s and R s are the reference incidence angle and distance, respectively.θ s and R s can be arbitrarily defined [19].By dividing Equations ( 6) and ( 7), we obtain [3]: Distance-corrected intensity I d and incidence angle-corrected intensity I a are expressed as [3]: where C 1 = f 3 (R s ) and C 2 = f 2 (cos θ s ) are two constants.Empirically, the distance and incidence angle effects were satisfactorily corrected by polynomials in [21] or by using reference targets in [19].
In this study, we aimed to eliminate the highlight phenomena that are merely related to the incidence angle effect.First of all, the original intensity is corrected by a polynomial, i.e., f 3 (R) = ∑ N 3 i=0 β i •R i , where β i and N 3 are polynomial coefficients, to eliminate the influence of the distance effect and obtain the distance corrected intensity . I d only depends on the reflectance and incidence angle, as shown by the right part of the first equation of Equation (9).

Proposed Method
Previous studies [19,21] have indicated that, even for pure Lambertian targets, Lambert's cosine law cannot be applied to correct the incidence angle effect because some commercial undisclosed instrumental effects may be mixed with the incidence angle effect.For natural surfaces, the incidence angle effect is, therefore, a combination of diffuse reflections, instrumental effects, and specular reflections.According to the first equation of Equation ( 9), the distance-corrected intensity I d can be expressed as: (10) where is a parameter related to target reflectance; I diffuse , I instrument , and I specular are intensity values caused by diffuse reflections, instrumental effects, and specular reflections, respectively.The diffuse reflections of a natural target are the same as that of a Lambertian target.Instrumental effects do not depend on target characteristics and are exactly the same for all scanned targets.Therefore, the former two parts of the right side of Equation ( 10) can be modeled by using a reference Lambertian target.Unfortunately, the specific instrumental effects are confidential due to commercial purposes and may differ significantly among different scanners.Regardless of what the specific instrumental effects are, a polynomial can be used to approximate the total effects of diffuse reflections and instrumental properties based on the Weierstrass approximation theorem [36].It should be noted that specular reflections should only be considered when incidence angles are smaller than 45 • .Thus, Equation ( 10) can be further written as: The polynomial parameters α i , β i , N 2 , and N 3 are the same for all targets and can be estimated by some Lambertian targets scanned at various incidence angles and distances.The estimation method is introduced in detail in [36].On the contrary, specular reflections depend on object surface smoothness and material.Consequently, k s and n should be determined individually.
The parameters k s and n were iteratively optimized within a homogeneous regions in [22] or empirically determined according to the data itself in [27].Considering that Equation ( 11) is a piecewise function and the first equation of Equation ( 11) is nonlinear, we proposed to estimate the polynomial parameters and specular reflection coefficients separately.Datasets of a certain surface with incidence angles both smaller and larger than 45 • should be sampled.According to the second equation of Equation (11), K 0 = I d /[∑ N 2 i=0 α i (cos θ) i ] can be calculated by using the sub-datasets with incidence angles larger than 45 • where no specular reflections exist in the intensity data.Then, the intensity value at incidence angles smaller than 45 • that is resulted from diffuse reflections and instrumental effects can be determined by the estimated polynomial.By subtracting the distance corrected intensity I d with the value of the polynomial, I specular can be estimated, as shown by Equation ( 12): To estimate the parameters of the Phong model, natural logarithms can be conducted on the two sides of Equation (12).Thus, Equation ( 12) can be written as: where e is the Euler's number and K = K 0 •k s .Equation ( 13) can be further written as: By conducting indirect adjustment of Equation ( 14) using the sub-datasets with incidence angles smaller than 45 • , the parameter X T = log e (K), n T = B T •B −1 •B T •L can be estimated where B and L are the matrix forms of 1 log e [cos(2θ)] and log e (M), respectively.Finally, the distance-corrected intensity can be further corrected for the incidence angle effect by the polynomial method according to the third part of Equation ( 8), as shown by Equation (15): The distance is calculated with the original 3D geometric coordinates of the scanned point and the scanner center.The surface normal is estimated by the best-fitting plane to a neighborhood of points surrounding the point of interest, and then the cosine of the incidence angle can be calculated by the incident radiation vector and the surface normal vector [21,36].

Experiments
Experiments on a wooden door in a laboratory were designed to validate the proposed method.The door, with a size of approximately 2 m × 1.4 m, is painted dark brown and surrounded by white lime walls.The surface of the door is smooth, whereas that of the wall is relatively rough.Experiments were conducted under indoor conditions through the use of the Faro Focus 3D 120, (FARO Technologies, Inc., Orlando, FL, USA) which delivers geometric information and returns intensity values recorded in 11 bits [0, 2048].Faro Focus 3D 120 is a continuous-wave 905 nm terrestrial laser scanner.The maximum distance is 120 m and the field-of-view is 360 • × 305 • .The beam divergence and the circular beam diameter at the exit are 0.009 • and 3.8 mm, respectively.
During the scanning period, all lights in the room were turned off and all curtains were pulled down to create a dark environment.This can avoid unnecessary errors from ambient light to make the results convincing and rigorous, although TLS intensity is insensitive to ambient light.The scan quality and scan resolution were set to 4 and 1/4, respectively.The scanner was moved from left to right with a total of five scans.The original intensity data were extracted in a point cloud image created by the standard software Faro SCENE 4.8.In our previous studies [21,36], the polynomial parameters of f 2 (cos θ) and f 3 (R) for the Faro Focus 3D 120 were calculated by some Lambertian targets scanned at various incidence angles and distances (Table 1).In this study, θ s and R s were chosen as 0 • and 5 m, respectively.

Results
The original intensity images of the five scans created by Faro SCENE are shown in Figure 3. Obviously, Scans 2 and 3 present a significant bright region (highlights) on the door because the incident ray is perpendicular to the bright region and the surface of the door is relatively smooth.On the contrary, Scans 1, 4, and 5 do not have apparent highlight regions because the incident radiation is perpendicular to the white lime wall whose surface is relatively rough.Datasets of the door with incidence angles from 0 • to 76 • were sampled.Firstly, the original intensity data were corrected for the distance effect by the polynomial method according to Equation ( 9) and the parameters in Table 1.As shown in Figure 4, the mean distance corrected intensity values for 0.5 • bin sizes, were used to estimate the parameters of the proposed method.The sub-datasets with incidence angles larger than 45 • were used to estimate K 0 .K 0 was estimated as 484.86.Then, the theoretical intensity values of the polynomial at incidence angles smaller than 45 • can be calculated.As shown in Figure 4a, the distance corrected intensity values of the door at incidence angles smaller than 45 • (particularly smaller than 15 • ) deviate significantly from the curve of the polynomial.By subtracting the polynomial values (red line in Figure 4a) with the distance corrected intensities, I specular can be estimated.Finally, according to Equation ( 14), log e (K), n was estimated as [215.06, 16.55] and k s = 0.44.distance corrected intensity values of the door at incidence angles smaller than 45° (particularly smaller than 15°) deviate significantly from the curve of the polynomial.By subtracting the polynomial values (red line in Figure 4a) with the distance corrected intensities, I can be estimated.Finally, according to Equation ( 14), log K , n was estimated as [215.06,16.55] and k = 0.44.With the obtained parameters of the Phong model, the curve of the first equation of Equation ( 11) (i.e., the curve of the sum of the polynomial and Phong model) can be drawn, as shown in Figure 4a.Apparently, the curve of the sum of the polynomial and Phong model fits optically with the distance corrected intensities at incidence angles smaller than 45°, whereas the curve of the polynomial can fit the intensities at incidence angles larger than 45°.This result proves the feasibility of our proposed method.Similarly, datasets of the white lime wall with incidence angles from 0° to 68° were sampled.The mean distance corrected intensity values for 0.5° bin sizes are shown in Figure 4b.K was estimated as 556.12 for the white lime wall.However, the distance-corrected intensity values are distributed closely around the curve of the polynomial at the entire scale of the incidence angles.This result indicates that the specular reflections of the wall are subtle and can be ignored; it also explains why there were no highlight regions when the incident ray was perpendicular to the white wall in Scans 1, 4, and 5.  Door data in Scans 2 and 3 were manually sampled using the polygon selector tool from the planar view in Faro SCENE to test the proposed method.The correction results are shown in Figure 5.The polynomial and reference targets methods were also used to correct the intensities, as shown in Figure 5b,c,f,g.However, the polynomial and reference targets methods failed to eliminate the specular reflection effects both in Scans 2 and 3. On the contrary, the highlight phenomenon was effectively eliminated by the proposed method, as shown in Figure 5d,h.The coefficient of variation (CV), namely, the ratio of the values of standard deviation (σ) and mean (u ) [22,26], can indicate the degree of dispersion of the intensity values from a homogenous target.The variation rate of CV, i.e., With the obtained parameters of the Phong model, the curve of the first equation of Equation ( 11) (i.e., the curve of the sum of the polynomial and Phong model) can be drawn, as shown in Figure 4a.Apparently, the curve of the sum of the polynomial and Phong model fits optically with the distance corrected intensities at incidence angles smaller than 45 • , whereas the curve of the polynomial can fit the intensities at incidence angles larger than 45 • .This result proves the feasibility of our proposed method.Similarly, datasets of the white lime wall with incidence angles from 0 • to 68 • were sampled.The mean distance corrected intensity values for 0.5 • bin sizes are shown in Figure 4b.K 0 was estimated as 556.12 for the white lime wall.However, the distance-corrected intensity values are distributed closely around the curve of the polynomial at the entire scale of the incidence angles.This result indicates that the specular reflections of the wall are subtle and can be ignored; it also explains why there were no highlight regions when the incident ray was perpendicular to the white wall in Scans 1, 4, and 5.
Remote Sens. 2017, 9, 853 8 of 15 distance corrected intensity values of the door at incidence angles smaller than 45° (particularly smaller than 15°) deviate significantly from the curve of the polynomial.By subtracting the polynomial values (red line in Figure 4a) with the distance corrected intensities, I can be estimated.Finally, according to Equation ( 14), log K , n was estimated as [215.06,16.55] and k = 0.44.With the obtained parameters of the Phong model, the curve of the first equation of Equation (11) (i.e., the curve of the sum of the polynomial and Phong model) can be drawn, as shown in Figure 4a.Apparently, the curve of the sum of the polynomial and Phong model fits optically with the distance corrected intensities at incidence angles smaller than 45°, whereas the curve of the polynomial can fit the intensities at incidence angles larger than 45°.This result proves the feasibility of our proposed method.Similarly, datasets of the white lime wall with incidence angles from 0° to 68° were sampled.The mean distance corrected intensity values for 0.5° bin sizes are shown in Figure 4b.K was estimated as 556.12 for the white lime wall.However, the distance-corrected intensity values are distributed closely around the curve of the polynomial at the entire scale of the incidence angles.This result indicates that the specular reflections of the wall are subtle and can be ignored; it also explains why there were no highlight regions when the incident ray was perpendicular to the white wall in Scans 1, 4, and 5.  Door data in Scans 2 and 3 were manually sampled using the polygon selector tool from the planar view in Faro SCENE to test the proposed method.The correction results are shown in Figure 5.The polynomial and reference targets methods were also used to correct the intensities, as shown in Figure 5b,c,f,g.However, the polynomial and reference targets methods failed to eliminate the specular reflection effects both in Scans 2 and 3. On the contrary, the highlight phenomenon was effectively eliminated by the proposed method, as shown in Figure 5d,h.The coefficient of variation (CV), namely, the ratio of the values of standard deviation (σ) and mean (u ) [22,26], can indicate the Door data in Scans 2 and 3 were manually sampled using the polygon selector tool from the planar view in Faro SCENE to test the proposed method.The correction results are shown in Figure 5.
The polynomial and reference targets methods were also used to correct the intensities, as shown in Figure 5b,c,f,g.However, the polynomial and reference targets methods failed to eliminate the specular reflection effects both in Scans 2 and 3. On the contrary, the highlight phenomenon was effectively eliminated by the proposed method, as shown in Figure 5d,h.The coefficient of variation (CV), namely, the ratio of the values of standard deviation (σ) and mean (u) [22,26], can indicate the degree of dispersion of the intensity values from a homogenous target.The variation rate of CV, i.e., δ = (CV ori − CV cor )/CV ori , where CV ori and CV cor are the coefficient of variations of the original and corrected intensities, respectively, was used to quantitatively evaluate the correction results.The statics of the intensities are shown in Table 2.In Scan 2, CV is 5.60% for the original intensities.It significantly decreases to 3.80% (δ = 32.14%)for the intensities corrected by the proposed method.However, CV is 5.35% (δ = 4.46%) and 5.40% (δ = 3.57%) for the intensities corrected by the polynomial and reference targets methods, respectively, which changes slightly compared with that of the original intensities.This result indicates that the polynomial and reference targets methods are invalid for the regions with incidence angles smaller than 45 • , i.e., these two methods cannot eliminate the highlight phenomenon caused by specular reflections.Similar results can be obtained for Scan 3, as shown in Table 2.  2. In Scan 2, CV is 5.60% for the original intensities.It significantly decreases to 3.80% (δ = 32.14%)for the intensities corrected by the proposed method.However, CV is 5.35% (δ = 4.46%) and 5.40% (δ = 3.57%) for the intensities corrected by the polynomial and reference targets methods, respectively, which changes slightly compared with that of the original intensities.This result indicates that the polynomial and reference targets methods are invalid for the regions with incidence angles smaller than 45°, i.e., these two methods cannot eliminate the highlight phenomenon caused by specular reflections.Similar results can be obtained for Scan 3, as shown in Table 2.

Discussion
Smooth surfaces (e.g., water [27], leaves [7], metal, and rubber products) widely exist in TLS scenarios and, thus, the elimination of specular reflection effects is highly important.To better validate the proposed method, a plastic curtain, a building facade decorated with small square ceramic tiles, a plywood door, a marble wall, an iron bookcase, and a rubber decorative board  The point clouds in the red dotted rectangles in Figure 6 were sampled to test the proposed method.The results are shown in Figure 8.The polynomial and reference targets methods are still invalid for these six surfaces.The mean values of δ for the polynomial and reference targets methods for the experimental targets are approximately 3.78% and 3.04%, respectively.Conversely, the proposed method can effectively correct the intensities of the highlight regions on these six surfaces.The values of δ when different parameters were used to eliminate the highlight phenomenon are shown in Table 4.The values of δ for the proposed method with optimal parameters are 4.88%, 18.64%, 11.00%, 61.29%, 69.98%, and 62.50% for the six surfaces.However, as expected the parameters of one surface are not suitable for another surface, i.e., the values on the diagonal are the largest among each column in Table 4.These results indicate that the parameters of the proposed method for different targets should be individually estimated based on the surface characteristics.

Conclusions
This study presents a new method to eliminate the specular reflection effects in the intensity data of smooth surfaces.By considering the effects of diffuse reflections, instrumental properties, and specular reflections, the incidence angle effect can be accurately corrected.The coefficients of variation of the intensity data are reduced by approximately 37.61%, 3.78%, and 3.04% for the proposed, polynomial, and reference targets methods, respectively.Since glossy or highly-reflective surfaces widely exist in TLS scenarios, elimination of specular reflection effects is inevitable and the proposed method can be used for applications (e.g., cultural relics modeling and preservation [20,37,38], point cloud supervision classification, and 3D visualization) based on TLS intensity data.To avoid the influence of specular reflection effects in the intensity data of smooth surfaces, a feasible way is to scan the target at large incidence angles, i.e., perpendicular scanning, is not preferred.However, large incidence angles may lead to poor quality point clouds [39].Therefore, this topic should be studied further.
Compared with existing methods which do not consider specular reflections, the proposed method is suitable for all natural surfaces.For rough surfaces whose specular reflections are relatively subtle, the second equation of Equation (15) with uniform parameters can be approximately applied to correct the incidence angle effect without considering the differences in surface characteristics.On the contrary, the parameters of the proposed method for glossy targets should be individually estimated based on the surface characteristics.By sampling datasets with incidence angles both smaller and larger than 45 • , the parameters of a certain target can be easily determined by two steps as stated in Section 2.
The Faro phase-based scanner was used to test the proposed method in this study.Theoretically, the incidence angle effect on intensity of all targets is a combination of diffuse reflections, specular reflections, and instrumental effects.For different scanners from other manufacturers, or using time-of-flight ranging methods, instrumental effects on intensity may be different.Therefore, the polynomial parameters N 2 and α i may be different.Regardless of the types and manufacturers of the scanners, the proposed method can be applied to other terrestrial laser scanners.However, this should be further tested and validated.Empirically, the scanned data of Faro Focus 3D 120 above 30 m are unreliable in practical applications because of the low accuracy, large amount of noise, and high level of uncertainty, although the theoretical maximum distance is 120 m [3].The parameters of f 3 (R) in our previous studies [21,36] were estimated by some Lambertian targets scanned at various distances below 30 m.Therefore, the polynomial parameters of f 3 (R) in Table 1 are suitable for targets scanned at distances below 30 m.However, it should be noted that the proposed method can be theoretically applied to arbitrary distances if the parameters of f 3 (R) are estimated by Lambertian targets scanned within the entire scanning range scale.

Figure 1 .
Figure 1.(a) At incidence angles larger than 45°, only diffuse reflections reach the receiver.(b) At incidence angles smaller than 45°, both diffuse and specular reflections can be received.The red dotted lines are perpendicular to the reflection directions.

Figure 1 .
Figure 1.(a) At incidence angles larger than 45 • , only diffuse reflections reach the receiver.(b) At incidence angles smaller than 45 • , both diffuse and specular reflections can be received.The red dotted lines are perpendicular to the reflection directions.

Figure 2 .
Figure 2. Phong model with different parameters.(a) Specular reflections are dominant at incidence angles smaller than 45° and must be considered (smooth surfaces).(b) Specular reflections are subtle (rough surfaces).

Figure 2 .
Figure 2. Phong model with different parameters.(a) Specular reflections are dominant at incidence angles smaller than 45 • and must be considered (smooth surfaces).(b) Specular reflections are subtle (rough surfaces).

Figure 3 .
Figure 3. Original intensity images created by Faro SCENE.(a) Scan 1.(b) Scan 2. (c) Scan 3. (d) Scan 4. (e) Scan 5. Highlights exist in Scans 2 and 3 because the surface of the door is smooth.Scans 1, 4, and 5 do not have highlight regions as the surface of the wall is relatively rough.

Figure 4 .
Figure 4. (a) Distance-corrected intensity values of the door and the curves of the fitting polynomial and Phong model.(b) Distance-corrected intensity values of the white lime wall and the curve of the fitting polynomial.

Figure 3 .
Figure 3. Original intensity images created by Faro SCENE.(a) Scan 1.(b) Scan 2. (c) Scan 3. (d) Scan 4. (e) Scan 5. Highlights exist in Scans 2 and 3 because the surface of the door is smooth.Scans 1, 4, and 5 do not have highlight regions as the surface of the wall is relatively rough.

Figure 3 .
Figure 3. Original intensity images created by Faro SCENE.(a) Scan 1.(b) Scan 2. (c) Scan 3. (d) Scan 4. (e) Scan 5. Highlights exist in Scans 2 and 3 because the surface of the door is smooth.Scans 1, 4, and 5 do not have highlight regions as the surface of the wall is relatively rough.

Figure 4 .
Figure 4. (a) Distance-corrected intensity values of the door and the curves of the fitting polynomial and Phong model.(b) Distance-corrected intensity values of the white lime wall and the curve of the fitting polynomial.

Figure 4 .
Figure 4. (a) Distance-corrected intensity values of the door and the curves of the fitting polynomial and Phong model.(b) Distance-corrected intensity values of the white lime wall and the curve of the fitting polynomial.
CV and CV are the coefficient of variations of the original and corrected intensities, respectively, was used to quantitatively evaluate the correction results.The statics of the intensities are shown in Table

Figure 5 .
Figure 5. Point cloud of the door colored by intensity.(a) Original intensities in Scan 2. (b) Intensities corrected by the polynomial method in Scan 2. (c) Intensities corrected by the reference targets method in Scan 2. (d) Intensities corrected by the proposed method in Scan 2. (e) Original intensity in Scan 3. (f) Intensities corrected by the polynomial method in Scan 3. (g) Intensities corrected by the reference targets method in Scan 3. (h) Intensities corrected by the proposed method in Scan 3.

Figure 5 .
Figure 5. Point cloud of the door colored by intensity.(a) Original intensities in Scan 2. (b) Intensities corrected by the polynomial method in Scan 2. (c) Intensities corrected by the reference targets method in Scan 2. (d) Intensities corrected by the proposed method in Scan 2. (e) Original intensity in Scan 3. (f) Intensities corrected by the polynomial method in Scan 3. (g) Intensities corrected by the reference targets method in Scan 3. (h) Intensities corrected by the proposed method in Scan 3.

Figure 8 .
Figure 8. Point cloud colored by intensity.(I) Original intensities.(II) Intensities corrected by the polynomial method.(III) Intensities corrected by the reference targets method.(IV,V) Intensities corrected by the proposed method with different parameters.(a) Curtain, (IV,V): parameters of the curtain and door.(b) Building facade, (IV,V): parameters of the building facade and curtain.(c) Plywood, (IV,V): parameters of the plywood and rubber.(d) Marble, (IV,V): parameters of the marble and bookcase.(e) Bookcase, (IV,V): parameters of the bookcase and rubber.(f) Rubber, (IV,V): parameters of the rubber and bookcase.

Figure 8 .
Figure 8. Point cloud colored by intensity.(I) Original intensities.(II) Intensities corrected by the polynomial method.(III) Intensities corrected by the reference targets method.(IV,V) Intensities corrected by the proposed method with different parameters.(a) Curtain, (IV,V): parameters of the curtain and door.(b) Building facade, (IV,V): parameters of the building facade and curtain.(c) Plywood, (IV,V): parameters of the plywood and rubber.(d) Marble, (IV,V): parameters of the marble and bookcase.(e) Bookcase, (IV,V): parameters of the bookcase and rubber.(f) Rubber, (IV,V): parameters of the rubber and bookcase.

Table 2 .
Statics of intensity values of the door.

Table 2 .
Statics of intensity values of the door.

Table 4 .
Values of δ (%) when parameters of one surface (row) were used to eliminate highlight phenomena of another surface (column).

Table 4 .
Values of δ (%) when parameters of one surface (row) were used to eliminate highlight phenomena of another surface (column).