# Measuring scale errors in a laser tracker's horizontal angle encoder through simple length measurement and two-face system tests.

1. Introduction

The spherical coordinates of a target (horizontal angle, vertical angle and range) reported by a laser tracker suffer from systematic errors due to various geometric and optical misalignments (1). Some misalignments such as encoder eccentricity, beam offset, beam tilt, etc., are modeled and collected for in the software. In this paper, we discuss an error source, the scale error in the horizontal angle encoder of the tracker, which is generally not considered during error modeling and compensation, but may potentially produce large errors in long length measurements made far away from the tracker.

The scale errors in the encoder are spacing errors in the gratings (or markings). In Coordinate Measuring Machines, scale errors are one source of linear displacement errors and can be calibrated using a laser interferometer. But an analogous technique to calibrate the angle encoder of a laser tracker is somewhat cumbersome and requires expensive and highly accurate instrumentation. For example, Gassner and Ruland (2) have performed a calibration of the angle encoder of their laser tracker using a calibrated rotary table whose accuracy exceeds that of the encoder under test. They report angle errors of about 20 [micro]rad peak-to-valley on a laser tracker. They further report that the second order harmonic is the dominant frequency in the error plot of their encoder. Ouyang et al., (3) have attempted to calibrate the angle errors in a tracker by mounting the tracker on a CMM and performing circle measurements in different planes.

As opposed to the direct methods for angle error estimation mentioned above, we propose an indirect method where we determine the parameters of a model for scale errors. The key advantage is the simplicity of the technique; we do not require expensive instrumentation such as a calibrated rotary table or even a calibrated reference length. On the other hand, our method does not provide a direct error map of the scale errors. But it is a method that can be implemented in any laboratory with reasonable degree of environmental control to obtain a fairly good estimate of the quality of the angle encoder in a tracker. Further, the method is very closely related to tests described in the ASME B89.4.19 Standard (4) and therefore manufacturers and users may already have necessary experience and materials to perform the test.

2. Approach

The errors in the scale e of the angle encoder may be decomposed into their Fourier components. The larger amplitudes are generally associated with low order harmonics, as the Gassner and Ruland (2) results suggest, and therefore higher order harmonics may be disregarded from consideration. Approximated as the sum of the first n harmonics, e is given by

e([empty set]) [approximately equal to] [n.summation over (i=1)][a.sub.i]cos(i[empty set]) + [n.summation over (i=1)][b.sub.i]sin(i[empty set]) (1)

where [a.sub.i] and [b.sub.i] are the amplitudes of the ith harmonic, and [empty set] is the measured horizontal (azimuth) angle. The true angle [[empty set].sub.t] assuming other sources of errors are absent, can then be given by: [[empty set].sub.t],=[empty set] - e ([empty set]).

Detecting Odd Orders From Two-Face Tests:

If the scale error contains only odd order harmonics, the error in the measured horizontal angle in the front face is equal in magnitude but opposite in sign to that measured in the back face. A two-face system test (the apparent distance between measurements of a stationary target in front face and then in back face) as described in (4) therefore produces a non-negative error (or apparent distance) in this case.

This sensitivity to front-face back-face measurements may, in principle, be exploited to determine the parameters ([a.sub.i],-, [b.sub.i]) for odd orders in Eq. (1). That is, half the difference between front-face azimuth measurements and back-face azimuth measurements to a stationary target placed at different azimuthal positions [empty set] yields errors e([empty set]) that may be used to determine parameters [a.sub.i] and [b.sub.i] for odd orders. It should be noted that in addition to being sensitive to odd order harmonics, such front-face back-face measurements are also sensitive to numerous other sources of misalignment such as collimation error, non-orthogonality between the standing axis and the transit axis, etc. The other sources generally produce errors that are constant (because vertical angle does not change between measurements performed at different azimuths); therefore residuals from the mean of half the difference between front face and back face measurements are to be considered when determining odd orders from Eq. (1).

It should be pointed out however that estimating odd order harmonics is not critical because their effect can be removed by averaging between front face and back face measurements. We therefore do not consider odd orders in detail in this paper.

Detecting Even Orders:

If the scale error contained purely even order terms (n = 2, 4, 6, ...), the two-face test (4) produces a result of zero because the error in the measured horizontal angle at [empty set] is identical in magnitude and sign to the error at [empty set] + 180[degrees]. Even order harmonics are therefore not sensitive to two-face tests.

It is therefore proposed to use simple length measurement system tests to determine parameters of the model in Eq. (1), for even order harmonics. The position and orientation of the length has to be chosen carefully so that it is sensitive primarily to encoder scale errors, and not to other geometric errors that may be a function of range or vertical angle.

Estimating low even order harmonics is critical because, if present, they cannot be removed by averaging front and back face measurements as mentioned earlier. The encoder's error map has to be compensated for even order harmonics to enhance the tracker's accuracy.

Experimental Setup:

A horizontal length L symmetrically placed with respect to the tracker and positioned at roughly tracker height (as described in the ASME B89.4.19 Standard) is ideally suited for this purpose (see Fig. 1). Because the range and the vertical angle to the two targets are equal, we eliminate any sources that may produce horizontal angle errors that are functions of nominal range and nominal vertical angle (for example, those arising from a squareness error, collimation error, beam offset along the transit axis, etc.). On the other hand, the test is sensitive to those sources that produce horizontal angle errors that are a function of the nominal horizontal angle, such as encoder scale errors and eccentricity of the encoder with respect to the vertical axis. Encoder eccentricity is a first order harmonic and is usually already compensated for in the software, leaving only higher order harmonics which are due to scale errors.

[FIGURE 1 OMITTED]

The test involves measuring a fixed length placed at different azimuthal positions of the tracker, in a manner similar to the horizontal length tests in the ASME B89.4.19 Standard. As a matter of practical convenience, we realize the length not by using a physical artifact, but as the distance between two targets on stands. After measuring the distance between the targets using the tracker, we rotate the tracker about its vertical axis and re-measure the distance between the targets. We repeat this process several times to cover a full circle (360[degrees] of the encoder). It is not necessary to take special care when performing this rotation. The mathematics involved in determining the scale error incorporates the actual range and azimuths as measured to the targets.

The error E in the measured length L can be calculated from some simple trigonometry as shown, see Fig. 2.

[FIGURE 2 OMITTED]

The angle [theta] is half the angle subtended by the reference length AB at the tracker, D is the distance from the tracker to the length AB, R is the range to the target, and [theta] is the azimuthal position of the tracker as measured to the center of the length AB (see Fig. 1 and Fig. 2). The azimuthal position of targets A is [empty set] + [theta]. The scale error at A is e([empty set] + [theta]). The range to target A is R. According to the adopted convention, a positive scale error implies the measured angle is larger than the true angle. We assume that e([empty set] + [theta]) is positive and therefore the measured angle is larger than the true angle. Therefore, the true position of point A is at [A.sub.1], where the azimuth is [empty set] + [theta] - e([empty set] + [theta]). [AA.sub.1] is the error vector. Its magnitude is Re([empty set] + [theta]). The component of the error along the length [AA.sub.2] is given by [AA.sub.1]cos([theta]) = R e([empty set] + [theta]) cos ([theta]). The difference in the errors at A and B (the azimuth at B is [empty set] - [theta]) is the net error in the length and is given by

