Highly Accurate GPS by averaging?

jarvist

New member
Hi Footleg,

If it's any use, my GPSBabel command to get info out was:
Code:
gpsbabel -T -i garmin -f /dev/ttyUSB0 -F testlog.dat -o xcsv,style=waypoint_data_log.style

The style file is here:
http://gist.github.com/615686

Output looks like this with NMEA input:
Code:
1277760772 51.xx6081666666664 -0.xx8856666666667 112.700000 11 5.400000 1.000000 5.500000
1277760774 51.xx6048333333336 -0.xx8883333333333 112.600000 11 2.900000 3.100000 4.300000
1277760776 51.xx6003333333336 -0.xx8911666666667 112.500000 11 2.900000 3.100000 4.300000
- truncation in the 7th decade of degrees, about 4cm.

And with garmin input:
Code:
 1277761219 51.xx9679420030139 -0.xx2414213831683 95.732574 -1 0.000000 0.000000 0.000000
1277761220 51.xx9553679179760 -0.xx2397661923045 95.686760 -1 0.000000 0.000000 0.000000
1277761221 51.xx9434218338981 -0.xx2381936684230 95.643234 -1 0.000000 0.000000 0.000000
- seemingly more numerical precision in the lat/long, but it could just be because it's transmitted in binary + so not obvious by inspection in decimal the degree of truncation.

(Obv. xx's are to anonymise my secret hideout)
 

footleg

New member
Thanks Jarvist, a really interesting and clear article. It appears the bottom line is that the errors in the calculated positions recorded at the same time by two different GPS receivers will not be the same, or even have much correlation to each other. So if you just subtract the position calculated by one receiver from another then you are as likely to get a less accurate relative position as a more accurate one. It appears that NMEA data from consumer units does not contain the error data for each satellite which is required to calculate a more accurate relative position. So using two consumer GPS receivers will not gain anything over just averaging the reading from a single receiver at the location and treating the reading as accurate to around 3m.

So given that differential GPS just does not appear possible from consumer receivers, lets just consider the accuracy of a single WAAS enabled GPS. This is quoted as accurate to 3m. This is presumably in ideal conditions. Can I assume that if I have a WAAS average being displayed on my receiver that the calculated position will be within 3m of the real position of the point I am measuring at? If this is the case, then does this mean that on one occasion I might get a reading 3m North of the real position, and on another occasion the reading might be 3m South of the real position? I am always within 3m of the real position, but might see up to 6m difference between readings taken on different days?

Finally, if I log data at my location over a long period (say several hour long sessions on different dates) and plot the positions on a scatter plot, I would expect to get something a bit like this:
before-sm.gif


Can I then correctly conclude that my real position is roughly in the centre of that splodge of points? Will this be more accurate than +/-3m ?
 

kdxn

New member
Nice thought. Here's some geeky stuff to ponder upon.......

Your cheapie handheld receivers use what we call Coarse Acquisition Code on the L1 frequency. The problem with this is that the signals get delayed by the ionosphere and the delay will depend upon where the satellite is with respect to you and where you are in the world. Low elevation satellite ranges tend to be noisier as they have more stuff to filter through, additional troposphere to deal with. The ionosphere noise varies with latitude, what time of day, what azimuth is has to you and also where we are in the 11 year solar cycle. What I am saying is that your simple handheld receiver has a lot to try and compensate for and it does this with an ionospheric model and a tropospheric model. Not to mention that the signal correlators on your typical handheld are good enough for a handheld - that is the ranges behind the position calculation do not have to be that good because most users are only expecting a few metres of accuracy - why increase the price of the product needlessly ? The WAAS and EGNOS augmentation satellite signals help improve upon these models but they are still models which have been forecast from prior data so your position is probably biased and may be biased for some time so averaging may look good but your average is probably biased. Comparing position fixes between two receivers is very dodgy because the sky view at both receivers is rarely the same so you can not guarantee that both receivers were using exactly the same satellites at the same time. Also you may have different receiver hardware, different firmware and my earlier comment about the fundamental accuracy of these handhelds. When your WAAS and EGNOS position indicates 3m acccuracy, that is based on the least squares fix of all the satellite ranges used with the augmentation signals. If you only have four satellite ranges but they all agree well with one another then the position statistic could be very good, however all four satellite ranges could be biased ! The 3m indication is a statistical amount - it does not mean that all your positions are within 3m - it suggests that a percentage of your position fixes are within 3m. Whats happening with the other percentage of positions ?......................could be anywhere. How big the percentage is depends upon the statistic being used, you often see 1,2...x sigma or CEP Circular Error Probability but here we are perhaps delving too deep.

The only way to get a really good improvement on your fix is to access the individual satellite ranges and preferably the carrier phase(s) too, some of the cheapies will provide low quality L1 phase. Get your data in a RINEX file and see the post processing online links that I posted earlier. These can process your raw data with respect to multiple ground stations equipped with very high quality geodetic receivers. If you wait a little while, you can also use very accurate post processed orbits, ionospheric and tropospheric data all of which can improve your handheld positions. However one last point, a cheapie can only get you so far because you run into the underlying limitations of the silicon and firmware inside. Do not forget you also have the issue of height - GPS gives you a spheroidal height - getting a mean sea level elevation requires the use of  a model - typically this is okay over short distances for relative height differences but there can be some dramatic differences between spheroid height and mean sea level. The UK has an east-west tip of a few metres - India has a 180m difference from south to north.

I am thinking that perhaps a presentation at HE 2011 could be useful. If there is enough interest, I could put one together and do a preliminary run through at the next Cave Electronics and Surveying Group meeting early 2011 with some example data from Mulu 2011 as I will be taking a professional GPS out there and should have the post processed results by then, subject to canopy coverage. GPS signals do not like wet foliage !
 

footleg

New member
Sounds like a great topic for a talk at the next cave tech symposium. I'd certainly find it very interesting.
 

kdxn

New member
OK, looks like I may be doing two talks (if time permits) at the next BCRA cave technology symposium to be held 9-10 April 2011 at Hulland Ward in Derbys. First will be Surface Mapping, second will hopefully be Laser Scanning.

Have dusted off my old GPS Garmin, DeLorme EarthMate and the big bad boy Leica SR530.

Get back from Mulu 31st March so will need to get a move on to include any results...............may have to be preliminary.............full results for HE 2011.
 

kdxn

New member
Thanks to those of you at the BCRA Tech symposium this last weekend.

Given the feedback received, I will propose a Surface Mapping talk for HE2011 and concentrate on practical methods for improving position using handheld (cheap) technology and also professional gear with some comparisons. Hopefully be of interest for UK and expeditions.
 

footleg

New member
I think that would be very interesting. I'll try not to take on too many jobs at Hidden Earth this year so I stand a small chance of being able to get to some talks!
 
Top