Results of a Surface Roughness Comparison between Stylus Instruments and Confocal Microscopes

This article presents the results of an LMM-R-2019 interlaboratory comparison. Such comparisons of different families of measuring instruments are one of the activities conducted among the calibration laboratories to maintain their ISO 17025 accreditation. Given that the study of surface roughness is becoming increasingly important in the field of dimensional metrology, the comparison focused on determining the Ra parameter on a pseudorandom metallic roughness standard using two types of measuring instruments: physical contact (stylus instruments) and optical (confocal microscopes). Among the aspects studied was whether the roughness measurements obtained using calibrated confocal microscopes could be compared with those using traditional methods since optical instruments obtain measurements more quickly and responsively than do stylus instruments. The results showed that roughness measurements using confocal microscopes are comparable with those from a traditional stylus instrument.


Accreditation
All products are launched commercially on the market after exhaustive quality control testing, which includes dimensional control. Measurements in general, and dimensional measurements in particular, play a very important role in everyday life. Being able to carry out correct measurements is very important from a financial point of view since it provides interchangeability for products and services on a global level. For a company to have sufficient credibility and confidence in its products or services, it needs an independent evaluation system, also known as accreditation. This internationally established tool helps to instil confidence that a company has performed a test or calibration correctly, among other conformity assessment activities. Examples of advantages that accreditation provides include • Reliability and recognition for a company's products and services, • Lower possibility of providing a service or making a product with defects, • Savings in the time and cost of re-evaluating products, • Opening business opportunities in international markets by reducing costs for manufacturers and exporters.
At the national level, each country has an entity that carries out these conformity assessment activities independently and rigorously. This may be a National Metrology Institute in accordance with the document, Quality Management Systems in the CIPM MRA-G-12 [1]. In Spain, the organization responsible for carrying out these tasks is the National

•
Evaluate random variation in the results of measurements conducted by different participating laboratories, • Determine systematic differences in laboratory measurement results, and • Determine the value of one or more physical properties for one or more artefacts.
The ISO/IEC 17025 standard also establishes that non-standard methods, methods developed in a specific laboratory and standard methods used outside this laboratory that are not within the scope of accreditation can be validated [3].

Evaluation of Surface Roughness
As commented by Leach et al. in [6], the measurement of the geometric properties of parts is a crucial factor in determining the functional properties of each component, as well as controlling the manufacturing parameters. As Pawlus et al. established in [7,8], the surface topography is the result of manufacturing, functioning conditions, and wear. One of the properties that is becoming increasingly important in the manufacturing industry is the evaluation and analysis of surface texture according to the study carried out by Todhunter et al. in [9]. The surface of a part separates it from what is not the part whether it be the environment or another object, as established by Leach et al. [6,10]. As Nieslony et al. established in [11], it is crucial to distinguish between a mechanical and an electromagnetic surface as defined in the standard ISO 25178-2 [12]: • Mechanical surface (Section 3.1.1.1 of [12]): "the boundary of the erosion, by spherical ball of radius r, of the locus of the centre of an ideal tactile sphere, also with radius r, rolled over the skin model of a workpiece". • Electromagnetic surface (Section 3.1.1.2 of [12]): "the surface obtained by the electromagnetic interaction with the skin model of a workpiece".
According to the ISO 25178-6 standard [13] and as stated by Leach in [14], there are three types of methods to analyse surface texture: • Line-profiling: These provide a 2-dimensional plot or profile of the surface, expressed mathematically as z(x). • Areal surface: These provide a 3-dimensional plot of the surface, expressed mathematically as z(x, y). Some measurement instruments that work with line profiling methods can work with areal surface methods, by extracting parallel profiles along the y-axis. • Area-integrating: These methods provide a numerical result dependent on the areaintegrating properties of the surface texture, by analysing the distribution of the scattered light. Figure 1 of the ISO 25178-6 standard [13] specifies the different measuring instruments used to analyse a surface by each of the methods described above. Caja et al. in [15] and Townsend et al. in [16] established another classification according to the surface analysis method whether by mechanical contact or optical methods. Pawlus et al. introduced in [7], the most used 2D roughness parameter level is the vertical parameter Ra (see Figure 1), which is calculated as As can be seen, for each evaluation length l p the part of the profile be becomes positive and the arithmetic mean of the profile | ( )| is calculate

