我正在使用带有空间参考 WSG84 的 gdal2tiles.py，它工作正常。但是，当我更改为空间参考 EPSG:3440 时。我收到以下错误：
首先，我运行命令：
gdal_translate -of VRT -a_srs EPSG:3440 -gcp 0 0 58.068451479718924 23.65512391903488 -gcp 21816 0 58.126966134442846 23.65512391903488 -gcp 21816 14871 58.126966134442846 23.6185834507829 myinputimage.png myoutput.vrt
比，我运行命令：
python gdal2tiles.py myoutput.vrt
我收到以下错误：
有没有办法修复这个错误或使用没有这个错误的类似空间参考。请注意，此错误仅在使用 WSG84 以外的空间参考时发生。
编辑：我想问题可能是空间参考 EPSG:3440 没有被包括在内，但问题仍然存在，是否可以将其添加到 gdal 参考数据库中。
EPSG:3440 是一个投影系统。精度 0.5 米的正确定义可能是：
+proj=utm +zone=40 +ellps=clrk80 +towgs84=-180.624,-225.516,173.919,-0.81,-1.898,8.336,16.71006 +units=m +no_defs
如果我看一下http://epsg.io/3440-1439它会告诉我：
Projected bounds:
-35283.71 1840661.51
789562.66 2942956.58
在上面的命令中，您分配看起来像大地测量（纬度经度数字）的地面控制点，例如 58.068451479718924。这很可能是错误的 - 因为您分配的数字不在预计的范围内。
您必须在 gdal_translate 数字的 -gcp 参数中使用更大的数字 - 已经预测并在上述范围内。
可以通过位于 http://epsg.io/3440-1439/map的在线界面从 WGS84 转换您的纬度/经度，并在 -gcp 参数中分配转换后的数字。
maptiler您已使用引用http://www.maptiler.com/的标签标记您的查询。该工具还可以直观地分配坐标，而无需创建 VRT - 因此您可以单击图像并将上述坐标复制并粘贴到投影系统中。MapTiler 生成一个查看器，在 OpenLayers 中显示 EPSG:3440 瓦片，或将这些瓦片转换为球形墨卡托瓦片，以便与 Google 或 OSM 叠加或在移动应用程序中使用。
顺便说一句，很高兴听到您使用我们开发的系统 gdal2tiles、maptiler 和http ://epsg.io/...