[gdal-dev] [Gdal-dev] python - reproject a shapefile

[Rutger]

Hey,

As you already mention in your post, you didn't specify a SRSfor the input,
that's probably the main issue. Inputs need to be georeferenced correctly. 

If you 'only' want to do a reprojection, you could consider using the
bindings to the command line utilities. This would simplify you code a lot,
reducing the possibility of bugs. Something like:

/dsvec = gdal.VectorTranslate('output.shp', 'input.shp' srcSRS=input_srs,
dstSRS="EPSG:32633")
dsvec = None /

Usually you don't have to specify the input SRS, but since yours is missing,
you can assign (overrule!) it like this.



--
Sent from: http://osgeo-org.1560.x6.nabble.com/GDAL-Dev-f3742093.html