Hi again, a brief update on my issue. I discovered that the operation works if I explicitly transform the coordinates in two steps. The first step I transform the coordinates from the projected CRS to the Geodetic CRS, then I can apply the selected Operation. Is this the correct approach?
Thank you Stefano Il giorno ven 5 mar 2021 alle ore 09:42 Stefano Iacovella < [email protected]> ha scritto: > Hi all, > > I am trying to use GDAL, via Python bindings, to convert coordinates. > I would like to select a specific transformation when converting among CRS > with different Datums. Reading the documentation my understanding is that I > need to create a osr.CoordinateTransformation with options parameter > valorized properly. > I am using binaries from OSGEO4WIN, the release are 3.1.4 for GDAL and > 6.3.2 for PROJ > This is a code snippet I am using to debug: > > from osgeo import ogr > from osgeo import osr > > source = osr.SpatialReference() > source.ImportFromEPSG(3003) > > target = osr.SpatialReference() > target.ImportFromEPSG(4326) > > options = osr.CoordinateTransformationOptions() > options.SetOperation("urn:ogc:def:coordinateOperation:EPSG::1169") > transform = osr.CoordinateTransformation(source, target, options) > > point = ogr.CreateGeometryFromWkt("POINT (1622987 4756952)") > point.Transform(transform) > > print point.ExportToWkt() > > I have no errors running it but the coordinates of the point, after the > transformation, are infinite. If I create the CoordinateTransformation > without options like this: > > transform = osr.CoordinateTransformation(source, target) > > The point coordinates after transformation are ok but of course the > conversion is performed using the Coordinate Operation automatically > selected by the system. > Am I doing something wrong? > > Thank you in advance for your help. > > All the best, > > Stefano > >
_______________________________________________ gdal-dev mailing list [email protected] https://lists.osgeo.org/mailman/listinfo/gdal-dev
