Aha! Thanks Even, setting PROJ_DEBUG=ON showed the problem. I had assumed when giving +geoidgrids=g2009conus.gtx, it would find the gtx in the current directory. But, now i see it looks in only one place, in the PROJ_LIB location. That's OK. I copied my .gtx files there, and now it works.
Yes, i'm using the SVN version of PROJ4. BTW, it would be really good to have a new release of PROJ4. The last was in 2009. And it is confusing, because the executables built from SVN report "4.7.1" when run, and Frank's page on http://trac.osgeo.org/proj/wiki/VerticalDatums says "from version 4.8.0 on there is support.." when we haven't reached that version yet. -Ben > -----Original Message----- > From: Even Rouault [mailto:[email protected]] > Sent: Wednesday, August 17, 2011 11:13 AM > To: [email protected] > > Le mercredi 17 aot 2011 03:24:58, Ben Discoe a crit : > > It is true that geocentric XYZ coordinates are not very useful, without > > being able to convert them to coordinate with orthometric ("geoid") ("sea > > level") height. > > > > Frank implemented 'preliminary' support for this, which i am preparing to > > test now. Before expecting it to work in GDAL/OGR, i thought i should try > > it down in PROJ.4 first, following the description on > > http://trac.osgeo.org/proj/wiki/VerticalDatums > > Not sure what you would be doing wrong, but works for me : > > $ PROJ_DEBUG=ON ~/gdal/proj/install-proj/bin/cs2cs +proj=latlong > +datum=WGS84 > +to +proj=utm +zone=10 +datum=WGS84 +geoidgrids=g2009conus.gtx > pj_open_lib(proj_def.dat): call fopen(/home/even/gdal/proj/install- > proj/share/proj/proj_def.dat) - succeeded > > -122.272778 37.871667 10 > pj_open_lib(g2009conus.gtx): call fopen(/home/even/gdal/proj/install- > proj/share/proj/g2009conus.gtx) - succeeded _______________________________________________ gdal-dev mailing list [email protected] http://lists.osgeo.org/mailman/listinfo/gdal-dev
