The Effect of Critical Distance in Digital Levelling

: Critical distance concerns precise digital levelling, which has inaccurate results at a certain sighting distance. The influence of critical distance on a measured height difference has been confirmed by calibrating certain digital levels and their appropriate code devices on a vertical comparator under laboratory conditions. The paper aims to explore the influence of critical distance on height differences obtained by precise digital levels of Leica NA3003 and DNA03 by experimental measurements realised in situ. The processing of the measurement results consisted of defining a random error on a station by using parameter estimation of an error model to specify a partial error on a station dependent on sighting distance. Then the processing phase continues with the finding of the relation between the sighting distance and the dispersion of height differences acquired by digital levelling under terrain conditions. The theoretical part involves the development of levelling accuracy theories that vary over time by view on random and systematic error propagation. The numerical and graphical solution of the experimental measurements involves ordering the height differences into sighting groups according to the sighting distance. The standard deviation computed in each sighting group represents a measure of the dispersion of height differences. Suppose the standard deviation in the sighting group in both independent experimental locations K1 and K2 exceeds twice the total standard deviation. In that case, it is most likely considered to be the influence of the critical distance, which is then compared with values obtained by laboratory calibration of the same digital levels.


Introduction
Digital levelling in engineering geodesy brings quick, precise and reliable measurements, often performed under real industrial or building discomfort conditions.The present experience from diagnostics of transport constructions realised by digital levels has brought increased occurrence of levelling errors in some cases.This problem and knowledge of the existence of critical sighting distance of specific digital levels demonstrated by Gassner and Rulland in [1,2] have started research aimed at finding the probable effect of this anomaly on the precision of height differences measured at various sighting distances.Moreover, mathematical modelling involved in the data processing should bring the view of the levelling error sources on a station and determine the ratio between errors dependent and independent of sighting distance.
Accuracy analysis of precise levelling is a highly developed theme by many researchers, who devoted their research to the calibration of the levelling system.For example, Baričevič and his colleagues [3] devoted their research to the automatic calibration of levelling staff, Kuchmister with his co-authors [4] describe the functionality of a vertical comparator based on an interferometric measuring system built in a laboratory at Wroclaw University, Takalo and Rouhiainen [5] from the Finnish Geodetic Institute describe a calibration system based on two concrete pillars for Zeiss DiNi instruments and a vertical laser system for rod testing and Woschitz and his colleagues describe the vertical comparator developed at Graz University in [6][7][8].Authors such as Rekus [9], Sjoberg [10], Yaprak [11], Zilkovski [12] and others have focused their research on accuracy performance and analysis of systematic and random levelling errors.Though the mentioned authors do not directly devote their research to the influence of critical sighting distance, their approach to the accuracy analysis was inspirative for the presented methodology.A short overview of the references concerning the analysis of digital levelling systems is in the second section of this paper.The difference between the presented research and those that are published in the cited references is in organising the special experimental measurements to find the effect of the critical sighting distance on levelling accuracy by using specific digital levels applied under various atmospheric conditions in various seasons.The reliability of the used methodology was verified by comparing the results from two independent localities.
The critical sighting distance should be caused by the systematic effect on measured height differences, therefore the research described in the paper consists of accuracy analysis of levelling data, which is based on classical theories of accuracy, which over time have adapted to new knowledge and technologies and brought a new representation of random and systematic components.The process of investigating the influence of the critical sighting distance consisted of an accuracy analysis of the height differences in the fixed points stabilised in situ on the concrete structure.The particular points have been measured by precise digital levels NA3003 and DNA3, each in four measurement cycles.The experimental measurements were organized into two independent levelling routes with a length of 0.500 km, while the length of intermediate sight on fixed points varied from 2 m to 35 m.The paper deals with a levelling error analysis to detect the partial errors that appeared as a sum of elementary errors accumulated on each station to differentiate errors dependent on the sighting distance.The total levelling error is composed of elementary errors accumulated in each backsight and foresight of a station.The main characteristics of accuracy respect this fact and involve the number of levelling stations, sighting distances and the measured height differences as main arguments in their formulas.Levelling error analysis is based on the basic formulas recommended by the International Association of Geodesy (IAG), which have been subjected to new views on the propagation of random and systematic errors over time.The current error analysis methods respect digital technologies, which brought new errors and changed the ratio of random and systematic error propagation [13].The well-established method of finding the influence of elementary errors in observed levelling data consists of parameter estimation of a mathematical model created from the components of levelling errors, which are (in)dependent on the sighting distance.The model estimations vary in each station, but the discrepancy in their numerical order should be the same.If not, suspicions of the influence of systematic error on a levelling station appear.Then the weighted standard deviations were computed for each meter of sighting distance, named sighting groups, to find the relationship between the height differences dispersion and the sighting distance to find the possible influence of critical distance for the specific digital level.Critical distance is a well-known term among geodesists and metrologists who know that a high probability of obtaining inaccurate results appears at a certain sighting distance during digital levelling.The existence of this critical distance on specific digital levels has been confirmed by calibrating levelling devices at the vertical comparator.As we have the values of critical distances of Leica levelling devices published in [1,2], we can compare the results of field measurements with theoretical values.

