The use of global navigation satellite systems for deformation analysis of the Dargovských hrdinov housing estate

Similar documents
Performance Testing of the Rockwell PLGR+ 96 P/Y Code GPS receiver

Uncertainty in measurements of power and energy on power networks

Calculation of the received voltage due to the radiation from multiple co-frequency sources

ESTIMATION OF DIVERGENCES IN PRECAST CONSTRUCTIONS USING GEODETIC CONTROL NETWORKS

MTBF PREDICTION REPORT

ANNUAL OF NAVIGATION 11/2006

To: Professor Avitabile Date: February 4, 2003 From: Mechanical Student Subject: Experiment #1 Numerical Methods Using Excel

PRACTICAL, COMPUTATION EFFICIENT HIGH-ORDER NEURAL NETWORK FOR ROTATION AND SHIFT INVARIANT PATTERN RECOGNITION. Evgeny Artyomov and Orly Yadid-Pecht

Comparative Analysis of Reuse 1 and 3 in Cellular Network Based On SIR Distribution and Rate

熊本大学学術リポジトリ. Kumamoto University Repositor

Time-frequency Analysis Based State Diagnosis of Transformers Windings under the Short-Circuit Shock

Webinar Series TMIP VISION

Dynamic Optimization. Assignment 1. Sasanka Nagavalli January 29, 2013 Robotics Institute Carnegie Mellon University

A Simple Satellite Exclusion Algorithm for Advanced RAIM

Discussion on How to Express a Regional GPS Solution in the ITRF

Control Chart. Control Chart - history. Process in control. Developed in 1920 s. By Dr. Walter A. Shewhart

Research of Dispatching Method in Elevator Group Control System Based on Fuzzy Neural Network. Yufeng Dai a, Yun Du b

Point Real-Time Kinematic Positioning

Evaluate the Effective of Annular Aperture on the OTF for Fractal Optical Modulator

A Comparison of Two Equivalent Real Formulations for Complex-Valued Linear Systems Part 2: Results

Multichannel Frequency Comparator VCH-315. User Guide

A Preliminary Study on Targets Association Algorithm of Radar and AIS Using BP Neural Network

Methods for Preventing Voltage Collapse

antenna antenna (4.139)

Application of Intelligent Voltage Control System to Korean Power Systems

A study of turbo codes for multilevel modulations in Gaussian and mobile channels

THE ARCHITECTURE OF THE BROADBAND AMPLIFIERS WITHOUT CLASSICAL STAGES WITH A COMMON BASE AND A COMMON EMITTER

A Current Differential Line Protection Using a Synchronous Reference Frame Approach

Generalized Incomplete Trojan-Type Designs with Unequal Cell Sizes

BaselineByCode: An Educational -Purpose Software Package for GPS Baseline Determination Using Code Measurements

POLYTECHNIC UNIVERSITY Electrical Engineering Department. EE SOPHOMORE LABORATORY Experiment 1 Laboratory Energy Sources

Comparison of Two Measurement Devices I. Fundamental Ideas.

Equivalent Circuit Model of Electromagnetic Behaviour of Wire Objects by the Matrix Pencil Method

High Speed ADC Sampling Transients

DETERMINATION OF WIND SPEED PROFILE PARAMETERS IN THE SURFACE LAYER USING A MINI-SODAR

Design of Shunt Active Filter for Harmonic Compensation in a 3 Phase 3 Wire Distribution Network

High Speed, Low Power And Area Efficient Carry-Select Adder

Study of the Improved Location Algorithm Based on Chan and Taylor

State Description of Wireless Channels Using Change-Point Statistical Tests

Guidelines for CCPR and RMO Bilateral Key Comparisons CCPR Working Group on Key Comparison CCPR-G5 October 10 th, 2014

A Novel GNSS Weak Signal Acquisition Using Wavelet Denoising Method

Chapter 29 GPS/GLONASS System Bias Estimation and Application in GPS/GLONASS Combined Positioning

Ultra Sonic Flow Rate Measurements in Large Pipes and open Channels

Digital Transmission

The Stability of Nodal Points in the Leveling Network in the Southwest Part of Israel

Rejection of PSK Interference in DS-SS/PSK System Using Adaptive Transversal Filter with Conditional Response Recalculation

A New Type of Weighted DV-Hop Algorithm Based on Correction Factor in WSNs

Simulation Analysis of GPS/GLONASS Absolute Positioning Performance in an Urban Canyon Environment

THEORY OF YARN STRUCTURE by Prof. Bohuslav Neckář, Textile Department, IIT Delhi, New Delhi. Compression of fibrous assemblies

GPS Attitude Determination Reliability Performance Improvement Using Low Cost Receivers

Technical Literature. SmartRTK: A Novel Method Of Processing Standardised RTCM Network RTK Information For High Precision Positioning

Interoperability Test Analysis between EGNOS and MSAS SBAS Systems

GLONASS Double Difference Ambiguity Resolution in Real-Time

RC Filters TEP Related Topics Principle Equipment

Modeling Power Angle Spectrum and Antenna Pattern Directions in Multipath Propagation Environment

1.0 INTRODUCTION 2.0 CELLULAR POSITIONING WITH DATABASE CORRELATION

Efficient Large Integers Arithmetic by Adopting Squaring and Complement Recoding Techniques

A Relative Positioning Technique with Spatial Constraints for Multiple Targets Based on Sparse Wireless Sensor Network

Strain Gauge Measuring Amplifier BA 660

Simple Models of EMI Filters for Low Frequency Range

NATIONAL RADIO ASTRONOMY OBSERVATORY Green Bank, West Virginia SPECTRAL PROCESSOR MEMO NO. 25. MEMORANDUM February 13, 1985

