<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=UTF-8">
</head>
<body>
<p class="MsoNormal" style="caret-color: rgb(0, 0,
0);font-variant-caps: normal;orphans:
auto;text-align:start;widows: auto;-webkit-text-size-adjust:
auto;-webkit-text-stroke-width: 0px;word-spacing:0px">Jennifer,<span
style="color:black"> <o:p></o:p></span></p>
<blockquote type="cite"
cite="mid:57B86F67-FD31-45B3-8D99-D0FBE9381029@nasa.gov">
<div class="WordSection1">
<p class="MsoNormal" style="caret-color: rgb(0, 0,
0);font-variant-caps: normal;orphans:
auto;text-align:start;widows: auto;-webkit-text-size-adjust:
auto;-webkit-text-stroke-width: 0px;word-spacing:0px">
<span style="color:black">Note that the projection parameters
are all the same (but with slightly different naming
conventions), but the CRS AUTHORITY line does not exist in
the GDAL version, is EPSG 102022 in the spatialreference.org
version, and is ESRI 102022 in the epsg.io version.</span></p>
</div>
</blockquote>
<p>you must use an outdated version of PROJ as with recent ones
SetFromUserInput("ESRI:102022") returns<br>
<br>
PROJCRS["Africa_Albers_Equal_Area_Conic",<br>
BASEGEOGCRS["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>
MEMBER["World Geodetic System 1984 (G2139)"],<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>
ID["EPSG",4326]],<br>
CONVERSION["Africa_Albers_Equal_Area_Conic",<br>
METHOD["Albers Equal Area",<br>
ID["EPSG",9822]],<br>
PARAMETER["Latitude of false origin",0,<br>
ANGLEUNIT["degree",0.0174532925199433],<br>
ID["EPSG",8821]],<br>
PARAMETER["Longitude of false origin",25,<br>
ANGLEUNIT["degree",0.0174532925199433],<br>
ID["EPSG",8822]],<br>
PARAMETER["Latitude of 1st standard parallel",20,<br>
ANGLEUNIT["degree",0.0174532925199433],<br>
ID["EPSG",8823]],<br>
PARAMETER["Latitude of 2nd standard parallel",-23,<br>
ANGLEUNIT["degree",0.0174532925199433],<br>
ID["EPSG",8824]],<br>
PARAMETER["Easting at false origin",0,<br>
LENGTHUNIT["metre",1],<br>
ID["EPSG",8826]],<br>
PARAMETER["Northing at false origin",0,<br>
LENGTHUNIT["metre",1],<br>
ID["EPSG",8827]]],<br>
CS[Cartesian,2],<br>
AXIS["(E)",east,<br>
ORDER[1],<br>
LENGTHUNIT["metre",1]],<br>
AXIS["(N)",north,<br>
ORDER[2],<br>
LENGTHUNIT["metre",1]],<br>
USAGE[<br>
SCOPE["Not known."],<br>
AREA["Africa"],<br>
BBOX[-35,-25,39,55]],<br>
ID["ESRI",102022]]<br>
<br>
</p>
<blockquote type="cite"
cite="mid:57B86F67-FD31-45B3-8D99-D0FBE9381029@nasa.gov">
<div class="WordSection1">
<p class="MsoNormal" style="caret-color: rgb(0, 0,
0);font-variant-caps: normal;orphans:
auto;text-align:start;widows: auto;-webkit-text-size-adjust:
auto;-webkit-text-stroke-width: 0px;word-spacing:0px"><span
style="color:black"><o:p></o:p></span></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">Is it best to represent the CRS with WKT
rather than importing using SetFromUserInput with EPSG:code or
ESRI:code? </p>
</div>
</blockquote>
If you use the WKT returned by GDAL/PROJ, then using it through
SetFromUserInput() will be equivalent as with auth:code<br>
<br>
<blockquote type="cite"
cite="mid:57B86F67-FD31-45B3-8D99-D0FBE9381029@nasa.gov">
<div class="WordSection1">
<p class="MsoNormal">If so, what source has the definitive WKT,
spatialreference.org or epsg.io?</p>
</div>
</blockquote>
<p>Well, none of them as they are generally outdated compared if you
use the latest of GDAL (actually PROJ). But there are a lot of
degree of freedoms in WKT (especially on the amount of metadata
you include), so there's no such thing as a canonical WKT. That
said I'd consider the one returned by PROJ to be reasonably a
reference one.<br>
</p>
<p><br>
</p>
<p>Just a slight correction to Howard's answer regarding:<br>
</p>
<pre class="moz-quote-pre" wrap="">
<a class="moz-txt-link-freetext" href="https://epsg.org/">https://epsg.org/</a> runs GDAL/PROJ too, but I don't know which version, and of course that isn't going to respond to ESRI codes.</pre>
<p></p>
<p><br>
</p>
<p>Actually, epsg.org does NOT run PROJ. This is their own software
stack, and I've actually reported to them bugs in their WKT output
for some CRS in the past (to my knowledge, it is all now fixed
now)</p>
<p>Even<br>
</p>
<blockquote type="cite"
cite="mid:57B86F67-FD31-45B3-8D99-D0FBE9381029@nasa.gov">
<div class="WordSection1">
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<br>
<fieldset class="moz-mime-attachment-header"></fieldset>
<pre class="moz-quote-pre" wrap="">_______________________________________________
gdal-dev mailing list
<a class="moz-txt-link-abbreviated" href="mailto:gdal-dev@lists.osgeo.org">gdal-dev@lists.osgeo.org</a>
<a class="moz-txt-link-freetext" href="https://lists.osgeo.org/mailman/listinfo/gdal-dev">https://lists.osgeo.org/mailman/listinfo/gdal-dev</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>