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 direct
transform but not the reverse. Any place where I can search what other then
Airy have the same issue?
echo 4.897 52.371 | mapproject -J+proj=airy
ERROR 1: PROJ: proj_create: Error -50: malformed pipeline
ERROR 6: Cannot find coordinate operations from
`PROJCRS["unknown",BASEGEOGCRS["unknown",DATUM["World Geodetic System
1984",ELLIPSOID["WGS
84",6378137,298.257223563,LENGTHUNIT["metre",1]],ID["EPSG",6326]],PRIMEM["Greenwich",0,ANGLEUNIT["degree",0.0174532925199433],ID["EPSG",8901]]],CONVERSION["unknown",METHOD["PROJ
airy"]],CS[Cartesian,2],AXIS["(E)",east,ORDER[1],LENGTHUNIT["metre",1,ID["EPSG",9001]]],AXIS["(N)",north,ORDER[2],LENGTHUNIT["metre",1,ID["EPSG",9001]]]]'
to `GEOGCRS["unknown",DATUM["World Geodetic System 1984",ELLIPSOID["WGS
84",6378137,298.257223563,LENGTHUNIT["metre",1]],ID["EPSG",6326]],PRIMEM["Greenwich",0,ANGLEUNIT["degree",0.0174532925199433],ID["EPSG",8901]],CS[ellipsoidal,2],AXIS["longitude",east,ORDER[1],ANGLEUNIT["degree",0.0174532925199433,ID["EPSG",9122]]],AXIS["latitude",north,ORDER[2],ANGLEUNIT["degree",0.0174532925199433,ID["EPSG",9122]]]]'
mapproject [ERROR]: Failed to create coordinate transformation between the
following
coordinate systems. This may be because they are not transformable,
or because projection services (PROJ.4 DLL/.so) could not be loaded.
mapproject [ERROR]: Source:
PROJCS["unknown",
GEOGCS["unknown",
DATUM["WGS_1984",
SPHEROID["WGS 84",6378137,298.257223563,
AUTHORITY["EPSG","7030"]],
AUTHORITY["EPSG","6326"]],
PRIMEM["Greenwich",0,
AUTHORITY["EPSG","8901"]],
UNIT["degree",0.0174532925199433,
AUTHORITY["EPSG","9122"]]],
PROJECTION["custom_proj4"],
UNIT["metre",1,
AUTHORITY["EPSG","9001"]],
AXIS["Easting",EAST],
AXIS["Northing",NORTH],
EXTENSION["PROJ4","+proj=airy +datum=WGS84 +wktext +no_defs"]]
GEOGCS["unknown",
DATUM["WGS_1984",
SPHEROID["WGS 84",6378137,298.257223563,
AUTHORITY["EPSG","7030"]],
AUTHORITY["EPSG","6326"]],
PRIMEM["Greenwich",0,
AUTHORITY["EPSG","8901"]],
UNIT["degree",0.0174532925199433,
AUTHORITY["EPSG","9122"]],
AXIS["Longitude",EAST],
AXIS["Latitude",NORTH]]
ERROR 10: Pointer 'hTransform' is NULL in 'OCTTransform'.
ERROR 10: Pointer 'hTransform' is NULL in 'OCTTransform'.
ERROR 10: Pointer 'hTransform' is NULL in 'OCTTransform'.
ERROR 10: Pointer 'hTransform' is NULL in 'OCTTransform'.
328249.003313 4987937.10145
_______________________________________________
gdal-dev mailing list
gdal-dev@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/gdal-dev
