Ошибка при использовании пространственной ссылки EPSG: 3440 в gdal2tiles.py

Question

Я использовал gdal2tiles.py с пространственной ссылкой WSG84, и он работал нормально. Однако, когда я изменил пространственную ссылку 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.

EDIT: Я думаю, проблема может заключаться в том, что пространственная ссылка EPSG: 3440 не используется, но вопрос остается, можно ли добавить его в базу данных ссылок gdal.

Answer 1

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. Это, скорее всего, неверно - поскольку присваиваемые вами номера не находятся в диапазоне прогнозируемых границ.

Вы должны очень использовать в параметре -gcp числа gdal_translate, которые больше - уже проецируются и находятся в указанном выше диапазоне.

Можно преобразовать ваш lat/lon из WGS84 через онлайн-интерфейс по адресу: http://epsg.io/3440-1439/map и назначить в аргументе -gcp преобразованные числа.

Вы пометили свой запрос тегом maptiler, который ссылается на http://www.maptiler.com/. Этот инструмент также сможет назначать координаты визуально, без необходимости создавать VRT - так что вы можете щелкнуть по своему изображению и скопировать & вставить указанные выше координаты в проецируемой системе. MapTiler создает зрителя, который отображает плитки EPSG: 3440 в OpenLayers или преобразует их в сферические плиты мерцателя для наложений с Google или OSM или для использования в мобильных приложениях.

BTW рад слышать, что вы использовать системы gdal2tiles, maptiler и http://epsg.io/, которые мы разработали ...