如何用一条GDAL命令将EPSG:4326的PNG瓦片转成EPSG:3857

当然可以！你完全不需要分成两步操作，直接用一条gdalwarp命令就能完成从无地理参考的EPSG:4326 PNG瓦片到EPSG:3857投影瓦片的转换，还能省去生成中间GeoTIFF文件的步骤。

合并后的命令 #

gdalwarp -of PNG -a_ullr 135.00000000000003 36.59788913307022 140.62500000000003 31.952162238024975 -a_srs EPSG:4326 -t_srs EPSG:3857 4326.png 3857.png

参数说明 #

• -of PNG：指定输出文件格式为PNG，对应你原来第二步的输出格式需求
• -a_ullr：给输入的PNG瓦片赋予左上角（Upper Left）和右下角（Lower Right）的地理坐标，参数值和你原来gdal_translate里的完全一致
• -a_srs EPSG:4326：明确指定输入文件的原始坐标系为EPSG:4326
• -t_srs EPSG:3857：设置目标坐标系为EPSG:3857，完成投影转换
• 末尾的4326.png是输入文件路径，3857.png是输出文件路径

这条命令相当于把gdal_translate赋予地理参考的操作和gdalwarp的投影转换操作合并到了一起，一步到位完成需求。

内容的提问来源于stack exchange，提问作者UsagiLabo