使用 Python 更改图像投影

【问题标题】：Changing image projection using Python使用 Python 更改图像投影
【发布时间】：2011-05-31 10:33:20
【问题描述】：

图像角坐标已知：

左上角：(70.20315, 20.92749)
右上角：(70.20315, 31.45240)
左下角：(59.70571, 20.92749)
右下：(59.70571, 31.45240)

（基本上是芬兰和瑞典/挪威/俄罗斯的一点点地图）

我找不到将图像投影更改为 Google 地图叠加图像的方法，最好使用 Python。

使用pyproj 和PIL 可以做到这一点，但我无法弄清楚如何做到这一点。我还找到了GDAL，更具体地说是gdalwarp，但找不到以正确方式更改投影的选项。

【问题讨论】：

您能否展示您当前用于将图像叠加层放入 Google 地图的代码？如果你的 Python 反过来调用 JavaScript，那应该是可能的。

标签： python google-maps map map-projections

【解决方案1】：

（从gis答案复制粘贴）

大量的反复试验表明它在没有指定源投影的情况下工作正常。仅使用 GDAL Python API 也可以实现相同的目的。重现步骤：

给TIF添加控制点（当然是先convert mapimage.png mapimage.tif）：

gdal_translate mapimage.tif mapimage_georeffed.tif -gcp 2300 2300 31.45240 59.70571 -gcp 0 0 20.92749 70.20315 -gcp 2300 0 31.45240 70.20315 -gcp 0 2300 20.92749 59.70571 -a_srs +proj=longlat

然后，GDAL python 绑定包括gdal2tiles.py，它会在地图的右投影上自动生成切片，以及示例地图页面。

为了记录，我并不是绝对肯定的预测是完全正确的，但在那个领域它似乎至少足够好。

【讨论】：