Преобразование координат из SVY21 в WGS84

Я пытаюсь преобразовать имеющиеся у меня координаты SVY21 в систему координат WGS84.

Я попытался использовать программное обеспечение GDAL и вставил следующее в команду после установки: ogr2ogr -f GeoJSON –t_srs WGS84 <jsonFileName>.json <shapefile>.shp

Это приводит к сбою:

Невозможно открыть источник данных 'WGS84' со следующими драйверами

И драйверы перечислены.

Также я попытался использовать пакет python utm для преобразования координаты utm N48 в WGS48, но мне не удалось преобразовать координаты SVY21 в эталонную точку ut48 N48.

Я нашел предложенное решение в этом старом посте запутанным с точки зрения терминологии, и я предпочел бы использовать какое-то существующее решение.

Есть ли "простое" решение, которое использует мои предыдущие подходы или быстро применимо к моей проблеме? Я также был бы рад, если бы кто-то мог дать мне более точное объяснение подхода, описанного в предыдущем посте.

Заранее спасибо.

Редактировать: я пытался написать скрипт Python, который реализует прямой расчет Vicenty из "старого поста", упомянутого ранее. Было бы также полезно, если бы кто-то мог подтвердить, что вычесть ложные координаты происхождения проекции (28001,642 мЕ,38744,572 мН) из необработанного svy21 (e,n) и использовать начало проекции (точка без отметки) (103° 50' 00", 1° 22' 00") правильно.

1 ответ

Решение

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

  1. Вы должны скопировать файлы *.shp (и *.shx) рядом с командой ogr2ogr (Источник: https://gis.stackexchange.com/questions/56652/how-to-move-and-open-shapefiles-with-ogr2ogr)

  2. gdal не сразу знает "EPSG:3414" - имя EPSG для SCY21, поэтому дайте ему подсказку, это в [УСТАНОВИТЬ РАСПОЛОЖЕНИЕ]/data (я не знаю, где он находится в установленной версии, особенно в Windows. Найдите файл с именем "gcs.csv" и укажите GDAL_LOCATION на содержащую его папку.

  3. определить t_srs и s_srs соответственно

  4. Определите формат вывода ("GEOJson")

Summa summarum:

GDAL_DATA="../data" ./ogr2ogr -f "GEOJson" "LaneMarking_wgs84.json" "LaneMarking.shp" -t_srs WGS84 -s_srs "EPSG:3414"

Повеселись:)

Другие вопросы по тегам