An Overview of Levelling Error Theories
One of the first official theories concerning levelling error propagation was presented by Ch.Lallemand, a French geophysicist and the first president of the International Union of Geodesy and Geophysics.He separated random and systematic influence in levelling measurements and defined international valid tolerances for comparing levelling accuracy in other countries (Bulletin géodésique 1930).In 1912, the Union accepted his theory and registered the official formulas for analysing levelling accuracy.Later, G. A. Rune tried to improve the Lallemand equations by estimating levelling standard deviation as a quadratic difference between the total and random error (Zeitschrift für Vermessungswesen 1930).
The head of the French Institut Géographique National in Paris J. Vignal has brought a new view on systematic errors in levelling that influence measurements up to a certain length, and after exceeding it, their variable influence is similar to the behaviour of random errors.After some objections from a member of the Finnish Academy of Sciences, T. J. Kukkamäki and J. Vignal modified his theories, which were subsequently accepted by the Geodetic Union in Oslo.J. Bohm and J. Sloboda summarised the existing theories of levelling error propagation in their book [14].
Over time, views on accuracy analysis and levelling error propagation have adapted to technological improvements in levelling instruments and their accessories.L. Sjöberg published a professional paper [10] concerning the influence of systematic and random errors on the precision of one-way levelling routes measured by a motorised technique.M. Craymer published a thesis [15] devoted to the analysis of Swiss levelling measurements from the point of systematic effects.P. Vaníček and M. Craymer devoted their research to the rod settlement effect on levelling precision [16].D. Zilkovski applied a "special adjustment" of levelling data based on using preliminary errors in evaluating levelling results that were published in [12].Rekus [9] and Saaranen [17] analysed robotized digital levelling in its practical applications.Yaprack [11] and Cvetkov [18] devoted their research to digital level applications in precise height systems.Besides field measurements, technological progress and automation also affected the calibration process of levelling systems.The authors H. Woschitz and F. Brunner described in [6][7][8] the calibration facilities and outcomes from testing certain digital levels at the Laboratory of Geodetic Metrology of the Graz University of Technology (TUG).The conception of the TUG vertical comparator inspired the construction of the calibration system at the Stanford Linear Accelerator Center (SLAC).The authors Gassner and Ruland presented in [1,2] the hardware and experimental results of the SLAC calibration system used to test levelling invar rods of the levelling devices Leica NA3003, DNA03, and Trimble DiNi12.Baričevič and his colleagues [3] also devoted their research to the calibration levelling staff by presenting a new calibration system.Interesting results concerning digital level calibration brought Takalo [5], Kuchmister [4] and Elhassan et al. [19] with their comparison of indoor and outdoor levelling accuracy of digital and optical levels.