E([theta], [empty set]) = Re([empty set] + [theta])cos[theta] - Re([empty set] - [theta]) cos [theta]. (2)

For the purposes of this experiment, the distance between the targets does not have to be calibrated by other means; it is only required that the distance not change during the duration of the experiments. We consider the residuals from the average of the multiple length measurements at different azimuths [empty set] as the errors due to scale, E, in Eq. (2) above. It should be noted that there may be other sources of systematic error such as, for instance, a ranging error in the tracker that will remain identical for all the lengths measured. These systematic sources produce a bias (from the true but unknown length between A and B) which is removed by discarding the mean of all the measured lengths and considering only the residuals. We can then fit Eq. (2) to the residuals to determine parameters ([a.sub.i], [b.sub.i]) for both even and odd harmonics. It should be noted that systematic errors in the ranging system are detected by a separate ranging test (4). We also point out that there may other sources of error, such as a cyclical error in range that may not be identical for all measured lengths, and must be accounted for in the uncertainty budget.

We note that the purpose of correcting scale errors is to improve the accuracy of the tracker, and that typically, high accuracy measurements involve averaging front face and back face measurements at every target position. Averaging front and back face measurements, as mentioned earlier, removes the influence of several geometric misalignments including odd order harmonics in the scale error. We therefore adopt this approach in our experiments also, hence averaging front and back face measurements means that we only attempt to fit even order harmonics to Eq. (2). As mentioned earlier, we can always assess the presence of odd order harmonics from difference of front face and back face measurements made at the different azimuths. Fig. 3 illustrates why the horizontal length test is sensitive to second order scale error in the encoder.

[FIGURE 3 OMITTED]

3. Sensitivity Analysis for Second Order Harmonic

Because the second order harmonic is expected to be the dominant error source, we take a closer look at optimizing the test variables (D and L) to obtain maximum sensitivity to the second order harmonic. We briefly address the issue of determining the length L between targets A and B, and distance D from the tracker to obtain maximum sensitivity to second order harmonic scale error in the encoder (see Fig. 1).

The error E in the measured length L due to a second order harmonic scale error can be given by

E([theta], [empty set]) = Rcos[theta][e([empty set] + [theta]) - e([empty set] - [theta])] = Rcos[theta][[a.sub.2]cos(2([empty set] + [theta])) + [b.sub.2]sin(2[[empty set] + [theta])) - [a.sub.2] cos(2([empty set] + [theta])) + [b.sub.2]sin(2([empty set] + [theta])) - [a.sub.2]cos(2([empty set] - [theta])) - [b.sub.2]sin(2([empty set] - [theta]))]. (3)

Eq. (3) can be simplified to

E([theta],[empty set]) = 2 Rcos[theta] sin2[theta] [square root of [a.sub.2.sup.2] + [b.sub.2.sup.2]] sin ([chi] - 2[empty set]) (4)

where

