I had to do the following modifications to get it work with libgdal-1.7.2 --- map.py.ori 2011-04-08 15:59:56.000000000 +0300 +++ map.py 2011-04-08 16:01:01.000000000 +0300 @@ -125,7 +125,7 @@ in_ds.SetGCPs( gcp_list, in_ds.GetProjection() ) out_srs = osr.SpatialReference() - out_srs.ImportFromEPSG(900913) + out_srs.ImportFromEPSG(3857) out_ds = gdal.AutoCreateWarpedVRT( in_ds, in_srs.ExportToWkt(), out_srs.ExportToWkt() ) gt = out_ds.GetGeoTransform()