<div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div>Hi all,</div><div><br></div><div>I figured I should ask here too if anyone know what the correct projection is for the Sandwell .img grids from <a href="http://topex.ucsd.edu">http://topex.ucsd.edu</a>.  I'm trying to keep the files in their original projections as I switch them to geotiffs.  I've asked David Sandwell directly too if he knows.</div><div><br></div><div>This works, but warps the data first:</div><div><br></div><div><div>gmt img2grd -V -R-180/180/-80/80 topo_18.1.img  -Gtest2.grd -T1 -D</div><div>gdal_translate -a_srs EPSG:4326 test2.grd test2.tif</div><div>gdalinfo test2.tif  # Results look believable</div><div>gdal_translate topo-18-1-epsg4326.tif topo-18-1-epsg4326-deflate.tif -co COMPRESS=DEFLATE -co PREDICTOR=3</div><div><br></div><div>Then imported into QGIS or Earth Engine as a normal user, things line up pretty well.</div><div><br></div><div>    <a href="https://code.earthengine.google.com/f43c8b13bafa26fd8c7f83ce9a919f6e">https://code.earthengine.google.com/f43c8b13bafa26fd8c7f83ce9a919f6e</a></div></div><div><br></div><div><br></div><div>I'd rather do it more like this:</div><div><br></div><div><div>gmt img2grd -V -R-180/180/-80/80 topo_18.1.img  -Gspherical-mercator-proj.grd -T1 -D -M</div><div>img2grd: Expects topo_18.1.img to be 21600 by 17280 pixels spanning 0/360.0/-80.738009/80.738009.</div><div>img2grd: To fit [averaged] input, your topo_18.1.img is adjusted to -R-180/180/-80.0023237126/80.0023237126.</div><div>img2grd: The output grid size will be 21600 by 16752 pixels.</div><div>img2grd: Created 21600 by 16752 Mercatorized grid file.  Min, Max values are -10914  8550</div><div><br></div><div>gives this which doesn't work as is:</div><div><br></div><div>gdalinfo spherical-mercator-proj.grd</div><div>Driver: netCDF/Network Common Data Format</div><div>Files: spherical-mercator-proj.grd</div><div>Size is 21600, 16752</div><div>Coordinate System is `'</div><div>Origin = (0.000000000000000,279.199999999999989)</div><div>Pixel Size = (0.016666666666667,-0.016666666666667)</div><div>Metadata:</div><div>  NC_GLOBAL#Conventions=COARDS, CF-1.5</div><div>  NC_GLOBAL#description=Spherical Mercator Projected with -Jm1 -R0/360/-80.0023237126/80.0023237126</div><div>  NC_GLOBAL#GMT_version=5.4.3 (r19528) [64-bit]</div></div><div><br></div><div>Should I be using one of 54004 or 41001?  e.g. <a href="https://epsg.io/?q=spherical+mercator">https://epsg.io/?q=spherical+mercator</a></div><div><br></div><div>Thanks!</div><div>-kurt</div><div><br></div><div dir="ltr" class="gmail_signature"></div></div></div></div></div>