Sorry to keep bugging you. Here is my code and the output it produces: import osr dst=osr.SpatialReference() print 'r1=',dst.ImportFromEPSG( 32119 ) print 'r2=',dst.ExportToWkt() print 'done' $ python trypy.py r1= 7 r2= done Do I need some EPSG database? Does it come with gdal?Thanks for your help gb