[chi] = [sin.sup.-1] ([[b.sub.2]/[[square root of [a.sub.2.sup.2] + [b.sub.2.sup.2]]]

For any given length L and distance D, E is maximum when sin ([chi] - 2 [empty set]) = 1. The azimuthal location [empty set] where E is maximum can therefore be determined from the above condition if [a.sub.2] and [b.sub.2] are known.

Our objective here is not only to determine where the maximum error occurs, but more importantly to determine optimal values for D and L so that the error at this angular position is as large as possible to achieve maximum sensitivity. Thus, the objective is to maximize the function [Florin] = 2R cos [theta] sin 2 [theta], which represents the sensitivity of the measurement, i.e., the length error in micrometers for 1 [micro]rad of second order scale error.

Expressing [Florin] in terms of the quantities of interest, we have

[Florin] = 2Rcos[theta]sin2[theta] = [[2[D.sup.2]L]/[[D.sup.2] + [[L.sup.2]/4]]] (5)

We plot [Florin] as a function of D and L in the contour plot shown in Fig. 4. Clearly, a large length L placed far away from the tracker provides maximum sensitivity. For example, the maximum error in the length for 1 [micro]rad of second order scale error is 5.6 [micro]m for a 3.5 m length placed 3.5 m away, but increases to 12.8 [micro]m for an 8 m length placed 8 m away.

[FIGURE 4 OMITTED]

The above sensitivity analysis does not however account for the decreasing accuracy in tracker measurements made farther away or for the effect of variations in the distance between the two targets due to short term thermal fluctuations. We address the issue of uncertainty in the calculated parameters [a.sub.2] and [b.sub.2], and in the resulting error map after we present the experimental results.

4. Experimental Results

We performed the length measurement system tests as described in Sec. 2 on six trackers from three different manufacturers and also at different combinations of L and D. Trackers A and B are older models and E is a newer model from one manufacturer. Tracker C is an older model and F is a newer model from a second manufacturer. Tracker D is an older model from a third manufacturer. The different combinations of L and D were chosen to assess the validity of the analysis discussed earlier. We present the results here.

Tracker A

As mentioned earlier, we measured each target in front and back face. We show a plot of the difference between front and back face azimuth measurements as a function of azimuthal position in Fig. 5. There is no evidence of any low order odd harmonics in the plot. We do not therefore consider mapping odd orders here. Also, as mentioned earlier, any odd order errors in Fig. 5 can be removed by averaging front and back face measurements. We therefore focus on low even orders next.

[FIGURE 5 OMITTED]

Figure 6 shows the length errors (residues) as a function of azimuth for two combinations of L and D. The length errors were calculated after averaging front face and back face measurements at each target position, thereby removing the influence of any odd order harmonics.

[FIGURE 6 OMITTED]

The dominant harmonic appears to be the second order from Fig. 6. A least-squares best-fit yielded an amplitude of 4.3 [micro]rad ([square root of.[a.sub.2.sup.2] + [b.sub.2.sup.2]]) for the second order harmonic for combination 1 (D = 7 m, L = 7.75 m) and 4.1 [micro]rad for combination 2 (D = 4.9 m, Z = 5.4 m). This agreement of the second order harmonic to within 0.2 [micro]rad is within its uncertainty as shown in the next section.

As a quick note, we point out that some simple calculations validate the analysis described in the previous section. The sensitivity [Florin] for combination 1 (D = 7 m, L = 7.75 m) is 11.8 [micro]m/[micro]rad. The sensitivity for combination 2(D = 4.9m, L = 5.4 m) is 8.3 [micro]m/[micro]rad. The errors (residuals from the mean) in the lengths are within [+ or -] 65 [micro]m for combination 1 from Fig. 6, and within [+ or -] 43 [micro]m for combination 2. The ratio of the maximum error to the sensitivity provides an estimate of the amplitude of second order harmonic scale error, and should be nearly the same regardless of the combination of L and D used. This of course assumes that the second order is the dominant term, as is the case here. For combination 1, that ratio is 65/11.8 = 5.5 [micro]rad. For combination 2, it is 43/8.3 = 5.2 [micro]rad. They agree fairly well, as expected. There is however some discrepancy between these estimates and the 4.3 [micro]rad (or the 4.1 [micro]rad) value calculated earlier because there are some other higher order even harmonics in the data as well.

A rigorous approach to the task of validating our claims would be to correct the error map of our tracker and perform the experiments again; we have in fact done so and the results of our experiments after correcting the tracker's azimuth encoder map are shown in Fig. 7. It is clear that there is a very small, if any, second order component in the length errors. There is substantial reduction in errors after collecting the error map; the errors are less than [+ or -] 25 [micro]m by removing just the second order alone.

[FIGURE 7 OMITTED]

The benefits of such second order correction are more clearly observed when measuring large lengths far away from the tracker. A further collection of other even order harmonics may also be performed, but we have not done so with our tracker because the amplitude of these higher orders are almost near the uncertainty achievable with our method. We should point out that downloading the error map from trackers, correcting it, and uploading the new map, while not very challenging, may require special support from the manufacturer of the tracker.

Tracker B

Tracker B is also an older model from the same manufacturer as tracker A. This tracker also displayed the second order error characteristic as described earlier for tracker A, but with smaller amplitude of about 2.5 [micro]rad.

Tracker C

Tracker C is an older model from a second manufacturer. Fig. 8 shows the difference between front and back face azimuth measurements to a target as a function of azimuth. Low odd orders, if present, are small in magnitude. As mentioned earlier, we do not compensate for odd orders; we simply average front and back face measurements and remove any odd orders.

[FIGURE 8 OMITTED]

Figure 9 shows the results from length tests performed for the two combinations discussed earlier. The dominant harmonic appears to be the second order. A least-squares best-fit yielded an amplitude of 5.0 [micro]rad ([square root of [a.sub.2.sup.2] + [b.sub.2.sup.2]]) for the second order harmonic for combination 1 (D = 7 m, L = 7.75 m) and 4.8 [micro]rad for combination 2(D = 4.9 m, L = 5.4 m). This agreement of the second order harmonic to within 0.2 [micro]rad is within its uncertainty as shown in the next section.

[FIGURE 9 OMITTED]

Fig. 10 shows the results from the length tests after correcting the tracker's error map for second order errors. There is substantial reduction in errors after correcting the error map; the errors are less than [+ or -]25 [micro]m by removing just the second order alone.

[FIGURE 10 OMITTED]

Tracker D

The amplitude of the second order was about 1 [micro]rad for this old tracker from a third manufacturer, much smaller than that seen for the other two trackers. The horizontal length tests however showed errors of about [+ or -] 40 [micro]m, possibly because of higher order harmonics.

5. Discussion

Low order errors can be removed in several different ways. Odd orders may be removed by having two read-heads positioned 180[degrees] apart on the encoder. Even orders may be removed by having two readheads 90[degrees] apart. The cost and complexity of having additional readheads may be mitigated by a calibration of the encoder through independent means such as by using a high accuracy rotary stage, a polygon with auto-collimator etc. The error mapping technique we have discussed in this paper is another approach for the same.

We have tested new tracker designs from two manufacturers, trackers E and F. The results of the horizontal length tests are shown in Fig. 11 for tracker E that intrinsically compensates for second order errors. The errors are much smaller than those in Fig. 6 and, as the manufacturer expected, no second order error is noticeable. We also tested a new tracker (Tracker F) from another manufacturer; the results are shown in Fig. 12. There is no evidence of second order errors.

[FIGURE 11 OMITTED]

[FIGURE 12 OMITTED]

6. A Note on Higher Order Harmonics

The preceding sections described low order harmonic errors in the scale and a simple technique to estimate those errors. There is an analogous effect at the high frequency end where periodic errors have been observed at a frequency representing the spacing between lines in an encoder. We do not directly address compensating high frequency encoder errors in this paper, but point out that a scheme similar to our proposed approach may be considered for testing purposes where an uncalibrated horizontal length can be measured by a tracker at closely spaced azimuthal positions so that multiple sampling positions are contained between two adjacent encoder lines. Two-face tests may also be performed so that multiple sampling positions are contained between two adjacent encoder lines. Such tests may reveal interpolation errors in the encoder.

7. Uncertainty Considerations

There is an uncertainty associated with each of the lengths that are measured as part of the test described in Sec. 2, and therefore an uncertainty in the parameters [a.sub.2] and [b.sub.2]. We address the issue of determining the uncertainty in these parameters here. Values are calculated for combination 1 (D = 7 m, L = 7.75 m; the range to targets R is 8 m and the half angle [theta] subtended by the length is 29[degrees]).

The tracker's ranging accuracy deteriorates farther away from the tracker. Manufacturers typically specify an MPE (maximum permissible error) for the range (for example, one of the trackers we tested had an MPE specification of 20 [micro]m + 0.8 [micro]m/m). Because the length measurements are always performed so that the target is at the same range (about 8 m) from the tracker, a large portion of the error in the range is expected to be due to systematic sources which produce the same errors in length for all horizontal length measurements. We only consider residuals from the mean of all length measurements. Therefore the error in the range due to systematic sources is not of any consequence, and using the MPE for estimating uncertainty is not appropriate for this type of measurement. On the other hand, the variation in repeated range measurements (i.e., the repeatability) is the noise in the measurement and therefore is the primary source of error that is relevant in our case. We experimentally determine the standard uncertainty due to repeatability in range as 1.5 [micro]m for one of the trackers we tested. Only the component of the ranging uncertainty along the length will actually contribute to a length error. Further, the length is the difference in measurements made at two targets and we average 2 measurements at each target location. Therefore the standard uncertainty in the measured length due to ranging error non-repeatability in the tracker is 1.5 sin[theta] = 0.7 [micro]m.

Manufacturers also sometimes specify an MPE for transverse measurements (for example, one of the trackers we tested had an MPE specification of 36 [micro]m + 6 [micro]m/m for angle measurements). As discussed in the previous paragraph, this MPE incorporates all sources that produce errors in the measured angle, including the second order harmonic scale error in the encoder. Further, this MPE is applicable for measurements made in the front face of the tracker, whereas all our measurements are performed in both faces and averaged. Therefore, using the MPE as a bound on the uncertainty is not appropriate for this type of measurement. Because it is the purpose of our experiments to determine the extent of one source of systematic error--namely the second order harmonic, it is necessary to primarily assess the repeatability of the angle measurement system which introduces noise into the measurement. We experimentally measured the standard uncertainty due to repeatability in the horizontal angle to be 0.0001[degrees] for one of the trackers we tested. At a range of 8 m, considering 2 targets, 2 repeat measurements at each target, and only the component along the length, the standard uncertainty in length due to repeatability in the measured horizontal angle is 8(0.0001[pi]/180)cos [theta] = 12.2 [micro]m.

As long lengths are considered in the experiments, stability of the length (realized between two targets on stands) during the duration of the experiment (about 1 hour) is important. To estimate the stability of the length, we mount an interferometer on one stand and a target on the other, and record the length in intervals of 5 seconds over a one hour period. The standard uncertainty in the measured lengths is about 1.2 [micro]m.

The terms discussed above may be summed in quadrature to yield a standard uncertainty of about 12.2 [micro]m in the length (for the 7.75 m length of combination 1). This uncertainty is the 'noise' and is to be contrasted against the 'signal' which is the error in the length (residual from mean length) due to the second order scale error of the encoder. The length error due to second order harmonic scale error is within [+ or -]65 [micro]m as pointed out earlier.

We perform 18 such length measurements from which we extract 2 parameters ([a.sub.2] and [b.sub.2]--the coefficients of the second order harmonic). A Monte Carlo Simulation yields a standard uncertainty of 0.4 [micro]rad on the parameters [a.sub.2] and [b.sub.2] resulting in a 0.4 [micro]rad uncertainty on an arbitrary point in the error map itself (the uncertainty in the scale error map e([empty set]) at any azimuthal position [empty set], where e([empty set]) = [a.sub.2]cos(2[empty set]) + [b.sub.2]sin(2[empty set]). Other simulations suggest that the standard uncertainty on the parameters remains 0.4 [micro]rad if an attempt is made to estimate both the 2nd and 4th orders from the 18 length measurements, and the standard uncertainty on an arbitrary point in the error map (uncertainty on e([empty set]) = [a.sub.2]cos(2[empty set])+[b.sub.2]sin(2[empty set])+[a.sub.4]cos(4[empty set])+[b.sub.4]sin(4[empty set]) is about 0.5 [micro]rad.

The standard uncertainty in the length for combination 2 (D = 4.9 m, L = 5.4 m; the range to targets R is 5.6 m and the half angle [theta] subtended by the length is 29[degrees]) can be calculated in a manner similar to that described above, and is 9 [micro]m. Again, we distinguish between this 'noise' and the 'signal', where the 'signal' is the error in the length (residual from mean) due to second order harmonic scale error, which is within [+ or -]43 [micro]m. The standard uncertainty in the two parameters calculated from 18 length measurements ([a.sub.2] and [b.sub.2]--the coefficients of the second order harmonic) is about 0.4 [micro]rad and there is also a 0.4 [micro]rad uncertainty on an arbitrary point in the error map itself.

The maximum signal-to-noise ratio for combination l (D = 7m, L = 7.75 m) is 65/12.2 = 5.3, and is marginally better than the maximum signal-to-noise ratio for combination 2 (D = 4.9 m, L = 5.4 m) which is 43/9 = 4.7. This does seem to suggest that a longer length placed farther away is more suitable for this test, but practical limitations such as space constraints will most likely limit the size of the length and its placement. A 7.75 m long length placed 7 m away does appear to be adequate in mapping out second order errors.

8. Summary

A tracker is designed to measure large lengths at long distances with a high degree of accuracy. The 7.75 m long length placed 7 m away that was used in this study is a fairly representative measurement for a tracker. Older trackers from two different manufacturers that we tested produced length errors as large as 65 [micro]m due to second order scale error. The maximum permissible error for this measurement is approximately 120 [micro]m based on a typical manufacturer's specifications, without considering temperature effects. The second order scale error is a significant fraction of the MPE.

The method described in this paper provides a simple technique for estimating scale errors in the horizontal angle encoder of a laser tracker. It does not require expensive instrumentation or even a calibrated artifact. It is quick and easy to perform and can be realized in an environment with reasonable degree of stability.

We did correct the horizontal angle encoder error map of two trackers which were from two different manufacturers. The collections were based on the horizontal length tests as described in this paper and, as expected, subsequent length measurement system tests revealed the absence of any second order terms. The length errors reduced from a maximum of about 65 [micro]m before correction to less than 25 [micro]m after collection. Newer trackers from two different manufacturers did not reveal this second order error; neither did an old tracker from the third manufacture.

Acknowledgements

The authors express their gratitude to Georg Gassner for detailed responses to queries about their technique and to Tyler Estler, John Palmateer, and Ted Vorburger for reviewing the manuscript.

Disclaimer: Commercial equipment and materials may be identified in order to adequately specify certain procedures. In no case does such identification imply recommendation or endorsement by the National Institute of Standards and Technology nor does it imply that the materials or equipment identified are necessarily the best available for the purpose.

9. References

(1) B. Muralikrishnan. D. Sawyer, C. Blackburn. S. Phillips. B. Borchardt. and W. T. Estler. ASME B89.4.19 Performance evaluation tests and geometric misalignments in laser trackers. Journal of Research of the National Institute of Standards and Technology. Volume 114. Number 1. January-February 2009. p. 21-35.

(2) Georg Gassner and Robert Ruland. Laser tracker calibration--Testing the angle measurement system. SLAC National Accelerator Laboratory Publication. SLAC-PUB-13476. December 2008.

(3) J. F. Ouyang. W. L. Liu. Y. G Yan. and D. X. Sun, Angular error calibration of laser tracker system. Advanced Laser Technologies 2005. Ivan A. Shcherbakov. Kexin Xu. Qingyue Wang. Alexander V. Priezzhev. and Vladimir I. Pustovoy, eds.. Proc. of SPIE Vol. 6344.

(4) ASME B89.4.19--2006 Performance Evaluation of Laser-Based Spherical Coordinate Measurement Systems. American Society of Mechanical Engineers, http://www.asme.org

About the authors: Bala Muralikrishnan is a guest researcher, Chris Blackburn is a physical science technician, Dan Sawyer is a mechanical engineer and Steve Phillips is a physicist; all are with the Precision Engineering Division of the NIST Manufacturing Engineering Laboratoiy. Robert Bridges is the Chief Scientist at FARO Technologies Inc.

The National Institute of Standards and Technology is an agency of the U.S. Department of Commerce.

B. Muralikrishnan, C. Blackburn, D. Sawyer, S. Phillips,

Precision Engineering Division, National Institute of Standards and Technology, Gaithersburg, MD 20899

and

R. Bridges

FARO Technologies, Kennett Square, PA 19348

bala.muralikrishnan@nist.gov

cluistopher.blackburn@nist.gov

daniel.sawyer@nist.gov

steven.phillips@nist.gov

robert.bridges@faro.com

The spherical coordinates of a target (horizontal angle, vertical angle and range) reported by a laser tracker suffer from systematic errors due to various geometric and optical misalignments (1). Some misalignments such as encoder eccentricity, beam offset, beam tilt, etc., are modeled and collected for in the software. In this paper, we discuss an error source, the scale error in the horizontal angle encoder of the tracker, which is generally not considered during error modeling and compensation, but may potentially produce large errors in long length measurements made far away from the tracker.

The scale errors in the encoder are spacing errors in the gratings (or markings). In Coordinate Measuring Machines, scale errors are one source of linear displacement errors and can be calibrated using a laser interferometer. But an analogous technique to calibrate the angle encoder of a laser tracker is somewhat cumbersome and requires expensive and highly accurate instrumentation. For example, Gassner and Ruland (2) have performed a calibration of the angle encoder of their laser tracker using a calibrated rotary table whose accuracy exceeds that of the encoder under test. They report angle errors of about 20 [micro]rad peak-to-valley on a laser tracker. They further report that the second order harmonic is the dominant frequency in the error plot of their encoder. Ouyang et al., (3) have attempted to calibrate the angle errors in a tracker by mounting the tracker on a CMM and performing circle measurements in different planes.

As opposed to the direct methods for angle error estimation mentioned above, we propose an indirect method where we determine the parameters of a model for scale errors. The key advantage is the simplicity of the technique; we do not require expensive instrumentation such as a calibrated rotary table or even a calibrated reference length. On the other hand, our method does not provide a direct error map of the scale errors. But it is a method that can be implemented in any laboratory with reasonable degree of environmental control to obtain a fairly good estimate of the quality of the angle encoder in a tracker. Further, the method is very closely related to tests described in the ASME B89.4.19 Standard (4) and therefore manufacturers and users may already have necessary experience and materials to perform the test.

2. Approach

The errors in the scale e of the angle encoder may be decomposed into their Fourier components. The larger amplitudes are generally associated with low order harmonics, as the Gassner and Ruland (2) results suggest, and therefore higher order harmonics may be disregarded from consideration. Approximated as the sum of the first n harmonics, e is given by

e([empty set]) [approximately equal to] [n.summation over (i=1)][a.sub.i]cos(i[empty set]) + [n.summation over (i=1)][b.sub.i]sin(i[empty set]) (1)

where [a.sub.i] and [b.sub.i] are the amplitudes of the ith harmonic, and [empty set] is the measured horizontal (azimuth) angle. The true angle [[empty set].sub.t] assuming other sources of errors are absent, can then be given by: [[empty set].sub.t],=[empty set] - e ([empty set]).

Detecting Odd Orders From Two-Face Tests:

If the scale error contains only odd order harmonics, the error in the measured horizontal angle in the front face is equal in magnitude but opposite in sign to that measured in the back face. A two-face system test (the apparent distance between measurements of a stationary target in front face and then in back face) as described in (4) therefore produces a non-negative error (or apparent distance) in this case.

This sensitivity to front-face back-face measurements may, in principle, be exploited to determine the parameters ([a.sub.i],-, [b.sub.i]) for odd orders in Eq. (1). That is, half the difference between front-face azimuth measurements and back-face azimuth measurements to a stationary target placed at different azimuthal positions [empty set] yields errors e([empty set]) that may be used to determine parameters [a.sub.i] and [b.sub.i] for odd orders. It should be noted that in addition to being sensitive to odd order harmonics, such front-face back-face measurements are also sensitive to numerous other sources of misalignment such as collimation error, non-orthogonality between the standing axis and the transit axis, etc. The other sources generally produce errors that are constant (because vertical angle does not change between measurements performed at different azimuths); therefore residuals from the mean of half the difference between front face and back face measurements are to be considered when determining odd orders from Eq. (1).

It should be pointed out however that estimating odd order harmonics is not critical because their effect can be removed by averaging between front face and back face measurements. We therefore do not consider odd orders in detail in this paper.

Detecting Even Orders:

If the scale error contained purely even order terms (n = 2, 4, 6, ...), the two-face test (4) produces a result of zero because the error in the measured horizontal angle at [empty set] is identical in magnitude and sign to the error at [empty set] + 180[degrees]. Even order harmonics are therefore not sensitive to two-face tests.

It is therefore proposed to use simple length measurement system tests to determine parameters of the model in Eq. (1), for even order harmonics. The position and orientation of the length has to be chosen carefully so that it is sensitive primarily to encoder scale errors, and not to other geometric errors that may be a function of range or vertical angle.

Estimating low even order harmonics is critical because, if present, they cannot be removed by averaging front and back face measurements as mentioned earlier. The encoder's error map has to be compensated for even order harmonics to enhance the tracker's accuracy.

Experimental Setup:

A horizontal length L symmetrically placed with respect to the tracker and positioned at roughly tracker height (as described in the ASME B89.4.19 Standard) is ideally suited for this purpose (see Fig. 1). Because the range and the vertical angle to the two targets are equal, we eliminate any sources that may produce horizontal angle errors that are functions of nominal range and nominal vertical angle (for example, those arising from a squareness error, collimation error, beam offset along the transit axis, etc.). On the other hand, the test is sensitive to those sources that produce horizontal angle errors that are a function of the nominal horizontal angle, such as encoder scale errors and eccentricity of the encoder with respect to the vertical axis. Encoder eccentricity is a first order harmonic and is usually already compensated for in the software, leaving only higher order harmonics which are due to scale errors.

[FIGURE 1 OMITTED]

The test involves measuring a fixed length placed at different azimuthal positions of the tracker, in a manner similar to the horizontal length tests in the ASME B89.4.19 Standard. As a matter of practical convenience, we realize the length not by using a physical artifact, but as the distance between two targets on stands. After measuring the distance between the targets using the tracker, we rotate the tracker about its vertical axis and re-measure the distance between the targets. We repeat this process several times to cover a full circle (360[degrees] of the encoder). It is not necessary to take special care when performing this rotation. The mathematics involved in determining the scale error incorporates the actual range and azimuths as measured to the targets.

The error E in the measured length L can be calculated from some simple trigonometry as shown, see Fig. 2.

[FIGURE 2 OMITTED]

The angle [theta] is half the angle subtended by the reference length AB at the tracker, D is the distance from the tracker to the length AB, R is the range to the target, and [theta] is the azimuthal position of the tracker as measured to the center of the length AB (see Fig. 1 and Fig. 2). The azimuthal position of targets A is [empty set] + [theta]. The scale error at A is e([empty set] + [theta]). The range to target A is R. According to the adopted convention, a positive scale error implies the measured angle is larger than the true angle. We assume that e([empty set] + [theta]) is positive and therefore the measured angle is larger than the true angle. Therefore, the true position of point A is at [A.sub.1], where the azimuth is [empty set] + [theta] - e([empty set] + [theta]). [AA.sub.1] is the error vector. Its magnitude is Re([empty set] + [theta]). The component of the error along the length [AA.sub.2] is given by [AA.sub.1]cos([theta]) = R e([empty set] + [theta]) cos ([theta]). The difference in the errors at A and B (the azimuth at B is [empty set] - [theta]) is the net error in the length and is given by

E([theta], [empty set]) = Re([empty set] + [theta])cos[theta] - Re([empty set] - [theta]) cos [theta]. (2)

For the purposes of this experiment, the distance between the targets does not have to be calibrated by other means; it is only required that the distance not change during the duration of the experiments. We consider the residuals from the average of the multiple length measurements at different azimuths [empty set] as the errors due to scale, E, in Eq. (2) above. It should be noted that there may be other sources of systematic error such as, for instance, a ranging error in the tracker that will remain identical for all the lengths measured. These systematic sources produce a bias (from the true but unknown length between A and B) which is removed by discarding the mean of all the measured lengths and considering only the residuals. We can then fit Eq. (2) to the residuals to determine parameters ([a.sub.i], [b.sub.i]) for both even and odd harmonics. It should be noted that systematic errors in the ranging system are detected by a separate ranging test (4). We also point out that there may other sources of error, such as a cyclical error in range that may not be identical for all measured lengths, and must be accounted for in the uncertainty budget.

We note that the purpose of correcting scale errors is to improve the accuracy of the tracker, and that typically, high accuracy measurements involve averaging front face and back face measurements at every target position. Averaging front and back face measurements, as mentioned earlier, removes the influence of several geometric misalignments including odd order harmonics in the scale error. We therefore adopt this approach in our experiments also, hence averaging front and back face measurements means that we only attempt to fit even order harmonics to Eq. (2). As mentioned earlier, we can always assess the presence of odd order harmonics from difference of front face and back face measurements made at the different azimuths. Fig. 3 illustrates why the horizontal length test is sensitive to second order scale error in the encoder.

[FIGURE 3 OMITTED]

3. Sensitivity Analysis for Second Order Harmonic

Because the second order harmonic is expected to be the dominant error source, we take a closer look at optimizing the test variables (D and L) to obtain maximum sensitivity to the second order harmonic. We briefly address the issue of determining the length L between targets A and B, and distance D from the tracker to obtain maximum sensitivity to second order harmonic scale error in the encoder (see Fig. 1).

The error E in the measured length L due to a second order harmonic scale error can be given by

E([theta], [empty set]) = Rcos[theta][e([empty set] + [theta]) - e([empty set] - [theta])] = Rcos[theta][[a.sub.2]cos(2([empty set] + [theta])) + [b.sub.2]sin(2[[empty set] + [theta])) - [a.sub.2] cos(2([empty set] + [theta])) + [b.sub.2]sin(2([empty set] + [theta])) - [a.sub.2]cos(2([empty set] - [theta])) - [b.sub.2]sin(2([empty set] - [theta]))]. (3)

Eq. (3) can be simplified to

E([theta],[empty set]) = 2 Rcos[theta] sin2[theta] [square root of [a.sub.2.sup.2] + [b.sub.2.sup.2]] sin ([chi] - 2[empty set]) (4)

where

[chi] = [sin.sup.-1] ([[b.sub.2]/[[square root of [a.sub.2.sup.2] + [b.sub.2.sup.2]]]

For any given length L and distance D, E is maximum when sin ([chi] - 2 [empty set]) = 1. The azimuthal location [empty set] where E is maximum can therefore be determined from the above condition if [a.sub.2] and [b.sub.2] are known.

Our objective here is not only to determine where the maximum error occurs, but more importantly to determine optimal values for D and L so that the error at this angular position is as large as possible to achieve maximum sensitivity. Thus, the objective is to maximize the function [Florin] = 2R cos [theta] sin 2 [theta], which represents the sensitivity of the measurement, i.e., the length error in micrometers for 1 [micro]rad of second order scale error.

Expressing [Florin] in terms of the quantities of interest, we have

[Florin] = 2Rcos[theta]sin2[theta] = [[2[D.sup.2]L]/[[D.sup.2] + [[L.sup.2]/4]]] (5)

We plot [Florin] as a function of D and L in the contour plot shown in Fig. 4. Clearly, a large length L placed far away from the tracker provides maximum sensitivity. For example, the maximum error in the length for 1 [micro]rad of second order scale error is 5.6 [micro]m for a 3.5 m length placed 3.5 m away, but increases to 12.8 [micro]m for an 8 m length placed 8 m away.

[FIGURE 4 OMITTED]

The above sensitivity analysis does not however account for the decreasing accuracy in tracker measurements made farther away or for the effect of variations in the distance between the two targets due to short term thermal fluctuations. We address the issue of uncertainty in the calculated parameters [a.sub.2] and [b.sub.2], and in the resulting error map after we present the experimental results.

4. Experimental Results

We performed the length measurement system tests as described in Sec. 2 on six trackers from three different manufacturers and also at different combinations of L and D. Trackers A and B are older models and E is a newer model from one manufacturer. Tracker C is an older model and F is a newer model from a second manufacturer. Tracker D is an older model from a third manufacturer. The different combinations of L and D were chosen to assess the validity of the analysis discussed earlier. We present the results here.

Tracker A

As mentioned earlier, we measured each target in front and back face. We show a plot of the difference between front and back face azimuth measurements as a function of azimuthal position in Fig. 5. There is no evidence of any low order odd harmonics in the plot. We do not therefore consider mapping odd orders here. Also, as mentioned earlier, any odd order errors in Fig. 5 can be removed by averaging front and back face measurements. We therefore focus on low even orders next.

[FIGURE 5 OMITTED]

Figure 6 shows the length errors (residues) as a function of azimuth for two combinations of L and D. The length errors were calculated after averaging front face and back face measurements at each target position, thereby removing the influence of any odd order harmonics.

[FIGURE 6 OMITTED]

The dominant harmonic appears to be the second order from Fig. 6. A least-squares best-fit yielded an amplitude of 4.3 [micro]rad ([square root of.[a.sub.2.sup.2] + [b.sub.2.sup.2]]) for the second order harmonic for combination 1 (D = 7 m, L = 7.75 m) and 4.1 [micro]rad for combination 2 (D = 4.9 m, Z = 5.4 m). This agreement of the second order harmonic to within 0.2 [micro]rad is within its uncertainty as shown in the next section.

As a quick note, we point out that some simple calculations validate the analysis described in the previous section. The sensitivity [Florin] for combination 1 (D = 7 m, L = 7.75 m) is 11.8 [micro]m/[micro]rad. The sensitivity for combination 2(D = 4.9m, L = 5.4 m) is 8.3 [micro]m/[micro]rad. The errors (residuals from the mean) in the lengths are within [+ or -] 65 [micro]m for combination 1 from Fig. 6, and within [+ or -] 43 [micro]m for combination 2. The ratio of the maximum error to the sensitivity provides an estimate of the amplitude of second order harmonic scale error, and should be nearly the same regardless of the combination of L and D used. This of course assumes that the second order is the dominant term, as is the case here. For combination 1, that ratio is 65/11.8 = 5.5 [micro]rad. For combination 2, it is 43/8.3 = 5.2 [micro]rad. They agree fairly well, as expected. There is however some discrepancy between these estimates and the 4.3 [micro]rad (or the 4.1 [micro]rad) value calculated earlier because there are some other higher order even harmonics in the data as well.

A rigorous approach to the task of validating our claims would be to correct the error map of our tracker and perform the experiments again; we have in fact done so and the results of our experiments after correcting the tracker's azimuth encoder map are shown in Fig. 7. It is clear that there is a very small, if any, second order component in the length errors. There is substantial reduction in errors after collecting the error map; the errors are less than [+ or -] 25 [micro]m by removing just the second order alone.

[FIGURE 7 OMITTED]

The benefits of such second order correction are more clearly observed when measuring large lengths far away from the tracker. A further collection of other even order harmonics may also be performed, but we have not done so with our tracker because the amplitude of these higher orders are almost near the uncertainty achievable with our method. We should point out that downloading the error map from trackers, correcting it, and uploading the new map, while not very challenging, may require special support from the manufacturer of the tracker.

Tracker B

Tracker B is also an older model from the same manufacturer as tracker A. This tracker also displayed the second order error characteristic as described earlier for tracker A, but with smaller amplitude of about 2.5 [micro]rad.

Tracker C

Tracker C is an older model from a second manufacturer. Fig. 8 shows the difference between front and back face azimuth measurements to a target as a function of azimuth. Low odd orders, if present, are small in magnitude. As mentioned earlier, we do not compensate for odd orders; we simply average front and back face measurements and remove any odd orders.

[FIGURE 8 OMITTED]

Figure 9 shows the results from length tests performed for the two combinations discussed earlier. The dominant harmonic appears to be the second order. A least-squares best-fit yielded an amplitude of 5.0 [micro]rad ([square root of [a.sub.2.sup.2] + [b.sub.2.sup.2]]) for the second order harmonic for combination 1 (D = 7 m, L = 7.75 m) and 4.8 [micro]rad for combination 2(D = 4.9 m, L = 5.4 m). This agreement of the second order harmonic to within 0.2 [micro]rad is within its uncertainty as shown in the next section.

[FIGURE 9 OMITTED]

Fig. 10 shows the results from the length tests after correcting the tracker's error map for second order errors. There is substantial reduction in errors after correcting the error map; the errors are less than [+ or -]25 [micro]m by removing just the second order alone.

[FIGURE 10 OMITTED]

Tracker D

The amplitude of the second order was about 1 [micro]rad for this old tracker from a third manufacturer, much smaller than that seen for the other two trackers. The horizontal length tests however showed errors of about [+ or -] 40 [micro]m, possibly because of higher order harmonics.

5. Discussion

Low order errors can be removed in several different ways. Odd orders may be removed by having two read-heads positioned 180[degrees] apart on the encoder. Even orders may be removed by having two readheads 90[degrees] apart. The cost and complexity of having additional readheads may be mitigated by a calibration of the encoder through independent means such as by using a high accuracy rotary stage, a polygon with auto-collimator etc. The error mapping technique we have discussed in this paper is another approach for the same.

We have tested new tracker designs from two manufacturers, trackers E and F. The results of the horizontal length tests are shown in Fig. 11 for tracker E that intrinsically compensates for second order errors. The errors are much smaller than those in Fig. 6 and, as the manufacturer expected, no second order error is noticeable. We also tested a new tracker (Tracker F) from another manufacturer; the results are shown in Fig. 12. There is no evidence of second order errors.

[FIGURE 11 OMITTED]

[FIGURE 12 OMITTED]

6. A Note on Higher Order Harmonics

The preceding sections described low order harmonic errors in the scale and a simple technique to estimate those errors. There is an analogous effect at the high frequency end where periodic errors have been observed at a frequency representing the spacing between lines in an encoder. We do not directly address compensating high frequency encoder errors in this paper, but point out that a scheme similar to our proposed approach may be considered for testing purposes where an uncalibrated horizontal length can be measured by a tracker at closely spaced azimuthal positions so that multiple sampling positions are contained between two adjacent encoder lines. Two-face tests may also be performed so that multiple sampling positions are contained between two adjacent encoder lines. Such tests may reveal interpolation errors in the encoder.

7. Uncertainty Considerations

There is an uncertainty associated with each of the lengths that are measured as part of the test described in Sec. 2, and therefore an uncertainty in the parameters [a.sub.2] and [b.sub.2]. We address the issue of determining the uncertainty in these parameters here. Values are calculated for combination 1 (D = 7 m, L = 7.75 m; the range to targets R is 8 m and the half angle [theta] subtended by the length is 29[degrees]).

The tracker's ranging accuracy deteriorates farther away from the tracker. Manufacturers typically specify an MPE (maximum permissible error) for the range (for example, one of the trackers we tested had an MPE specification of 20 [micro]m + 0.8 [micro]m/m). Because the length measurements are always performed so that the target is at the same range (about 8 m) from the tracker, a large portion of the error in the range is expected to be due to systematic sources which produce the same errors in length for all horizontal length measurements. We only consider residuals from the mean of all length measurements. Therefore the error in the range due to systematic sources is not of any consequence, and using the MPE for estimating uncertainty is not appropriate for this type of measurement. On the other hand, the variation in repeated range measurements (i.e., the repeatability) is the noise in the measurement and therefore is the primary source of error that is relevant in our case. We experimentally determine the standard uncertainty due to repeatability in range as 1.5 [micro]m for one of the trackers we tested. Only the component of the ranging uncertainty along the length will actually contribute to a length error. Further, the length is the difference in measurements made at two targets and we average 2 measurements at each target location. Therefore the standard uncertainty in the measured length due to ranging error non-repeatability in the tracker is 1.5 sin[theta] = 0.7 [micro]m.

Manufacturers also sometimes specify an MPE for transverse measurements (for example, one of the trackers we tested had an MPE specification of 36 [micro]m + 6 [micro]m/m for angle measurements). As discussed in the previous paragraph, this MPE incorporates all sources that produce errors in the measured angle, including the second order harmonic scale error in the encoder. Further, this MPE is applicable for measurements made in the front face of the tracker, whereas all our measurements are performed in both faces and averaged. Therefore, using the MPE as a bound on the uncertainty is not appropriate for this type of measurement. Because it is the purpose of our experiments to determine the extent of one source of systematic error--namely the second order harmonic, it is necessary to primarily assess the repeatability of the angle measurement system which introduces noise into the measurement. We experimentally measured the standard uncertainty due to repeatability in the horizontal angle to be 0.0001[degrees] for one of the trackers we tested. At a range of 8 m, considering 2 targets, 2 repeat measurements at each target, and only the component along the length, the standard uncertainty in length due to repeatability in the measured horizontal angle is 8(0.0001[pi]/180)cos [theta] = 12.2 [micro]m.

As long lengths are considered in the experiments, stability of the length (realized between two targets on stands) during the duration of the experiment (about 1 hour) is important. To estimate the stability of the length, we mount an interferometer on one stand and a target on the other, and record the length in intervals of 5 seconds over a one hour period. The standard uncertainty in the measured lengths is about 1.2 [micro]m.

The terms discussed above may be summed in quadrature to yield a standard uncertainty of about 12.2 [micro]m in the length (for the 7.75 m length of combination 1). This uncertainty is the 'noise' and is to be contrasted against the 'signal' which is the error in the length (residual from mean length) due to the second order scale error of the encoder. The length error due to second order harmonic scale error is within [+ or -]65 [micro]m as pointed out earlier.

We perform 18 such length measurements from which we extract 2 parameters ([a.sub.2] and [b.sub.2]--the coefficients of the second order harmonic). A Monte Carlo Simulation yields a standard uncertainty of 0.4 [micro]rad on the parameters [a.sub.2] and [b.sub.2] resulting in a 0.4 [micro]rad uncertainty on an arbitrary point in the error map itself (the uncertainty in the scale error map e([empty set]) at any azimuthal position [empty set], where e([empty set]) = [a.sub.2]cos(2[empty set]) + [b.sub.2]sin(2[empty set]). Other simulations suggest that the standard uncertainty on the parameters remains 0.4 [micro]rad if an attempt is made to estimate both the 2nd and 4th orders from the 18 length measurements, and the standard uncertainty on an arbitrary point in the error map (uncertainty on e([empty set]) = [a.sub.2]cos(2[empty set])+[b.sub.2]sin(2[empty set])+[a.sub.4]cos(4[empty set])+[b.sub.4]sin(4[empty set]) is about 0.5 [micro]rad.

The standard uncertainty in the length for combination 2 (D = 4.9 m, L = 5.4 m; the range to targets R is 5.6 m and the half angle [theta] subtended by the length is 29[degrees]) can be calculated in a manner similar to that described above, and is 9 [micro]m. Again, we distinguish between this 'noise' and the 'signal', where the 'signal' is the error in the length (residual from mean) due to second order harmonic scale error, which is within [+ or -]43 [micro]m. The standard uncertainty in the two parameters calculated from 18 length measurements ([a.sub.2] and [b.sub.2]--the coefficients of the second order harmonic) is about 0.4 [micro]rad and there is also a 0.4 [micro]rad uncertainty on an arbitrary point in the error map itself.

The maximum signal-to-noise ratio for combination l (D = 7m, L = 7.75 m) is 65/12.2 = 5.3, and is marginally better than the maximum signal-to-noise ratio for combination 2 (D = 4.9 m, L = 5.4 m) which is 43/9 = 4.7. This does seem to suggest that a longer length placed farther away is more suitable for this test, but practical limitations such as space constraints will most likely limit the size of the length and its placement. A 7.75 m long length placed 7 m away does appear to be adequate in mapping out second order errors.

8. Summary

A tracker is designed to measure large lengths at long distances with a high degree of accuracy. The 7.75 m long length placed 7 m away that was used in this study is a fairly representative measurement for a tracker. Older trackers from two different manufacturers that we tested produced length errors as large as 65 [micro]m due to second order scale error. The maximum permissible error for this measurement is approximately 120 [micro]m based on a typical manufacturer's specifications, without considering temperature effects. The second order scale error is a significant fraction of the MPE.

The method described in this paper provides a simple technique for estimating scale errors in the horizontal angle encoder of a laser tracker. It does not require expensive instrumentation or even a calibrated artifact. It is quick and easy to perform and can be realized in an environment with reasonable degree of stability.

We did correct the horizontal angle encoder error map of two trackers which were from two different manufacturers. The collections were based on the horizontal length tests as described in this paper and, as expected, subsequent length measurement system tests revealed the absence of any second order terms. The length errors reduced from a maximum of about 65 [micro]m before correction to less than 25 [micro]m after collection. Newer trackers from two different manufacturers did not reveal this second order error; neither did an old tracker from the third manufacture.

Acknowledgements

The authors express their gratitude to Georg Gassner for detailed responses to queries about their technique and to Tyler Estler, John Palmateer, and Ted Vorburger for reviewing the manuscript.

Disclaimer: Commercial equipment and materials may be identified in order to adequately specify certain procedures. In no case does such identification imply recommendation or endorsement by the National Institute of Standards and Technology nor does it imply that the materials or equipment identified are necessarily the best available for the purpose.

9. References

(1) B. Muralikrishnan. D. Sawyer, C. Blackburn. S. Phillips. B. Borchardt. and W. T. Estler. ASME B89.4.19 Performance evaluation tests and geometric misalignments in laser trackers. Journal of Research of the National Institute of Standards and Technology. Volume 114. Number 1. January-February 2009. p. 21-35.

(2) Georg Gassner and Robert Ruland. Laser tracker calibration--Testing the angle measurement system. SLAC National Accelerator Laboratory Publication. SLAC-PUB-13476. December 2008.

(3) J. F. Ouyang. W. L. Liu. Y. G Yan. and D. X. Sun, Angular error calibration of laser tracker system. Advanced Laser Technologies 2005. Ivan A. Shcherbakov. Kexin Xu. Qingyue Wang. Alexander V. Priezzhev. and Vladimir I. Pustovoy, eds.. Proc. of SPIE Vol. 6344.

(4) ASME B89.4.19--2006 Performance Evaluation of Laser-Based Spherical Coordinate Measurement Systems. American Society of Mechanical Engineers, http://www.asme.org

About the authors: Bala Muralikrishnan is a guest researcher, Chris Blackburn is a physical science technician, Dan Sawyer is a mechanical engineer and Steve Phillips is a physicist; all are with the Precision Engineering Division of the NIST Manufacturing Engineering Laboratoiy. Robert Bridges is the Chief Scientist at FARO Technologies Inc.

The National Institute of Standards and Technology is an agency of the U.S. Department of Commerce.

B. Muralikrishnan, C. Blackburn, D. Sawyer, S. Phillips,

Precision Engineering Division, National Institute of Standards and Technology, Gaithersburg, MD 20899

and

R. Bridges

FARO Technologies, Kennett Square, PA 19348

bala.muralikrishnan@nist.gov

cluistopher.blackburn@nist.gov

daniel.sawyer@nist.gov

steven.phillips@nist.gov

robert.bridges@faro.com

Printer friendly Cite/link Email Feedback | |

Author: | Muralikrishnan, B.; Blackburn, C.; Sawyer, D.; Phillips, S.; Bridges, R. |
---|---|

Publication: | Journal of Research of the National Institute of Standards and Technology |

Article Type: | Report |

Geographic Code: | 1USA |

Date: | Sep 1, 2010 |

Words: | 5451 |

Previous Article: | Software assurance using structured assurance case models. |

Next Article: | MEMS Young's modulus and step height measurements with round robin results. |

Topics: |

## Reader Opinion