Comment puis-je convertir un fichier de formes d'AGD66 en GDA94 à partir d'un script, en utilisant des bibliothèques open source

Question

Comment convertir un fichier de formes d'AGD66 en GDA94 par programmation, en utilisant des bibliothèques open source? Je ne veux pas utiliser arcgisscripting parce que je n'ai pas de licence.

Ceci doit être automatisable. Un script bash ou python serait acceptable.

Ceci est un peu plus compliqué qu'une reprojection normale, car un ellipsoïde différent est utilisé entre ces systèmes de référence de coordonnées, et donc une grille de distorsion doit être utilisée.

Answer 1

Désolé de répondre à votre question, mais ici, c'est pour la postérité. Utilisation de la instructions here, d'abord télécharger la grille de distorsion:

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 

Ensuite, utilisez ogr2ogr reprojeter:

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

Utilisez les bibliothèques OpenSource GDAL.

Pour convertir les coordonnées de leur propre chef nous le GGAL transformons utilitaire: http://www.gdal.org/gdaltransform.html

Pour convertir une utilisation toute shapefile: http://www.gdal.org/ogr2ogr.html

Il est un exemple de ligne de commande pour un shapefile au bas de la page .