Re: [gdal-dev] GDAL Python Projection Question

2010-04-12 Thread Kyle Shannon
I use AutoCreateWarpedVRT() to create a virtual dataset, then you CreateCopy() to write your data using the proper Driver: vrt_ds = AutoCreateWarpedVRT( ds, ... ) dst_ds = driver.CreateCopy('dst.out', vrt_ds, ...) dst_ds = None kss # Kyle Shannon Physical Science Te

Re: [gdal-dev] ogr2ogr BeTA2007 grid

2010-04-12 Thread Iván Sánchez Ortega
El día Monday 12 April 2010 13:47:44, Tilo Villwock dijo: > ogr2ogr -s_srs "EPSG:3396 +nadgrids=BETA2007.gsb +wktext" t_srs EPSG:25832 > target_shapefile.shp source_shapefile.shp H. Try the following: forget about epsg:3396, and instead specify all the individual options from epsg:3396.

[gdal-dev] GDAL Python Projection Question

2010-04-12 Thread Discourse Maps
Hi, I am attempting to use GDAL with Python to reproject a raster (called 'deforest') from Geographic Coordinate System to NAD83_UTM10N. Most of the projection tutorials seem to be in C... Can anyone offer assistance for Python? My current code: os.chdir(gp.workspace) ds = gdal.Ope

[gdal-dev] ogr2ogr BeTA2007 grid

2010-04-12 Thread Tilo Villwock
Hi, I am trying to reproject some shapefiles from GKPD83 to ETRS89/UTM with ogr2ogr. I am using the nadgrids option to specify the custom transformation grid (BeTA2007.gsd), but for some reason it is being ignored during the process. Now I know from previous posts on this list that I am suppose