Beam quality measurements with Shack-Hartmann wavefront sensor and M2-sensor: comparison of two methods

Passive Filters. References: Barbow (pp ), Hayes & Horowitz (pp 32-60), Rizzoni (Chap. 6)

@IJMTER-2015, All rights Reserved 383

Learning Ensembles of Convolutional Neural Networks

Transformer winding modal parameter identification based on poly-reference least-square complex frequency domain method

Walsh Function Based Synthesis Method of PWM Pattern for Full-Bridge Inverter

Priority based Dynamic Multiple Robot Path Planning

Fiber length of pulp and paper by automated optical analyzer using polarized light (Five-year review of T 271 om-12) (no changes since Draft 1)

Sensors for Motion and Position Measurement

A GBAS Testbed to Support New Monitoring Algorithms Development for CAT III Precision Approach

TECHNICAL NOTE TERMINATION FOR POINT- TO-POINT SYSTEMS TN TERMINATON FOR POINT-TO-POINT SYSTEMS. Zo = L C. ω - angular frequency = 2πf

ECE315 / ECE515 Lecture 5 Date:

Satellite Attitude Determination Using GPS Receiver Based on Wahba Cost Function

AFV-P 2U/4U. AC + DC Power Solutions. series. Transient Generation for Disturbance Tests. only. High Performance Programmable AC Power Source

Low Sampling Rate Technology for UHF Partial Discharge Signals Based on Sparse Vector Recovery

PPS Positioning in Weak Signal GPS Environments using a TIDGET Sensor

Performance Analysis of Multi User MIMO System with Block-Diagonalization Precoding Scheme

RECOMMENDATION ITU-R P Multipath propagation and parameterization of its characteristics

Fast Code Detection Using High Speed Time Delay Neural Networks

Impact of Interference Model on Capacity in CDMA Cellular Networks. Robert Akl, D.Sc. Asad Parvez University of North Texas

2142-Presentation. Advanced Conference on Seismic Risk Mitigation and Sustainable Development May Asi-Sisma Presentations

Topology Control for C-RAN Architecture Based on Complex Network

Graph Method for Solving Switched Capacitors Circuits

Performance Analysis of the Weighted Window CFAR Algorithms

1. REVIEW 2. DELIVERY SET

Mesh Equations. Introduction

Parameter Free Iterative Decoding Metrics for Non-Coherent Orthogonal Modulation

Fine Time Aiding in Unsynchronised Cellular Systems: the Benefits for GPS Integration in Mobile Phones

AOA Cooperative Position Localization

Development of an UWB Rescue Radar System - Detection of Survivors Using Fuzzy Reasoning -

Research on Peak-detection Algorithm for High-precision Demodulation System of Fiber Bragg Grating

An Adaptive Over-current Protection Scheme for MV Distribution Networks Including DG

Source Localization by TDOA with Random Sensor Position Errors - Part II: Mobile sensors

IEE Electronics Letters, vol 34, no 17, August 1998, pp ESTIMATING STARTING POINT OF CONDUCTION OF CMOS GATES

RESULTS OF FOUR YEARS MONITORING OF DISPLACEMENTS OF ROCK BLOCKS IN THE KRKONOŠE MTS.

Detection and Mitigation of GPS Spoofing Based on Antenna Array Processing

Total Electron Content (TEC) and Estimation of Positioning Error Using Malaysia Data

STUDY ON LINK-LEVEL SIMULATION IN MULTI- CELL LTE DOWNLINK SYSTEM

MASTER TIMING AND TOF MODULE-

Transcription:

Acta Montanstca Slovaca Ročník 16 (011), číslo 4, 319-37 The use of global navgaton satellte systems for deformaton analyss of the Dargovských hrdnov housng estate Pavel Kukučka 1, Slavomír abant and Gabrel Wess In ths paper, the deformaton analyss of the Dargovských hrdnov housng estate usng global navgaton satellte systems s descrbed. Theoretcally, t deals wth the reasons of formaton of slope deformatons, the ways of ther measurement, and wth the recent development and avalablty of GNSS. Partcular attenton s pad to the descrpton of the deformaton analyss technque, begnnng wth workng out the project, locatng by the GNSS, up to processng and testng of the deformaton network. The practcal part s based on the theoretcal knowledge and descrbes the measurement procedure, the equpment used, workng methodology and testng of the bondng network. Based on the results obtaned by the deformaton analyss, the paper offers nformaton about the recent condton of the nvestgated area. Key words: deformatonal analyss, global navgaton satellte systems, GMM Introducton The global navgaton satellte systems (GNSS) are ncreasngly used to montor landslde areas for the needs of deformaton analyss. The advantage of these measurements s manly n ther hgh relablty and accuracy, as well as ther ndependence from the season and tme of measurements. The man advantage s n elmnaton of the necessty of drect vsblty between the ponts whch s,however, replaced by the necessty of satsfactory recepton of the sgnal from GNSS satelltes. The only condton s drect access of satellte sgnals at the measured ponts. The hghest accuracy of usng GNSS s acheved by relatve measurements of the phase of a carrer wave. Statc relatve-phase measurement s a basc and most often used method n geodesy. Over the last decade, GNSS have essentally ncreased ther accuracy that s a precondton for ther use also n the feld of engneerng geodesy for determnng shfts of buldng structures and slope deformatons. Slope deformatons are processes n whch, for dfferent reasons, the stablty of rock slopes s affected resultng n the sldng movement of rock masses. They pose a danger to all the planned and accomplshed structures whether we are talkng about people s dwellngs or buldngs whch make our lfe on ths planet easer. A slope s prone to slde due to ts geologcal structure, propertes of rocks, hydrogeologcal condtons, etc. The cause of slope movements s gravtaton. The most mportant passve force s the compactness of rock. At present 4 GNSS systems are n use or beng developed (www1): Global Postonng System (GPS NAVSTAR), Global navgaton Satellte System (GONASS), Galleo - a system beng bult by the EU wth global coverage by the year 014, Compass - an ndependent system beng developed by Chna wth global coverage by the year 017. In Slovaka, there are networks n operaton whch are used for spatal determnng of poston: SKPOS and eca SmartNet. SKPOS s a permanent servce consstng of the network of reference statons that processes and n real tme provdes geocentrc coordnates for precse locaton of objects and phenomena. The reference statons recevng GNSS sgnals are postoned n geodetc ponts. The ponts are manly placed on the roofs of the buldngs of and regstres whch are connected to the department vrtual prvate network (VPS) to provde hgh qualty electronc communcaton. In the year 010, other 3 reference statons SKVT (Vranov nad Topľou), BREZ (Brezno), JABO (Jaslovské Bohunce) became avalable, thus ncreasng the total number of reference statons to 6. (www, www3) eca SmartNet s an ntegrated software package that enables to centrally control and operate ndvdual reference statons or the whole networks of reference statons. Ths network s a modular system wth new advanced solutons for hgh-accuracy RTK networks wth extensve coverage, centralzed dstrbuton of data and admnstraton of data access. At present, the network ncludes 0 reference statons. In comparson wth SKPOS that uses the technology of Vrtual reference statons (VRS), the eca SmartNet network uses MAC 1 Ing. Pavel Kukučka, Ing. Slavomír abant, PhD., prof. Ing. Gabrel Wess, PhD., Techncal unverzty of Košce, Faculty of mnng, ecology, procces control, and geotechnology, Insttute of Geodesy, Cartography and Geographc Informaton Systems, Park Komenského 19, 040 01 Košce, Slovak Republc, pavel.kukucka@tuke.sk, slavomr.labant@tuke.sk, gabrel.wess@tuke.sk. 319

