Как я могу преобразовать шейп-файл из AGD66 в GDA94 из скрипта, используя библиотеки с открытым исходным кодом - PullRequest
Новая
       34

Как я могу преобразовать шейп-файл из AGD66 в GDA94 из скрипта, используя библиотеки с открытым исходным кодом

0 голосов
/ 04 июня 2010

Как я могу преобразовать шейп-файл из AGD66 в GDA94 программно, используя библиотеки с открытым исходным кодом?Я не хочу использовать arcgisscripting, потому что у меня нет лицензии.

Это должно быть автоматизировано.Сценарий bash или python был бы приемлем.

Это немного сложнее, чем обычное перепроецирование, потому что между этими системами координат используется другой эллипсоид, и поэтому необходимо использовать сетку искажения.

Ответы [ 2 ]

2 голосов
/ 04 июня 2010

Используйте библиотеки OpenSource GDAL.

Для преобразования координат самостоятельно используется утилита преобразования GGAL: http://www.gdal.org/gdaltransform.html

Для преобразования всего шейп-файла используйте: http://www.gdal.org/ogr2ogr.html

В нижней части страницы приведен пример командной строки для шейп-файла.

0 голосов
/ 04 июня 2010

Извините, что отвечаю на свой вопрос, но здесь речь идет о потомстве. Используя инструкции здесь , сначала загрузите сетку искажения: 

wget http://www.icsm.gov.au/icsm/gda/gdatm/national66.zip
unzip national66.zip
mkdir -p ~/bin
mv "A66 National (13.09.01).gsb" ~/bin/a66_national.gsb
rm national66.zip

Затем используйте ogr2ogr для перепроектирования: 

ogr2ogr -f "ESRI Shapefile" -s_srs "+proj=longlat +ellps=aust_SA +nadgrids=~/bin/a66_national.gsb +wktext" -t_srs EPSG:4283 outputgda94.shp inputagd66.shp
...