¿Incorporar una proyección ArcGIS compatible con ESRI mediante GDAL?

Question

¿Existe alguna forma de incrustar una proyección ArcGIS compatible con ESRI utilizando las herramientas GDAL?

Actualmente estoy usando gdalwarp para crear archivos raster que traducen raster de WGS84 (EPSG:4326) a Spherical Mercator (EPSG:3857). También crea Geotiff con etiquetas y tfw -file.

gdalwarp -q -overwrite -s_srs '+init=epsg:4326' -t_srs '+init=epsg:3857' -r bilinear -of GTiff -co COMPRESS=LZW -co TILED=YES -co TFW=YES...

Todavía no se abre correctamente al abrir en ESRI ArcGIS Desktop 10.4 pero complining falta proyección. ¿Es posible crear un archivo raster que se abra directa y correctamente en el entorno ArcGIS y que se pueda producir utilizando las herramientas GDAL?

Answer 1

Si el archivo .tfw no existe, ArcGIS se ve obligado a leer la georreferenciación a partir de las etiquetas GeoTIFF y luego, con suerte, lee e interpreta también las etiquetas de proyección.

Si ArcGIS no puede interpretar correctamente la proyección a partir de las etiquetas GeoTIFF, puede crear un archivo ESRI .prj con gdalsrsinfo http://www.gdal.org/gdalsrsinfo.html

gdalsrsinfo -o wkt_esri epsg:3857

PROJCS["WGS_1984_Web_Mercator_Auxiliary_Sphere",GEOGCS["GCS_WGS_1984",DATUM["D_WGS_1984",SPHEROID["WGS_1984",6378137.0,298.257223563]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Mercator_Auxiliary_Sphere"],PARAMETER["False_Easting",0.0],PARAMETER["False_Northing",0.0],PARAMETER["Central_Meridian",0.0],PARAMETER["Standard_Parallel_1",0.0],PARAMETER["Auxiliary_Sphere_Type",0.0],UNIT["Meter",1.0]]