使用 Python 更改图像投影

发布于 2024-11-10 18:56:42 字数 770 浏览 4 评论 0原文

图像角点坐标已知:

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

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

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

可以使用 pyproj 和 PIL,但我不知道该怎么做。我还找到了 GDAL,更具体地说 gdalwarp,但找不到以正确方式更改投影的选项。

示例图片

Image corner coordinates are known:

  • Top left: (70.20315, 20.92749)
  • Top right: (70.20315, 31.45240)
  • Bottom left: (59.70571, 20.92749)
  • Bottom right: (59.70571, 31.45240)

(Basically map over Finland and a little bit of Sweden/Norway/Russia)

I couldn't find a way to change image projection to Google maps overlay image, preferably using Python.

It's possible to do that with pyproj and PIL, but I wasn't able to figure out how to do that. I also found GDAL, and more specifically gdalwarp, but couldn't find options for changing projection in a correct way.

Example image

如果你对这篇内容有疑问,欢迎到本站社区发帖提问 参与讨论,获取更多帮助,或者扫码二维码加入 Web 技术交流群。

扫码二维码加入Web技术交流群

发布评论

需要 登录 才能够评论, 你可以免费 注册 一个本站的账号。

评论(1

草莓酥 2024-11-17 18:56:42

(从 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,它会自动生成图块地图的右投影以及示例地图页面。

根据记录,我并不绝对肯定的预测完全正确,但在该领域,它似乎至少足够好。

(Copy-paste from gis answer)

Extensive trial-and-error showed it worked fine without specifying source projection. Same can be achieved using only GDAL Python API. Steps to reproduce:

Adding control points to TIF (of course, first 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

Then, GDAL python bindings include gdal2tiles.py, which autogenerates tiles on right projection for maps, and also sample map pages.

For the record, I'm not absolutely positive projection is exactly correct, but over that area it seems to be at least good enough.

~没有更多了~
我们使用 Cookies 和其他技术来定制您的体验包括您的登录状态等。通过阅读我们的 隐私政策 了解更多相关信息。 单击 接受 或继续使用网站,即表示您同意使用 Cookies 和您的相关数据。
原文