Α complete processing methodology for 3D monitoring using GNSS receivers

Similar documents
ST 524 NCSU - Fall 2008 One way Analysis of variance Variances not homogeneous

8-node quadrilateral element. Numerical integration

Jones vector & matrices

From Structural Analysis to FEM. Dhiman Basu

A Note on Estimability in Linear Models

The Hyperelastic material is examined in this section.

Grand Canonical Ensemble

Soft k-means Clustering. Comp 135 Machine Learning Computer Science Tufts University. Mixture Models. Mixture of Normals in 1D

External Equivalent. EE 521 Analysis of Power Systems. Chen-Ching Liu, Boeing Distinguished Professor Washington State University

COMPLEX NUMBER PAIRWISE COMPARISON AND COMPLEX NUMBER AHP

CHAPTER 7d. DIFFERENTIATION AND INTEGRATION

September 27, Introduction to Ordinary Differential Equations. ME 501A Seminar in Engineering Analysis Page 1. Outline

Linear Algebra Provides a Basis for Elasticity without Stress or Strain

10/7/14. Mixture Models. Comp 135 Introduction to Machine Learning and Data Mining. Maximum likelihood estimation. Mixture of Normals in 1D

CHAPTER 33: PARTICLE PHYSICS

THREE DIMENSIONAL GEOMETRY MAINTENANCE FOR FORMATION FLYING ON ELLIPTIC ORBITS

Heisenberg Model. Sayed Mohammad Mahdi Sadrnezhaad. Supervisor: Prof. Abdollah Langari

GPC From PeakSimple Data Acquisition

Review - Probabilistic Classification

Radial Cataphoresis in Hg-Ar Fluorescent Lamp Discharges at High Power Density

Fakultät III Univ.-Prof. Dr. Jan Franke-Viebach

Outlier-tolerant parameter estimation

Lecture 14. Relic neutrinos Temperature at neutrino decoupling and today Effective degeneracy factor Neutrino mass limits Saha equation

The Fourier Transform

EEC 686/785 Modeling & Performance Evaluation of Computer Systems. Lecture 12

Discrete Shells Simulation

Analyzing Frequencies

Lecture 23 APPLICATIONS OF FINITE ELEMENT METHOD TO SCALAR TRANSPORT PROBLEMS

Abstract Interpretation: concrete and abstract semantics

FEFF and Related Codes

A NON-LINEAR MODEL FOR STUDYING THE MOTION OF A HUMAN BODY. Piteşti, , Romania 2 Department of Automotive, University of Piteşti

1. Basic Concepts. 1.1 What is Kinematics 1.2 Solution Methods 1.3 Choice of Coordinates 1.4 2D versus 3D. Professor Wan-Suk Yoo

Static/Dynamic Deformation with Finite Element Method. Graphics & Media Lab Seoul National University

EDGE PEDESTAL STRUCTURE AND TRANSPORT INTERPRETATION (In the absence of or in between ELMs)

Chapter 6 Student Lecture Notes 6-1

Phy213: General Physics III 4/10/2008 Chapter 22 Worksheet 1. d = 0.1 m

Observer Bias and Reliability By Xunchi Pu

Electrochemical Equilibrium Electromotive Force. Relation between chemical and electric driving forces

Journal of Theoretical and Applied Information Technology 10 th January Vol. 47 No JATIT & LLS. All rights reserved.

Application of Local Influence Diagnostics to the Linear Logistic Regression Models

JEE-2017 : Advanced Paper 2 Answers and Explanations

From Structural Analysis to Finite Element Method


167 T componnt oftforc on atom B can b drvd as: F B =, E =,K (, ) (.2) wr w av usd 2 = ( ) =2 (.3) T scond drvatv: 2 E = K (, ) = K (1, ) + 3 (.4).2.2

Three-Node Euler-Bernoulli Beam Element Based on Positional FEM

Fakultät III Wirtschaftswissenschaften Univ.-Prof. Dr. Jan Franke-Viebach

Heating of a solid cylinder immersed in an insulated bath. Thermal diffusivity and heat capacity experimental evaluation.

A general N-dimensional vector consists of N values. They can be arranged as a column or a row and can be real or complex.

Group Codes Define Over Dihedral Groups of Small Order

HORIZONTAL IMPEDANCE FUNCTION OF SINGLE PILE IN SOIL LAYER WITH VARIABLE PROPERTIES

Econ107 Applied Econometrics Topic 10: Dummy Dependent Variable (Studenmund, Chapter 13)

Linear Algebra. Definition The inverse of an n by n matrix A is an n by n matrix B where, Properties of Matrix Inverse. Minors and cofactors

A NEW GENERALISATION OF SAM-SOLAI S MULTIVARIATE ADDITIVE GAMMA DISTRIBUTION*

