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;
[hidden email]
_______________________________________________
UK mailing list
[hidden email]
http://lists.osgeo.org/mailman/listinfo/uk