Pavel Kukučka, Slavomír abant and Gabrel Wess: The use of global navgaton satellte systems for deformaton analyss of the Dargovských hrdnov housng estate (Master Auxlary Concept) whch s a newer technology. Moreover, MAC network soluton provdes homogenous accuracy ndependent from the dstance from the reference staton. (www4). ocaton of deformaton The Dargovských hrdnov housng estate, also known as Furča s part of the town of Košce located n the dstrct Košce III and s an ndependent cadastral area wth acreage of 1109 hectares. The housng estate stretches n the Košce basn on the top of a major morphologcal elevaton between the valley of the Rver Hornád n the west and the Torysa Rver n the east. The area n queston (Fg. 1) s demarcated by the roads Prešovská cesta, Sečovská cesta, the streets Treda armádneho generála. Svobodu and Herlanska. It s stuated n the Košce basn n the sub-unt of the Toryska hlly land and les 10 80 m above the sea level. ocated n the area there can be found the Opátsky fault (of the Zemplnsko- Berehovsky system) that runs along the foothll of the western front slopes as well as of the Hornad system,.e. western Kosce fault and Kostoľanský fault. The area les n the warm clmatc zone T3 wth average annual temperatures varyng from + 8.5 to 8.9 C and wth on average 118 frost days per year, and average depth of sol freezng 111 cm. Accordng to SHMU measurements, the average long-term total amount of ranfall s 580 660 mm. Fg. 1. Area map. Measurement of the network and evaluaton of sgnals Measurements were performed usng two types of doublefrequency recevers eca GPS900CS, eca GPS100 and a dgtal levellng devce Topcon D-101C. eca GPS100 (Fg. ) and eca GPS900CS (Fg. 3) are double-frequency geodetc RTK recevers that allow to use statc, fast statc, knematc and On the fly methods for measurements. Topcon D-101C allows to perform the levellng measurement usng the method of techncal, precse and very accurate levellng as well as layng-out. Pror to the measurement tself, thorough reconnassance of the terran was carred out n order to nvestgate the condton of the pont feld and select sutable observed ponts (herenafter OP) n the rsky areas. After the reconnassance, t was found out that n the area there were a lot of ponts servng for assessment of slope movements, however, some of these ponts were reconstructed and some of them were destroyed. All the ponts were stablzed as nterconnected geodetc ponts (lned drlled wells ftted wth a metal ppe and flled wth concrete). On the surface of each pllar there was cemented a metal plate wth a drlled hole 16 mm n dameter whch provded dependent centrng and the benchmark was ftted at the lower part of the pllars (Fg. 3). The choce of sutable OP n the rsky areas was dependent on the densty of the bult-up area and vegetaton cover. 5 OP marked B -6, B-10, C-0, P-III-1, P-IV-1 (Fg. 4) were selected for montorng slope deformatons stablzed wth the lned drlled wells where there was an ncreased antcpaton of movement. Selecton of these ponts was performed takng nto consderaton the opton of the best possble recepton of sgnals from the satelltes. Harngeš, Fg.. Measurng at the pont Harngeš. Fg. 3. Measurng at the pont P-IV-1. 30

