The Analysis of Scientific and Commercial Softwares Accuracy in GPS Observation Processing

The engineers use various softwares for observation processing and GPS networks, generally everyone who uses a device, he uses softwares with the same brand for the processing of it GPS processing networks by defined default parameters and models. The most outstanding issue is that, it should be understood, the application of which kinds of softwares will help us to gain the desirable accuracy and whether the accuracy of coordinate estimation is acceptable in all commercial softwares or not. In this article, we compared a series of softwares which engineers use for their data processing with powerful scientific software that meanwhile it is too hard and time consuming to learn it. We have considered accuracy and validity of scientific software, which gives acceptable results, as a baseline to see how much errors and difference is between this software and the commercial softwares, so we selected Bernese as a scientific software and LGO, Ashtech Solutions, CGO softwares as commercial softwares and we analyzed them and compared the applied models of the softwares with each other. To do so, we choose a GPS network from Tehran city and discussed it, in order to evaluate the measured accuracy for baselines and also coordinate system that after data processing, it was understood that if the polygon of a network is small, commercial softwares can be utilized with very high acceptable accuracies and even by Ephemerides Precise data, their accuracy will increase more, but if our network is more extended and bigger, there is no other choice except using scientific softwares to calculate coordinate that performs a processing by many parameters.


Introduction
In the middle of 1990s, GPS inventors introduced Real-Time Kinematic (RTK) method, that by means of it, we can achieve the desirable accuracies in surveying.In this method, the user usually is in need of a reference station which is located at maximum 10 km distance, in order to determine the position with centimeter accuracy.The applications of this method account for: positioning with high accuracy for navigation systems, guiding auto-run machines, unmanned aerial vehicles, marine missions like: dredging and large-scale map's update [1].
In general, RTK is a process in which GPS signal corrections are communicated from a reference station's receiver with known position to one or more rover receivers in a real-time strategy [2].It means that, by short time observations, the determination of dynamic position is possibly prepared and by carrier code signals, the position of rover station is determined with centimeter accuracy immediately [3].
By increasingly development of GNSS satellite navigation and emergence of new companies for production of related receivers, numerous softwares have been developed for processing of these data that in the most general form, these softwares can be divided into scientific and commercial groups.Scientific softwares that mainly are prepared by universities or reliable scientific centers are programmed for utilization in precise and professional tasks.In addition, commercial softwares which are produced by those companies that produce GNSS receivers are applied for data processing of general users and engineers.Commercial softwares not only don't have the complexity of professional and scientific softwares, they have a plain graphical and user friendly environment that it has simplified and quicken dealing with these softwares.While scientific softwares are much more complex and have specific and professional settings that make it too hard to handle them [4].Nowadays, we are all witness for emergence of many commercial softwares that, it is very crucial to choose the best and the most precise software out of these softwares.Considering the ideas of geomatics engineers, these commercial softwares: CGO, LGO and Ashtech Solutions, are the most used softwares in Iran that we compare this commercial software with a scientific software named Bernese to receive the best answers from comparisons.In this study, we don't aim to describe special software, but the goal is to compare softwares with default methods and an algorithm of softwares.The reason for using the softwares defaults is that, geomantic engineers presented that they don't change any default except the zone of region while working with the softwares.Our theory is that the results and accuracy of scientific softwares is more precise than the results mined from commercial softwares.So our main goal is to compare calculated baselines and estimated coordinate of a process from a network of GPS stations of GNSS servers by scientific and commercial software with each other and conduction of a comparison among commercial softwares, regarding their accuracy in baselines, coordinate and Standard Deviation.Of course, it should be mentioned that, the present softwares or definable models for each software will be totally different from each other.

Models of Bernese Scientific Softwares
In this part, we will introduce various models of this scientific software and we will explain each of them briefly.

Offset Model for Antenna Phase Center of Satellite and Receiver
Offset model of satellite antenna phase and variation of electrical phase of receiver antenna and offset center of receiver's antenna have many effects on PPP processing.So models with ANTEX format have been applied to delete these effects.In Bernese software, ".rel" format is used to omit these effects and in fact, rel. is the ANTEX format of Bernese software.To achieve this format, in addition of using the products of CODE, PHCCNV software can be applied that in Bernese software is prepared for conversion of ANTEX format to ".rel" format [5].

Error of Satellite Clock and Error of Satellite Position
Like other scientific softwares, Bernese software also, uses IGS products to decrease the orbital errors of satellite and the errors of satellite clock, so Bernese software owns this ability to decrease the errors of satellite clocks and errors of satellite clock by using ".clk" and ".sp3" formats.Furthermore, to decrease the orbital errors of satellites, it can use ".sp3" and ".pre" and ".eph" formats too.

Model of Relocation from Ocean Upload
In order to omit the effect of relocation caused by ocean upload, Bernese software uses ".blq" standard suffix of these models.One of the most famous centers that gives these models is Onsala observatory.The users can refer to the website of this observatory and fill the related form to receive these models [6].

Earth Rotation Parameters
Earth rotation parameters are those parameters that are presented in ".eop" standard format.These parameters both are presented by IGS processing center weekly and also they are shown in two (A) and (B) Bultons by IERS.
Bulton (A) includes fast earth rotation parameters that are presented daily and Bulton (B) are earth rotation parameters which are published monthly.

Satellite State Space Model
Satellite state space model are models which are released annually by AIUB.These models reveal satellites' safety at different episodes and Bernese software uses these models for proving the satellites' safety at observing time and if satellites' safety is not proved, that satellites' observation would be deleted from observation list [7].

Satellite Ephemeris Model
Satellite Ephemeris models are binary files with ".eph" format.These models are presented by solar dynamics observatory, Jet Propulsion Laboratory and NASA.
Up to now, three series of these models are presented in these numbers: DE100, DE200 and DE400.Bernese software uses these models to calculate the position of moon, sun and other planets.
Since, the binary format of these models are not free, Bernese software uses a tool to convert numeral format of these models into binary format of Bernese [8].

Ocean Tide Model
The ocean tide models have ".tid"Suffix that most of them are prepared by space research center of Texas University.Bernese software exploits these models for calculation of variations in gravity caused by oceanic tides.
Note: in scientific softwares, we don't have model for ionosphere and troposphere and the software itself calculates them scientifically [10].

Applied Models in Commercial Softwares
As mentioned before, in scientific softwares, we don't have any ionospheric and tropospheric errors and the software calculates them scientifically but in commercial softwares, we decrease or omit them by models.In commercial software, we tune the ionosphere error model in automatic state that if the length of network is large, L3 frequency will be used that in this condition, the ionosphere error would be deleted thoroughly but the problem of this work is that the noise is very much.But if the length of network is short, L1 + L2 frequency is used that reduces the errors.

Applied Models by LGO Commercial Software for Processing
Parameters and cases that should be set and investigated before processing by LGO software are mentioned in Table 1.
Table 1.Applied models by LGO commercial software for processing.

Applied models by Ashtech Solutions Commercial Software for Processing
Parameters and cases which should be set and investigated before processing in Ashtech Solutions software are mentioned in Table 2.

Applied Models by CGO Commercial Software for Processing
Some of changeable parameters in CGO software are mentioned in Table 3.

Research Zone
Selected zones for GPS processing network are from RTK network in Tehran city.Tehran's real-time positioning stations are designed in a way that can cover all 22 regions of Tehran city.In this case study network, we have used 6 regions that (TEHN) is considered as a reference point and others (M804, M318, M022, M020 and TMIC) as rover points.These points of Tehran's City council which are used as permanent stations are chosen.These places are monitored structurally and technically and it is good to say, the points are in a complete safety.
Furthermore, to 20 km distance from each station out of network, it is possible to have access to centimeter accuracy in Single Based format.RTK system is like differential corrections with a difference, this method applies a local reference Base that messages corrections with high accuracy to a rover receiver by radio communication.Thus, the distance and direct vision of base to rover station is an outstanding parameter in this method that causes some limitations while dealing with it, this implies that by increasing the distance of rover station from base one, the accuracy of positioning would face with reduction and the reason is that, ionospheric and tropospheric error models are not equal in different regions and also, it is because of systematic effects on Ephemerids.In fact, limitations such as distance between base and rover stations and mortality of radio   A solution that nowadays, is applied in many countries to solve these kinds of problems is to exploit some reference stations and in fact it is formation of a national or local network that sends RTK information and in hence, just by exploiting one GPS receiver, centimeter accuracy is achieved in real-time state.In this method, that is called Network RTK (opposite to the first state that is named Single Base RTK); an unknown station can use one or some known stations (reference) for positioning.The increase in number of reference stations leads to adjustment condition's optimization of observational errors and getting access to more precise results.In fact, the reason for insisting on RTK networking method is that, it can achieve accuracy more than 3 centimeters just only by exploiting one human force.In Figure 1, you can see an outline of the points in urban plan and also satellite images are shown.

The location of zones along with Baseline on city map
The location of zones on satellite images

Numeral Results of Comparison between Commercial Software and Scientific Software
After entering the information and data to each of softwares, we started to process GPS network that the results are given in Tables 4-6.The results of each commercial software have been compared with scientific software.In these tables, the names of stations and their 3-dimentional positions are mentioned that easily, we can investigate and compare the accuracy of each calculated cases [11].While typing the data of positions to LGO software, the same mentioned parameters in part (3-1) which are the defaults of the software, have been applied but for outputs of Processed data and the final adjustment, the software gives the three dimensional default and extracts X,Y coordinates out of reference ellipsoid (WGS84) and reports by modeling of ellipsoid removal from geoid and regarding to the region's zone, it transfers reference ellipsoid to geoid and gives the correct height of the point compared to the geoid [12].
For Ashtech Solutions software like (2-3) part, to enter data, the same default parameters of software were used and for taking final output from this software, the same strategy of LGO software was applied.
In CGO software, based on the introduced default for the software in part , the data were entered.But while adjustment for coordinate and achieving the final report of software, it gives two options to user for adjustment: a) 2D adjustment and b) 3D adjustment.When adjustment is carried out two-dimensionally, (X,Y) coordinate estimation is performed by reference ellipsoid (WGS84) and gives very acceptable results that in Table 6, part CGO and North and East columns, the presented coordinate of two-dimensional adjustment options of software are depicted, but when we select the option of three-dimensional adjustment option, the results both for planimetric coordinate (X,Y) and height coordinate (Z), they face with marked and significant variations, while in three-dimensional adjustment option, height coordinate is calculated just based on reference ellipsoid and does not give a model for ellipsoid removal from geoid.This means that, (Z) coordinate that should be estimated based on geoid, is calculated by ellipsoid.This subject can be observed by calculated (Z) coordinate that is shown in Table 6 for (TEHN) control point.Of course, this numeral amount is ellipsoid removal in this region that almost, it is a fixed number for all (Z) coordinates in all points.
In following, we can mention to Tables 7-10 that depict the results of calculations of baseline lengths by each software among main stations.Now, we compare the commercial softwares with scientific softwares that have been as the baseline of research in picture format and by showing the difference amount of 2D coordinate that the main goal of this paper is to investigate the difference between 2D coordinate and the difference between baseline lengths.
The comparison of Bernese scientific software and LGO commercial software for baseline length is shown in Figure 2 and the results of medium difference are presented in Table 11.Also, the comparison of Bernese scientific software and Ashtech Solutions commercial software are presented for baseline lengths in Figure 3 and the result of dimension difference is written in Table 12.
And finally, the comparison of Bernese scientific software and CGO commercial software for baseline lengths are depicted in Figure 4 and the results of dimension difference are mentioned in Table 13.
And at last, it should be mentioned that in planimetric coordinate estimation (X, Y), among commercial softwares, CGO has the best function and its results      are very similar and close to Bernese scientific software, but in estimation of height coordinate (Z), Ashtech Solutions software has higher accuracy than other commercial softwares.

Conclusion
In this research, the goal is not to bold any specific commercial software but it is to compare these softwares with reliable scientific software in order to help the engineers to compare their applied softwares with real results and to know the accuracy of their work.According to the results, it was obvious that commercial softwares just can present acceptable accuracy in small polygons considering their capacity and few numbers of parameters for processing of GPS data and if there is any need for higher accuracy in small polygons amounts, we can achieve higher accuracy (the difference is in mm scale) by entering Ephemerids Precise data after some weeks and for bigger polygons, scientific softwares can be applied because these kinds of softwares are more functional and have numerous parameters for GPS data processing because if commercial softwares are utilized, the results will be unacceptable.
waves lead to undesirable accuracy in farther distances and long periods of time and high fees.

Figure 1 .
Figure 1.Location of these zones in satellite images.

Figure 2 .
Figure 2. Graphic presentation of comparison between Bernese scientific software and LGO software.

Figure 3 .
Figure 3. Graphic presentation of comparison between Bernese scientific software and Ashtech Solutions commercial software.

Figure 4 .
Figure 4. Graphic presentation of comparison between Bernese scientific software and CGO commercial software.

Table 2 .
Applied models by Ashtech Solutions commercial software for processing.

Table 3 .
Applied models by CGO commercial software for processing.

Table 4 .
Numeral results of comparison between LGO software and Bernese software.

Table 5 .
Numeral results of comparison between Ashtech Solutions software and Bernese software.

Table 6 .
Numeral results of comparison between CGO software and Bernese software.

Table 7 .
Baseline length in Bernese scientific software.

Table 8 .
Baseline length in LGO commercial software.

Table 9 .
Baseline length in Ashtech Solutions commercial software.

Table 10 .
Baseline length in CGO commercial software.

Table 11 .
Comparison of Bernese scientific software results and LGO commercial software at stations dimension.

Table 12 .
Comparison of Bernese scientific software results and Ashtech Solutions commercial software at stations dimension.

Table 13 .
Comparison of Bernese scientific software results and CGO commercial software at stations dimension.