Basic Formulas for an Analysis of Levelling Accuracy
The accuracy analysis of levelling data consists of applying well-known formulas, which have been adapted to new levelling technologies to bring a new perspective on the representation of the random and systematic components.The basic formulas of classical accuracy theories are based on defining the sources of levelling errors, the representation of the components of random and systematic influences and the theory of error propagation.Lallemand differentiated the terms random and systematic error in levelling measurement and estimated the variable systematic errors from the cumulative frequency of two levelling routes.The estimation process was based on a hypothesis that assumes an upward or downward tendency of linear regression, possibly caused by the systematic effect.His theory of accuracy is based on the height differences in bidirectional levelling adjusted to the length of one kilometre of the levelling route.The reliable estimation of the accuracy of a one-kilometre-long levelling route is derived by using the law of error propagation of levelling errors.It is given by the quadratic mean of double levelling, known as unit standard deviation.Lallemand recalculated height differences in two-way levelling on kilometre differences by dividing by root squared of the quantity R, representing the length of the levelling route in km.Then the most used accuracy estimation of a one-kilometrelong levelling route follows from the law of error propagation and is given by the quadratic mean of double levelling.This accuracy characteristic is known as the unit standard deviation σ 0 and is given by the following formula [14]: where ρ is the difference between height differences computed in each direct and reverse levelling section with length R given in kilometres.Over time, the restrictive conditions for the universal use of Lallemand formulas appear.His main opponent was Vignal, whose theory consisted of a determination of a limited length of the levelling route, within which unit standard deviation presented in Formula (1) includes only the effect of random errors and short-periodical variable systematic ones.Other levelling lines that exceed the limited length contain also long-periodic systematic errors, which influence levelling observations with the same sign, so they need not be eliminated by bidirectional levelling.The Vignal theory for accuracy analysis of precise levelling is based on the following statements: • Lallemand's deviations of the cumulative height differences in bidirectional levelling are caused not only by systematic errors but also by random noise.• The influence of systematic errors defined by Lallemand is valid only within a certain distance (limited length) of the levelling route, beyond which they behave as variable systematic errors dispersed around the mean systematic error.

•
Levelling variance is expressed as a root square of the total variance: The quantity represents a random and systematic component with variable behaviour in the frame of limit length and beyond, and it behaves as a constant.The systematic influence can be calculated according to the following equation [20]: where K is a constant equal to 2 or 3, Z is the limit length and R m is the mean distance of the levelling sections.The value ζ can be determined from the levelling closures or using differences in the endpoints of the regression line, estimated from the cumulative differences in the height difference in bidirectional levelling in the levelling sections.
Assuming n stations in the levelling route, the total variance is equal to the following relationship [9,14,21]: Random errors η in the levelling route increase in proportion to the number of stations, and the occurrence of systematic error ξ depends on the square of the number of stations.The behaviour of the random and systematic component varies in the frame of limit length, and beyond, it behaves as a constant.The systematic components can be determined from the levelling closures or using differences in the endpoints of the regression line, estimated from the cumulative height differences in bidirectional levelling.

Analysis of Experimental Data
The experimental part of the research involved field observations organised in eight cycles over two years.Digital levels of Leica NA3003 and Leica DNA03 have been used, each in four measurement cycles under various atmospheric conditions.Data intended for accuracy analysis were generated from the intermediate sights acquired on each station in two independent 500 m long levelling routes named K1 and K2.The levelling process was organised by the method of differential levelling using 10 benchmarks fixed on the concrete footings of power line poles.The particular points were stabilized by iron pins fixed to the concrete surface and the sighting distance varied from cca 2 m to approximately 40 m.The data file consisted of the set of measured height differences obtained at the appropriate sighting distances.The amount of the measured intermediate sights was 1490 in locality K1 and 1509 in locality K2, whereby the number of intermediate sights from one station varied in dependence on the visibility to the particular point.The accuracy analysis adapted to this fact by involving weights in the error model and the apriori standard deviations.

