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
_______________________________________________ gdal-dev mailing list [email protected] https://lists.osgeo.org/mailman/listinfo/gdal-dev
