Take the 2-minute tour ×
Geographic Information Systems Stack Exchange is a question and answer site for cartographers, geographers and GIS professionals. It's 100% free, no registration required.

I am using the GDAL API for Python for a script that takes only EPSG codes for inputs to specify source and target coordinate systems.

While the cs2cs command seems to handle quite well transformation between NAD27 (EPSG:4267) and NAD83 (EPSG:4326) coordinates in North America, I would like to have similar capabilities for AGD66 (EPSG:4202) to GDA94 (EPSG:4283) and AGD84 (EPSG:4203) to GDA94 (EPSG:4283) in Australia.

When using the -v option with cs2cs to look at the proj string being used i get the following:

    >cs2cs -v +init=epsg:4267
    ># ---- From Coordinate System ----
    >#Lat/long (Geodetic alias)
    >#
    ># +init=epsg:4267 +proj=longlat +datum=NAD27 +no_defs +ellps=clrk66
    ># +nadgrids=@conus,@alaska,@ntv2_0.gsb,@ntv1_can.dat
    ># ---- To Coordinate System ----
    >#Lat/long (Geodetic alias)
    >#
    ># +proj=latlong +datum=NAD27 +ellps=clrk66
    ># +nadgrids=@conus,@alaska,@ntv2_0.gsb,@ntv1_can.dat 

The nadgrids parameters tells cs2cs to use the best fitting datum grid shift file from conus, alaska, ntv2_0.gsb, and ntv1_can.dat.

For AGD66 (EPSG:4202) cs2cs returns:

    >cs2cs -v +init=epsg:4202
    ># ---- From Coordinate System ----
    >#Lat/long (Geodetic alias)
    >#
    ># +init=epsg:4202 +proj=longlat +ellps=aust_SA
    ># +towgs84=-117.808,-51.536,137.784,0.303,0.446,0.234,-0.29 +no_defs
    ># ---- To Coordinate System ----
    >#Lat/long (Geodetic alias)
    >#
    ># +proj=latlong +ellps=aust_SA
    ># +towgs84=-117.808,-51.536,137.784,0.303,0.446,0.234,-0.29

There is no nadgrids parameter but a towgs84 instead. A grid shift file from the ICSM; "A66 National (13.09.01).gsb" contains the shift between AGD66 and GDA94.

How can I make GDAL/cs2cs/proj identify this file when using EPSG codes?

I tried to reinstall GDAL/cs2cs/proj using the osgeo4w installer with A66 National (13.09.01).gsb within the directory containing the other datum shift files prior to the reinstallation but without success.

share|improve this question
    
Just for the record: We discussed the topic on the gdal-dev mailing list: osgeo-org.1560.x6.nabble.com/… –  AndreJ Feb 25 at 12:04

Your Answer

 
discard

By posting your answer, you agree to the privacy policy and terms of service.

Browse other questions tagged or ask your own question.