<div dir="ltr"><div>Both examples look like what I was looking for</div><div><br></div><div>Thanks a lot Even :)</div><div><br></div><div>Stefano<br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">Il giorno mar 19 apr 2022 alle ore 12:17 Even Rouault <<a href="mailto:even.rouault@spatialys.com">even.rouault@spatialys.com</a>> ha scritto:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Hi,<br>
<br>
looking at <br>
<a href="https://github.com/OSGeo/PROJ/blob/master/data/sql/other_transformation_custom.sql" rel="noreferrer" target="_blank">https://github.com/OSGeo/PROJ/blob/master/data/sql/other_transformation_custom.sql</a> <br>
should give you good hints on how to add custom transformations, here <br>
under a PROJ authority<br>
<br>
If you want to create your own authority, you'll also need to do similar <br>
steps as the ones dones for the NKG authority in <br>
<a href="https://github.com/OSGeo/PROJ/blob/master/data/sql/nkg_post_customizations.sql" rel="noreferrer" target="_blank">https://github.com/OSGeo/PROJ/blob/master/data/sql/nkg_post_customizations.sql</a><br>
<br>
Even<br>
<br>
Le 19/04/2022 à 08:55, Stefano Iacovella a écrit :<br>
> Dear all,<br>
><br>
> I am wondering which is the best approach to use a custom operation to <br>
> transform coordinates from a datum to another.<br>
> I have managed to use a custom operation by creating a specific <br>
> pipeline and set it as:<br>
><br>
> options = osr.CoordinateTransformationOptions()<br>
><br>
> options.SetOperation('''+proj=pipeline<br>
>                         +step +proj=axisswap +order=2,1<br>
>                         +step +proj=unitconvert +xy_in=deg +xy_out=rad <br>
> +step +proj=push +v_3<br>
>                         +step +proj=cart +ellps=intl<br>
>                         +step +proj=helmert +x=-104.1 +y=-49.1 +z=-9.9 <br>
> +rx=0.971 +ry=-2.917 +rz=0.714 +s=-11.68  +convention=position_vector<br>
>                         +step +inv +proj=cart +ellps=WGS84<br>
>                         +step +proj=pop +v_3<br>
>                         +step +proj=unitconvert<br>
>                         +xy_in=rad<br>
>                         +xy_out=deg<br>
>                         +step +proj=axisswap +order=2,1''')<br>
><br>
> transform = osr.CoordinateTransformation(source_crs, target_crs, options)<br>
><br>
> I have a few custom operation, then I would like to understand if it <br>
> is possible to include their definitions in the PROJ configuration, <br>
> with a different autorithy then EPSG. So that I can use it in <br>
> something similar to:<br>
><br>
> options.SetOperation('urn:ogc:def:coordinateOperation:my_authority::900001')<br>
><br>
> I tried to read documentation of GDAl and PROJ but didn't find, <br>
> probably my fault, any clear explanation on how to customize <br>
> configuration to include custom operation.<br>
><br>
> Thank you all in advance for any hint<br>
><br>
> Stefano<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>
<br>
-- <br>
<a href="http://www.spatialys.com" rel="noreferrer" target="_blank">http://www.spatialys.com</a><br>
My software is free, but my time generally not.<br>
<br>
</blockquote></div>