Hi Anthony<br><br>Thanks for the tips.<br><br>I checked the epsg file and also read the man page for cs2cs which I should have done first.<br><br>Using cs2cs -lp and -ld showed that nzgd49 and nzmg were available but nzgd2000 was conspicuous by it's absence.
<br><br>The epsg file had the following for nzgd2000:<br><br>+proj=longlat +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +no_defs<br><br>So I then tried the following:<br><br>cs2cs +proj=longlat +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +no_defs +to +proj=nzmg +datum=nzgd49 +nadgrids=nzgd2kgrid0005
<br><br>Which seemed to work (the first set of coordinates entered get droppped, but I think I saw a bug logged for that already).<br><br>Using the following coordinates:<br><br>175.232739997 -37.3476467750<br><br>Gave me:
<br><br>2708093.01 6425845.96<br><br>Whereas the LINZ site (<a href="http://www.linz.govt.nz/rcs/linz/pub/web/root/core/SurveySystem/GeodeticInfo/CoordinateConversions/coordinateconversions/index.jsp">http://www.linz.govt.nz/rcs/linz/pub/web/root/core/SurveySystem/GeodeticInfo/CoordinateConversions/coordinateconversions/index.jsp
</a><br><br>Gave me the following for the conversion:<br><br>175.232739997
-37.347646775
2708091.7091
6425845.7365<br><br>Which is 2 meters different in the x coordinate. <br><br>I then set PROJ_DEBUG to see what was going on and it said it couldn't open /usr/local/share/proj/nzgd2kgrid0005, so I changed the command to add in the .gsb on the end and got the correct answer (saw somehwre on the net with the .gsb dropped off so incorrectly thought it was not needed).
<br><br>cs2cs +proj=longlat +ellps=GRS80 +a=6378137 +rf=298.257222101 +to +proj=nzmg +datum=nzgd49 +nadgrids=nzgd2kgrid0005.gsb<br>pj_open_lib(proj_def.dat): call fopen(/usr/local/share/proj/proj_def.dat) - succeeded<br>pj_open_lib(proj_def.dat): call fopen(/usr/local/share/proj/proj_def.dat) - succeeded
<br><br>pj_open_lib(nzgd2kgrid0005.gsb): call fopen(/usr/local/share/proj/nzgd2kgrid0005.gsb) - succeeded<br>NTv2 NZNAT 141x141: LL=(166,-48) UR=(180,-34)<br>pj_apply_gridshift(): failed to find a grid shift table for<br>
location (0.0000000dW,0.0000000dN)<br> tried: nzgd2kgrid0005.gsb<br>* * 0.00<br>175.232739997 -37.3476467750<br>NTv2 - loading grid NZNAT<br>pj_open_lib(nzgd2kgrid0005.gsb): call fopen(/usr/local/share/proj/nzgd2kgrid0005.gsb) - succeeded
<br>pj_apply_gridshift(): used NZNAT<br>2708091.71 6425845.74 -0.00<br><br>Success, thank you.<br><br>I mainly use postgres and have data stored in NZGD2000/latlong(SRID=4167), which I need to convert to NZMG (SRID=27200)
<br><br>How do I set up proj so that it uses the grid shift for this transform. Do I have to add something somewhere to make it use it?<br>Postgres doco states if it is there it will use it but how???? cs2cs had a lot of trouble, and using transform in postgres gives the same result as cs2cs not finding the gsb file.
<br><br>Simon<br><br><div><span class="gmail_quote">On 5/17/06, <b class="gmail_sendername">Anthony Dunk</b> <<a href="mailto:AnthonyD@rocketmail.com">AnthonyD@rocketmail.com</a>> wrote:</span><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
Hi Simon,<br><br>I think the problem is that NZ datums are not defined in PROJ (I found<br>the same was true of Australian datums). You will need to provide the<br>ellipsoid parameters for these datums instead. i.e. +a=6378388.0
<br>+rf=297.0 for NZGD49 and +a=6378137.0 +rf=298.257222 for NZGD2000.<br><br>So, to convert from NZGD200 Lat/Lon to NZGD49 Lat/Lon try this:<br><br>cs2cs +proj=latlong +a=6378137.0 +rf=298.257222 +to +proj=latlong<br>+a=
6378388.0 +rf=297.0 +nadgrids=nzgd2kgrid0005<br><br>Anthony.<br><br>--- Simon O'Malley <<a href="mailto:simonom@gmail.com">simonom@gmail.com</a>> wrote:<br><br>> Hi<br>><br>> I am trying to test the datum shift (using
nzgd2kgrid0005.gsb file)<br>> for<br>> conversions between NZGD2000 and NZGD49(NZMG).<br>><br>> I have downloaded the proj-datumgrid-1.3.zip file, and extracted<br>> everything<br>> under the nad directory.
<br>><br>> configured, make, and make install. However the gsb file does not get<br>> copied<br>> to /usr/local/share/proj<br>><br>> I copied it across but postgis does not appear to use it i.e. my<br>> transformed
<br>> coordinates are not correct, as verified by the linz website and<br>> other<br>> software that uses the grid file.<br>><br>> I then tried to use cs2cs to validate it but all attempts at<br>> specifying the
<br>> projections and datums gave me errors.<br>><br>> I have tried multiple versions of:<br>><br>> cs2cs +proj=latlong +datum=nzgd2000 +nadgrids=nzgd2kgrid0005 +to<br>> +proj=nzmg<br>> +datum=nzgd49<br>
><br>> After setting PROJ_DEBUG I get:<br>><br>> pj_open_lib(proj_def.dat): call<br>> fopen(/usr/local/share/proj/proj_def.dat) -<br>> succeeded<br>> Using from definition: proj=latlong datum=nzgd2000
<br>> nadgrids=nzgd2kgrid0005<br>> Rel. 4.4.9, 29 Oct 2004<br>> <cs2cs>:<br>> projection initialization failure<br>> cause: unknown elliptical parameter name<br>> program abnormally terminated<br>>
<br>> Tried upper and lower case NZGD2000, NZGD49, NZMG, latlong, longlat,<br>> without<br>> nadgrids???<br>> Do I have to add an ellipsoid to the projection definitions for<br>> cs2cs?<br>><br>> The grid shift file does not appear to mentioned in any of the files
<br>> in<br>> /usr/local/share/proj is that normal?<br>><br>> Some pointers would be appreciated.<br>><br>> Simon<br>> > _______________________________________________<br>> Proj mailing list<br>
> <a href="mailto:Proj@lists.maptools.org">Proj@lists.maptools.org</a><br>> <a href="http://lists.maptools.org/mailman/listinfo/proj">http://lists.maptools.org/mailman/listinfo/proj</a><br><br><br><br><br>__________________________________________________
<br>Do You Yahoo!?<br>Tired of spam? Yahoo! Mail has the best spam protection around<br><a href="http://mail.yahoo.com">http://mail.yahoo.com</a><br>_______________________________________________<br>Proj mailing list<br>
<a href="mailto:Proj@lists.maptools.org">Proj@lists.maptools.org</a><br><a href="http://lists.maptools.org/mailman/listinfo/proj">http://lists.maptools.org/mailman/listinfo/proj</a><br></blockquote></div><br>