<div dir="ltr">Hi guys,<div><br></div><div>I'm having weird behavior using ogr and osr through python while reprojecting point geometries. </div><div><br></div><div>Following code:</div><div>wgs = osr.SpatialReference()<br>wgs.ImportFromEPSG(4326)<br><br>etrs = osr.SpatialReference()<br>etrs.ImportFromEPSG(3763)<br><br>geom = ogr.Geometry(ogr.wkbPoint)<br>geom.AddPoint(-9.10165,38.73981)<br>#geom.AddPoint(38.73981,-9.10165)</div><div>geom.AssignSpatialReference(wgs)<br><br>reproj = osr.CoordinateTransformation(wgs, etrs)<br>geom.Transform(reproj)</div><div><br></div><div>If I run it as it is I get wrong coordinates in the target SRS. However if I switch X,Y when building the geometry (the commented line) I will get the correct coordinates. What is weird is that I have another script reading from geonames dataset and there this behavior doesn't happen.</div><div><br></div><div>I'm using GDAL 3.4.1 installed through conda on Ubuntu 20.04.</div><div><br></div><div>Any hints?</div><div>thanks in advance<br clear="all"><div><br></div>-- <br><div dir="ltr" class="gmail_signature" data-smartmail="gmail_signature">Hugo Martins<br></div></div></div>