[Proj] NZGD2000 to NZGD49

[Simon O'Malley]

Hi

I am trying to test the datum shift (using nzgd2kgrid0005.gsb file) for
conversions between NZGD2000 and NZGD49(NZMG).

I have downloaded the proj-datumgrid-1.3.zip file, and extracted everything
under the nad directory.

configured, make, and make install. However the gsb file does not get copied
to /usr/local/share/proj

I copied it across but postgis does not appear to use it i.e. my transformed
coordinates are not correct, as verified by the linz website and other
software that uses the grid file.

I then tried to use cs2cs to validate it but all attempts at specifying the
projections and datums gave me errors.

I have tried multiple versions of:

cs2cs +proj=latlong +datum=nzgd2000 +nadgrids=nzgd2kgrid0005 +to +proj=nzmg
+datum=nzgd49

After setting PROJ_DEBUG I get:

pj_open_lib(proj_def.dat): call fopen(/usr/local/share/proj/proj_def.dat) -
succeeded
Using from definition: proj=latlong datum=nzgd2000 nadgrids=nzgd2kgrid0005
Rel. 4.4.9, 29 Oct 2004
<cs2cs>:
projection initialization failure
cause: unknown elliptical parameter name
program abnormally terminated

Tried upper and lower case NZGD2000, NZGD49, NZMG, latlong, longlat, without
nadgrids???
Do I have to add an ellipsoid to the projection definitions for cs2cs?

The grid shift file does not appear to mentioned in any of the files in
/usr/local/share/proj is that normal?

Some pointers would be appreciated.

Simon
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://lists.maptools.org/pipermail/proj/attachments/20060517/f02e1457/attachment.html