[gdal-dev] ERROR 6: Failed to initialize PROJ.4 with "..."

Chris Marsh chris.marsh at usask.ca
Thu Oct 19 12:19:52 PDT 2017


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 <even.rouault at spatialys.com>
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 <chris.marsh at usask.ca> 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
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/gdal-dev/attachments/20171019/062364d8/attachment.html>


More information about the gdal-dev mailing list