如何使用开源库从脚本将 shapefile 从 AGD66 转换为 GDA94

Question

如何使用开源库以编程方式将 shapefile 从 AGD66 转换为 GDA94？我不想使用 arcgisscripting 因为我没有执照。

这需要是自动化的。bash 或 python 脚本是可以接受的。

这比普通重投影稍微复杂一些，因为这些坐标参考系之间使用不同的椭球体，因此需要使用畸变网格。

Answer 1

很抱歉回答自己的问题，但这适用于后代。使用 此处说明, ，首先下载畸变网格：

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

然后使用ogr2​​ogr重新投影：

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

转换整个ShapeFile的使用：http://www.gdal.org/ogr2ogr.html

页面底部有一个命令行示例。