I've recently run into some issues when trying to reproject rasters and shapefiles using GDAL to the ESRI Web Mercator (auxiliary sphere) projection (EPSG:3857), but have come up with a workaround that I wanted to share.  <div>
<br></div><div>After specifying "EPSG:3857" as the target projection to gdalwarp or ogr2ogr, the resulting dataset does not overlay nicely in ArcMap/ArcGIS Server with ESRI-projected datasets or ArcGIS Online basemaps.  If you instead use the code "EPSG:102113", it results in the data overlaying correctly, but requires an on-the-fly datum transformation when drawing in EPSG:3857.  I've tried specifying an ESRI prj file along with the "ESRI::" prefix to gdalwarp, but that too failed to reproject properly.</div>
<div><br></div><div>To avoid a datum transformation, I've merged the information from an ESRI 3857 .prj file with the EPSG:102113 definition found in $GDAL_DATA/esri_extra.wkt (specifically modifying the PROJ4 extension definition to use a semiminor axis of 6356752.314245179, which is what ESRI expects), and created a .prf file, which I've been able to use successfully with gdalwarp and ogr2ogr.</div>
<div><br></div><div>The WKT of the .prf file is as follows:</div><div><br></div><div>---------------------------------------------------------------------------<br><div><div><div>PROJCS['WGS_1984_Web_Mercator_Auxiliary_Sphere',</div>
<div>    GEOGCS['GCS_WGS_1984',</div><div>        DATUM['D_WGS_1984',</div><div>            SPHEROID['WGS_1984',6378137.0,298.257223563]</div><div>        ],</div><div>        PRIMEM['Greenwich',0.0],</div>
<div>        UNIT['Degree',0.0174532925199433]</div><div>    ],</div><div>    PROJECTION['Mercator_Auxiliary_Sphere'],</div><div>    PARAMETER['False_Easting',0.0],</div><div>    PARAMETER['False_Northing',0.0],</div>
<div>    PARAMETER['Central_Meridian',0.0],</div><div>    PARAMETER['latitude_of_origin',0.0],</div><div>    PARAMETER['Standard_Parallel_1',0.0],</div><div>    PARAMETER['Auxiliary_Sphere_Type',0.0],</div>
<div>    UNIT['Meter',1.0],</div><div>    EXTENSION["PROJ4",</div><div>        "+proj=merc +a=6378137 +b=6356752.314245179 +lat_ts=0.0 +lon_0=0.0 +x_0=0.0 +y_0=0.0 +k=1.0 +units=m +nadgrids=@null +wktext +no_defs"</div>
<div>    ],</div><div>    AUTHORITY["EPSG","3857"]</div><div>]</div></div><div>---------------------------------------------------------------------------</div><div><br></div><div>Hope this solution saves someone else some time.</div>
<div>Jason</div><div><br clear="all"><div>--</div>Jason Greenlaw<div>Contractor, ERT, Inc.<br><div>NOAA/NOS/OCS/CSDL<div><a href="http://nowcoast.noaa.gov" target="_blank">http://nowcoast.noaa.gov</a></div><div><a href="mailto:Jason.Greenlaw@noaa.gov" target="_blank">Jason.Greenlaw@noaa.gov</a></div>
</div></div><br>
</div></div></div>