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

Question

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

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

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

Answer 1

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

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 

Answer 2

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

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

Чтобы преобразовать целое использование шейпа файла: http://www.gdal.org/ogr2ogr.html

Там есть пример командной строки для шейпа в нижней части страницы ,