Sometimes (eg. doing it for a bulk of geodatasets) converting and transforming geodatasets from their Source-CRS to another with command line tools like ogr2ogr (http://www.gdal.org/ ) can be helpful.
With ogr2ogr it’s also an easy task to use NTv2-based transformations like for Austria the AT_GIS_GRID.
The following syntax did a good job in my case – the results have been identical using ArcMap or QGIS.
Before you start download the AT_GIS_GRID from the BEV (http://www.bev.gv.at).
The basic syntax for a geodataset based on the Austrian MGI (in this for a dataset in GK/West – EPSG:31254):
ogr2ogr -f "ESRI Shapefile" -s_srs "+proj=tmerc +lat_0=0 +lon_0=10.33333333333333 +k=1 +x_0=0 +y_0=-5000000 +ellps=bessel +towgs84=577.326,90.129,463.919,5.137,1.474,5.297,2.4232 +units=m +no_defs +nadgrids=AT_GIS_GRID.gsb" -t_srs "EPSG:4258" Targetfile.shp Sourcefile.shp
In my case using ogr2ogr worked really good and FAST ! The results did not differ from a reference dataset transformed with ArcMap.
Depending on the Source-SRS replace the “-s_srs”-Parameters. At the end of the article the most common are listed.
If you don’t like working with the command line use AT_GIS_GRID Transformations with…
- AT_GIS_GRID with ArcMap 10.2.x
- AT_GIS_GRID with QGIS 2.x and ArcGIS 10.0
- AT_GIS_GRID with the QGIS-NTv2 Plugin
The most common “s_srs”-definitions for Austria (for those who like Copy&Paste, but don’t forget to integrate the AT_GIS_GRID in the syntax)
EPSG:31254 (Austria GK West)
+proj=tmerc +lat_0=0 +lon_0=10.33333333333333 +k=1 +x_0=0 +y_0=-5000000 +ellps=bessel +towgs84=577.326,90.129,463.919,5.137,1.474,5.297,2.4232 +units=m +no_defs
EPSG:31255 (Austria GK Central)
+proj=tmerc +lat_0=0 +lon_0=13.33333333333333 +k=1 +x_0=0 +y_0=-5000000 +ellps=bessel +towgs84=577.326,90.129,463.919,5.137,1.474,5.297,2.4232 +units=m +no_defs
EPSG:31256 (Austria GK East)
+proj=tmerc +lat_0=0 +lon_0=16.33333333333333 +k=1 +x_0=0 +y_0=-5000000 +ellps=bessel +towgs84=577.326,90.129,463.919,5.137,1.474,5.297,2.4232 +units=m +no_defs
EPSG:31257 (Austria GK M28)
+proj=tmerc +lat_0=0 +lon_0=10.33333333333333 +k=1 +x_0=150000 +y_0=-5000000 +ellps=bessel +towgs84=577.326,90.129,463.919,5.137,1.474,5.297,2.4232 +units=m +no_defs
EPSG:31258 (Austria GK M31)
+proj=tmerc +lat_0=0 +lon_0=13.33333333333333 +k=1 +x_0=450000 +y_0=-5000000 +ellps=bessel +towgs84=577.326,90.129,463.919,5.137,1.474,5.297,2.4232 +units=m +no_defs
EPSG:31259 (Austria GK M34)
+proj=tmerc +lat_0=0 +lon_0=16.33333333333333 +k=1 +x_0=750000 +y_0=-5000000 +ellps=bessel +towgs84=577.326,90.129,463.919,5.137,1.474,5.297,2.4232 +units=m +no_defs