<div dir="ltr"><div dir="ltr">> I'm a bit extreme when I say "don't use PROJ strings", but otherwise people <br>> don't get it :-) So yes if you don't care about datum transformations, don't <br>> use unusual axis orders, you might still use them. They won't disappear in a <br>> foreseable future.<br>> The GDAL OGRSpatialReference class has a number of helper methods to help you <br>> build your own custom CRS, like SetLAEA(), SetLCC(), etc. and SetGeogCS()<br></div><div><br></div>Oh excellent, that's less scary than I'd understood. <div><br><div> I'll look at exposing those custom builders in R. <div><br></div><div>Cheers, Mike. <br></div><div><br></div><div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Wed, Oct 30, 2019 at 9:33 AM Even Rouault <<a href="mailto:even.rouault@spatialys.com">even.rouault@spatialys.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">On mercredi 30 octobre 2019 09:22:11 CET Michael Sumner wrote:<br>
> Great discussion, thanks to both of you.  I have a question.<br>
> <br>
> I wonder how we ought to specify a projection that we desire, I work in<br>
> contexts where 1.5m or even 1.5 km is not an unreasonable level of<br>
> precision, and we get a lot of value out of specifying local projections<br>
> (at ocean basin scale) in LAEA, LCC, STERE and others  (OMERC is especially<br>
> good for hermisphere-scale migrating birds, very custom and not related to<br>
> human activities) - PROJ strings are obviously an easy to specify the<br>
> family and centre point of the projection, and we have a lot of workflows<br>
> (computer and human) for this. The datum is a bit irrelevant a lot of the<br>
> time but we always include it and try to explain it to users.<br>
> <br>
> Are we expected to express a full WKT or WKT2 string for this?   I know we<br>
> can write wrappers in some language, but there doesn't seem to be a culture<br>
> around crafting custom projections that I can latch onto. Usually I see<br>
> people looking for an EPSG code for their region, and as often as not they<br>
> end up using UTM because "it's a projection" and no authority has ever done<br>
> a survey in their part of the Southern Ocean and published it with EPSG.<br>
> <br>
> Is something like the gdaltransform pipelining syntax the most concise we<br>
> can aim for, or do we somehow cache our own local EPSG-like codes?<br>
<br>
I'm a bit extreme when I say "don't use PROJ strings", but otherwise people <br>
don't get it :-) So yes if you don't care about datum transformations, don't <br>
use unusual axis orders, you might still use them. They won't disappear in a <br>
foreseable future.<br>
The GDAL OGRSpatialReference class has a number of helper methods to help you <br>
build your own custom CRS, like SetLAEA(), SetLCC(), etc. and SetGeogCS()<br>
<br>
The PROJJSON notation might also be slightly easier to work with than WKT2 (it <br>
is completely equivalent to it, just using a more standard encoding):<br>
example at <a href="https://github.com/OSGeo/PROJ/pull/1547" rel="noreferrer" target="_blank">https://github.com/OSGeo/PROJ/pull/1547</a><br>
Schema at <a href="https://proj.org/schemas/v0.1/projjson.schema.json" rel="noreferrer" target="_blank">https://proj.org/schemas/v0.1/projjson.schema.json</a><br>
<br>
Even<br>
<br>
-- <br>
Spatialys - Geospatial professional services<br>
<a href="http://www.spatialys.com" rel="noreferrer" target="_blank">http://www.spatialys.com</a><br>
</blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr" class="gmail_signature">Michael Sumner<br>Software and Database Engineer<br>Australian Antarctic Division<br>Hobart, Australia<br>e-mail: <a href="mailto:mdsumner@gmail.com" target="_blank">mdsumner@gmail.com</a></div></div></div></div></div>