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 wrote a Python script reprojecting points from a source coordinate reference system (SCRS) to a target one (TCRS) both identified by an EPSG code. I would like my python script to retrieve details about the grid shift file it uses when doing reprojections from NAD83 to NAD27.

Such information is given by CS2CS in a terminal window when the variable PROJ_DEBUG is set to ON.

Can my python script retrieve the same information? Any ideas?


My script is:

from osgeo import osr

#Input coordinates
xin = 436898.0
yin = 5577554.0

# Transformation setup from NAD83 UTM Z15 to NAD27 UTM Z15
srs = osr.SpatialReference()
trs = osr.SpatialReference()
rsTrans = osr.CoordinateTransformation(srs,trs)

# Doing transformation
(xout, yout, zout) = rsTrans.TransformPoint(xin,yin)
share|improve this question
add comment

Know someone who can answer? Share a link to this question via email, Google+, Twitter, or Facebook.

Your Answer


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

Browse other questions tagged or ask your own question.