使用GDAL将Robinson转换为EPSG:3857
使用GDAL将Robinson转换为EPSG:3857
提问于 2021-10-12 17:57:41
我在试着把罗宾逊投影转换成墨卡托。例如，我使用这个image。
首先，我应用地理标记：
gdal_translate -a_ullr -180 90 180 -90 -a_srs ESRI:54030 source.jpg source_tagged.tif
并最终将其传递给墨卡托：
gdalwarp -t_srs ESRI:54030 -s_srs EPSG:3857 source_tagged.tif target.tif
结果在垂直方向上略有延伸，但与适当的墨卡托相去甚远。我做错了什么？
回答 1
Stack Overflow用户
回答已采纳
发布于 2021-10-13 22:02:20
您的命令中有几个问题。第一个是-180 90不在全局Robinson投影GeoTIFF的左上角像素中。左上角将类似于-338.2187147689 90，右下角将是338.2187147689 -90。但是，您还将srs指定为ESRI:54030，因此这些边界需要使用投影坐标，而不是经度/经度。从您的映像生成GeoTIFF的命令为：
gdal_translate -a_ullr -17005833.3305252 8625154.47184994 17005833.3305252 -8625154.47184994 -a_srs ESRI:54030 source.jpg source_tagged.tif您的第二个命令切换了-t_srs和-s_srs。考虑到您要投影到EPSG:3857，您还需要提供界限，因为墨卡托在两极是无穷大的。因此，更新后的命令将如下所示：
gdalwarp -s_srs ESRI:54030 -t_srs EPSG:3857 -te -180 -81 180 81 -te_srs EPSG:4326 source_tagged.tif target.tif我使用这些命令从您的示例中生成了以下投影图像(在运行它们之前裁剪图像中的空格)。
页面原文内容由Stack Overflow提供。腾讯云小微IT领域专用引擎提供翻译支持
原文链接：
https://stackoverflow.com/questions/69545017
复制相关文章
相似问题
腾讯云开发者