Identification of the Type of Elementary Errors
The process of finding the influence of critical sighting distance starts with the identification of the levelling errors depending on the sighting distance by parameter estimation in the error model and testing the homogeneity of estimated variances by involving statistical hypothesis testing that can help validate the mathematical model used.
Böhm [14] divided levelling errors into four groups.The first group consists only of random influences with different sizes and signs at each station.These random errors can be divided into those that depend on the length of the levelling sight, e.g., reading error, error of spirit level, the inaccurate position of the instrument, etc., and those that depend on the number of stations in a levelling route, e.g., errors arising due to movement of an instrument or levelling staff or scale errors.The second group of levelling errors involves only systematic effects depending on the levelling route direction and is mainly influenced by the oscillation of the Earth's crust.The influence of external conditions belongs to the third group of levelling errors such as the intensity of sunlight and the gravity of the Moon and Sun acting systematically and mainly in long-period measurements.The fourth group of levelling errors depends on the measured height difference.A typical representative of this group is levelling refraction with its systematic effect [11,14].
The current levelling technologies have brought easier manipulation and a reduction in the subjective errors represented mainly by targeting and reading errors.On the other hand, new sources of errors appear resulting from their coding and demodulation principles.Ingensand classifies these errors into four groups according to their impact on levelling [22,23]:

•
Bad illumination caused by various intensities of natural light or inhomogeneous light intensity caused by shadows at the levelling bar.
• Atmospheric influences such as turbulences cause blurred images, and refraction, which causes deviation of the line of sight.

•
Mechanical influences such as vibrations (deviation of the line of sight), settlement of the instrument and bar and bar centring and inclination.

•
Instrumental behaviour such as thermal effects (deviation of the line of sight), interference of code element size and pixels (wrong results at certain distances) and bad compensator function.
In addition to the instrumental errors, natural and personal errors appear during levelling such as curvature, refraction, ground settlement and instability, effects of heat on the instrument, parallax, staff out of plumb, etc. Böhm in [14] summarises information about systematic and random error propagation and creates an error model according to the relation of its arguments to the sighting distance: Components of the partial error δ represent the sum of elementary errors on a station [14], which are independent of the sighting distance (e.g., instrument movement), increase with the root of the sighting distance (e.g., targeting error), increase in proportion to the sighting distance (e.g., non-horizontal sightline) and grow with the square of the sighting distance (e.g., refraction error).Provided these components are random and uncorrelated, the partial variance of a levelling sight is equal to the mean squared error value calculated by the law of error propagation as follows [24]: In case of the existence of mathematical or physical correlations of the error components, it is recommended to use the general law of propagation of levelling errors [25,26]: The difference between both Formulas ( 6) and ( 7) is the addition of the component with mixed covariances σ ij generally considered a measure of the linear relationship of the relevant error components and the corresponding differential components f i , f j .In practice, Formulas ( 6) and ( 7) are simplified to the following form: where σ ∆H 2 represents the accuracy of the measured height difference while unit variance σ 0 2 is computed from Lallemand's Formula (1) and R is the length of a levelling section.If the length of the levelling route is only a few kilometres long, the square root of the variance (8) gives a reliable estimation of the levelling precision.
The estimation process consists of finding out the representation of the types of errors, which can numerically vary from sight to sight, but the discrepancy in their numerical order may indicate a possible occurrence of measurement inhomogeneity caused by internal or external effects.The arguments of the error model (5) were estimated in each observation cycle by applying the least-squares method in the nonlinear regression.The solution of the estimation process consists of a vector of unknown parameters and an appropriate covariance matrix [24,25]: with A as the design matrix, l as a vector of input values, and P as a symmetrical weighted matrix with the particular weights laying on the diagonal of the matrix and equal to the number of measured height differences in a group.The covariance matrix was estimated from the Gauss-Markov model by using the law of error propagation where σ 0 2 is an a-posteriori unit variance and Q(x) is a cofactor matrix of unknown estimations.The estimated coefficients of the error model are shown in Table 1.The discrepancy in the numerical order of estimated parameters led to the verifying and validation of the model by method of variance analysis.The variances estimated in each sighting group were subjected to statistical hypothesis testing to validate the mathematical model.The null hypothesis assumes the homogeneity of measurements expressed by the equality of the variances [27].The Cochran variance test was used for this purpose to verify the tested value C computed as the ratio between the largest variance and the sum of variances according to the equation: The critical value C u was computed for the significance level α = 0.05, the number of data series N and the number of sighting groups n, according to the following formula [27]: Symbol F α represents the quantile of F-distribution with three arguments α/N and the degrees of freedom of used datafiles (N − 1) and sighting groups (n − 1) in a file.The Cochran test is used to identify an outlier file with a variance for which test value C exceeds a critical value C > C u .The arguments of the hypothesis testing are published in Table 2 and point to confirmation of the null hypothesis, which proves the homogeneity of tested data files and refers to verifying the mathematical model used.

