Context We are trying to reproject a GeoJSON polygon to EPSG:4326. In some cases, depending on whether we use the ogr2ogr command line or the equivalent in python, the resulting polygon is different. In those situations, the command line gives us the appropriate polygon, whereas the python code gives us a distorted polygon.
Question Any idea why there is a difference between the command line and the python results? Is there a mistake in our interpretation of the code/how we converted it in python? Sample GeoJSON: https://pastebin.com/hUiVVjFv Proj4: "+proj=stere +lat_0=90 +lat_ts=60 +lon_0=264 +x_0=0 +y_0=0 +R=6371229 +units=m +no_defs" Command line that works: ogr2ogr -s_srs "+proj=stere +lat_0=90 +lat_ts=60 +lon_0=264 +x_0=0 +y_0=0 +R=6371229 +units=m +no_defs" -t_srs EPSG:4326 result_command_line.json geojson.json Result: https://pasteboard.co/cTpZBRamEYv8.jpg Python code (we would expect the same output as the command line): https://pastebin.com/QgWDBL1W Result: https://pasteboard.co/8iSv0vbnQlHw.jpg Thank you! Nicolas D.
_______________________________________________ gdal-dev mailing list gdal-dev@lists.osgeo.org https://lists.osgeo.org/mailman/listinfo/gdal-dev
