如何解决将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 举报,一经查实,本站将立刻删除。