将GDALwarp与非比例尺地图一起使用

问题描述

我的输入地图是不按比例缩放地图的JPG扫描。我想将其变形为EPSG：3857，以便可以强制扩展。结果不像我期望的那样翘曲。

采用的方法是将9个分散的GCP（地面控制点）添加到新的TIF文件中，但否投影。 GCP已添加为X，Y，Longitide，Latitude，Z = 0。 gdal_edit（从Python）执行此步骤没有错误，并且gdalinfo in.tif现在显示：

    Size is 3421,1925
Coordinate System is `'
GCP Projection =
GCP[  0]: Id=1,Info=
          (462,111) -> (-1.541,50.759,0)
GCP[  1]: Id=2,Info=
          (732,657) -> (-1.496,50.704,0)
GCP[  2]: Id=3,Info=
          (1998,1764) -> (-1.285,50.588,0)
GCP[  3]: Id=4,Info=
          (1911,687) -> (-1.291,50.701,0)
GCP[  4]: Id=5,Info=
          (1896,108) -> (-1.3,50.762,0)
GCP[  5]: Id=6,Info=
          (2691,360) -> (-1.158,50.732,0)
GCP[  6]: Id=7,Info=
          (2739,1098) -> (-1.149,50.659,0)
GCP[  7]: Id=8,Info=
          (3108,756) -> (-1.084,50.689,0)
GCP[  8]: Id=9,Info=
          (2163,1320) -> (-1.254,50.635,0)

后续的warp命令为：

gdalwarp -of GTiff -t_srs EPSG:3857 in.tif out-warped-3857.tif

它运行，但是生成的地图被错误地拉​​伸，从3421x1925扩展到3852x1381。将输出投影更改为EPSG：4326，结果与我所期望的完全相同。

有三幅图像可以说明：

参考图

输入地图

输出图

问题：这是添加投影并强制地图缩放的最佳方法吗？

解决方法

暂无找到可以解决该程序问题的有效方法，小编努力寻找整理中！

如果你已经找到好的解决方法，欢迎将解决方案带上本链接一起发送给小编。

小编邮箱:dio#foxmail.com (将#修改为@）