Acta Montanstca Slovaca Ročník 16 (011), číslo 4, 319-37 Varkapa and Šroká ponts were used for processng as reference ponts (herenafter marked as H, V, S) stablzed by naled sgns wth holes n stone prsms of the sze 0 x 0 x 70 cm protected by concrete rngs and protectve rods (Fg. ). (Sabova, 007, AMS) Before the expermental measurement, whch was performed on 13 June 009, a sutable method of measurement was selected to ensure the hghest accuracy of observaton. The selected method approprate for the ntended measurement was a statc method wth 3 peces of eca GPS100 equpment and 5 peces of eca GPS900CS equpment (+1 back- up) avalable to carry out the mea- these types of equpment, apart from ther techncal parameters, les manly n the fact that eca GPS100 receves only GPS sgnals and eca GPS900CS allows to recevee sgnals from the GPS satellte system as well as surement. The dfference between from GONASS (GNSS recever). That s because of ths very dfference that the posto- was nng of the equpment reassessed wth the concluson that GPS (eca GPS100) recevers should be nstalled on the reference ponts (here- was not blocked and a smaller number of observaton satelltes would be requred. GNSS (eca GPS900CS) recevers were pos- toned on the OP where a clear vew of the sky s sometmes lmted and t was necessary to ncrease the number of obser- vaton satelltes. In ths case these were GONASS satelltes that were used. nafter RP) H, V, S, where the sky Fg. 4. Stuaton of the deformaton networks. Before the measurement t was necessary to set the same mode of measurement for every pece of equpment to avod dstorton of the results of the measurements due to the wrong settng. eca GPS900CSS equpment was postoned on all of the OP usng a Zess pad and two types of fxng bolts. A clamp bolt was used for grppng of the Zess pad to the pllar thus provdng vertcal centrng; the second one was used for grppng of the antenna n the Zess pad. evellng was performed usng pre-rectfed optcal centrng devce mbedded nto the Zess pad. Recevers eca GPS100 were postoned on the RP on trpods wth three telescopc legs usng the Zess pad and a specal bolt. A rectfed optcal centrng devce was used for levellng the devces. Durng the observaton, the heght dfference of the recevers antennas and benchmarks of the ponts was determned by the levellng method usng a dgtal levellng devce Topcon D 101C. For each measurement the horzontal poston was adjusted so that the horzontal adjustment lne n the telescope of the levellng devce was the same as the heght determned on the antenna (Fg.5) and consequently the heght was electroncally read off several tmes for more precse adjustment on the code levellng lath placed on the benchmark of the pont. The resultng heght was calculated on the bass of the arthmetcal average. The control measurement of the slant dstance was performed by applyng an addtonal measurng tape to compare or eventually detectt gross errors. 4-hour observatons were recorded at all ponts Fg. 5. Feld of vew by the measured heght of wth the tme of measurement from 11:00-15:00. antenna. 31

Pavel Kukučka, Slavomír abant and Gabrel Wess: The use of global navgaton satellte systems for deformaton analyss of the Dargovských hrdnov housng estate Processng of the measured data was carred out usng the company software eca Geo Offce 5.0 (GO) where the fact was taken nto account that n the area of the town of Košce there were RP of two ndependent from each other permanent GNSS servces, namely KUKE reference staton of SKPOS servve placed on the buldng of the and regstry n Košce and KOSI reference staton of eca SmartNet network placed on Delluss pavlon on the campus of the Techncal Unversty. The data from the reference statons for the gven tme of measurement were receved from Ms Šalátová (data from KUKE) and from Mr E. Frohmanna, PhD (data from KOSI) upon wrtten requests n the exchange format RINEX. In order to compare the pont poston determned by SKPOS and eca SmartNet servce, n the software processng t was always only one pont was chosen as a RP from the par of KUKE and KOSI ponts and the other was not used so that two ndependent sets of the coordnates of the OP n the coordnate system Unfed Trgonometrc Cadastral Network 03 (UTCN 03) were obtaned. By comparson of the ndvdual sets of coordnates t was proved that the data provded by the state permanent servce of SKPOS observaton statons and the prvate eca SmartNet network were at the same accuracy level (Tab. 1). For the purposes of the deformaton analyss performed ndependent processng of the bondng network was performed whch corresponded to processng n the year 003 when the poston of the ponts was carred out usng the GNSS method and the heght was determned usng the hght precson levellng (HP) method. In the year 003 the ponts H, V, S were used as RP and because of that for the year 009 the KOS and KUKE ponts were excluded from processng. In the year 003, the HP method was used to determne the heghts of the ponts usng the levellng devce Zess N 007 wth an average klometer error 0,5 mm/km and n the year 009 determnaton of the heghts was carred out usng the GNSS statc method wth the vertcal poston determnaton 10 mm + 0.5 D ppm. As the yearly measurements dd not have the same precson level, and the ntal comparson of heghts by the GNSS method wth the year 003 showed consderable dfferences, processng of the 3D network was cancelled, n whch sgnfcant heght dfferences negatvely affected the overall results, and only postonng processng of the bondng network was carred out. Intal pre-processng of the observaton data was performed n GO 5.0 software envronment where, for the sake of updatng, transformaton of the coordnates of the OP and RP from the coordnate system ETRS 89 nto the local coordnate system UTCN 03 (Unfed Trgonometrc Cadastral Network 03) was carred out usng the transformaton parameters UTCN 03. Transformaton of the year 003 used for the deformaton analyss was performed usng smlar transformaton (Šütt, 1997), where on the bass on the respectve formulas the accuracy of the transformaton of ndvdual coordnates of the OP was determned wth an average value of the mean error of the transformaton 17.68 mm. Gauss-Markov model wth full rank Postonng coordnates of 5 OP and 3 RP as well coordnate dfferences between them (Fg. 5) were used for D bondng processng of the estmates of unknown parameters usng Gauss-Markov model (GMM) gven by: v = AdCˆ d = A( Cˆ C ) ( ), Σ = s Q, 0 Tab. 1. Dfferences of coordnates of ponts. KUKE KOSI SKPOS SmartNet [m] [m] [mm] number of ponts - functonal part, - stochastc part, coordnate dfferences The structure of ndvdual values n ths network s defned by m = 5 GNSS vectors obtaned through observatons, n = m = 50 observaton components, b = 8 OP and k = b = 16 determned by parameters, wth k = 10 beng only for the OP. (Gašnec et al., 000; Pukanská et al., 007, Sabová et al., 007; Wess et al., 008): Vector of observatons (50,1) s created by poston coordnate dfferences between ndvdual network ponts whch are taken from the processng of vector n GO 5.0 software. The sequence of coordnate dfferences s gven n alphabetcal order of vector startng ponts based frstly on observed ponts and then on RP wth coordnate dfferences X, Y. Ths order s mantaned for the entre adjustment. B-6 B-10 C-0 P-III- 1 P-IV- 1 X 138170.736 138170.736-0.01 Y 61135.93 61135.93 0.04 Z 45.954 45.954 0.00 X 13886.019 13886.019 0.05 Y 60850.333 60850.333-0.0 Z 80.35 80.35 0.00 X 138054.799 138054.799-0.01 Y 61450.470 61450.470 0.04 Z 18.45 18.45 0.00 X 137837.38 137837.38-0.03 Y 6138.307 6138.307-0.03 Z 60.160 60.160 0.00 X 13841.35 13841.35 0.01 Y 61173.400 61173.400-0.0 Z 65.198 65.198 0.00 (1) 3

