<div dir="ltr">Even,<div><br></div><div>Thank you for catching the quote problem. I've fixed my WKT accordingly.</div><div><br></div><div>Unfortunately altering the SRS to GDAL's default EPSG:3857 definition results in the same alignment problem when viewed in ArcGIS Desktop/Server, so I altered it back to my custom WKT (using double quotes this time) - now it's aligning properly again, but I'm still seeing the same error:</div>
<div><br></div><div> ERROR 6: No translation for Mercator_Auxiliary_Sphere to PROJ.4 format is known.</div><div><br><div class="gmail_extra">It still produces valid output with 1.10.0, but with 1.10.1 it still exits prematurely with status 1 and produces an all-zero image.<br clear="all">
<div><div dir="ltr"><br></div></div><div>It seems that something was changed between 1.10.0 and 1.10.1 that made this a fatal error.</div><div><br></div><div>Thanks,</div><div>Jason</div>
<br><br><div class="gmail_quote">On Wed, Jan 22, 2014 at 3:07 PM, Even Rouault <span dir="ltr"><<a href="mailto:even.rouault@mines-paris.org" target="_blank">even.rouault@mines-paris.org</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex">
Le mercredi 22 janvier 2014 20:27:19, Jason Greenlaw - NOAA Affiliate a écrit :<br>
<div>> When attempting to use a polygon shapefile (in custom web mercator<br>
> projection) with the gdalwarp cutline operation, operating on a GeoTIFF<br>
> (with same projection), the following error is output in both versions<br>
> 1.10.0 and 1.10.1:<br>
><br>
> ERROR 6: No translation for 'Mercator_Auxiliary_Sphere' to PROJ.4<br>
> format is known.<br>
<br>
</div>The issue is that esri-3857.prf contains single quote characters instead of<br>
double quotes, and oceans_glakes_50m_3857.prj contain single quotes inside<br>
double quotes. Both are not valid AFAIK<br>
<br>
You should be able to overcome this with :<br>
<br>
ogr2ogr -a_srs EPSG:3857 oceans_glakes_50m_3857.shp<br>
oceans_glakes_50m_3857_fixed.shp<br>
<br>
and<br>
<br>
gdalwarp -s_srs EPSG:4326 -t_srs EPSG:3857 -r bilinear -srcnodata<br>
<div><div> 9999.0 -dstnodata 9999.0 -tr 4000 9000 -overwrite -co COMPRESS=LZW -co<br>
TILED=yes sst.flt sst_projected.tif<br>
<br>
><br>
> With GDAL 1.10.0, the operation still completes successfully with a zero<br>
> exit status, but with GDAL 1.10.1 the operation fails with a nonzero exit<br>
> status, and the resulting GeoTIFF has zeros for all pixels.<br>
><br>
> Here is the command I'm using:<br>
><br>
> gdalwarp -co COMPRESS=LZW -cutline oceans_glakes_50m_3857.shp<br>
> -srcnodata 9999.0 -dstnodata 9999.0 -overwrite sst_global.tif<br>
> sst_global_cut.tif<br>
><br>
><br>
> The input shapefile was produced by reprojecting with ogr2ogr from<br>
> EPSG:4326 to my custom projection (WKT is included below), which is stored<br>
> in a *.prf file. The source GeoTIFF is single-band with Float32 pixel<br>
> type.<br>
><br>
> The custom WKT can be found here: <a href="http://jgreenlaw.org/gdal/esri-3857.prf" target="_blank">http://jgreenlaw.org/gdal/esri-3857.prf</a><br>
><br>
> The cutline shapefile can be downloaded at:<br>
> <a href="http://jgreenlaw.org/gdal/oceans_glakes_50m_3857.tar.gz" target="_blank">http://jgreenlaw.org/gdal/oceans_glakes_50m_3857.tar.gz</a> (900kb)<br>
><br>
> The source GeoTIFF can be generated from this FLT/HDR raster:<br>
> <a href="http://jgreenlaw.org/gdal/sst_flt.tar.gz" target="_blank">http://jgreenlaw.org/gdal/sst_flt.tar.gz</a> (11mb) by running these commands:<br>
><br>
> gdalwarp -s_srs EPSG:4326 -t_srs esri-3857.prf -r bilinear -srcnodata<br>
> 9999.0 -dstnodata 9999.0 -tr 4000 9000 -overwrite -co COMPRESS=LZW -co<br>
> TILED=yes sst.flt sst_projected.tif<br>
><br>
> gdalwarp -co COMPRESS=LZW -te -20037507.0671618 -19971868.8804086<br>
> 20037507.0671618 19971868.8804086 -srcnodata 9999.0 -dstnodata 9999.0 -tr<br>
> 4000 9000 sst_projected.tif sst_global.tif<br>
><br>
> Does anyone know if this behavior is to be expected when the above error is<br>
> output? Is there a better method for reprojecting to Web Mercator that<br>
> will align properly in ArcGIS?<br>
><br>
> Thanks<br>
> -Jason<br>
><br>
><br>
> Background:<br>
> The reason for using custom WKT for reprojecting to Web Mercator is because<br>
> I have experienced alignment problems displaying GDAL-produced GeoTIFFs in<br>
> ArcGIS when reprojecting using "EPSG:3857" as the target SRS (this seems to<br>
> be related to <a href="http://trac.osgeo.org/gdal/ticket/3962" target="_blank">http://trac.osgeo.org/gdal/ticket/3962</a>).<br>
><br>
> I generated the WKT below by copying the 102113 definition in<br>
> $GDAL_DATA/esri_extra.wkt and merging information from an ESRI Web Mercator<br>
> *.prj file. Using this custom WKT as the target SRS with gdalwarp normally<br>
> results in rasters that align properly with ArcGIS Web Mercator data<br>
> (without on-the-fly datum transformation).<br>
><br>
><br>
> ---------------------------------------------------------------------------<br>
</div></div>> ------------------- Custom WKT for reprojecting to ESRI Web Mercator using<br>
<div><div>> GDAL (<br>
> <a href="http://jgreenlaw.org/gdal/esri-3857.prf" target="_blank">http://jgreenlaw.org/gdal/esri-3857.prf</a>):<br>
> ---------------------------------------------------------------------------<br>
> -------------------<br>
><br>
> PROJCS['WGS_1984_Web_Mercator_Auxiliary_Sphere',<br>
> GEOGCS['GCS_WGS_1984',<br>
> DATUM['D_WGS_1984',<br>
> SPHEROID['WGS_1984',6378137.0,298.257223563]<br>
> ],<br>
> PRIMEM['Greenwich',0.0],<br>
> UNIT['Degree',0.0174532925199433]<br>
> ],<br>
> PROJECTION['Mercator_Auxiliary_Sphere'],<br>
> PARAMETER['False_Easting',0.0],<br>
> PARAMETER['False_Northing',0.0],<br>
> PARAMETER['Central_Meridian',0.0],<br>
> PARAMETER['latitude_of_origin',0.0],<br>
> PARAMETER['Standard_Parallel_1',0.0],<br>
> PARAMETER['Auxiliary_Sphere_Type',0.0],<br>
> UNIT['Meter',1.0],<br>
> EXTENSION["PROJ4",<br>
> "+proj=merc +a=6378137 +b=6356752.314245179 +lat_ts=0.0 +lon_0=0.0<br>
> +x_0=0.0 +y_0=0.0 +k=1.0 +units=m +nadgrids=@null +wktext +no_defs"<br>
> ],<br>
> AUTHORITY["EPSG","3857"]<br>
> ]<br>
<br>
</div></div><span><font color="#888888">--<br>
Geospatial professional services<br>
<a href="http://even.rouault.free.fr/services.html" target="_blank">http://even.rouault.free.fr/services.html</a><br>
</font></span></blockquote></div><br></div></div></div>