Hi, I build GDAL by specifying --with-static-proj4 on the configure command line. This does not seem to have resolved the problem.
echo 2 49 | build-release/lib/gdal/bin/gdaltransform -s_srs EPSG:4326 > -t_srs EPSG:32632 > -11818.9526229417 5451106.95020125 0 Cheers Chris On 18 October 2017 at 12:58, Even Rouault <[email protected]> wrote: > On mercredi 18 octobre 2017 10:20:41 CEST Chris Marsh wrote: > > I was wondering if anyone had thoughts on this? I posted to the proj.4 > > list, but seems like this may be a gdal problem and not proj.4 as the > > proj.4 command line has no problem with this string. > > I have also tried the newest gdal without success. > > Try building --with-static-proj4 to avoid the issues with the dlopen()'ing > of > libproj.so > > Check if the following works > > echo 2 49 | gdaltransform -s_srs EPSG:4326 -t_srs EPSG:32632 > > > > Cheers > > Chris > > > > On 12 October 2017 at 11:41, Chris Marsh <[email protected]> wrote: > > > Hi, > > > > > > I'm using GDAL (and by extension proj.4) to reproject a set of WGS84 > > > coordinates to North America Albers Conic Conformal. > > > > > > On all but 1 computer I run my code on, this works as expected. > However, > > > when trying to run on one of Compute Canada's newest machines, I get > the > > > error > > > > > > ERROR 6: Failed to initialize PROJ.4 with `+proj=aea +lat_1=20 > +lat_2=60 > > > +lat_0=40 +lon_0=-96 +x_0=0 +y_0=0 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 > > > +units=m +no_defs'. > > > > > > I have tried building and linking against various proj.4 versions, > none of > > > which solves this problem. > > > > > > The code is straightfoward > > > OGRSpatialReference insrs; > > > > > > insrs.SetWellKnownGeogCS("WGS84"); > > > > > > OGRSpatialReference outsrs; > > > outsrs.importFromProj4("+proj=aea +lat_1=20 +lat_2=60 +lat_0=40 > > > > > > +lon_0=-96 +x_0=0 +y_0=0 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m > > > +no_defs"); > > > > > > OGRCoordinateTransformation* coordTrans = NULL; > > > coordTrans = OGRCreateCoordinateTransformation(&insrs, &outsrs); > > > > > > coordTrans->Transform(1, &longitude, &latitude) > > > > > > with the error appearing on the Transform call. The proj.4 string is > from > > > using the ExportToWkt of GDAL for the WKT of Albers Conic Conformal. > (long > > > story as to why i'm doing this). I have tried using the WKT string and > > > that > > > results in the same error -- seems like it is just converted to proj.4 > > > behind the scenes. > > > > > > I am using gdal 2.2.1 and proj4 4.8 (4.9 has other issues I'm running > > > into) > > > > > > $gcc -v > > > Using built-in specs. > > > COLLECT_GCC=gcc > > > COLLECT_LTO_WRAPPER=/cvmfs/soft.computecanada.ca/nix/store/ > > > lvjwgn6hpngyy6k4xqcqa9h2cxy3fl30-gfortran-5.4.0/libexec/gcc/ > > > x86_64-unknown-linux-gnu/5.4.0/lto-wrapper > > > Target: x86_64-unknown-linux-gnu > > > Configured with: > > > Thread model: posix > > > gcc version 5.4.0 (GCC) > > > > > > $uname -a > > > Linux gra-login4 3.10.0-514.21.1.el7.x86_64 #1 SMP Thu May 25 17:04:51 > UTC > > > 2017 x86_64 GNU/Linux > > > > > > Cheers > > > Chris > > > -- > Spatialys - Geospatial professional services > http://www.spatialys.com >
_______________________________________________ gdal-dev mailing list [email protected] https://lists.osgeo.org/mailman/listinfo/gdal-dev