Acta Montanstca Slovaca Ročník 16 (011), číslo 4, 319-37 A vector of approxmate coordnates of OP C (10,1) s created by poston coordnates of OP, whch were taken from ntal processng n GO 5.0 software. A vector of approxmate observatons (50,1) s gven by relatve coordnate dfferences of approxmate coordnates of gven ponts = f ( C ). The order of ndvdual dfferences s the same as n vector of observatons. The vector of observaton addtves d (50,1) s gven by the dfference of the elements of the observaton vector and the vector of approxmate observatons : d = and s determned n mlmetres. Cofactor matrx Q (50,50) s a dagonal matrx wth cofactors on the man dagonal. Indvdual cofactors σ were calculated from the relatonshp Σ X, Y q =, where 0 = σ σ, σ = ( 5mm + 0.5 ppm) σ 0 n and the accuracy of the determnaton of the length of ndvdual vectors usng the statc method quoted by the manufacturer s 5 mm + 0,5 mm/km of the length. Confguraton matrx A (partal dervatons, desgn) characterses network geometry (confguraton) of the network. As there were 3 RP n the gven network, the matrx A (50,16) was dvded nto actve part A (50,10), whch was further processed and passve part A (50,6), whch belonged to 3 RP and was defned as: f ( C ) A =. If vectors and C are related to each other, the coeffcents wll have values { 1,1} C C = Cˆ and f not, they wll have value {} 0. Vector of estmates of auxlares of determned coordnates dc ˆ (10,1 n mm was determned usng a matrx multplcaton: dcˆ 1 1 = ( A Q A) A Q ( ) = N A Q d. () The vector of adjusted coordnates of the OP Ĉ (10,1) s determned by addton of the vector o of approxmate coordnates C and the vector of estmates of coordnate addtves for the calculated ponts dc ˆ : Cˆ = C + dcˆ. (3) Vector of correctons v (50,1) of observed values: v = AdCˆ d. (4) Vector of adjusted measured values ˆ (50,1): = + v (5) v Q v Estmated varance faktor s 0 (1,1) : s0 =. ( n k) (6) The covarance matrx Σ Ĉ (10,10) of adjusted coordnates X, Y Σ s Q Cˆ = 0 C ˆ (7) where on the man dagonal there are varances showng the accuracy of the estmates of the calculated coordnates. Covarance matrx of the adjusted values of the observed quanttes s Σ ) (50,50) : Σ ) = s Q) 0 (8) Redundancy matrx R (50,50) s the product of cofactor matrx 1 1 R Q Q = ( Q AN A ) Q. (9) = V Q and cofactor matrx of correctons Q : Tab. 3. Coordnates of OP and ther accuracy. V Table shows measured, approxmate and adjusted numercal data that relate to the measured values havng vector sze (50,1) gvng unts of quanttes. The total of the elements on the man dagonal of the redundancy matrx fulflls the equaton gven by the relatonshp to the redundances of the network: n tr( R) = r + r +... + r = r = r = n k = 50 10 = 40. 1 n = 1 1 Table 3 contans approxmate and adjusted numercal data that relate to the determned coordnates of the sze of vectors (10,1) gven wth the unts of ndvdual quanttes. Ponts coordnates o C (10,1) dc ˆ Ĉ σ (10,1) (10,1) (10,1) Ĉ [ m ] [ mm ] [ m ] [ mm ] X B-6 138170.73-0.39 138170.73 0.37 Y B-6 61135.97 0.38 61135.97 0.38 X B-10 13886.015 0.1 13886.015 0.40 Y B-10 60850.3 0.66 60850.33 0.4 X C-0 138054.795-0.14 138054.795 0.38 Y C-0 61450.467-0.39 61450.467 0.38 X P-III-1 137837.31-0.61 137837.30 0.38 Y P-III-1 6138.30 0.11 6138.30 0.38 X P-IV-1 13841.354-1.44 13841.353 0.38 Y P-IV-1 61173.393-0.0 61173.393 0.39 33