Detection of the Influence of Critical Sighting Distance
Calibrations of certain levelling devices performed on a vertical comparator are described in [1,2,8].They have brought results concerning the scale determination, the detection of the possible errors of height deviations measured in the end sections of the staff, the determination of the influence of damaged code elements on the height readings, and the identification of the critical sighting distance of specific digital levels.The scale determination is based on using two separate runs of the levelling staff to detect and compare the edges of all code elements of the staff.The scale value of staff is then determined from a linear regression model.The second published calibration procedure is based on finding the height deviations, which rise by sighting at the end sections of levelling code staffs.It starts with defining the useable area of the staff, which varies from approximately 2.80-2.98m to a 3 m long staff, of both Leica and Trimble systems.According to Woschitz [8], the height reading beyond the useable staff area might be wrong by more than 0.5 mm at the sighting distance of 30 m.The reason might be in an asymmetric pixel image on both ends of a staff or the refraction effect, which appears mainly by sighting the lower parts of the levelling staff.Experimental measurements at the Stanford SLAC have produced the results, published by Gassner et al. [1], which contain formulas for computing the ideal sighting section on the end of the levelling staff to avoid corrupted height differences for Leica DNA03.For the lower end, Gassner suggests the following formula: H lower [mm] = visible code on the staff + 20 + 7 × sighting distance [m] (13) and for the upper end of the levelling staff: where H lower and H upper are the ideal sighting sections on the lower and upper ends of levelling staff equal to 10 millimetres at a sighting distance of 10 m.Both Equations ( 13) and ( 14) are valid for sighting distances of up to 15 m.If the sighting distance is less than 3 m, the correct results are only possible in a usable area of 0.078-1.899m using a two-meter-long levelling staff.The third published calibration procedure concerns damaged code elements' influence on the height readings.Reference measurements realised at the TUG have brought knowledge concerning errors arising from damaged code elements, which depend on the used correlation function, which varies for the different levels.Investigation results of the critical sighting distance obtained on the SLAC vertical comparator described in [1] confirmed the assumptions of Woschitz and Brunner [6] that define the critical distance of Leica NA3000, which is 15 m.According to the obtained results, the critical distance occurs when the size of code lines projected onto the CCD array is exactly the size of one pixel or if multiple code lines are mapped to a whole number of pixels.This means that for Leica DNA03, one code element sized 2.025 mm is projected onto the CCD array with the size of one pixel at a distance of 26.7 m, and for Trimble DiNi12, the code element with a width of 20 mm corresponds to the critical distance of 10.98 m or its multiple.These described calibration outputs inspire our research based on experimental measurements to find the influence of the critical sighting distance of digital levels NA3003 and DNA03.This influence was searched by comparing height differences obtained by certain sighting distances, which were arranged into sighting groups with computed standard deviations, which represent a measure of dispersions of height differences.Standard deviations twice greater than the total standard deviation point out suspicious measurement dispersion, which indicates the existence of gross error possibly caused by critical distance influence.Table 3 shows the measured height differences arranged into sighting groups and the appropriate standard deviations that point to their dispersion in the frame of a sighting group.The total value of standard deviation was computed in the reference system separately for each digital level and each experimental locality and is displayed at the end of Table 3.The mean standard deviations are graphically displayed separately for the levelling device Leica NA3003 in Figure 1 and for the DNA03 in Figure 2.