Quantum Mechanics I - Session 4

Transfer Functions. Convenient representation of a linear, dynamic model. A transfer function (TF) relates one input and one output: ( ) system

SCITECH Volume 5, Issue 1 RESEARCH ORGANISATION November 17, 2015

History. One breakthrough paper Cousot & Cousot 77 (?) Abstract Interpretation. Inspired by. Lecture 16. Enthusiastically embraced by the community

Economics 600: August, 2007 Dynamic Part: Problem Set 5. Problems on Differential Equations and Continuous Time Optimization

Properties of ferromagnetic materials, magnetic circuits principle of calculation

??? Dynamic Causal Modelling for M/EEG. Electroencephalography (EEG) Dynamic Causal Modelling. M/EEG analysis at sensor level. time.

COHORT MBA. Exponential function. MATH review (part2) by Lucian Mitroiu. The LOG and EXP functions. Properties: e e. lim.

Questions k 10k 100k 1M Speaker. output

MATCHED FILTER BOUND OPTIMIZATION FOR MULTIUSER DOWNLINK TRANSMIT BEAMFORMING

Physics 256: Lecture 2. Physics

Dealing with quantitative data and problem solving life is a story problem! Attacking Quantitative Problems

Text: WMM, Chapter 5. Sections , ,

Lucas Test is based on Euler s theorem which states that if n is any integer and a is coprime to n, then a φ(n) 1modn.

A Packet Buffer Evaluation Method Exploiting Queueing Theory for Wireless Sensor Networks

Decision-making with Distance-based Operators in Fuzzy Logic Control

EE 570: Location and Navigation: Theory & Practice

Journal of Chemical and Pharmaceutical Research, 2014, 6(7): Research Article

Basic Electrical Engineering for Welding [ ] --- Introduction ---

SENSOR MODEL FOR AIRBORNE CCD LINEAR SCANNERS

Study of Dynamic Aperture for PETRA III Ring K. Balewski, W. Brefeld, W. Decking, Y. Li DESY

Ερωτήσεις και ασκησεις Κεφ. 10 (για μόρια) ΠΑΡΑΔΟΣΗ 29/11/2016. (d)

Folding of Regular CW-Complexes

MCB137: Physical Biology of the Cell Spring 2017 Homework 6: Ligand binding and the MWC model of allostery (Due 3/23/17)

Slide 1. Slide 2. Slide 3 DIGITAL SIGNAL PROCESSING CLASSIFICATION OF SIGNALS

Higher-Order Discrete Calculus Methods

Add sodium hydroxide solution

Relate p and T at equilibrium between two phases. An open system where a new phase may form or a new component can be added

cycle that does not cross any edges (including its own), then it has at least

On Selection of Best Sensitive Logistic Estimator in the Presence of Collinearity

On-Line Safe Path Planning in Unknown Environments

Naresuan University Journal: Science and Technology 2018; (26)1

Function Spaces. a x 3. (Letting x = 1 =)) a(0) + b + c (1) = 0. Row reducing the matrix. b 1. e 4 3. e 9. >: (x = 1 =)) a(0) + b + c (1) = 0

Section 11.6: Directional Derivatives and the Gradient Vector

UNIT 8 TWO-WAY ANOVA WITH m OBSERVATIONS PER CELL

OPTIMAL TOPOLOGY SELECTION OF CONTINUUM STRUCTURES WITH STRESS AND DISPLACEMENT CONSTRAINTS

Study interaction between intensive circularly polarized laser and hydrogen atom using a matrix method

ACOUSTIC WAVE EQUATION. Contents INTRODUCTION BULK MODULUS AND LAMÉ S PARAMETERS

FINITE ELEMENT METHOD II Autumn 2015

[ ] 1+ lim G( s) 1+ s + s G s s G s Kacc SYSTEM PERFORMANCE. Since. Lecture 10: Steady-state Errors. Steady-state Errors. Then

Lecture 08 Multiple View Geometry 2. Prof. Dr. Davide Scaramuzza

Shape sensing of aerospace structures by coupling of isogeometric analysis and inverse finite element method

A Probabilistic Characterization of Simulation Model Uncertainties

Logistic Regression I. HRP 261 2/10/ am

A New Refinement of Jacobi Method for Solution of Linear System Equations AX=b

Math 656 March 10, 2011 Midterm Examination Solutions

APPLICABILITY OF LINEARIZED DUSTY GAS MODEL FOR MULTICOMPONENT DIFFUSION OF GAS MIXTURES IN POROUS SOLIDS. Jelena Markovi and Radovan Omorjan

Transcription:

7-5-5 NATIONA TECHNICA UNIVERSITY OF ATHENS SCHOO OF RURA AND SURVEYING ENGINEERING DEPARTMENT OF TOPOGRAPHY AORATORY OF GENERA GEODESY Α complt procssng mthodology for D montorng usng GNSS rcvrs Gorg Pantazs Assocat Profssor of NTUA gpanta@cntralntuagr Sofa, ulgara 7- May 5 5 6 9 8 Two Masurmnt campagns (I, II) aslns Masurmnt 5 4 7 rlatv statc postonng mthod 4 Only non-trval baslns ar usd for th adustmnt n ordr to nsur dffrnt condtons and to avod bas at th calculatons FIG Workng Wk 5

7-5-5 Mthodology soluton of ach basln For ach campagn Dtrmnaton of th componnts ΔΧ, ΔΥ, ΔΖ nth Cartsan gocntrc rfrnc systm btwn th occupaton ponts and rlabl stmaton of th wghts of ΔX, ΔY andδz lnar quatons ar formd and th unknown gocntrc coordnats X, Y, Z ar calculatd by a last squars adustmnt th full VCV matrx s provdd for ach campagn us of th Excl or Matlab softwar Mthodology Th dffrncs of th coordnats Χ, Υ, Ζ, btwn succssv campagns ar calculatd transformaton of Χ, n a local procton plan Υ, E V δe,n,up Ζ, th full VCV matrx s calculatd accordng to th low ofpropagaton rrors by usng th full VCV matrx of th dsplacmnts V δx,y,z th rror llpss or llpsods of th control ponts ar calculatd for th dtcton of th pont s dsplacmnt N Up FIG Workng Wk 5

7-5-5 st mthod Wghts stmaton Αn obctv stmaton of th rrors Χ Y a prlmnary last squar qual wght adustmnt could b appld accordng to thr th mthod of ndrct obsrvatons or th mthod of condton quatons ndrct obsrvatons,qual wght adustmnt Χ = Χ Χ Υ = Υ Υ Ζ = Ζ Ζ Th numbr of quatons of ach systm s qual to th masurd baslns Z nd mthod Th condton quatons ar formd by ntwork s loops closur by usng th masurd ΔΧ, ΔΥ, ΔΖ as follows Χ + Χ + Χ Y + Y + Y Z + Z + Z k k = k k = k k = Th numbr of quatons n vry systm s qual to th numbr of th unary loops Th obctv rrors of th unknown componnts for ach pont of th ntwork s prsntd by th root of thvaranc for ach on n VCV matrcs Χ Y Z nd mthodmprcal Wghts stmaton Thmsclosurofthunaryloopsofthntwork(mc)sth rror, whch th loop contans for thr partcpatng baslns So a dcnt stmaton of ths rror for ach componnt ΔΧ,ΔΥ,ΔΖ s gvn by th followng quaton Χ mc loopx Y mc Thn th man rrors of th baslns componnts dtrmnaton ar calculatd as follows sthnumbrofthloops Y Z Χ t = Χ m Consdrng that for ach basln, th followng quatons ar vald = + Χ Χm Assumng that = Χ Χ Χ thn Χ Y Y loopy Z t= Y m = + Y Ym = Z = Z Y Χ Χ Y Z mc loop Z Zm Z m t = Z = + Y Y Z Z Z FIG Workng Wk 5

7-5-5 Absolut dsplacmnts calculaton Accordng to th ndrct obsrvatons mthod th followng systm of rgular quatons s formd Χ = Χ Χ Υ = Υ Υ Ζ = Ζ Ζ ) T - Ô Ô x = (A Ρ A) Á Ρ l = N Á Ρ l Th absolut poston changs of ntwork s ponts n btwn two squntal masurmnt campagns (I and II) ar calculatd δ IΙ-Ι Χ Y Z Χ+ Y + = Z+ Χ n Yn Z n I II Th varancs and covarancs of ths changs V δx, Y,Z = VX,Y,Z + VX,Y,Z Th changs ΙΙ Χ Ι,,, of ach pont must b convrtd n Υ, Ζ, an orntd local plan procton, δeast, δnorth and δup n ordr to b mor prcptbl and to dfn thr drctons and thr trnds n rlaton to th arth s surfac Absolut dsplacmnts calculaton th total rotaton matrx S A for th (n-) unknown ponts of th ntwork S A sn ë sn ö cos ë cos ö cos ë = cos ë sn ö sn ë cos ö sn ë cosö sn ö sn ö sn ë + + + cos ë cos ö cos ë + + sn ö + + cos ë + sn ë cosö sn ë + + cos ö sn ö + + sn ë n- n- n- sn ö cosë cosö cos ë n- n- cos ë sn ö n- n- n- sn ë cos ö sn ë n- n- cos ö n- sn ö n- th poston changs of ach pont (,, ) n a local procton plan E N Up Th CV matrx for th componnts E, N, Up ar calculatd accordng to th law of propagaton rrors by usng th approprat J matrx as T δe, N,Up δx,y,z J matrx s formd by th partal drvaton of th prvous quaton J=S A V = J V J FIG Workng Wk 5 4