Pavel Kukučka, Slavomír abant and Gabrel Wess: The use of global navgaton satellte systems for deformaton analyss of the Dargovských hrdnov housng estate After processng the measured quanttes t was necessary to verfy by statstcal testng f n the vector of the measured values there were any measurements contanng gross errors. Tab.. Measured data and ther accuracy. GNSS vector o (50,1) d man dag. v (50,1) (50,1) (50,1) ˆ σ (50,1) (50,1) ˆ man dag. Q (50,1) from pont to pont X/ Y [ m ] [ m ] [ mm ] [ mm ] [ m ] [ mm ] R (50,1) B-6 B-10 X 691.83 691.83 0.00 0.63 0.51 691.84 0.44 0.71 Y -85.605-85.605 0.50 0.546-0. -85.605 0.43 0.67 B-6 C-0 X -115.936-115.937 1.00 0.511-0.75-115.937 0.41 0.68 Y 314.539 314.540-0.50 0.55-0.7 314.539 0.41 0.70 B-6 P-III-1 X -333.41-333.411-0.50 0.555 0.7-333.411 0.4 0.69 Y 10.375 10.375 0.00 0.508-0.6 10.375 0.40 0.69 B-6 P-IV-1 X 41.6 41.6-0.50 0.536-0.55 41.61 0.4 0.69 Y 37.465 37.466-1.00 0.496 0.61 37.466 0.41 0.67 B-6 H X 655.899 655.899 0.00 1.145 0.39 655.899 0.37 0.7 Y -687.3-687. -1.00 0.63 0.6-687. 0.38 0.70 B-6 S X -5566.74-5566.74 0.00.180 0.39-5566.74 0.37 0.73 Y 4468.873 4468.871.00 1.751 -.38 4468.871 0.38 0.68 B-6 V X 395.189 395.189 0.00 0.568 0.39 395.189 0.37 0.68 Y -568.787-568.786-1.00 1.119 0.6-568.786 0.38 0.67 B-10 C-0 X -807.1-807.0-1.00 0.659 0.74-807.0 0.44 0.68 Y 600.144 600.145-1.00 0.61-0.05 600.144 0.44 0.69 B-10 P-III-1 X -104.694-104.694 0.00 0.709-0.73-104.695 0.45 0.70 Y 387.981 387.980 0.50 0.567-1.04 387.979 0.43 0.69 B-10 P-IV-1 X -449.663-449.661 -.00 0.580 0.44-449.663 0.44 0.71 Y 33.069 33.071 -.00 0.553 1.33 33.070 0.43 0.68 B-10 S X -658.05-658.05-0.50.475 0.38-658.05 0.40 0.88 Y 4754.479 4754.476.50 1.858-3.16 4754.475 0.4 0.78 B-10 V X -96.094-96.094 0.00 0.548-0.1-96.094 0.40 0.94 Y -83.183-83.181-1.50 1.036 0.84-83.18 0.4 0.9 C-0 P-III-1 X -17.474-17.474 0.00 0.53-0.47-17.474 0.4 0.76 Y -1.164-1.165 1.00 0.530-0.50-1.164 0.41 0.88 C-0 P-IV-1 X 357.557 357.559-1.50 0.560 0.0 357.558 0.4 0.94 Y -77.073-77.074 0.50 0.544-0.13-77.074 0.4 0.91 C-0 H X 771.836 771.836 0.00 1.179 0.14 771.836 0.38 0.71 Y -1001.76-1001.76 0.00 0.703 0.39-1001.76 0.38 0.83 C-0 S X -5450.804-5450.805 0.50.133-0.36-5450.805 0.38 0.88 Y 4154.333 4154.331 1.50 1.636-1.11 4154.331 0.38 0.80 C-0 V X 511.16 511.16 0.00 0.593 0.14 511.16 0.38 0.94 Y -883.37-883.36-0.50 1.14 0.89-883.36 0.38 0.91 P-III-1 P-IV-1 X -575.033-575.033 0.50 0.607-1.3-575.034 0.43 0.77 Y 64.910 64.909 1.50 0.501-1.63 64.909 0.41 0.88 P-III-1 H X 989.309 989.310-1.00 1.46 1.61 989.311 0.38 0.89 Y -789.597-789.597 0.00 0.655-0.11-789.597 0.38 0.79 P-III-1 S X -533.331-533.331 0.50.045 0.11-533.330 0.38 0.93 Y 4366.497 4366.496 1.00 1.713-1.11 4366.496 0.38 0.9 P-III-1 V X 78.601 78.600 1.00 0.641-0.39 78.601 0.38 0.78 Y -671.16-671.161-1.00 1.149 0.89-671.161 0.38 0.88 P-IV-1 H X 414.79 414.77.50 1.073-1.06 414.78 0.38 0.87 Y -74.688-74.688 0.50 0.640-0.48-74.688 0.39 0.77 P-IV-1 S X 153.569 153.567.00 0.519-0.56 153.568 0.38 0.73 Y -606.53-606.5-1.00 1.130 1.0-606.5 0.39 0.87 H V X -60.709-60.710 0.50 1.09-0.50-60.710 0.00 1.00 Y -1881.564-1881.564 0.00 0.95 0.00-1881.564 0.00 1.00 S V X 5961.931 5961.931-0.50.346 0.50 5961.931 0.00 1.00 Y -7037.659-7037.657 -.50.830.50-7037.657 0.00 1.00 Testng of the network The adjusted network structure was tested usng varous methods. Sutablty of selecton of GMM used for adjustment was verfed by a global test of an estmatng model and the presence of remote measurements was tested by a Pope Tau method, Krűger's test and testng of dverged values. 34

