Top menu

модул

Search

GeoStatistica

СНПЧ А7 Курск, обзоры принтеров и МФУ

Field GIS (GeoStatistics)

Accuracy measurement of different point area with two receiver Leica-1200 (Real Time Kinematic) and receiver Trimble-Juno

 

 

Table of Contents

Introduction

  • Create .txt file and Downloading data to RStudio
  • Plotting points compare with RTK(3D)
  • An accuracy of averaging and not with RTK point 
  • Finding better code for present data and analyze
  • Boxplot
  • Find the difference between each point and RTK

Conclusion 

 

 

Introduction

 

  RStudio is a free and open-source integrated development environment (IDE) for R, a programming language for statistical computing and graphics. RStudio was founded by JJ Allaire, creator of the programming language ColdFusion. Hadley Wickham is the Chief Scientist at RStudio. RStudio is available in two editions: RStudio Desktop, where the program is run locally as a regular desktop application; and RStudio Server, which allows accessing RStudio using a web browser while it is running on a remote Linux server. Prepackaged distributions of RStudio Desktop are available for Windows, OS X, and Linux.

RStudio is available in open source and commercial editions and runs on the desktop (Windows, OS X, and Linux) or in a browser connected to RStudio Server or RStudio Server Pro (Debian, Ubuntu, Red Hat Linux, CentOS, open SUSE and SLES).

RStudio is written in the C++ programming language and uses the Qt framework for its graphical user interface.

Work on RStudio started at around December 2010 and the first public beta version (v0.92) was officially announced in February 2011

The R language is widely used among statisticians and data miners for

developing statistical software and data analysis. Polls and surveys of data miners are showing R's popularity has increased substantially in recent years. 

 The aim is to get familiar systematic analysis and statistic calculation between two receiver Leica 1200 series with geodetic accuracy and GIS-level GNSS receiver Trimble Juno for measuring point in different environment open area, vegetated are and build up area. Measuring point with GNSS receiver Trimble Juno five days continuously two times on each point because of averaging and non-averaging time between them around 10 seconds first non-averaging then averaging 30 seconds.
After that, we measure with receiver Leica 1200 series the same point one time in real time kinematic observation. The date is with GNSS receiver Trimble Juno from 04.05 till 04.09.2016 and with Leica 1200 04.10.2016.
The main point is that we need by visualized and analyze with RStudio tools the accuracy between that points measuring by two receivers the same place and find root square error such as called
Euclidean distance. In mathematics, the Euclidean distance or Euclidean metric is the "ordinary" (i.e. straight-line) distance between two points in Euclidean space. With this distance, Euclidean space becomes a metric space
The Euclidean distance between points p and q is the length of the line segment connecting them ({p}{q}). 

In Cartesian coordinates, if p = (p1, p2,..., pn) and q = (q1, q2,..., qn) are two points in Euclidean n-space, then the distance (d) from p to q, or from q to p is given by the Pythagorean formula:

1 stat

After correcting my data from two receivers in Excel I had two tables:

Table.1 Points measured by receiver Trimble Juno five days

ID Date Environmen Mean East North Elevation Nsatell
1 04.05.2016 open 0 464622.661 262740.598 242.96 5
2 04.05.2016 open 1 464622.736 262739.155 239.18 5
3 04.05.2016 vegar 1 464374.325 262597.769 255.84 9
4 04.05.2016 vegar 0 464376.801 262599.131 253.56 9
5 04.05.2016 Urb 0 464735.843 262911.108 250.06 4
6 04.05.2016 Urb 1 464731.031 262912.539 257.9 5
7 04.06.2016 vegar 0 464378.648 262598.145 259.76 6
8 04.06.2016 vegar 1 464375.537 262597.033 249.42 5
9 04.06.2016 open 0 464616.794 262741.154 228.86 7
10 04.06.2016 open 1 464617.863 262737.332 234.21 7
11 04.06.2016 Urb 0 464724.057 262922.419 231.06 9
12 04.06.2016 Urb 1 464730.624 262917.898 237.99 9
13 04.07.2016 vegar 0 464374.407 262598.65 246.76 5
14 04.07.2016 vegar 1 464372.627 262603.706 246.78 4
15 04.07.2016 open 1 464642.031 262733.663 265.05 5
16 04.07.2016 open 0 464626.442 262741.405 244.56 5
17 04.07.2016 Urb 1 464725.964 262923.625 235.21 5
18 04.07.2016 Urb 0 464731.299 262921.819 243.96 7
19 04.08.2016 vegar 0 464373.244 262601.468 246.76 7
20 04.08.2016 vegar 1 464368.593 262598.488 241.94 6
21 04.08.2016 open 0 464619.266 262740.149 227.86 8
22 04.08.2016 open 1 464619.015 262737.363 228.35 7
23 04.08.2016 Urb 1 464728.081 262913.763 239.69 8
24 04.08.2016 Urb 0 464733.053 262906.005 245.46 7
25 04.09.2016 Urb 0 464715.624 262921.016 245.86 8
26 04.09.2016 Urb 1 464729.588 262912.962 264.32 7
27 04.09.2016 open 1 464618.285 262738.431 231.79 7
28 04.09.2016 open 0 464619.065 262737.745 231.86 8
29 04.09.2016 vegar 0 464374.706 262596.231 256.46 8
30 04.09.2016 vegar 1 464373.493 262584.28 278.68 7
               

 Table2: Points measured by receiver Leica 1200(RTK)

