[gdal-dev] Which Proj.4 transforms are available in GDAL?
even.rouault at spatialys.com
Fri Aug 25 04:52:14 PDT 2017
> I'm implementing Proj.4 in GMT via GDAL and now, in the the testing stage,
> I'm using data from
> but to my surprise lots of projections are not implemented in GDAL. For
GDAL SRS model is based on WKT, so only projections that have a mapping to WKT 1 (and for which a developer cared enough to write the needed code !) are nominally supported
> gdaltransform -s_srs EPSG:4326 -t_srs +proj=airy +ellps=WGS84 +units=m
For unsupported projections, add +wktext. A dummy WKT representation is then built that captures the full proj.4 string in a EXTENSION node. Can also serve for supported projections if you want to overrite the WKT -> proj.4 string conversion done by GDAL
$ gdalsrsinfo "+proj=airy +ellps=WGS84 +units=m +wktext"
PROJ.4 : +proj=airy +ellps=WGS84 +units=m +wktext
OGC WKT :
EXTENSION["PROJ4","+proj=airy +ellps=WGS84 +units=m +wktext"]]
$ gdaltransform -s_srs EPSG:4326 -t_srs "+proj=airy +ellps=WGS84 +units=m +wktext"
> ERROR 1: Translating source or target SRS failed:
> gdaltransform -s_srs EPSG:4326 -t_srs +proj=bonne
> ERROR 6: Failed to initialize PROJ.4 with `+proj=bonne +lon_0=0 +lat_1=0
> +x_0=0 +y_0=0 +ellps=WGS84 +units=m +no_defs'.
I thing the issue is that lat_1=0 is invalid for bonne. Any non zero value is OK
> Is there a place where I can find the list of available proj.4 transforms?
You mean projections supported by GDAL (as far as I know, most if not all projections supported by GDAL translate then to proj.4) ?
Look at the various OSRSetXXXXX() methods in
Or starting there
> And also found some that do no agree between GDAL and proj. For example
> echo 4.897 52.371 | proj +proj=vandg +ellps=WGS84 +units=m
> 489005.93 6431581.02
> gdaltransform -s_srs EPSG:4326 -t_srs "+proj=vandg +ellps=WGS84 +units=m"
> 4.897 52.371
> 488459.293282675 6424391.46347316 0
If you add --debug on, you'll see that GDAL adds a +R_A parameter. I haven't dug up into history to know the reason for that.
If you add +wktext, the original string will be passed unmodified to proj.4
Spatialys - Geospatial professional services
-------------- next part --------------
An HTML attachment was scrubbed...
More information about the gdal-dev