7-5-5 Absolut dsplacmnts calculaton Contnuously th chang vctor r and ts barng b, n rgard to north, r = (δe ) + (δn ) δε b = arctan I, II δν Th horzontal dsplacmntscould b chckd by applyng at a glanc gnral on dmnson chck f r < ó II z and r I, < ó z thn thr s no horzontal dsplacmnt I, II δe δn f r > ó II λ δe I, and r > ó λ I, II thn thr s horzontal dsplacmnt δn th full chck procdur th absolut rror llps s drawn for ach pont for a spcfc confdnc lvl and th dsplacmnt vctor of ach pont s ovr dsgnd th vrtcal dsplacmnts dtcton f δup < σ I, II z thn s no vrtcal dsplacmnt of th pont, othrws δup pont has a vrtcal dsplacmnt for th slctd confdnc lvl w a total approach of th absolut dsplacmnt s chck could b don by th calculaton of th rror llpsod s axs for ach pont ów ów óu u óu óv óv v postonng vctors Rlatv dsplacmnts calculaton In ordr to calculat th rlatv dsplacmnts btwn two pontsandofthntworkthchang svctors, ÄÍ ÄUp, btwn two squntal masurmnt campagns I and II ar calculatd by usng th quatons Ε,, Ε ÄÍ II I = ( Ε Ε, ) ( Ε Ε) = Ε II É = ( Ν Ν, ) ( Ν Ν ) = Ν Ε Ν δup Ι, ΙΙ II I Ι, ΙΙ Ι, ΙΙ = (Up Up, ) (Up Up ) = δup δup Samchcks FIG Workng Wk 5 5

7-5-5 Dscusson Th nflunc of th full CV matrx us s vry mportant for th dsplacmnts dtrmnaton as t maks dffrnc not only to th magntud of th rror llpss axs but also to th orntaton of ts man axs a δr =llps wth Full CV matrx =llps only wth varancs mm mm N N δr E b Th mscalculaton of th llps lads to wrong conclusons about th dsplacmnts of th control pont, as th dsplacmnt vctor may l accdntally outsd or nsd th llps δr δr E Conclusons ThlackofthfullCVmatrxasoutput th ovrstmatd standard rrors of th baslns soluton as wll as th black box followd procdur, ar th man dsadvantags of th maorty commrcal GNSS softwars whn usd n th D montorng In th advantags of th proposd procssng mthodology ar rgstrd th lnar quatons, whch ar formd, rlas th procdur from approxmatons Th ntr procdur can b carrd out n an asy Excl or Matlab nvronmnt as smpl quatons systms ar solvd thus no spcal softwar dvlopmnt s rqurd FIG Workng Wk 5 6

7-5-5 Conclusons Th wght dfnton proposd tchnqu avods th unralstcally optmstc standard rrors calculaton du tothgnssabltytocollctplthoraofdata Thrby, t nsurs th rlablty of th adustmnt as t llustrats th obctv achvd standard rrors n th orgnal capturd data Th us of spcfc rotaton matrx for ach pont n ordr to calculat thr th absolut and rlatv dsplacmnts accordng to th law of propagaton s rror nsur th corrctnss of th rsults Th full CV matrx formaton allows th accurat rror llps or rror llpsod calculaton, th rght valuaton of th dsplacmnts Conclusons Thcomparsonofthszandthrotatonofthrror llpsswhcharformdbyusngthfullcvmatrxor th dfcnt on prov that thr s a strong possblty to xtract dffrnt conclusons for a pont s dsplacmnt as manly th llps s orntaton s compltly dffrnt Th proposd procssng mthodology allows th total survllanc of th adustmnt s stps, th obctv wghts dfnton and th full CV matrx formaton t s valuatd as ffcnt and rlabl for such a trustabl and srous actvty as th D montorng by usng GNSS rcvrs FIG Workng Wk 5 7

7-5-5 NATIONA TECHNICA UNIVERSITY OF ATHENS SCHOO OF RURA AND SURVEYING ENGINEERING DEPARTMENT OF TOPOGRAPHY AORATORY OF GENERA GEODESY Α complt procssng mthodology for D montorng usng GNSS rcvrs Gorg Pantazs Assocat Profssor of NTUA gpanta@cntralntuagr Sofa, ulgara 7- May 5 FIG Workng Wk 5 8