A Direct Georeferencing Method for Terrestrial Laser Scanning Using GNSS Data and the Vertical Deflection from Global Earth Gravity Models

Total Page:16

File Type:pdf, Size:1020Kb

A Direct Georeferencing Method for Terrestrial Laser Scanning Using GNSS Data and the Vertical Deflection from Global Earth Gravity Models sensors Article A Direct Georeferencing Method for Terrestrial Laser Scanning Using GNSS Data and the Vertical Deflection from Global Earth Gravity Models Edward Osada , Krzysztof So´snica , Andrzej Borkowski *, Magdalena Owczarek-Wesołowska and Anna Gromczak Institute of Geodesy and Geoinformatics, Wrocław University of Environmental and Life Sciences, ul, Grunwaldzka 53, 50-357 Wrocław, Poland; [email protected] (E.O.); [email protected] (K.S.); [email protected] (M.O.-W.); [email protected] (A.G.) * Correspondence: [email protected]; Tel.: +48-71-230-5609 Received: 8 May 2017; Accepted: 21 June 2017; Published: 24 June 2017 Abstract: Terrestrial laser scanning is an efficient technique in providing highly accurate point clouds for various geoscience applications. The point clouds have to be transformed to a well-defined reference frame, such as the global Geodetic Reference System 1980. The transformation to the geocentric coordinate frame is based on estimating seven Helmert parameters using several GNSS (Global Navigation Satellite System) referencing points. This paper proposes a method for direct point cloud georeferencing that provides coordinates in the geocentric frame. The proposed method employs the vertical deflection from an external global Earth gravity model and thus demands a minimum number of GNSS measurements. The proposed method can be helpful when the number of georeferencing GNSS points is limited, for instance in city corridors. It needs only two georeferencing points. The validation of the method in a field test reveals that the differences between the classical georefencing and the proposed method amount at maximum to 7 mm with the standard deviation of 8 mm for all of three coordinate components. The proposed method may serve as an alternative for the laser scanning data georeferencing, especially when the number of GNSS points is insufficient for classical methods. Keywords: terrestrial laser scanning; georeferencing; reference frame transformation; GNSS positioning; spatial directional intersection; Earth gravity model; deflection of the vertical 1. Introduction Terrestrial laser scanning is one of the most efficient methods for providing the highly accurate and dense point clouds that can subsequently be used for measuring Earth surface processes [1], monitoring the deformations caused by natural hazards [2,3] or generating three-dimensional spatial models, e.g., [4–8]. The point clouds collected by a terrestrial laser scanner require georeferencing [9,10], i.e., the collected points have to be transformed to a reference frame that has a well-defined origin, orientation and scale. The reference ellipsoids of the Geodetic Reference System 1980 (GRS80) [11] or the World Geodetic System 1984 (WGS84) are commonly used as the geocentric coordinate frames for georeferencing of spatial data, because, e.g., the GPS coordinates are directly measured and estimated in WGS84. In terrestrial laser scanning, the point cloud is typically transformed into the geocentric coordinate frame using a 7-parameter Helmert transformation [6]. Three coordinate components of the laser scanner’s origin, the scale, and three orientation angles with respect to the geocentric reference frame are needed to perform such a transformation. The set of seven transformation parameters Sensors 2017, 17, 1489; doi:10.3390/s17071489 www.mdpi.com/journal/sensors Sensors 2017, 17, 1489 2 of 12 is determined using the least squares method on a basis of several georeferencing points whose coordinates are measured both by the laser scanner and a global navigation satellite system (GNSS) receiver. Alternatively, the optical head of a laser scanner may be equipped with several mounted GPS antennas [12–14]. Other methods of georeferencing of terrestrial scanner data comprise back-sighting, sensor-driven, and data-driven approaches [10,15] which also rely on control point targets. In some cases, for instance in city corridors, the availability of control points determined by GNSS can be limited due to deteriorated GNSS signal. Therefore, additional terrestrial measurements may be needed in order to ensure the sufficient number of control points in such situations. To avoid this possible additional effort, we propose a method for direct georeferencing of point clouds obtained using terrestrial laser scanning that allows for providing coordinates in the geocentric frame using a minimum number of GNSS measurements. In the proposed method, three coordinate components of the laser scanner’s measurement origin, two components describing the vertical deflection of the scanner’s rotational axis and one angle of the horizontal orientation of the scanner’s measurement reference frame provide a minimum set of parameters that are needed for georeferencing and are entirely sufficient for a transformation to the geocentric reference frame. To apply the proposed method, the scanner must be levelled. The transformation parameters are estimated using the least squares method on a basis of the geocentric coordinates of the laser scanner’s origin, the geocentric coordinates of one GNSS point and the components of the vertical deflection of the scanner’s rotational axis determined in the gravitational force field of the Earth Gravity Model 2008 (EGM2008) [16] which gives the components of axis deflection with respect to the reference ellipsoid GRS80. Furthermore, we show results from a field test that proves correctness of the results obtained using proposed method. The paper is structured as follows: First, we introduce laser scanner geocentric orientation parameters, then an adjustment procedure for determination of the orientation parameters is described. In the following section, a field test of the proposed method is provided. Finally the method is discussed in an extended context and conclusions are drawn. 2. Methods 2.1. The Laser Scanner Geocentric Orientation Parameters The external geocentric orientation parameters of the laser scanner’s measurement frame, assuming the equality of the scale in the both reference frames, reads as follows (see Figure1): • X0, Y0, Z0 geocentric GRS80 coordinates of the origin P of the laser scanner’s measurement frame (x, y, z), • S, x, h orientation angles of the measuring frame (x, y, z) with respect to the external reference frame (X, Y, Z). The x, h orientation angles constitute the components of the laser scanner vertical axis deflection with respect to the normal to the GRS80 reference ellipsoid (Figure1), which is caused by the inhomogeneity of the local gravity field due to the local mass distribution. The directional horizontal angle S is called the instrument horizontal orientation constant. The coordinates x, y, z of the point Q measured by a laser scanner can be obtained by (Figure1): 0 1 0 1 x s cos a sin b B C B C x = @ y A = @ s sin a sin b A (1) z s cos b + i − j where s—spatial distance, a—horizontal direction, b—vertical angle (corrected due to refraction), i—height of the laser scanner above ground point P, j—height of the reflector above measured point Q, j = 0 for the point clouds and j > 0 for the georeferencing points. The origin of the laser scanner coordinate system (x, y, z) is usually defined as a point of the electro-optical center of the scanner or a point of the intersection of the horizontal and vertical rotation axes of the scanner or a zero Sensors 2017, 17, 1489 3 of 12 Sensors 2017, 17, 1489 3 of 12 intersection of the horizontal and vertical rotation axes of the scanner or a zero distance measurement point of the scanner [17,18]. The connection of the computed vertical coordinate z scos j with the distance measurement point of the scanner [17,18]. The connection of the computed vertical coordinate z coordinate according to Formula (1) yields as follows: z z i (see Figure 1). z = s cos b − j with the z coordinate according to Formula (1) yields as follows: z = z + i (see Figure1). z – laser scanner rotation axis Z Q(X, Y, Z) plumb line EGM direction j North direction x s in the z – ellipsoid normal horizontal y (x, y) - plane i x y P(X0, Y0, Z0) The WGS84 Z0 reference φ Y ellipsoid λ X0 Y0 X Figure 1. Six parameters of the laser scanner geocentric orientation: X0, Y0, Z0, S, x, h. Figure 1. Six parameters of the laser scanner geocentric orientation: X0, Y0, Z0, Σ, ξ, η. The coordinates (x, y, z) of a particular point Q measured by a laser scanner are converted into the The coordinates (x, y, z) of a particular point Q measured by a laser scanner are converted into (X, Y, Z) external reference frame according to the formula (e.g., [19,20]): the (X, Y, Z) external reference frame according to the formula (e.g., [19,20]): 0 1 0 1 0 1 X X X0X x x 0 T B Y C = B Y0 C + (R(S) · Q(x, h, j) · P(j,Tl)) B y C (2) @ A Y@ Y0A RΣQ ,,P, y @ A (2) Z Z0 z Z Z0 z or equivalently: or equivalently: T X = X0 + (R(S) · Q(x, h, j) · P(j, l)) x (3) T X X0 RΣQ,,P, x (3) where: 0 1 where: − sin j cos l − sin j sin l cos j B C P(j, l) = @ − sin l cos l 0 A (4) sin cos sin sin cos cos j cos l cos j sin l sin j P, 0 sin cos 01 (4) 1 −h tan j −x cosB cos cos sin sinC Q(x, h, j) = @ h tan j 1 −h A (5) x1 tan h 1 0 1 Q ,, costanS sin1S 0 (5) R(S) = B − sin S cosS 0 C1 (6) @ A 0 0 1 cosΣ sinΣ 0 with j, l being the geodetic latitude and longitude of the laser scanner’s or total station’s reference RΣ sinΣ cosΣ 0 (6) frame origin P.
Recommended publications
  • Fitting Gravimetric Geoid Models to Vertical Deflections
    Journal of Geodesy (submitted) 1 Fitting gravimetric geoid models to vertical deflections 2 3 W.E. Featherstone ( ) 4 Western Australian Centre for Geodesy & The Institute for Geoscience Research, 5 Curtin University of Technology, GPO Box U1987, Perth, WA 6845, Australia 6 Tel: +61 8 9266-2734; Fax: +61 8 9266-2703; Email: [email protected] 7 8 D.D. Lichti 9 Department of Geomatics Engineering & The Centre for Bioengineering Research and 10 Education, Schulich School of Engineering, The University of Calgary, 11 2500 University Drive NW, Calgary, Alberta, T2N 1N4, Canada 12 Tel: +1 403 210-9495; Fax: +1 403 284 1980; Email: [email protected] 13 14 Abstract . Regional gravimetric geoid and quasigeoid models are now commonly fitted 15 to GPS-levelling data, which simultaneously absorbs levelling, GPS and quasi/geoid 16 errors due to their inseparability. We propose that independent vertical deflections are 17 used instead, which are not affected by this inseparability problem. The formulation is 18 set out for geoid slopes and changes in slopes. Application to 1080 astrogeodetic 19 deflections over Australia for the AUSGeoid98 model shows that it is feasible, but the 20 poor quality of the historical astrogeodetic deflections led to some unrealistic values. 21 22 Keywords : Gravimetric geoid errors, vertical deflections, vertical datum errors 23 24 1. Introduction 25 Fitting regional gravimetric geoid or quasigeoid models to GPS-levelling data has 26 become a widespread practice. A principal objection to this is the inseparability of 27 errors among the levelling and local vertical datum (LVD), GPS and gravimetric 28 quasi/geoid model (cf.
    [Show full text]
  • The Deflection of the Vertical, from Bouguer to Vening-Meinesz, and Beyond – the Unsung Hero of Geodesy and Geophysics
    EGU21-596 https://doi.org/10.5194/egusphere-egu21-596 EGU General Assembly 2021 © Author(s) 2021. This work is distributed under the Creative Commons Attribution 4.0 License. The Deflection of the Vertical, from Bouguer to Vening-Meinesz, and Beyond – the unsung hero of geodesy and geophysics Christopher Jekeli Ohio State University, School of Earth Sciences, Division of Geodetic Science, United States of America ([email protected]) When thinking of gravity in geodesy and geophysics, one usually thinks of its magnitude, often referred to a reference field, the normal gravity. It is, after all, the free-air gravity anomaly that plays the significant role in terrestrial data bases that lead to Earth Gravitational Models (such as EGM96 or EGM2008) for a multitude of geodetic and geophysical applications. It is the Bouguer anomaly that geologists and exploration geophysicists use to infer deep crustal density anomalies. Yet, it was also Pierre Bouguer (1698-1758) who, using the measured direction of gravity, was the first to endeavor a determination of Earth’s mean density (to “weigh the Earth”), that is, by observing the deflection of the vertical due to Mount Chimborazo in Ecuador. Bouguer’s results, moreover, sowed initial seeds for the theories of isostasy. With these auspicious beginnings, the deflection of the vertical has been an important, if not illustrious, player in geodetic history that continues to the present day. Neglecting the vertical deflection in fundamental surveying campaigns in the mid to late 18th century (e.g., Lacaille in South Africa and Méchain and Delambre in France) led to errors in the perceived shape of the Earth, as well as its scale that influenced the definition of the length of a meter.
    [Show full text]
  • Toolbox for Analysis of Flexural Isostasy (TAFI)—A MATLAB Toolbox for Modeling Flexural Deformation of the Lithosphere GEOSPHERE; V
    Software Contribution GEOSPHERE Toolbox for Analysis of Flexural Isostasy (TAFI)—A MATLAB toolbox for modeling flexural deformation of the lithosphere GEOSPHERE; v. 13, no. 5 Sumant Jha, Dennis L. Harry, and Derek L. Schutt Department of Geosciences, Colorado State University, 1482 Campus Delivery, Fort Collins, Colorado 80523, USA doi:10.1130/GES01421.1 7 figures; 5 tables; 2 supplemental files ABSTRACT typically involving subsidence in the vicinity of the load (assuming the load represents an addition of weight to the lithosphere) and uplift in surrounding CORRESPONDENCE: sumant .jha@colostate .edu This paper describes a new interactive toolbox for the MATLAB computing areas (the peripheral “bulge”) (Fig. 1). The magnitude and wavelength of the platform called Toolbox for Analysis of Flexural Isostasy (TAFI). TAFI supports flexural deflection depends on the magnitude of the load and the strength of CITATION: Jha, S., Harry, D.L., and Schutt, D.L., two- and three-dimensional (2-D and 3-D) modeling of flexural subsidence the lithosphere, with greater strength resulting in more broadly distributed 2017, Toolbox for Analysis of Flexural Isostasy (TAFI)—A MATLAB toolbox for modeling flexural and uplift of the lithosphere in response to vertical tectonic loading. Flexural (larger-wavelength) deformation. Flexure results in variations in the regional defor ma tion of the lithosphere: Geosphere, v. 13, deformation is approximated as bending of a thin elastic plate overlying an gravity field due to deformation of density interfaces within the lithosphere. no. 5, p. 1555–1565, doi:10.1130/GES01421.1. inviscid fluid asthenosphere. The associated gravity anomaly is calculated by Such density interfaces may include the crust-mantle boundary, the interface summing the anomalies produced by flexure of each density interface within between the sedimentary basin fill and underlying rocks, the seafloor in marine Received 26 August 2016 the lithosphere, using Parker’s algorithm.
    [Show full text]
  • An Analysis of Vertical Deflections Derived from High-Degree Spherical
    Journal of Geodesy (1999) 73: 10 ± 22 An analysis of vertical de¯ections derived from high-degree spherical harmonic models C. Jekeli Department of Civil and Environmental Engineering and Geodetic Science, The Ohio State University, 2070 Neil Ave, Columbus, OH 43210, USA e-mail: [email protected]; Tel.: +1 614 292 7117; Fax: +1 614 292 2957 Received: 9 December 1997 / Accepted: 21 August 1998 Abstract. The theoretical dierences between the Hel- from the obvious limited resolution imposed by the mert de¯ection of the vertical and that computed from a degree of truncation, vertical de¯ections from spherical truncated spherical harmonic series of the gravity ®eld, harmonic models dier, in principle, from their geomet- aside from the limited spectral content in the latter, ric counterpart mostly because of the curvature of the include the curvature of the normal plumb line, the normal plumb line. This is a systematic eect that is permanent tidal eect, and datum origin and orientation usually small, but may be signi®cant at very high ¯ight osets. A numerical comparison between de¯ections altitudes: for example, if the spherical harmonic models derived from spherical harmonic model EGM96 and are used to compensate inertial navigation systems astronomic de¯ections in the conterminous United (INS) for the eects of gravitation. A rigorous compar- States (CONUS) shows that correcting these systematic ison is presented here that identi®es this and all other eects reduces the mean dierences in some areas. dierences, including the question of permanent tide Overall, the mean dierence in CONUS is reduced from eects.
    [Show full text]
  • Catherine L. Mayes, David T. Sandwell, and Lawrence A. Lawver
    PALEOCEANOGRAPHIC MAPPING PROJECT L i t;~- ...1 ''\1 ! ~f; ,.- F PROGRESS REPORT No. 43 - 0888 '··· ; (.c:- F-;.-:.... Tectonic History and New Isochron Chart of the South Pacific by Catherine L. Mayes, David T. Sandwell, and Lawrence A. Lawver University of Texas Institute for Geophysics Technical Report No. 95 TECfONIC HISTORY AND NEW ISOCHR.ON CHART OF THE SOUTH PACIFIC Catherine L. Mayes 1,2 David T. Sandwell2 and Lawrence A. Lawver 2 1 The Department of Geological Sciences and Instititue for Geophysics, P.O. Box 7456, University of Texas at Austin, Austin, Texas 78713-7456 2 Institute for Geophysics, University of Texas at Austin, 8701 Mopac Blvd., Austin, Texas 78751 Submitted to: Journal of Geophysical Research May, 1988 1 ABS1RACf We have developed an internally consistent isochron chart and a tectonic history of the South Pacific using a combination of new satellite altimeter data and shipboard magnetic and bathymetric data. Highly accurate, vertical deflection profiles (1-2 Jlrad), derived from 22 repeat cycles of Geosat altimetry, reveal subtle lineations in the gravity field associated with the South Pacific fracture zones. These fracture zone lineations are correlated with sparse shipboard bathymetric identifications of fracture zones and thus can be used to determine paleo-spreading directions in uncharted areas. The high density of Geosat altimeter profiles reveals previously unknown details in paleo-spreading directions on all of the major plates. Magnetic anomaly identifications and magnetic lineation interpretations from published sources were combined with these fracture zone lineations to produce a tectonic fabric map. Tne tectonic fabric was then used to derive new plate reconstructions for twelve selected times in the Late Cretaceous and Cenozoic.
    [Show full text]
  • Influence of the Earth's Topographic Masses On
    M. Rezo, D. Markovinović, M. Šljivarić Utjecaj Zemljinih topografskih masa na otklone vertikala ISSN 1330-3651 (Print), ISSN 1848-6339 (Online) UDC/UDK 528.42.024:551.43 INFLUENCE OF THE EARTH’S TOPOGRAPHIC MASSES ON VERTICAL DEFLECTION Milan Rezo, Danko Markovinović, Marko Šljivarić Original scientific paper Modern scientific research is largely focused on interdisciplinary aspect of individual technical disciplines in research and scientific approach in discovering new insights. In the field of geo-sciences, there is a wide range of mutual methods of measuring physical parameters needed for the determination of geologic and structural composition of the underground area, the density of surface layers, the Earth’s gravity field, i.e. local gravity changes. This paper explains the influences of topography of the Earth masses in the determination of vertical deflections needed for the transformation of the heights by means of trigonometric levelling methods in the real gravity field acceleration. Keywords: relief correction and residual terrain modelling, topographic isostatic reduction, topographic reduction, vertical deflection Utjecaj Zemljinih topografskih masa na otklone vertikala Izvorni znanstveni članak Suvremena znanstvena istraživanja u velikoj su mjeri usmjerena k interdisciplinarnosti pojednih tehničkih grana u istraživačkom i znanstvenom pristupu iznalaženja novi spoznaja. U području geo-znanosti postoji široki spektar zajedničkih metoda mjerenja različitih fizikalnih parametara radi određivanja geološko-strukturnog sastava
    [Show full text]
  • Hirt C. (2010) Prediction of Vertical Deflections from High-Degree Spherical Harmonic Synthesis and Residual Terrain Model Data
    Citation: Hirt C. (2010) Prediction of vertical deflections from high-degree spherical harmonic synthesis and residual terrain model data. Journal of Geodesy 84(3): 179-190. DOI 10.1007/s00190-009-0354-x. Journal of Geodesy manuscript No. (will be inserted by the editor) Christian Hirt Prediction of vertical deflections from high-degree spherical harmonic synthesis and residual terrain model data Received: date / Accepted: date Abstract This study demonstrates that in mountainous areas the use of residual terrain model (RTM) data significantly improves the accuracy of vertical deflections obtained from high-degree spherical harmonic synthesis. The new Earth gravitational model EGM2008 is used to compute vertical deflections up to a spherical harmonic degree of 2160. RTM data can be constructed as difference between high-resolution SRTM (Shuttle Radar Topography Mission) elevation data and the terrain model DTM2006.0 (a spherical harmonic terrain model that complements EGM2008) providing the long-wavelength reference surface. Because these RTM elevations imply most of C. Hirt Western Australian Centre for Geodesy & The Institute for Geoscience Research Curtin University of Technology GPO Box U 1987 Perth, WA 6845, Australia E-mail: [email protected] E-mail: [email protected] 2 Christian Hirt the gravity field signal beyond spherical harmonic degree of 2160, they can be used to augment EGM2008 vertical deflection predictions in the very high spherical harmonic degrees. In two mountainous test areas - the German and the Swiss Alps - the combined use of EGM2008 and RTM data was successfully tested at 223 stations with high-precision astrogeodetic vertical deflections from recent zenith camera observations (accuracy of about 0.1 arc seconds) available.
    [Show full text]
  • Lunar Deflections of the Vertical and Their Distribution
    Vol.3, No.5, 339-343 (2011) Natural Science http://dx.doi.org/10.4236/ns.2011.35045 Lunar deflections of the vertical and their distribution Jinyun Guo1,2*, Yu Sun1,2, Xiaotao Chang2,3, Fanlin Yang1,2 1College of Geodesy and Geomatics, Shandong University of Science and Technology, Qingdao, China; *Corresponding Author: [email protected] 2Key Laboratory of Surveying and Mapping Technology on Island and Reef of SBSM, Qingdao, China; 3Satellite Surveying and Mapping Application Center of SBSM, Beijing, China. Received 4 February 2011; revised 20 March 2011; accepted 6 April 2011. ABSTRACT interior, from crust to core and further understanding of the thermal evolution of the Moon [1]. Lunar gravity The deflection of the vertical reflects the mass field, which contains abundant information related to distribution and density anomaly of celestial lunar internal mass distribution, plays a key role in orbit bodies. Lunar deflections of the vertical include determination and controlling lunar spacecrafts. directional information of the Moon’s gravity Abundant lunar gravity information can be extracted field. SGM90d, recovered from SELENE mission, from the ground tracking data of lunar satellite and ob- revealed the lunar far side gravity field for the servations from lunar explorers, which are scientific data first time in history owes to 4-way Doppler data. used to refine the lunar gravity field model [2]. Study of Lunar deflections of the vertical and their me- lunar gravity with ground tracking data began with the ridional and prime vertical components are cal- Soviet Union spacecraft Luna 10 lunched in 1966. Many culated from SGM90d, and then their global lunar probing missions have been implemented since distributions are also given in the paper.
    [Show full text]
  • Determination of Astro-Geodetic Vertical Deflections Using Digital Zenith Camera System in Istanbul, Turkey
    Determination of Astro-Geodetic Vertical Deflections using Digital Zenith Camera System in Istanbul, Turkey Kerem HALICIOGLU, Rasim DEN İZ, and Haluk OZENER, Turkey Key words : Digital Zenith cameras, Astro-Geodetic Geoid, Vertical Deflections SUMMARY In this study we introduce the first observations of a Digital Zenith Camera System (DZCS) used in Istanbul, Turkey. In the scope of this study we discuss the astro-geodetic observations performed on a test station using a CCD camera, Schmidt-Cassegrain type telescope with 14 inches aperture, two inclinometers, and a dual frequency geodetic GPS receiver which is used for time measurements and ellipsoidal coordinates of the station with a computer control unit for data capturing and system control. The observations performed with the equipment stated above are based in Bogazici University (BOUN) Kandilli Observatory and Earthquake Research Institute (KOERI), and Istanbul Technical University. The results of the first test observations show the limits of the system designed and they help us to plan future studies. We performed the test observations using a larger aperture telescope and astronomic CCD camera in TÜB İTAK National Observatory (TUG) in order to compare the performance of different type of instruments for astro-geodetic observations. This article presents the preliminary research we conducted in 2011, which guides us to develop a new DZCS for our current astro-geodetic studies that is supported by The Scientific and Technological Research Council of Turkey. TS04B - Heights, Geoid and Gravity, 6017 1/10 Kerem Halicioglu, Rasim Deniz, and Haluk Ozener Determination of Astro-Geodetic Vertical Deflections using Digital Zenith Camera System in Istanbul, Turkey FIG Working Week 2012 Knowing to manage the territory, protect the environment, evaluate the cultural heritage Rome, Italy, 6-10 May 2012 Determination of Astro -Geodetic Vertical Deflections using Digital Zenith Camera System in Istanbul, Turkey Kerem HALICIOGLU, Rasim DEN İZ, and Haluk OZENER, Turkey 1.
    [Show full text]
  • Calculation of Deflection of the Vertical Components: Analyzing the Gps, Levelling Measurements and Their Distribution Geometry
    CALCULATION OF DEFLECTION OF THE VERTICAL COMPONENTS: ANALYZING THE GPS, LEVELLING MEASUREMENTS AND THEIR DISTRIBUTION GEOMETRY HGMP Wijesiriwardhana1, RMSP Rathnayake1, HMI Prasanna1, and DR Welikanna1 1Department of Surveying and Geodesy, Faculty of Geomatics, Sabaragamuwa University of SriLanka Email (Corresponding author): [email protected] KEYWORDS: Deflection of the vertical, Least-squares, GPS-Levelling, A priori variance, EGM2008 ABSTRACT: Precise calculation of survey networks, such as reduction of directions, azimuths, zenith angles, slope distances onto the ellipsoid and geoid modeling processes are topically demanded applications of deflection of vertical. Astro-geodetic and gravimetric techniques, incorporated with GPS data and terrestrial data, are complex and time consuming methods of determining this value. This study resolved the deflection of the vertical components in north-south (휉 ) and east-west (휂 ) of a fixed point by using GPS and levelling measurements, and investigated the effect of distribution geometry of the ancillary stations around the selected station for the calculated results. This process was enhanced by considering both spreaded distribution geometry of the stations and skewed distribution geometry of the stations. The deflection of vertical components were calculated for a fixed point in the area of Sabaragamuwa University of Sri Lanka, using two test networks for both spreaded and skewed distribution geometries with 6 ancillary stations and a test station ranging from 300 to 1200 meters. Geoidal heights and ellipsoidal heights were obtained by using levelling and GPS measurements. General least-squares solution applied to mathematical model that relates the unknown parameters (휉, 휂) with GPS and levelling observables. Deflection of vertical components 휂 and 휉 of the test station in the spreaded distribution network were 휂 = -2.49086695" ± 4.29708426" and 휉 = -25.08863327"± 3.74482408".
    [Show full text]
  • Arxiv:1210.1616V3 [Astro-Ph.IM] 3 Mar 2013
    Technical Report for CTIO, July 2012. Last updated 1 March 2013 Preprint typeset using LATEX style emulateapj v. 08/22/09 ACCURATE GEODETIC COORDINATES FOR OBSERVATORIES ON CERRO TOLOLO AND CERRO PACHON´ Eric E. Mamajek1,2 Technical Report for CTIO, July 2012. Last updated 1 March 2013 ABSTRACT As the 50th anniversary of the Cerro Tololo Inter-American Observatory (CTIO) draws near, the author was surprised to learn that the published latitude and longitude for CTIO in the Astronomical Almanac and iraf observatory database appears to differ from modern GPS-measured geodetic posi- tions by nearly a kilometer. Surely, the position for CTIO could not be in error after five decades? The source of the discrepancy appears to be due to the ∼30" difference between the astronomical and geodetic positions { a systematic effect due to vertical deflection first reported by Harrington, Mintz Blanco, & Blanco (1972). Since the astronomical position is not necessarily the desired quantity for some calculations, and since the number of facilities on Cerro Tololo and neighboring Cerro Pach´on has grown considerably over the years, I decided to measure accurate geodetic positions for all of the observatories and some select landmarks on the two peaks using GPS and Google Earth. Both sets of measurements were inter-compared, and externally compared to a high accuracy geodetic position for a NASA Space Geodesy Program survey monument on Tololo. I conclude that Google Earth can currently be used to determine absolute geodetic positions (i.e. compared to GPS) accurate to roughly ±0.15" (±5 m) in latitude and longitude without correction, or approximately ±0".10 (±3 m) with correction.
    [Show full text]
  • Critical Vertical Deflection of Buried HDPE Pipes
    Critical Vertical Deflection of Buried HDPE Pipes A dissertation presented to the faculty of the Russ College of Engineering and Technology of Ohio University In partial fulfillment of the requirements for the degree Doctor of Philosophy Xiao Han April 2017 © 2017 Xiao Han. All Rights Reserved 2 This dissertation titled Critical Vertical Deflection of Buried HDPE Pipes by XIAO HAN has been approved for the Department of Civil Engineering and the Russ College of Engineering and Technology by Teruhisa Masada Professor of Civil Engineering Dennis Irwin Dean, Russ College of Engineering and Technology 3 ABSTRACT HAN, XIAO, Ph.D., April 2017, Civil Engineering Critical Vertical Deflection of Buried HDPE Pipes Director of Dissertation: Teruhisa Masada Profile-wall high-density polyethylene (HDPE) pipes are used at increasing rates to convey surface drainage under roadways. This is because these products are cost effective, easy to handle and resistant to environmental elements. The HDPE pipes are derived from a thermoplastic material and much more flexible compared to the conventional pipes made of concrete and metals. Because of this inherent flexibility and other concerns that exist with thermoplastic materials, many highway agencies in the U.S. specify a limit on the HDPE pipe’s vertical deflection. For example, the Ohio Department of Transportation (ODOT) currently imposes a threshold vertical deflection of 7.5% on all HDPE drainage pipes installed in the ground. Many other state DOTs have similar vertical deflection limits for HDPE pipes. The origin of this vertical deflection limit is not well understood in the engineering community. This study aims at locating the foundation of the threshold vertical deflection enforced by the highway agencies.
    [Show full text]