<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=UTF-8">
</head>
<body>
<p>Daniel,</p>
<p>It is likely one of the ESRI specific things at
<a class="moz-txt-link-freetext" href="https://raw.githubusercontent.com/Esri/projection-engine-db-doc/master/csv/pe_list_geogtran.csv">https://raw.githubusercontent.com/Esri/projection-engine-db-doc/master/csv/pe_list_geogtran.csv</a><br>
</p>
<p>The situation is that the import script of the ESRI database does
import them in the PROJ database with some ad-hoc logic, but
proj_create() has no support for those strings themselves.
However, as they are in the database, you can use them by name
(provided yours is a registered one, and not a custom one), with
proj_create_from_name() (or if you find its code in the above
list, with proj_create_from_database())</p>
<p>e.g.<br>
</p>
<p>projinfo -k operation Palestine_1923_To_WGS_1984_1X<br>
PROJ string:<br>
+proj=pipeline<br>
+step +proj=axisswap +order=2,1<br>
+step +proj=unitconvert +xy_in=deg +xy_out=rad<br>
+step +inv +proj=longlat +a=6378300.789 +b=6356566.435<br>
+step +proj=push +v_3<br>
+step +proj=cart +a=6378300.789 +b=6356566.435<br>
+step +proj=helmert +x=-181 +y=-122 +z=225<br>
+step +inv +proj=cart +ellps=WGS84<br>
+step +proj=pop +v_3<br>
+step +proj=unitconvert +xy_in=rad +xy_out=deg<br>
+step +proj=axisswap +order=2,1<br>
<br>
WKT2:2019 string:<br>
COORDINATEOPERATION["Palestine_1923_To_WGS_1984_1X",<br>
SOURCECRS[<br>
GEOGCRS["Palestine 1923",<br>
DATUM["Palestine 1923",<br>
ELLIPSOID["Clarke 1880
(Benoit)",6378300.789,293.466315538981,<br>
LENGTHUNIT["metre",1]]],<br>
PRIMEM["Greenwich",0,<br>
ANGLEUNIT["degree",0.0174532925199433]],<br>
CS[ellipsoidal,2],<br>
AXIS["geodetic latitude (Lat)",north,<br>
ORDER[1],<br>
ANGLEUNIT["degree",0.0174532925199433]],<br>
AXIS["geodetic longitude (Lon)",east,<br>
ORDER[2],<br>
ANGLEUNIT["degree",0.0174532925199433]],<br>
ID["EPSG",4281]]],<br>
TARGETCRS[<br>
GEOGCRS["WGS 84",<br>
ENSEMBLE["World Geodetic System 1984 ensemble",<br>
MEMBER["World Geodetic System 1984 (Transit)"],<br>
MEMBER["World Geodetic System 1984 (G730)"],<br>
MEMBER["World Geodetic System 1984 (G873)"],<br>
MEMBER["World Geodetic System 1984 (G1150)"],<br>
MEMBER["World Geodetic System 1984 (G1674)"],<br>
MEMBER["World Geodetic System 1984 (G1762)"],<br>
ELLIPSOID["WGS 84",6378137,298.257223563,<br>
LENGTHUNIT["metre",1]],<br>
ENSEMBLEACCURACY[2.0]],<br>
PRIMEM["Greenwich",0,<br>
ANGLEUNIT["degree",0.0174532925199433]],<br>
CS[ellipsoidal,2],<br>
AXIS["geodetic latitude (Lat)",north,<br>
ORDER[1],<br>
ANGLEUNIT["degree",0.0174532925199433]],<br>
AXIS["geodetic longitude (Lon)",east,<br>
ORDER[2],<br>
ANGLEUNIT["degree",0.0174532925199433]],<br>
ID["EPSG",4326]]],<br>
METHOD["Geocentric translations (geog2D domain)",<br>
ID["EPSG",9603]],<br>
PARAMETER["X-axis translation",-181,<br>
LENGTHUNIT["metre",1],<br>
ID["EPSG",8605]],<br>
PARAMETER["Y-axis translation",-122,<br>
LENGTHUNIT["metre",1],<br>
ID["EPSG",8606]],<br>
PARAMETER["Z-axis translation",225,<br>
LENGTHUNIT["metre",1],<br>
ID["EPSG",8607]],<br>
OPERATIONACCURACY[1.5],<br>
USAGE[<br>
SCOPE["Not known."],<br>
AREA["Israel - onshore and offshore."],<br>
BBOX[29.45,32.99,33.53,35.69]],<br>
ID["ESRI",108220]]<br>
</p>
<p><br>
</p>
<p>If it is a custom one, the logic at
<a class="moz-txt-link-freetext" href="https://github.com/OSGeo/PROJ/blob/master/scripts/build_db_from_esri.py#L1480">https://github.com/OSGeo/PROJ/blob/master/scripts/build_db_from_esri.py#L1480</a>
might help <br>
</p>
<p>Even<br>
</p>
<p><br>
</p>
<div class="moz-cite-prefix">Le 29/05/2021 à 00:00, Daniel Schilling
a écrit :<br>
</div>
<blockquote type="cite"
cite="mid:CALCct-xkwk0e=wsQUWBeNi+=mit-KXaLkBbqrcG1-rPVgnVfRQ@mail.gmail.com">
<meta http-equiv="content-type" content="text/html; charset=UTF-8">
<div dir="ltr">
<div>Greetings! </div>
<div>I am trying to create a bound CRS using
`proj_crs_create_bound_crs`, but I only have a WKT string
(`GEOGTRAN`...) for the transformation, not an EPSG code. I
cannot figure out how to do it.</div>
<div><br>
</div>
<div>Neither `proj_create` nor `proj_create_from_wkt` seem able
to create transformations. The only way I have been able to
create a transformation so far is by using
`proj_create_from_database`, but again - I don't have an EPSG
code for this. Is there a way to do this? I am using Proj 8.</div>
<div><br>
</div>
<div>Thank you in advance for your help!<br>
</div>
<div>Daniel Schilling<br>
</div>
</div>
<br>
<fieldset class="mimeAttachmentHeader"></fieldset>
<pre class="moz-quote-pre" wrap="">_______________________________________________
PROJ mailing list
<a class="moz-txt-link-abbreviated" href="mailto:PROJ@lists.osgeo.org">PROJ@lists.osgeo.org</a>
<a class="moz-txt-link-freetext" href="https://lists.osgeo.org/mailman/listinfo/proj">https://lists.osgeo.org/mailman/listinfo/proj</a>
</pre>
</blockquote>
<pre class="moz-signature" cols="72">--
<a class="moz-txt-link-freetext" href="http://www.spatialys.com">http://www.spatialys.com</a>
My software is free, but my time generally not.</pre>
</body>
</html>