使用GDAL将罗宾逊投影转换为EPSG:3857投影问题求助

核心错误：投影转换参数颠倒 #

你在gdalwarp步骤中错误调换了源投影和目标投影的参数：

• 第一步gdal_translate操作已经给源影像标记了投影为ESRI:54030（罗宾逊投影），因此投影转换时-s_srs（指定源数据投影）参数应取值ESRI:54030，-t_srs（指定目标投影）参数应取值你需要的EPSG:3857
• 你之前的命令相当于把罗宾逊投影的源文件误判为墨卡托投影，再次转换为罗宾逊投影，自然只有轻微拉伸，和预期结果不符

补充注意事项 #

EPSG:3857（Web墨卡托）仅支持南北纬85.0511°之间的区域，转换后南北极超出有效范围的部分会被自动裁切，属于正常现象。

正确操作命令 #

1. 地理标记步骤（你当前的操作正确，无需修改）

gdal_translate -a_ullr -180 90 180 -90 -a_srs ESRI:54030 source.jpg source_tagged.tif

2. 投影转换步骤（修正参数顺序）

gdalwarp -s_srs ESRI:54030 -t_srs EPSG:3857 source_tagged.tif target.tif

如果需要自定义输出分辨率，可增加-tr参数（单位为米，符合EPSG:3857的投影单位），示例：

gdalwarp -s_srs ESRI:54030 -t_srs EPSG:3857 -tr 1000 1000 source_tagged.tif target.tif

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