[Proj] OSTN02 question: GRS80 instead of Airy.

Roger Oberholtzer roger at opq.se
Mon Oct 30 16:25:08 EST 2006


On Mon, 2006-10-30 at 15:04 -0500, Frank Warmerdam wrote:
> Roger Oberholtzer wrote:
> > No takers on this one?
> > 
> > On Wed, 2006-10-25 at 12:58 +0200, Roger Oberholtzer wrote:
> >> I use the C API to do transformations. I am converting from GPS values
> >> in WGS84 to OSTN02 for the United Kingdom. My current code is:
> >>
> >>        fromProj = pj_init_plus(
> >>                 "+proj=latlong "
> >>                 "+ellps=WGS84 "
> >>                 "+towgs84=0,0,0 "
> >>                 "+nodefs")
> >>
> >>         toProj   = pj_init_plus(
> >>                 "+proj=tmerc "
> >>                 "+lat_0=49 "
> >>                 "+lon_0=-2 "
> >>                 "+k=0.9996012717 "
> >>                 "+x_0=400000 "
> >>                 "+y_0=-100000 "
> >>                 "+ellps=airy "
> >>                 "+nadgrids=gb.gsb "
> >>                 "+units=m "
> >>                 "+no_defs")
> >>
> >> Note that I do a grid shift, using data provided by a kind member of
> >> this list.
> >>
> >> My question is, how can I use the GRS80 ellipsoid instead of Airy? 
> 
> Roger,
> 
> Well, you would just change +ellps=airy to +ellps=GRS80.
> 
> But perhaps your question was more subtle?

Or not. If I change the +ellps=airy to +ellps=grs80, leaving all else
the same, I see no difference at all.

Ad an example, for lat=50.34180000 lng=-4.82595000, I get
Easting=199020.37000 Northing=52918.58000 for both airy and grs80.

It could be that these should be the same. But I thought grs80 would
result in some small difference.

--
Roger Oberholtzer





More information about the Proj mailing list