coordinate systemgdalqgis
I have tried different open source tools to transform a world data shapefile from EPSG:4326 to EPSG:3857 but I always have an error. I'm not sure, but is it because of Antarctica?
I have tried with GDAL and QGIS and both are unable to do it. Why is it so complicated and how to do it?
The .prj is:
GEOGCS["GCS_WGS_1984",DATUM["D_WGS_1984",SPHEROID["WGS_1984",6378137.0,298.257223563]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]]
you might try exporting the data or possibly run it through a fix geometry (I use et geotools).
It looks to me lie with the new meridian you are splitting some polygons and that is causing the fill to "leak" out.
This was a proposed feature for QGIS which seems to have slipped down the priorities of the developers for the time being. However the feature request page outlines the syntax of the aux.xml file so it should be possible to write a simple python script to get all of the required parameters from the raster file and SpatialReference.org and write them to an xml file with the required filename.
XML parsing is a massive pain and generally it is best to use a pre-written parser such as beautifulsoup, however since this will write the same file structure every time it should be possible to manage in pure python.
Best Answer
The Google Mercator projection is usually bound for aereas between 85.0511° North and South.
If your data includes 90° North or South, the reprojection is mathematically not possible.
See also http://wiki.openstreetmap.org/wiki/Slippy_map_tilenames for all kinds of lat/lon to Google Mercator conversions.