Acta Montanstca Slovaca Ročník 16 (011), číslo 4, 319-37 s0 In the global test, test statstcs T gven by the calculaton T = ( n k) = 0, 806 was compared wth σ 0 the crtcal value wth dvson χ (ch-square) at the level of sgnfcance α = 0, 05 and grades of freedom f = ( n k) : χ α ( n k) = 1, 484 where n = m = 50 represents the number of observaton components and k = b = 10 represents the number of determned parameters. As the crtcal value was hgher than the testng crteron χ α > T, the test dd not confrm any dfferences between the mathematcal model of processng and the observatons used, therefore t can be consdered as not dstorted and observatons as not contanng gross errors. (Caspary, 1987) Deformaton analyss Assessment of changes n the area n queston n the md perod t and t + 1 was carred out based on changes n the coordnates of the ponts n the gven area. Snce, as the result of carryng out bondng processng of the network the coordnates of the OP - B-6, B-10, C-0, P-III-1 and P-IV-1 n the UTCN system for the year 009 were adjusted, the deformaton analyss was performed. The stablty of the gven ponts n the perod between the years 003 and 009 was nvestgated. Through compa-rson of the adjusted coordnates of the OP n the years 003 and 009 the deformaton vector was determned gven by: Tab. 4. Testng of the networks wth local tests. testng of Pope Tau method Krűger's test from pont to pont dverged values TKRIT TX TY TKRIT TX TY TKRIT TX TY B-6 B-10 0.744 0.358 0.740 0.354 0.17 0.064 B-6 C-0 1.56 0.43 1.66 0.418 0.36 0.077 B-6 P-III-1 0.435 0.440 0.430 0.436 0.078 0.083 B-6 P-IV-1 0.889 1.033 0.887 1.034 0.163 0.196 B-6 H 0.380 0.871 0.376 0.868 0.054 0.157 B-6 S 0.67 1.84 0.64 1.901 0.08 0.4 B-6 V 0.580 0.618 0.575 0.613 0.108 0.088 B-10 C-0 1.063 0.070 1.065 0.069 0.180 0.01 B-10 P-III-1 1.003 1.650 1.003 1.688 0.165 0.99 B-10 P-IV-1 0.696.14 0.69.47 0.1 0.395 B-10 S 0.47.39 0.44.55 0.04 0.88 B-10 V 0.186 0.89 0.184 0.890 0.034 0.130 C-0 P-III-1 3.80 0.776 0.809 4.091 0.77 0.805.010 0.14 0.149 C-0 P-IV-1 0.319 0.0 0.316 0.00 0.057 0.037 C-0 H 0.135 0.51 0.133 0.507 0.019 0.088 C-0 S 0.51 0.894 0.48 0.89 0.07 0.108 C-0 V 0.04 0.846 0.0 0.843 0.037 0.117 P-III-1 P-IV-1 1.983.753.06 3.00 0.358 0.551 P-III-1 H 1.511 0.154 1.537 0.15 0.10 0.07 P-III-1 S 0.081 0.873 0.080 0.870 0.009 0.104 P-III-1 V 0.537 0.868 0.533 0.866 0.096 0.13 P-IV-1 H 1.08 0.679 1.085 0.674 0.159 0.10 P-IV-1 S 0.899 1.008 0.897 1.008 0.173 0.144 H V 0.485 0.000 0.480 0.000 0.077 0.000 S V 0.31 1.463 0.318 1.484 0.034 0.144 t009 003 dcˆ = t009 003 t009 003 t009 003 dxˆ dyˆ ˆ dz = t003 t003 t003 ˆ t X ˆ t Y ˆ t Z 009 009 009 Xˆ Yˆ Zˆ, = 1,,..,5. (10) Table 5 gves the dfferences of the coordnates for all the OP. The values of these dfferences were up to 6 mm only for the pont B-6 and P-VI-1 where the dfference of coordnates n the drecton of the y-axs was 16.84 mm and n the drecton of the x-axs 1.9 mm. Postonal shft (Fg. 6) of these ponts was the bggest for the ponts B-6 and P-VI-1 whch was 17.37 mm and 14.30 mm whereas for other ponts t was up to 4.74 mm. Through subsequent testng by a global test and a local test t was nvestgated whether coordnate dfferences occured through movememts of ponts or arose based on the errors n the measurements and processng. Tab. 5. Coordnates of OP n the perod of 003 and 009. Epoch 009 Epoch 003 coordnate Postonal Bod dfferences shft [m] [m] [mm] [mm] B-6 X 138170.73 138170.736-4.5 Y 61135.97 61135.911 16.84 17.37 B-10 X 13886.015 13886.019-4.10 Y 60850.33 60850.3 0.8 4.19 C-0 X 138054.795 138054.798 -.75 Y 61450.467 61450.463 3.86 4.74 P-III-1 X 137837.30 137837.317 3.43 Y 6138.30 6138.30 0.36 3.45 P-IV-1 X 13841.353 13841.365-1.9 6117 6117-6.13 14.30 35