ID Date Environmen rtkEst rtkN Elev
1 04.05.2016 open 464618.030 262737.426 222.444
2 04.05.2016 open 464618.030 262737.426 222.444
3 04.05.2016 vegar 464373.619 262596.806 231.446
4 04.05.2016 vegar 464373.619 262596.806 231.446
5 04.05.2016 Urb 464723.499 262922.801 219.273
6 04.05.2016 Urb 464723.499 262922.801 219.273
7 04.06.2016 vegar 464373.619 262596.806 231.446
8 04.06.2016 vegar 464373.619 262596.806 231.446
9 04.06.2016 open 464618.030 262737.426 222.444
10 04.06.2016 open 464618.030 262737.426 222.444
11 04.06.2016 Urb 464723.499 262922.801 219.273
12 04.06.2016 Urb 464723.499 262922.801 219.273
13 04.07.2016 vegar 464373.619 262596.806 231.446
14 04.07.2016 vegar 464373.619 262596.806 231.446
15 04.07.2016 open 464618.030 262737.426 222.444
16 04.07.2016 open 464618.030 262737.426 222.444
17 04.07.2016 Urb 464723.499 262922.801 219.273
18 04.07.2016 Urb 464723.499 262922.801 219.273
19 04.08.2016 vegar 464373.619 262596.806 231.446
20 04.08.2016 vegar 464373.619 262596.806 231.446
21 04.08.2016 open 464618.030 262737.426 222.444
22 04.08.2016 open 464618.030 262737.426 222.444
23 04.08.2016 Urb 464723.499 262922.801 219.273
24 04.08.2016 Urb 464723.499 262922.801 219.273
25 04.09.2016 Urb 464723.499 262922.801 219.273
26 04.09.2016 Urb 464723.499 262922.801 219.273
27 04.09.2016 open 464618.030 262737.426 222.444
28 04.09.2016 open 464618.030 262737.426 222.444
29 04.09.2016 vegar 464373.619 262596.806 231.446
30 04.09.2016 vegar 464373.619 262596.806 231.446

Position of points taken by Trimble Juno in 3D

 2 stat

Different height points Trimble Juno

3 stat

Plotting real RTK points taken by Leica 1200

4 stat

 

5 stat

 

 Distribution around RTK point in vegetated  area

6 stat

Number of tracking satellite and PDOP dependency

7 stat

 

Euclidian 3D distance
Environment Open area Buildup Vegetated
Mean(0:1) Aver NonAv Aver NonAv Aver NonAv
04.05.2016 17,47 21,27 35,17 40,67 24,42 22,46
04.06.2016 7,52 11,77 11,9 20,62 28,8 18,08
04.07.2016 49,1 23,99 16,2 25,91 15,4 16,84
04.08.2016 5,99 6,19 32,5 22,79 16 11,76
04.09.2016 9,4 9,48 25 48,87 27,8 46,51

 

Euclidian 2D distance without high
Environment Open area Buildup Vegetated
Mean(0:1) Aver NonAv Aver NonAv Aver NonAv
04,05,2016 5,01 5,61 17,00 12,73 1,19 3,94
04,06,2016 3,93 0,19 0,68 8,65 5,2 1,93
04,07,2016 24,3 9,31 2,6 7,86 2,01 6,97
04,08,2016 0,99 2,99 19,3 10,13 4,68 5,3
04,09,2016 1,04 1,08 8,07 11,57 1,23 12,53

 

: