[OSGeo-UK] OSGB36 WGS84 conversion problem -SOLVED

[John A Stevenson]

Thanks to those that got back to me.  I have solved it now, and attached 
some information sent to the GRASS list in case it it useful to anyone else.

Cheers

John

-----------------

Hi,

I had a problem importing GPS data and features digitised from Google 
Earth (both WGS84 datum) into locations with a British National Grid 
(OSBG36) projection.  I was working in a region in northwest Scotland, 
and I found that everything that I imported was plotting about 100 ENE 
of where it should do.

I have finally worked out why, but after a lot of searching, so I am 
posting a summary of what I found here so that google will pick it up in 
the archives if anyone else has the same problem in future.

The short answer is that the projection parameters that proj4 was using 
were wrong, but they are fixed now.  Noone using an up-to-date version 
of proj4 should have this problem, but the Dag rpm repository still has 
an old one.

The problem was that going to OSGB36 requires both a projection and a 
datum conversion {1}.  Prior to October 2007 {2}, proj was only carrying 
out the projection, thus resulting in the large offset.  You can check 
if you have the new version by running 'proj -v' or by looking at your 
epsg file:

cat /usr/share/proj/epsg | grep -A 1 "British National Grid"

# OSGB 1936 / British National Grid
<27700> +proj=tmerc +lat_0=49 +lon_0=-2 +k=0.9996012717 +x_0=400000 
+y_0=-100000 +ellps=airy +datum=OSGB36 +units=m +no_defs  <>

The new versions have +datum=OSGB36.

If you have an old version, you can correct it by replacing the line with:

+proj=tmerc +lat_0=49 +lon_0=-2 +k=0.999601 +x_0=400000 +y_0=-100000 
+ellps=airy 
+towgs84=446.448,-125.157,542.060,0.1502,0.2470,0.8421,-20.4894 +units=m 
+no_defs <>

A complication is that OSGB36 is slightly distorted with respect to GPS 
projections (such as WGS84 and ETRS89) {3}.  This offset is small, and 
is only important for higher precision surveying.  Many searches about 
OSGB36 offsets bring up pages relating to this.  If you want to 
compensate for this too,
you can download a nadgrid file and use it {4}.  For my data, this moved 
the points by about 1 m.

I hope that this helps anyone with the same problem,

Cheers

John



-- 


Dr John Stevenson
Postdoctoral Research Associate
School of Earth, Atmospheric and Environmental Sciences
Williamson Building (Room 2.42)
University of Manchester
Manchester M13 9PL, UK
tel. +44(0)161 306 6585; fax. +44(0)161 306 9361;
john.stevenson at manchester.ac.uk