<div dir="ltr">You guys are the best - thank you. I thought I was going crazy. Brad - your code worked great on version 3.6.0 but when I was running this locally on 3.0.4. I get the wrong answer at the end. No problem, though, because this is fixed now. Thank you!!<div><br></div><div>>>> transform_func = osr.CoordinateTransformation(srs2, srs)<br>>>> u = ogr.Geometry(ogr.wkbPoint)<br>>>> u.AddPoint(0,0)<br>>>> u.Transform(transform_func)<br>0<br>>>> print(u.ExportToWkt())<br>POINT (0 0 0)<br></div><div>>>> import osgeo<br>>>> osgeo.__version__<br>'3.0.4'<br></div><div><br></div><div><br></div><div>Meanwhile, works great in 3.6.0</div><div><br></div><div>>>> import osgeo.osr as osr<br>>>> import osgeo.ogr as ogr<br>>>> p = 'PROJCS["WGS_1984_Web_Mercator_Auxiliary_Sphere",GEOGCS["GCS_WGS_1984",DATUM["D_WGS_1984",SPHEROID["WGS_1984",6378137.0,298.257223563]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Mercator_Auxiliary_Sphere"],PARAMETER["False_Easting",0.0],PARAMETER["False_Northing",0.0],PARAMETER["Central_Meridian",-41.34804059746178],PARAMETER["Standard_Parallel_1",0.0],PARAMETER["Auxiliary_Sphere_Type",0.0],UNIT["Meter",1.0]]'<br>>>> p2 = """PROJCS["WGS_1984_Web_Mercator_Auxiliary_Sphere",GEOGCS["GCS_WGS_1984",DATUM["D_WGS_1984",SPHEROID["WGS_1984",6378137.0,298.257223563]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Mercator_Auxiliary_Sphere"],PARAMETER["False_Easting",0.0],PARAMETER["False_Northing",0.0],PARAMETER["Central_Meridian",0.0],PARAMETER["Standard_Parallel_1",0.0],PARAMETER["Auxiliary_Sphere_Type",0.0],UNIT["Meter",1.0]]"""<br>>>> srs2 = osr.SpatialReference(wkt=p2)<br>>>> srs = osr.SpatialReference(wkt=p)<br></div><div>>>> transform_func = osr.CoordinateTransformation(srs2, srs)<br>>>> u = ogr.Geometry(ogr.wkbPoint)<br>>>> u.AddPoint(0,0)<br>>>> u.Transform(transform_func)<br>0<br>>>> print(u.ExportToWkt())<br>POINT (4602842.82460905 0.0 0)<br></div><div>>>> osgeo.__version__<br>'3.6.0'<br></div><div><br></div><div><br></div><div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Thu, Mar 16, 2023 at 3:14 AM Brad Hards <<a href="mailto:bradh@frogmouth.net">bradh@frogmouth.net</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">On Thursday, 16 March 2023 3:13:59 PM AEDT Tom Hayden wrote:<br>
> >>> u.Transfrom(transform_func)<br>
Did this line not give you an error?<br>
<br>
In any case, I used:<br>
<br>
<br>
import osgeo.osr as osr<br>
import osgeo.ogr as ogr<br>
p = 'PROJCS["WGS_1984_Web_Mercator_Auxiliary_Sphere",GEOGCS["GCS_WGS_1984",DATUM["D_WGS_1984",SPHEROID["WGS_1984",6378137.0,298.257223563]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Mercator_Auxiliary_Sphere"],PARAMETER["False_Easting",0.0],PARAMETER["False_Northing",0.0],PARAMETER["Central_Meridian",-41.34804059746178],PARAMETER["Standard_Parallel_1",0.0],PARAMETER["Auxiliary_Sphere_Type",0.0],UNIT["Meter",1.0]]'<br>
p2 = """PROJCS["WGS_1984_Web_Mercator_Auxiliary_Sphere",GEOGCS["GCS_WGS_1984",DATUM["D_WGS_1984",SPHEROID["WGS_1984",6378137.0,298.257223563]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Mercator_Auxiliary_Sphere"],PARAMETER["False_Easting",0.0],PARAMETER["False_Northing",0.0],PARAMETER["Central_Meridian",0.0],PARAMETER["Standard_Parallel_1",0.0],PARAMETER["Auxiliary_Sphere_Type",0.0],UNIT["Meter",1.0]]"""<br>
srs2 = osr.SpatialReference(wkt=p2)<br>
srs = osr.SpatialReference(wkt=p)<br>
transform_func = osr.CoordinateTransformation(srs2, srs)<br>
u = ogr.Geometry(ogr.wkbPoint)<br>
u.AddPoint(0,0)<br>
u.Transform(transform_func)<br>
print(u.ExportToWkt())<br>
<br>
which resulted in:<br>
POINT (4602842.82460905 0.0 0)<br>
<br>
<br>
Brad<br>
<br>
<br>
_______________________________________________<br>
gdal-dev mailing list<br>
<a href="mailto:gdal-dev@lists.osgeo.org" target="_blank">gdal-dev@lists.osgeo.org</a><br>
<a href="https://lists.osgeo.org/mailman/listinfo/gdal-dev" rel="noreferrer" target="_blank">https://lists.osgeo.org/mailman/listinfo/gdal-dev</a><br>
</blockquote></div>