Discussion
The first part of the paper aimed to verify the dependence of the accuracy of measured height differences on the sighting distance.The measurement and processing methodology of digital levelling was based on the assumption of eliminating the subjective effects on the measurement accuracy and appearance of the new kinds of instrumental errors resulting from their digital technology as was mentioned in [23].This means that the ratio between random and systematic influence should change.The current levelling accuracy characteristics are based on the cumulation of partial random and systematic effects on a station, and therefore the experimental measurements consisted of verifying the representation of elementary errors, divided into four groups according to their dependence on sighting distance [14].The main aim of the paper was to detect the influence of critical sighting distance of specific digital levels on measured height differences and compare it with the results found by Gassner, Ruland and Dix in [1,2] via calibration on a vertical comparator.The experiment was based on measuring height differences in fixed points from different sighting distances with digital levels of Leica NA3000 and Leica DNA03.For this purpose, the standard deviations were calculated in each sighting group created around the average sighting distance and were compared with the total standard deviation that represents tolerance.The double-increased standard deviation that appeared in the same sighting group in independent observational localities K1 and K2 can show the probable influence of critical sighting distance.
The results in Table 3 and Figures 1 and 2 show the increased standard deviations for Leica NA3003 around the sighting distances 12 m, 22 m, 27 m and 30 m. Twice higher standard deviations of the Leica DNA03 level were observed at distances of 21 m, 29, 30 m and 32 m while the values with sighting distances over 35 m were not involved in the analysis because of the increased dispersion of measured values (see Figures 1 and 2).The sighting distances with increased values of standard deviations of both digital levels do not correspond to the critical ones published by Gassner and colleagues [1,2].However, the experimental results require further investigation.The uncertainties corresponding to detected sighting distances either showed a different systematic source, critical sighting distance varies from instrument to instrument or it is necessary to revise the research methodology.The second part of the experiment consisted of an estimation of particular errors on a station represented by error model arguments to find the ratio between dependent and independent errors on a sighting distance.Though the representation of the levelling errors varied from station to station, the common characteristic in their numerical order can be found in the whole levelling route in both experimental localities.The results from the estimation of the mathematical model that is published in Table 1 and

Discussion
The first part of the paper aimed to verify the dependence of the accuracy of measured height differences on the sighting distance.The measurement and processing methodology of digital levelling was based on the assumption of eliminating the subjective effects on the measurement accuracy and appearance of the new kinds of instrumental errors resulting from their digital technology as was mentioned in [23].This means that the ratio between random and systematic influence should change.The current levelling accuracy characteristics are based on the cumulation of partial random and systematic effects on a station, and therefore the experimental measurements consisted of verifying the representation of elementary errors, divided into four groups according to their dependence on sighting distance [14].The main aim of the paper was to detect the influence of critical sighting distance of specific digital levels on measured height differences and compare it with the results found by Gassner, Ruland and Dix in [1,2] via calibration on a vertical comparator.The experiment was based on measuring height differences in fixed points from different sighting distances with digital levels of Leica NA3000 and Leica DNA03.For this purpose, the standard deviations were calculated in each sighting group created around the average sighting distance and were compared with the total standard deviation that represents tolerance.The double-increased standard deviation that appeared in the same sighting group in independent observational localities K1 and K2 can show the probable influence of critical sighting distance.
The results in Table 3 and Figures 1 and 2 show the increased standard deviations for Leica NA3003 around the sighting distances 12 m, 22 m, 27 m and 30 m. Twice higher standard deviations of the Leica DNA03 level were observed at distances of 21 m, 29, 30 m and 32 m while the values with sighting distances over 35 m were not involved in the analysis because of the increased dispersion of measured values (see Figures 1 and 2).The sighting distances with increased values of standard deviations of both digital levels do not correspond to the critical ones published by Gassner and colleagues [1,2].However, the experimental results require further investigation.The uncertainties corresponding to detected sighting distances either showed a different systematic source, critical sighting distance varies from instrument to instrument or it is necessary to revise the research methodology.The second part of the experiment consisted of an estimation of particular errors on a station represented by error model arguments to find the ratio between dependent and independent errors on a sighting distance.Though the representation of the levelling errors varied from station to station, the common characteristic in their numerical order can be found in the whole levelling route in both experimental localities.The results from the estimation of the mathematical model that is published in Table 1 and were validated by an analysis of variances show that digital level NA3003 has approximately 49.9% levelling errors that are independent of the sighting distance and DNA03 has approximately 65% dependent errors dependent on the sighting distance.Estimations in Table 1 confirm the well-known expectations, which are also presented in [14], that error arguments c and d have insignificant importance.

