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

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

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

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

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

И драйвера указаны.

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

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

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

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

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


person MNue    schedule 26.10.2016    source источник
comment
Можете ли вы опубликовать точную команду ogr2ogr, которую вы используете? Я думаю, что у вас есть синтаксис, который отличается от того, что вы опубликовали. В ошибке говорится, что вы пытаетесь открыть набор данных с именем WGS84, но опубликованная вами команда не пытается это сделать.   -  person Logan Byers    schedule 26.10.2016
comment
Я загрузил скриншот моей команды. Я действительно не вижу разницы.   -  person MNue    schedule 26.10.2016
comment
Вы уверены, что между файлами .json и .shp есть пробел?   -  person Logan Byers    schedule 26.10.2016
comment
...filename1.jsonxfilename2.shp В команде есть пробел, где x находится в этой его части.   -  person MNue    schedule 26.10.2016


Ответы (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»)

Подведение итогов:

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

Развлекайся :)

person EGOrecords    schedule 27.10.2016