<div dir="ltr"><br><br>On Fri, Mar 8, 2019 at 12:41 PM Even Rouault <<a href="mailto:even.rouault@spatialys.com">even.rouault@spatialys.com</a>> wrote:<br>><br>> On vendredi 8 mars 2019 12:27:37 CET Markus Metz wrote:<br>> > On Fri, Mar 8, 2019 at 10:49 AM Roger Bivand <<a href="mailto:Roger.Bivand@nhh.no">Roger.Bivand@nhh.no</a>> wrote:<br>> > > Since rgdal::make_EPSG() is facing the same problems of listing tabulated<br>> > > EPSG fields as g.proj -l, I was very happy to see Markus' code in<br>> > > g.proj/main.c mentioned in this thread, and have used this approach in<br>> ><br>> > <a href="https://r-forge.r-project.org/scm/viewvc.php/pkg/src/proj_info6.cpp?view=mar">https://r-forge.r-project.org/scm/viewvc.php/pkg/src/proj_info6.cpp?view=mar</a><br>> > kup&root=rgdal<br>> > > However, there are plenty of messages such as: "proj_as_proj_string:<br>> > > Unsupported conversion method: Lambert Conic Conformal (West<br>> ><br>> > Orientated)". I<br>> ><br>> > > haven't installed GRASS trunk with PROJ6, so I can't see whether g.proj -l<br>> > > also sees the same messages.<br>> ><br>> > I get the same messages with g.proj, but no error code from PROJ error<br>> > handlers.<br>> ><br>> > > If it does, maybe we could ask on the proj list<br>> ><br>> > yes, we should ask. For now g.proj is ignoring those codes for which no<br>> > proj string can be retrieved.<br>><br><div>> I'm not sure what the question is exactly :-) <br></div><div><br></div><div>My question is, what should GRASS do with a SRS without proj string? The best answer is probably: use WKT instead of a proj string. But most of the GRASS code uses (for historical reasons) some form of proj syntax, also to store projection info for a GRASS location. If there is no proj string, GRASS thinks there is no SRS and falls back to a generic xy system which is not always correct.</div><div><br></div><div>Markus M<br></div></div>