python - 使用 Python 更改图像投影

Question

图像角坐标已知：

• 左上角：（70.20315、20.92749）
• 右上角：（70.20315、31.45240）
• 左下：（59.70571、20.92749）
• 右下角：（59.70571、31.45240）

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

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

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

Answer 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自动生成地图右侧投影的切片，以及示例地图页面。

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