3

I have some TIFs in EPSG4326 the GeoTiff info looks like:

GEOGCS["WGS 84",
    DATUM["WGS_1984",
        SPHEROID["WGS 84",6378137,298.257223563,
            AUTHORITY["EPSG","7030"]],
        AUTHORITY["EPSG","6326"]],
    PRIMEM["Greenwich",0],
    UNIT["degree",0.0174532925199433],
    AUTHORITY["EPSG","4326"]]

How do convert them into EPSG:3857 so that the GeoTiff info reads something like:

PROJCS["Google Maps Global Mercator",
    GEOGCS["WGS 84",
        DATUM["WGS_1984",
            SPHEROID["WGS 84",6378137,298.257223563,
                AUTHORITY["EPSG","7030"]],
            AUTHORITY["EPSG","6326"]],
        PRIMEM["Greenwich",0],
        UNIT["degree",0.0174532925199433],
        AUTHORITY["EPSG","4326"]],
    PROJECTION["Mercator_1SP"],
    PARAMETER["central_meridian",0],
    PARAMETER["scale_factor",1],
    PARAMETER["false_easting",0],
    PARAMETER["false_northing",0],
    UNIT["metre",1,
        AUTHORITY["EPSG","9001"]]]
  • gdal_translate -of GTiff -s_srs epsg:4326 -t_srs epsg:3857 input.tif output.tif – user30184 Jul 2 '14 at 20:19
  • Sorry, of course gdalwarp. – user30184 Jul 3 '14 at 5:57
7

Use the gdalwarp utility.

gdalwarp -s_srs epsg:4326 -t_srs epsg:3857 input.tif output.tif
1

Use EPSG:3857. That is the official EPSG code number for Google mercator projection.

900913 was incorporated in GDAL some time ago, but is now dropped in favour of 3857.

  • Ok thanks, but how do I do the actual conversion to EPSG:3857 – nickponline Jul 2 '14 at 18:57
  • See the answer by Luke for the command line syntax. – AndreJ Sep 15 '16 at 18:20

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Not the answer you're looking for? Browse other questions tagged or ask your own question.