[gdal-dev] PROJ4 Miller +R_A enforced

[Martin.Franzke at telekom.de]

Hi,

trying to transform from a "Miller WGS84" coordinate system, whose exact definition I don't know I found that the following two transformations yield the same result, regardless of whether +R_A is present or not.

$ echo "781046.71 6190010.02" | gdaltransform -s_srs "+proj=mill +lat_0=0 +lon_0=0 +x_0=0 +y_0=0 +datum=WGS84 +units=m +R_A" -t_srs EPSG:31466
2571974.08458298 5649095.03802563 5685.34965490364

$ echo "781046.71 6190010.02" | gdaltransform -s_srs "+proj=mill +lat_0=0 +lon_0=0 +x_0=0 +y_0=0 +datum=WGS84 +units=m" -t_srs EPSG:31466
2571974.08458298 5649095.03802563 5685.34965490364

This is due to https://github.com/OSGeo/gdal/blob/master/gdal/ogr/ogr_srs_proj4.cpp line 1852ff. I wonder, whether this is on purpose or not. At least using QGIS I get 2571781.3732362 5622898.28359721 in the latter case.

Kind regards,
Martin


-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/gdal-dev/attachments/20190121/1cb9f76e/attachment.html>