Pavel Kukučka, Slavomír abant and Gabrel Wess: The use of global navgaton satellte systems for deformaton analyss of the Dargovských hrdnov housng estate The stablty or nstablty of all the ponts was tested by the global test wth t003 o hypothess : ˆ t009 H ˆ 0 C C = 0. The testng statctcs T = 0, 419 was compared wth the crtcal value F - random varable T at the sgnfcance level α = 0,05 and degrees of freedom f = k 10, f = n k = 50 10 : F ( α, f, f ),065. α 1 1 = 1 = It s gven that T < Fα, therefore the dfferences of coordnates can be consdered as a statstcally nsgnfcant stochastc change of the poston of a pont, and because of that all the OP can be consdered as stable. Though the global test dd not prove the self-movement of the ponts, localzaton testng was conducted to examne the stablty of ndvdual ponts. For each pont of the deformaton network the localzaton testng statstcs T was used gven by the relatonshp: ) ) dc Q ) dc dc T = F ( f = k1, f = n k), α s0 k1 where k 1 = s the number of the tested parameters of the partcular pont whch was compared wth the crtcal value T KRIT ( F α ) (Tab. 6) at the sgnfcance level α = 0, 05 and degrees of freedom f = k 1 =, f = n k = 50 10. (Wess et al., 007). From Table 6 t s obvous that n the local testng of ponts n all the cases the crtcal value was bgger than the testng crteron, hence the OP - B-6, B-10, C-0, P-III-1 a P-IV-1 as well the area defned by these ponts could be consdered as stable. Graphc testng usng confdence ellpses (Sabová et al., 00) confrmed the conclu-sons drawn from global and local testng of the stablty of the OP located n the Dargovských hrdnov housng estate (Fg. 6). As the posto-nal vector of the pont shft n the perod of the years 003 and 009 was located nsde the confdence ellpse (Tab. 7), there was no shft and therefore t was about the accumulaton of the survey errors. = (1,,...,5) Tab. 6 ocal testng of ponts ponts T T KRIT B-6 B-10 C-0 P-III-1 P-IV-1 X 1.094 8.11 Y 0.634 8.081 X 0.14 8.05 Y 0.070 8.03 X 0.057 7.997 Y 0.061 7.971 X 0.045 7.947 Y 0.09 7.93 X 0.00 7.901 Y 0.001 7.880 (11) Tab.7. Parameters of the relatve confdence ellpses a [mm] b [mm] σ α [ g ] B-6 0.49 0.45 0.00 B-10 1.36 1.09 0.00 C-0 0.6 0.46 0.00 P-III-1 0.57 0.48 0.00 P-IV-1 0.73 0.58 0.00 Fg. 6 Graphc testng wth confdence ellpses n the perod between 003 and 009. If a pont n the year 009 was outsde the ellpse, then we could state that there was a shft whch would be supported by the statstcal testng. From the conclusons of the deformaton analyss t appears that n the perod between the years 00 and 009 there occurred no shfts of ponts whatsoever and that coordnate dfferences that arose were only due to survey errors whch was proven by the dentfcaton test. 36

Acta Montanstca Slovaca Ročník 16 (011), číslo 4, 319-37 Based on these conclusons, we can consder the area represented by the OP - B-6, B-10, C-0, P-III-1 and P-IV-1 as stable, however t s necessary to further montor ths area. Concluson Montorng of landslde areas s necessary to ensure safety and dentfy rsk factors. When montorng t s crucal to pay partcular attenton to developng the project and the choce of the method of work. It s always essental to consder economc, tme and accuracy crtera. The use of GNSS for the needs of deformaton nvestgaton appeared to be most useful both n terms of tme, relablty and precson. The advantage of ths choce was that determnaton of the vertcal poston of the ponts was less accurate what n case of processng the 3D deformaton network would affect the overall results. For that reason t s more advsable to determne the poston of the surveyed ponts usng the method of very precse levellng. As follows from the results of the deformaton analyss, coordnate dfferences of the surveyed ponts are statstcally sgnfcant whch s proved by the stablty of the area. References Caspary, W. F.: Concepts of network and deformaton analyss. 1 st edton. Kensngthon: School of surveyng The Unversty of New South Wales, 1987. 187 p., ISBN 0-85839-044-. Gašnec, J., Gašncová, S., Rajnak, M.: To the problem of the adjustment of D geodetc networks. In: Current ssues of land surveyng and engneerng geodesy : proceedngs from the nterenatonal scentfc semnar : Herľany, 9 10 October 000, FBERG, TU Košce, pp. 38-4, 000, ISBN 80-7099-595-5. Pukanská, K. and Wess, G.: The accuracy of locaton ponts usng GPS technology. Coal - Ores - Geologcal Survey, 14, 9, pp. 30-35, 007, ISSN 110-7697. Sabová, J., Gašncová, S.: Graphc testng n deformaton montorng. Acta Montanstca Slovaca. Vol. 7, 00, no., pp. 17-130, ISBN 1335-1788. Sabová, J., Jakub, V.: Geodetc deformaton montorng. 1 st edton. Košce: Edtor centre and edtoral offce AMS, F BERG, Techncal Unversty of Košce, 007. 18p., ISBN 978-80-8073-788-7. Sabová, J., Pukanská, K.: Projekt der Deformatonsuntersuchungen. Acta Montanstca Slovaca. Vol. 1, 007, Specal Issue no. 3, pp 516-519, ISBN 1335-1788. Šütt, J.: Devatonless geodetc transformatons. Acta Montanstca Slovaca, Vol., 1997, no. 1, pp. 1-10, ISBN 1335-1788. Tometz,.: Košce - Dargovských hrdnov housng estate, horzontal wells. Regstered n the Archeve of the geologcal survey, TUKE, 004. Wess, G., Gašnec, J., Engel, J., abant, S., Rákay, ml., Š.: Effect ncorrect ponts of the ocal Geodetc Network at results of the adjustment. Acta Montanstca Slovaca, Vol. 13, 008, 4, pp. 485-490, ISBN 1335-1788. Wess, G., Jakub, V.: The test verfcaton of 3D geodetc ponts and ther changes. Acta Montanstca Slovaca. 1, Specal Issue 3, 007, pp. 61-616. ISSN 1335-1788. www1: Global navgaton satelte system [onlne]. [ct. 010.3.5]. Avalable at: http://en.wkpeda.org/wk/global_navgaton_satellte_system www: Geodetc and Cartographc Insttute Bratslava [onlne]. [ct. 010.3.3]. Avalable at: http://www.gku.sk/predmet-cnnost/geodetcke-zaklady www3: Web server SKPOS [onlne]. [ct. 010.7.15]. Avalable at: http://www.skpos.gku.sk/ www4: SmartNet Slovaka [onlne]. [ct. 010.3.0]. Avalable at: http://www.geotech.sk/srv_gtch.htm 37