[gdal-dev] OGR/PROJ with +geoidgrids

[Ben Discoe]

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:even.rouault at mines-paris.org]
> Sent: Wednesday, August 17, 2011 11:13 AM
> To: gdal-dev at lists.osgeo.org
> 
> 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