Coordinate transformations with pipelines

--000000000000063443064b8cb4b9
Content-Type: text/plain; charset="UTF-8"

I am usign PostGIS 3.5 with PostgreSQL 17.0

Here are the details for my installation:

"postgis_full_version"
"POSTGIS=""3.5.3 3.5.3"" [EXTENSION] PGSQL=""170""
GEOS=""3.13.1-CAPI-1.19.2"" PROJ=""8.2.1 NETWORK_ENABLED=OFF
URL_ENDPOINT=https://cdn.proj.org
USER_WRITABLE_DIRECTORY=C:\Windows\ServiceProfiles\NetworkService\AppData\Local/proj

DATABASE_PATH=C:\Program
Files\PostgreSQL\17\share\contrib\postgis-3.5\proj\proj.db"" (compiled
against PROJ 8.2.1)
LIBXML=""2.12.5"" LIBJSON=""0.12"" LIBPROTOBUF=""1.2.1"" WAGYU=""0.5.0
(Internal)"" (core procs from ""3.5.2 3.5.2"" need upgrade)"


I have done the following coordinate transformation
from  ED50 / UTM zone 33N  (espg:23033)
to  ETRS89 / UTM zone 33N  (epsg:25833

SELECT
  ST_X(ST_Transform(ST_Point(541000,4516000,23033),25833)) as X,
  ST_Y(ST_Transform(ST_Point(541000,4516000,23033),25833)) as Y;

-- result
"x" "y"
540998.1383329581 4515921.565374701

I have read that the best way to utilize PROJ.6 capabilities is by allowing
the library to freely choose the best possible transformation paths
between the origin and the destination CRSes.

Does someone know how to find which specific transformation pipeline has
been chosen by the PROJ.6 library in this specific transformation ?

Antonio

--000000000000063443064b8cb4b9
Content-Type: text/html; charset="UTF-8"
Content-Transfer-Encoding: quoted-printable

<div dir=3D"ltr"><div>I am usign PostGIS 3.5 with PostgreSQL 17.0</div><div=
><br></div><div>Here are the details for my installation:</div><div><br></d=
iv><div>&quot;postgis_full_version&quot;<br>&quot;POSTGIS=3D&quot;&quot;3.5=
.3 3.5.3&quot;&quot; [EXTENSION] PGSQL=3D&quot;&quot;170&quot;&quot; GEOS=
=3D&quot;&quot;3.13.1-CAPI-1.19.2&quot;&quot; PROJ=3D&quot;&quot;8.2.1 NETW=
ORK_ENABLED=3DOFF <br>URL_ENDPOINT=3D<a href=3D"https://cdn.proj.org">https=
://cdn.proj.org</a> USER_WRITABLE_DIRECTORY=3DC:\Windows\ServiceProfiles\Ne=
tworkService\AppData\Local/proj <br>DATABASE_PATH=3DC:\Program Files\Postgr=
eSQL\17\share\contrib\postgis-3.5\proj\proj.db&quot;&quot; (compiled agains=
t PROJ 8.2.1) <br>LIBXML=3D&quot;&quot;2.12.5&quot;&quot; LIBJSON=3D&quot;&=
quot;0.12&quot;&quot; LIBPROTOBUF=3D&quot;&quot;1.2.1&quot;&quot; WAGYU=3D&=
quot;&quot;0.5.0 (Internal)&quot;&quot; (core procs from &quot;&quot;3.5.2 =
3.5.2&quot;&quot; need upgrade)&quot;</div><div><br></div><div><br></div><d=
iv>I have done the following coordinate transformation<br>from=C2=A0 ED50 /=
 UTM zone 33N

=C2=A0(espg:23033)<br></div><div>to=C2=A0 ETRS89 / UTM zone 33N

=C2=A0(epsg:25833<br></div><div><br>SELECT <br>=C2=A0 ST_X(ST_Transform(ST_=
Point(541000,4516000,23033),25833)) as X,<br>=C2=A0 ST_Y(ST_Transform(ST_Po=
int(541000,4516000,23033),25833)) as Y;<br><br>-- result<br>&quot;x&quot;	&=
quot;y&quot;<br>540998.1383329581	4515921.565374701<br><br>I have read that=
 the best way to utilize PROJ.6 capabilities is by allowing the library to =
freely choose the best possible transformation paths <br>between the origin=
 and the destination CRSes.<br><br>Does someone know how to find which spec=
ific transformation pipeline has been chosen by the PROJ.6 library in this =
specific transformation ?</div><div><br></div><div>Antonio</div><div><br><b=
r></div></div>

--000000000000063443064b8cb4b9--