<div dir="ltr"><div class="gmail_extra"><div><div class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr">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.</div></div><div>I have also tried the  newest gdal without success.</div><div>Cheers</div><div>Chris</div></div></div></div></div></div></div></div>
<br><div class="gmail_quote">On 12 October 2017 at 11:41, Chris Marsh <span dir="ltr"><<a href="mailto:chris.marsh@usask.ca" target="_blank">chris.marsh@usask.ca</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><div class="m_-2742599217339033247gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr">Hi,<br></div></div></div></div></div></div></div></div></div><div dir="ltr"><div><br></div><div>I'm using GDAL (and by extension proj.4) to reproject a set of WGS84 coordinates to North America Albers Conic Conformal.</div><div><br></div><div>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</div><div><br></div><div>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'.<br></div><div><br></div><div>I have tried building and linking against various proj.4 versions, none of which solves this problem.</div><div><br></div><div>The code  is straightfoward</div><div><div>OGRSpatialReference insrs;</div><div>    insrs.SetWellKnownGeogCS("<wbr>WGS84");</div><div><br></div><div>    OGRSpatialReference outsrs;</div><div>    outsrs.importFromProj4("+proj=<wbr>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");</div><div><br></div><div>    OGRCoordinateTransformation* coordTrans =  NULL;</div><div>    coordTrans = OGRCreateCoordinateTransformat<wbr>ion(&insrs, &outsrs);</div><div><br></div><div>    coordTrans->Transform(1, &longitude, &latitude)</div></div><div><br></div><div>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.</div><div><br></div><div>I am using gdal 2.2.1 and proj4 4.8 (4.9 has other issues I'm running into)</div><div><br></div><div><div>$gcc -v</div><div>Using built-in specs.</div><div>COLLECT_GCC=gcc</div><div>COLLECT_LTO_WRAPPER=/cvmfs/<a href="http://soft.computecanada.ca/nix/store/lvjwgn6hpngyy6k4xqcqa9h2cxy3fl30-gfortran-5.4.0/libexec/gcc/x86_64-unknown-linux-gnu/5.4.0/lto-wrapper" target="_blank">sof<wbr>t.computecanada.ca/nix/store/<wbr>lvjwgn6hpngyy6k4xqcqa9h2cxy3fl<wbr>30-gfortran-5.4.0/libexec/gcc/<wbr>x86_64-unknown-linux-gnu/5.4.<wbr>0/lto-wrapper</a></div><div>Target: x86_64-unknown-linux-gnu</div><div>Configured with:</div><div>Thread model: posix</div><div>gcc version 5.4.0 (GCC)</div></div><div><br></div><div><div>$uname -a</div><div>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</div></div><div><br></div><div>Cheers</div><span class="HOEnZb"><font color="#888888"><div>Chris</div><div><br></div></font></span></div></div>
</blockquote></div><br></div></div>