Conclusions
The calibration tests of digital levels Leica NA3003 and DNA03 published by researchers Gassner [1,2] introduced the concrete values of critical sighting distances, which caused the increased occurrence of levelling errors when the size of code lines responds to the size of one pixel of the CCD camera or its multiple.The described experiment, which consisted of many intermediate sights measured in two independent levelling routes, did not confirm the influence of a critical distance of 15 m for level NA3003 and 26.7 m for level DNA03 on the precision of measured height differences.However, digital levels showed systematic influences at other sighting distances: 12 m, 22 m, 27 m and 30 m for NA3003 and 21 m, 29 m, 30 m and 32 m for DNA03.The results did not correspond to the multiple code sizes of both digital levels.The reason may be the influence of another systematic error that is assumed to be verified by an independent methodology in future.When analysing the results, it is necessary to take into account that they are compared with findings obtained under significantly more accurate laboratory conditions using highly accurate and specialized calibration equipment.Even if there is an influence of the critical sighting distance on digital levelling, it may not have been detected by the proposed methodology, based on the detection of measurement uncertainty that often appears with an increased value of the standard deviation.In the future, it might be useful to compare results from precise digital and analogue levels in places where the influence of critical sighting distance is expected or to compare precise digital devices of different manufacturers.The analysis of elementary error propagation on a station leads to finding the ratio between errors dependent on and independent of the sighting distance through parameter estimation in the mathematical "error" model.The variances (6) computed from the estimates of error parameters displayed in Table 1 show that digital level NA3003 has approximately 50% of levelling errors that are dependent on the sighting distance and DNA03 has approximately 35%.Because of the precise levelling method, the accuracy analysis concerned only height differences, which do not exceed a sighting distance of 35 m. Results above this limited sighting distance detected much more increased dispersions.
Since the research aimed to compare the results of laboratory tests of specific level-ling devices, obtained by the authors in [1,2], which are naturally much more accurate than field measurements, the methodology of the measurements analysis was based on the assumption that the dispersions of height differences obtained at the critical sighting distance increase much more.The research output pointed to the occurrence of levelling errors in places that do not match the assumptions.Therefore, the task of future experiments will be to compare the results of several digital levelling systems and to focus the methodology of measurement analysis on the evaluation of the height differences found in these suspicious measurement places.The estimation process of finding the proportion of levelling errors depending on the sighting distance has produced an interesting result.The unknown estimations of the mathematical model showed that DNA03 has a larger proportion of errors depending on the sighting distance (65%) than the older instrument NA3003 (50%).

Figure 1 .
Figure 1.The dependence of the levelling accuracy on the mean sighting distance of the level Leica NA3003.

Figure 1 . 13 Figure 2 .
Figure 1.The dependence of the levelling accuracy on the mean sighting distance of the level Leica NA3003.

Figure 2 .
Figure 2. The dependence of levelling accuracy on the mean sighting distance of the level Leica DNA03.

Table 1 .
The regression coefficients estimated from the error model.

Table 2 .
Mathematical model validation by Cochran variance test.

Table 3 .
Mean standard deviations arranged into the sighting groups.

Table 3 .
Mean standard deviations arranged into the sighting groups.