I am trying to reproject a netCDF file which is in Lat Lon (EPSG:4326) projection to EPSG:3857 (Google mercator) (on GDAL 3.4.3, released 2022/04/22). If I reproject to a tif file, it is reprojected correctly. If I reproject to a nc file it reproject incorrectly/different location that tif file. I am giving the commands I ran below. Correctness is determined by viewing them in QGIS: (All tif files show up at exactly same location as a_stamped_latlon.nc. agoog.nc show up at a different location ).
I will be happy to post the sample files [The only way I know to send files to open a bug, but very likely its some details I don't know about different formats].: a.nc is only 42KB. Here is how to reproduce it. gdalwarp -s_srs EPSG:4326 -t_srs EPSG:3857 -co COMPRESS=DEFLATE -r cubic a.nc agoog.nc gdalwarp -s_srs EPSG:4326 -t_srs EPSG:3857 -co COMPRESS=DEFLATE -r cubic -co "GEOTIFF_KEYS_FLAVOR=STANDARD" agoog_std.tif gdalwarp -s_srs EPSG:4326 -t_srs EPSG:3857 -co COMPRESS=DEFLATE -r cubic -co "GEOTIFF_KEYS_FLAVOR=STANDARD" a.nc agoog_std.tif gdalwarp -s_srs EPSG:4326 -t_srs EPSG:3857 -co COMPRESS=DEFLATE -r cubic -co "GEOTIFF_KEYS_FLAVOR=ESRI_PE" a.nc agoog_esri.tif gdalwarp -s_srs EPSG:4326 -t_srs EPSG:3857 -co COMPRESS=DEFLATE -r cubic -co "GEOTIFF_KEYS_FLAVOR=STANDARD" a.nc agoog_std.tif gdal_translate -a_srs EPSG:4326 a.nc a_stamped_latlon.nc GDAL 3.4.3, released 2022/04/22
_______________________________________________ gdal-dev mailing list gdal-dev@lists.osgeo.org https://lists.osgeo.org/mailman/listinfo/gdal-dev
