Right. I changed the signature and removed the CPLFree() and now in GMT it also
... half works. The middle warnings are print from GMT but the other come from
GDAL and are like those you suppressed in your gdaltransform commit.
So, I should try to catch this type of cases where one can do the
> but raise the error when I change the order of the source and target (i.e.
> create the inverse mapping)
Sure, as that would trigger the non existing reverse operation.
>
> pSrcSRS = "+proj=airy +ellps=WGS84 +units=m"
> pDstSRS = "+proj=latlong +datum=WGS84"
>
> In this case, it even raises
I used the gdaltranslate example because I knew that it would much simpler for
you to test but, as I said first, this issue was raised by a call to GDAL from
GMT.
Looked again into it and what happens is that the
gmt_OGRCoordinateTransformation() call runs fine when
pSrcSRS = "+proj=latlong
On lundi 17 juin 2019 17:14:21 CEST Joaquim Manuel Freire Luís wrote:
> Well, I call it via GMT but that test is in GMT for a while and it used to
> pass.
>
> Anyway, I'll have to investigate this further because I rebuilt GDAL and now
> the gdaltransform works, but when I call it from GMT is
rom: Even Rouault
|>Sent: Monday, June 17, 2019 6:00 PM
|>To: Joaquim Manuel Freire Luís
|>Cc: gdal-dev@lists.osgeo.org
|>Subject: Re: [gdal-dev] Failure to convert to Airy
|>
|>> used to work before the GDAL3+PROJ6 combo.
|>
|>That would be surprising: with GDAL 2.4
> used to work before the GDAL3+PROJ6 combo.
That would be surprising: with GDAL 2.4 and PROJ < 6, you also get an error,
but with a different message:
$ gdalinfo --version
GDAL 2.4.1, released 2019/03/15
$ echo 4.897 52.371 | \
gdaltransform -s_srs +proj=longlat -t_srs "+proj=airy
geo.org
|>Cc: Joaquim Manuel Freire Luís
|>Subject: Re: [gdal-dev] Failure to convert to Airy
|>
|>On lundi 17 juin 2019 16:08:36 CEST Joaquim Manuel Freire Luís wrote:
|>> Even
|>>
|>> I'm seeing this failure with the airy conversion
|>>
|>> echo 4.897 52.371 | gd
On lundi 17 juin 2019 16:08:36 CEST Joaquim Manuel Freire Luís wrote:
> Even
>
> I'm seeing this failure with the airy conversion
>
> echo 4.897 52.371 | gdaltransform -s_srs +proj=longlat -t_srs +proj=airy
>
> ERROR 1: PROJ: proj_create: Error -50: malformed pipeline
> ERROR 6: Cannot find
Even
I'm seeing this failure with the airy conversion
echo 4.897 52.371 | gdaltransform -s_srs +proj=longlat -t_srs +proj=airy
ERROR 1: PROJ: proj_create: Error -50: malformed pipeline
ERROR 6: Cannot find coordinate operations from
`PROJCRS["unknown",BASEGEOGCRS["unknown",DATUM["World