•
Contact stylus instruments: These mechanically contact the surface being scanned using a tip, the movement of which is converted into a signal as the function z(x).
The z-coordinate of the contact point is determined by the transverse trajectory of the instrument axes, the shape of the tip, the deviation of the probe and the local inclination of the surface [6,13,14]. These instruments provide measurements of the mechanical surface of a sample.

•
Optical instruments: These emit electromagnetic radiation reflected from the surface and detected by a sensor. The intensity of the detected radiation is used to calculate the spatial coordinates [6,14]. There are two types: those that measure the real surface of the study piece (using line-profiling or areal surface methods) and those that analyse scattered light (area-integrating methods). These instruments provide measurements of the electromagnetic surface of a sample. Among the different types of optical instruments are point autofocus, focus variation, interferometric (phase shifting and coherence), holographic digital microscopy, and confocal microscopes.
Mathia et al. also highlighted in [17] the importance of scanning tunnelling microscopes (STMs) and atomic force microscopes (AFMs) in the analysis of surface roughness. These instruments are usually used as scanning probe microscopes. There are other hybrid measuring instruments that allow the measurement of a surface with both optical and contact techniques. An example is the Infinite Focus Measurement (IFM) machine used by Nieslony et al. in [11]. Surface texture is quantified by a series of parameters established by mathematical analysis. As Królczyk et al. established in [18], the emerging aim in science is to determine functional surface parameters of the samples. Determination of the 2D roughness parameters is included in the ISO 4287 [19], ISO 12085 [20] and ISO 13565-2 [21] specification standards, while the 3D roughness parameters are included in the ISO 25178-2 standard [12]. Among the different standards, it is possible to find almost 300 parameters to characterize the surface, according to Mathia et al. in [17]. They could be classified into height, spatial, hybrid, functional, end and other, as Pawlus et al. commented in [8]. In this article, we focus on 2D roughness parameters.
According to the results of the survey conducted by Todhunter et al. in [9] and as Pawlus et al. introduced in [7], the most used 2D roughness parameter at the industrial level is the vertical parameter Ra (see Figure 1), which is calculated as As can be seen, for each evaluation length l p the part of the profile below the midline becomes positive and the arithmetic mean of the profile |z(x)| is calculated for that length. The advantage of this parameter is that it is affected less than other parameters (such as Rz, Rp, Rv, Rc and Rt) by scratches that the standard may have on its surface, foreign particles present in the groove, or by noise introduced during measurement.
Another widely used parameter in industrial environments is Rz (see Figure 2). It is calculated as the sum of height of the maximum profile peak height Z p and the maximum profile valley depth Z v [19], as detailed in Equation (2): Materials 2022, 15,5495 The advantage of this parameter is that it is affected less than other paramete Rz, Rp, Rv, Rc and Rt) by scratches that the standard may have on its surfa particles present in the groove, or by noise introduced during measurement. Another widely used parameter in industrial environments is Rz (see Fig  calculated as the sum of height of the maximum profile peak height and the profile valley depth [19], as detailed in Equation (2) Note that this parameter is very sensitive to scratches, dirt, and other surfa As Caja et al. pointed out in [15], there are substantial differences between the a surface by physical contact with stylus instruments and by optical measur ments, due to the difference in the physical measurement principles of each. I surface is measured using different types of measuring instruments, the result parable, however different they may seem at first. According to Leach et al. in [ measuring instruments present some potential advantages such as the absence of damage to the surface of the sample and a shorter time spent performing th ment. However, optical measuring instruments historically have not always pro ues equivalent to those obtained by contact stylus instruments. Leach et al. als in [22] that the problems presented by optical measuring instruments make ne a priori understanding of the nature of the surface. ILCs are very useful for the tasks. In this study, we show the ILC LMM-R-2019 surface roughness results in authors in 2019 participated with other nationally accredited laboratories. Am things, the intention was to compare the measurements of the parameter Ra optical measuring instruments (confocal microscopes) with those obtained u instruments from several ENAC accredited laboratories.

Interlaboratory Comparison LMM-R-2019
The following objectives were pursued in the ILC, LMM-R-2019: • Study the performance of the participating laboratories (accredited if p determining the Ra and Rz parameters for a roughness standard. These ar used and studied roughness parameters according to the study of Todhun [9]. • Determine the degree of compatibility among the different laboratory resu 2.4 details the procedure.

•
Study the behaviour of optical measuring instruments with respect to con ods. According to Leach et al. in [22], there are no widely accepted metho brating optical measuring instrument. For this reason, this ILC could be u Note that this parameter is very sensitive to scratches, dirt, and other surface defects. As Caja et al. pointed out in [15], there are substantial differences between the analysis of a surface by physical contact with stylus instruments and by optical measuring instruments, due to the difference in the physical measurement principles of each. If the same surface is measured using different types of measuring instruments, the results are comparable, however different they may seem at first. According to Leach et al. in [22], optical measuring instruments present some potential advantages such as the absence of the risk of damage to the surface of the sample and a shorter time spent performing the measurement. However, optical measuring instruments historically have not always provided values equivalent to those obtained by contact stylus instruments. Leach et al. also point out in [22] that the problems presented by optical measuring instruments make necessary an a priori understanding of the nature of the surface. ILCs are very useful for these types of tasks. In this study, we show the ILC LMM-R-2019 surface roughness results in which the authors in 2019 participated with other nationally accredited laboratories. Among other things, the intention was to compare the measurements of the parameter Ra using two optical measuring instruments (confocal microscopes) with those obtained using stylus instruments from several ENAC accredited laboratories.

Interlaboratory Comparison LMM-R-2019
The following objectives were pursued in the ILC, LMM-R-2019: • Study the performance of the participating laboratories (accredited if possible) by determining the Ra and Rz parameters for a roughness standard. These are the most used and studied roughness parameters according to the study of Todhunter et al. in [9]. • Determine the degree of compatibility among the different laboratory results. Section 2.4 details the procedure.

•
Study the behaviour of optical measuring instruments with respect to contact methods. According to Leach et al. in [22], there are no widely accepted methods for calibrating optical measuring instrument. For this reason, this ILC could be used to validate the calibration procedures of the different confocal microscopes.

Participants
Two types of laboratories participated.
• Accredited Calibration Laboratories (ACLs): A total of 7 ACLs participated from the public and industrial sectors and university centre laboratories. Only six laboratories in Spain are accredited by ENAC (the Spanish Accreditation Body) under ISO 17025 [3] to perform calibrations in roughness, and all of them participated in this ILC. All six are accredited to perform calibrations of roughness standards by determining Ra and Rz roughness parameters (amplitude parameters [19]). For this reason, the ILC was limited only to those two amplitude parameters. Only two of these six ACLs are accredited for the determination of Sm (spacing parameter [19]). The seventh ACL that participated in the ILC is accredited by ENAC in dimensional metrology but not roughness. Among them was the organiser of the ILC, the Metrology and Metrotechnics Laboratory (Laboratorio de Metrología y Metrotecnia, LMM) at Madrid Polytechnic University (Universidad Politécnica de Madrid, UPM).

•
Research centres: The two participating research centres were the Dimensional Metrology Laboratory (LMD) and the Laser Centre (CL), both from UPM. It should be noted that the CL and LMM are inter-related, as they share staff. However, unlike the LMD, they do not share facilities.
All laboratories were traceable to the Spanish Metrology Centre (CEM). It was also remarkable that all of them performed roughness measurements in industrial environments. For this reason, the calibration and measuring procedures should be as simple and robust as possible.

Measuring Instruments
Two types of measuring instruments were used.

•
Stylus. All laboratories except CL used a stylus. These are classic measuring instruments on which all traceable and precision roughness measurements have been based so far. Their main advantage is that they are accepted worldwide for roughness measurement, and there is a lot of know-how based on its use. There are many wellestablished manufacturing processes that have quality assurance based on roughness measurements with these instruments. However, they have several limitations. As Leach et al. describe in [22], one of the most critical parameters is the shape of the stylus. Typically, it has a 2 µm radius tip, which may damage the surface of the samples and provide a too-coarse filter for certain applications. By taking images at different heights, a 3D reconstruction of the sample surface can be created using specialised software. These are optical measuring instruments the use of which is expanding as they make possible the carrying out of non-contact measurements of 3D roughness parameters. However, as stated by Leach et al. in [22], large differences have been observed in the results compared with those provided by traditional stylus instruments. In addition, all the acquired know-how and the procedures and techniques developed for stylus instruments may not work as well with these measuring instruments.
The instruments used in the two research centres were LEICA DCM3D confocal microscopes. The CL one was calibrated according to the procedure proposed by Mínguez-Martínez et al. in [23], while the procedure proposed by Wang et al. in [24] was used for the one in LMD. It is also necessary to say that the CL confocal microscope was not able to provide the Rz roughness parameter, so this measuring instrument only participated in the determination of Ra roughness parameter. Measurements were made of the object in Figure 3. This is a metallic, pseudorandom roughness standard, HALLE brand, model the roughness parameter Ra was to be determined on its standard face. The nom ues of the standard are = 0.8 μm and = 4.5 μm. The measuring face is as i in the figure.

ILC Description
The organising laboratory, LMM, made an initial measurement and sent the s to the second laboratory. From then on, each laboratory sent the standard to the n according to a previously established schedule. After receiving the standards, e ticipating laboratory had to send an email to the LMM with a report on the stat standard, especially of its measurement surface. Finally, the standard reached t which performed a second measurement.
Each measurement was made in the area indicated according to the followi ciples: • The scan had to be located as close as possible to the measurement face axis line collinear with the arrow engraved on the face).

•
The evaluation length could not exceed the area enclosed by the dark m graved on the right and left of the measurement face.

•
The evaluation length had to be centred with respect to the measurement fa In accordance with the reference standard ISO 4288 [19], the following meas conditions were established.
As confocal microscopes cannot generally reach these and values, the were required to be as close as possible. In addition, if the temperature could no trolled, the LMM had to be informed to take this into account.
After the laboratories had performed the measurements, the results were r on a calibration certificate. This included the measurement result (the paramete This is a metallic, pseudorandom roughness standard, HALLE brand, model A1, and the roughness parameter Ra was to be determined on its standard face. The nominal values of the standard are Ra = 0.8 µm and Rz = 4.5 µm. The measuring face is as indicated in the figure.

ILC Description
The organising laboratory, LMM, made an initial measurement and sent the standard to the second laboratory. From then on, each laboratory sent the standard to the next one, according to a previously established schedule. After receiving the standards, each participating laboratory had to send an email to the LMM with a report on the status of the standard, especially of its measurement surface. Finally, the standard reached the LMM which performed a second measurement.
Each measurement was made in the area indicated according to the following principles: • The scan had to be located as close as possible to the measurement face axis (straight line collinear with the arrow engraved on the face).

•
The evaluation length could not exceed the area enclosed by the dark marks engraved on the right and left of the measurement face.

•
The evaluation length had to be centred with respect to the measurement face.
In accordance with the reference standard ISO 4288 [19], the following measurement conditions were established.
As confocal microscopes cannot generally reach these l p and l r values, the lengths were required to be as close as possible. In addition, if the temperature could not be controlled, the LMM had to be informed to take this into account.
After the laboratories had performed the measurements, the results were recorded on a calibration certificate. This included the measurement result (the parameter Ra expressed in µm, to the appropriate number of significant figures), the expanded uncertainty (with coverage factor k = 2, which corresponds to a coverage probability of~95%), the values of l p and l r used and the measurement temperature. This was all done according to International Metrology Vocabulary (VIM) terminology [25]. All measurement uncertainties were calculated according to the ISO-GUM guide [26] or document EA-04/02 M:2013 [27], the standard procedures used by ACLs.
When all measurements were done, each participant emailed the others a protected zip file with the calibration certificate. When all participants confirmed receipt of the email to the MM, it sent another email to all the participants who sent the code to be able to decompress the zip file. This meant all participants could simultaneously access the calibration results.
The LMM prepared a final report for the ILC, LMM-R-2019, summarising the results from all the participants with different statistical parameters, which made it possible to compare the performance of each laboratory.

Estimation of the Standard Value and Statistical Parameters
This section establishes the procedure for estimating the measured standard value, which was assumed to be stable over time. The one used during the ILC can be described as a real number for which the estimate x 0 (with measurement uncertainty u(x 0 )) had to be evaluated. In general, if the number of laboratories is n, the problem can be described by means of the vector equation described by Cox in [28]: As all the laboratories were traceable to the CEM, the correlation among them was not zero, so the covariance matrix of the independent term → x had to be calculated. Equation (4) describes how to calculate this matrix: It should be noted that • Each term u(x i ) corresponds to the uncertainty obtained by each of the laboratories. • Each term r ij corresponds to the correlation between laboratories i and j.
As all ACLs are traceable to the CEM, the correlation among independent laboratories was considered to be r ij = 0.1.

•
As two sets of measurements were performed in the LMM (at the beginning and end of the ILC) using the same equipment under similar conditions and with the same technician, the correlation between them must be the highest of all. Therefore, the value r = 0.6 was assigned.

•
As two measurements were made in the LMD (one with a stylus instrument and the other with a confocal microscope), the correlation between these two measurements must be slightly higher than those with other laboratories but not as high as the two measurements made by the LMM. Therefore, the value r = 0.5 was assigned. • For measurements made in the CL, as it was related to the LMM via its technicians, the correlation among them and those made in the LMM must be higher than between the CL and other laboratories but less than those of the LMD. Therefore, the value r = 0.4 was assigned.
It can be shown that the solution to the matrix problem can be reduced to the following expressions: Once the measured value of the standard was estimated, the statistical parameters were calculated to evaluate the performance of each laboratory, which is included in the reference standard ISO 17043 [5]. Please note that due to the correlation between x i and reference value x 0 in Equation (9) there is a negative sign in the denominator instead of the usual positive sign when x i and x 0 are statistically independent.
Standard uncertainty of differences : Normalized error or Compatibility Index : 3. Results Table 1 shows the Ra values (expressed as x i ) and the associated expanded uncertainty (with coverage factor k = 2, that provides a coverage probability of 95%) obtained by each of the participating laboratories. All units are in µm. As can be seen, the results are very close to the nominal values of the material standard (Ra = 0.8 µm and Rz = 4.5 µm). Of all the results, we considered that the measurements made with the confocal microscopes and the estimation of the uncertainty associated with these measurements should be highlighted.
Attending to the Ra roughness parameter results, confocal microscopes did not seem to present significant deviations from the laboratories that used stylus instruments, considering the uncertainty. This meant that the deviations were lower than the uncertainty. For example, the maximum deviation was between the confocal microscope of CL and Lab 2, Lab 4 and Lab 5 (0.09 µm), which was lower than the uncertainty of the CL and very similar to that of Lab2, Lab 4 and Lab 5. It is also important to highlight that the confocal microscope's uncertainties (0.12 and 0.16 µm) are clearly higher than the uncertainties of the stylus instruments (by approximately 0.08 µm on average), that is, between 50 and 100% higher.
Attending to Rz roughness parameter results, the maximum variability observed was 0.44 µm. It was important to consider that this parameter was more statistically unstable, which is the reason that greater variability is observed in relative terms. The maximum difference between stylus instrument and confocal microscope measurements is 0.34 µm, which is clearly lower than the uncertainty of the confocal microscope (0.48 µm) and is on the order of most uncertainties obtained with the stylus instruments. For this reason, it could be said that the uncertainty for confocal microscopes is similar to those for the stylus instruments.
The measurements for the CL confocal microscope were made with the 10× objective, the same with which the confocal microscope was calibrated. This objective provided a distance between planes (voxel height) of 2 µm. Measurements were taken so that the roughness parameter Ra could be measured both along the instrument xand y-axes (see Figure 4).  The confocal microscope provides 3D images like the one in Figure 5.

Ra Roughness Parameter
The measurement uncertainty of the equipment used by the CL was estimated from the formula described by Mínguez-Martínez et al. in [23]: where ̅ is the mean Ra measurement; ( ) is the uncertainty associated with the z-axis scale correction factor ( ); ( ̅ ) is the repeatability; is the number of measurements made; ( ̅ ) is the uncertainty associated with the positive bias of the measurement instrument; ( ) is the uncertainty associated with the variation in the positive bias of the measurement instrument; and is the uncertainty component due to instrument noise. In this case, the value = 0 was assigned as it is small enough to be considered negligible compared to the height of the voxel.
Similarly, the estimation of uncertainties of the equipment used by the LMD-Conf considers the following sources of uncertainty [24]: Measurement noise (evaluated when measuring at standard flatness), The confocal microscope provides 3D images like the one in Figure 5.  The confocal microscope provides 3D images like the one in Figure 5.

Ra Roughness Parameter
The measurement uncertainty of the equipment used by the CL was estimated from the formula described by Mínguez-Martínez et al. in [23]: where ̅ is the mean Ra measurement; ( ) is the uncertainty associated with the z-axis scale correction factor ( ); ( ̅ ) is the repeatability; is the number of measurements made; ( ̅ ) is the uncertainty associated with the positive bias of the measurement instrument; ( ) is the uncertainty associated with the variation in the positive bias of the measurement instrument; and is the uncertainty component due to instrument noise. In this case, the value = 0 was assigned as it is small enough to be considered negligible compared to the height of the voxel.
Similarly, the estimation of uncertainties of the equipment used by the LMD-Conf considers the following sources of uncertainty [24]: Measurement noise (evaluated when measuring at standard flatness),

Ra Roughness Parameter
The measurement uncertainty of the equipment used by the CL was estimated from the formula described by Mínguez-Martínez et al. in [23]: where R is the mean Ra measurement; u(c z ) is the uncertainty associated with the z-axis scale correction factor (c z ); s R is the repeatability; m is the number of measurements made; u b is the uncertainty associated with the positive bias of the measurement instrument; u(δb) is the uncertainty associated with the variation in the positive bias of the measurement instrument; and u noise is the uncertainty component due to instrument noise. In this case, the value u noise = 0 was assigned as it is small enough to be considered negligible compared to the height of the voxel. Similarly, the estimation of uncertainties of the equipment used by the LMD-Conf considers the following sources of uncertainty [24]: • Measurement repeatability, • Z-axis calibration coefficient, • Measurement noise (evaluated when measuring at standard flatness), • Variation in light intensity used by the equipment, and • Variation in the sample orientation (lack of perpendicularity of the sample with respect to the equipment optical axis).
Considering the results in Table 1, the correlation coefficients r ij and Equation (4), the covariance matrix is calculated as: In this case, since all r ij = 0, the matrix has all terms. The values are generally less than 1. Exceptions are for those laboratories where the correlation r ij > 0.1.
Using this matrix and Equations (5) and (6), x 0 and U(x 0 ), respectively, can be calculated: The value of x 0 is slightly higher than the certified value; however, it is contained in the interval (x 0 − U(x 0 ), x 0 + U(x 0 )). Based on the values of Equations (12) and (13) and those of Table 1, the graph in Figure 6 was obtained. The blue point represents the value of each x i from Table 1; the vertical bars represent the uncertainty U(x i ) from Table 1; the continuous orange line represents the x 0 value calculated with Equation (12); and the discontinuous orange lines represent the maximum and minimum values of U(x 0 ) calculated with Equation (13).
In this case, since all ≠ 0, the matrix has all terms. The values are generally less than 1. Exceptions are for those laboratories where the correlation > 0.1. Using this matrix and Equations (5) and (6), 0 and ( 0 ), respectively, can be calculated: The value of 0 is slightly higher than the certified value; however, it is contained in the interval ( 0 − ( 0 ), 0 + ( 0 )). Based on the values of Equations (12) and (13) and those of Table 1, the graph in Figure 6 was obtained. The blue point represents the value of each from Table 1; the vertical bars represent the uncertainty ( ) from Table 1; the continuous orange line represents the 0 value calculated with Equation (12); and the discontinuous orange lines represent the maximum and minimum values of ( 0 ) calculated with Equation (13).

Rz Roughness Parameter
Analogous to what has been done in Section 4.1, the covariance matrix was calculated: In addition, since all r ij = 0, the matrix has all terms. The values are generally less than 1 except for those laboratories where the correlation r ij > 0.1. Using this matrix and Equations (5) and (6), x 0 and U(x 0 ), respectively, can be calculated: The value of x 0 is slightly higher than the certified value; however, it is contained in the interval (x 0 − U(x 0 ), x 0 + U(x 0 )). Taking into account the values of Equations (12) and (13) and those of Table 1, the graph in Figure 7 is obtained. The blue point represents the value of each x i , from Table 1, the vertical bars represent the uncertainty U(x i ) from Table 1; the continuous orange line represents the x 0 value, calculated with Equation (15); and the discontinuous orange lines represent the maximum and minimum values of U(x 0 ) calculated with Equation (16).  (14) In addition, since all ≠ 0, the matrix has all terms. The values are generally less than 1 except for those laboratories where the correlation > 0.1. Using this matrix and Equations (5) and (6), 0 and ( 0 ), respectively, can be calculated: The value of 0 is slightly higher than the certified value; however, it is contained in the interval ( 0 − ( 0 ), 0 + ( 0 )). Taking into account the values of Equations (12) and (13) and those of Table 1, the graph in Figure 7 is obtained. The blue point represents the value of each , from Table 1, the vertical bars represent the uncertainty ( ) from Table  1; the continuous orange line represents the 0 value, calculated with Equation (15); and the discontinuous orange lines represent the maximum and minimum values of ( 0 ) calculated with Equation (16). As can be seen, in all cases, the uncertainty bars cover the interval of the average values of the measurand ( 0 − ( 0 ), 0 + ( 0 )). Therefore, the measurement results from the different participating laboratories can be said to be compatible with each other. The result that deviates the most is the one obtained using the CL confocal microscope although this is still compatible with the other laboratories. As can be seen, in all cases, the uncertainty bars cover the interval of the average values of the measurand (x 0 − U(x 0 ), x 0 + U(x 0 )). Therefore, the measurement results from the different participating laboratories can be said to be compatible with each other. The result that deviates the most is the one obtained using the CL confocal microscope although this is still compatible with the other laboratories.
The calculation results for the statistical parameters from Equations (7)-(9) are shown in Table 2. According to the results in Table 2, it is necessary to highlight that the maximum difference between x 0 and the different x i is 9.68% for Ra and 6.37% for Rz. If we consider the averages of the different x i , the difference with x 0 is 2.4% for Ra and 3.1% for Rz. As can be seen, the maximum deviation between the average value and that obtained by each of the participating laboratories was 0.08 µm, which is that of the CL confocal microscope. However, the average value x 0 is in the range (x i − U(x i ), x i + U(x i )) for this laboratory. It also can be seen that the measurement uncertainties obtained with confocal microscopes are higher than those for stylus instruments. Theoretically, this result seems to be negative. However, the most important reading is that comparable results were found when using two measurement instruments having completely different measurement principles. This result is important since the use of optical measurement instruments is simpler than that of stylus instruments and can lead to savings for ACLs in money, time and effort.
In the end, the most important analysis corresponded to the E n number (sometimes known as "normalized error" or "compatibility index"). According to ISO 17043 [5], when |E n | ≤ 1, the result can be considered satisfactory. Therefore, all results from Table 1 for Ra or for Rz are satisfactory, especially those obtained with confocal microscopes.

Conclusions
The procedure followed for this ILC roughness study, LMM-R-2019, is submitted in this report. It asked participating accredited laboratories research centres to determine the surface roughness of a measurement standard from the vertical parameters Ra and Rz. In this ILC, the participating laboratories used both mechanical contact (stylus) and optical (confocal microscopes) measurement instruments.
The results were very good, and showed that the results from the different laboratories were compatible with each other. The largest deviation was between a confocal microscope and one of the laboratories that used a stylus instrument. The mean square of the deviation E n turned out to be 0.33 in the case of the Ra roughness parameter and 0.32 in the case of the Rz roughness parameter. If the laboratories had accurately estimated their uncertainties, the root mean square value would have been close to 0.5.
This allowed us to say that they seemed to be working correctly within a wide safety margin. In other words, the uncertainties of these laboratories were probably on the order of 0.06 µm while they were providing clearly higher uncertainties.
Finally, special mention should be made of the results from a confocal microscope. It is known that the measurement of roughness with optical instruments have significant biases compared to measurements made with traditional contact stylus instruments. However, in this comparison, the measurements made with confocal microscopes were fully compatible with those made with contact stylus instruments. The problem with the microscopes, may have been due in many cases to an incorrect calibration or one that was not focused on carrying out roughness measurements. In this case, it seems that the confocal instruments used were correctly calibrated and their uncertainties reasonably well estimated, resulting in compatibility with the results from contact stylus instruments.
This ILC was designed to measure one kind of material standard. In the future, we will repeat the measurements using another material standard (e.g., isotropic and anisotropic/periodic and random) with different roughness values.

Conflicts of Interest:
The authors declare no conflict of interest.