Geographic Information Systems Stack Exchange is a question and answer site for cartographers, geographers and GIS professionals. It's 100% free, no registration required.

Sign up
Here's how it works:
  1. Anybody can ask a question
  2. Anybody can answer
  3. The best answers are voted up and rise to the top

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

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.