<div dir="ltr"><div>Several different methods might be available to reproject from one CRS to another CRS. PROJ6 can select the most appropriate method if a bounding box is provided. This means that the selected method depends on the bounding box and that results of a reprojection can differ depending on the provided bounding box.</div><div><br></div><div>In GRASS, this bounding box can be obtained from the current region. That means if the current region changes, the method and thus the results of coordinate reprojection might change. This effect can not be underestimated.</div><div><br></div><div>With the current pull request #118 the current region is used to help PROJ6 select the best method, and [r|v].import should work as before. Results can be different compared to PROJ5 or earlier.</div><div><br></div><div>IMHO, there is no way around that users become more familiar with details of coordinate reprojection and need to read the output of [r|v].proj carefully regarding different methods known to PROJ6.<br></div><div><br></div><div>PROJ6 will not use the best method if any required datum transformation grid is not available. Users will need to obtain the corresponding grid themselves, which raises another question: where to put this grid? On Linux, this would be e.g. /usr/share/proj, but if the user can not write to /usr/share/proj, the grids need to be saved somewhere in $HOME, and PROJ6 must look at that place. This problem, i.e. the place where PROJ6 should look for in $HOME is not yet solved AFAICT (proj-6.2.0). <br></div><div><br></div><div>Markus M<br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Wed, Sep 4, 2019 at 9:46 PM Markus Metz <<a href="mailto:markus.metz.giswork@gmail.com">markus.metz.giswork@gmail.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"><div dir="ltr"><br><br>On Wed, Sep 4, 2019 at 6:01 PM Markus Neteler <<a href="mailto:neteler@osgeo.org" target="_blank">neteler@osgeo.org</a>> wrote:<br>><br>> On Wed, Sep 4, 2019 at 4:31 AM Anna Petrášová <<a href="mailto:kratochanna@gmail.com" target="_blank">kratochanna@gmail.com</a>> wrote:<br>> > On Tue, Sep 3, 2019 at 4:22 AM Markus Metz <<a href="mailto:markus.metz.giswork@gmail.com" target="_blank">markus.metz.giswork@gmail.com</a>> wrote:<br>> >><br>> >> Hi all,<br>> >><br>> >> there is a new PR for PROJ6 support in GRASS<br>> >> <a href="https://github.com/OSGeo/grass/pull/118" target="_blank">https://github.com/OSGeo/grass/pull/118</a><br>> >><br>> >> There are two important changes when using PROJ6:<br>> >><br>> >> First, reprojection with v.proj and r.proj is no longer always possible without the user making informed decisions. The reason is that there can be several different operations available to reproject coordinates from one CRS to another CRS. These different operations are listed and the user has to provide the appropriate operation with the pipeline option, taking care of any axisswap.<br>> ><br>> > first, thanks for all the work. Second, I don't see how most users are supposed to know what to pick. Is there perhaps a way to pick a good default? I just can't imagine not having r.import/v.import...<br>><br>> I see the same problem: users won't know what to select if defaults<br>> are gone with PROJ 6.<br><div><br></div><div>We can provide information that enables users to make an informed decision. The options listed by PROJ6 are ordered by guessed best choice, i.e. the first one is, given the provided information, the best choice. But users need to review the list of options.</div><div><br></div><div>Reprojection from one CRS to another CRS was never easy. For both raster and vector data, some preparation was always needed to decide on appropriate settings. PROJ6 provides methods to improve the accuracy of reprojected coordinates, depending on the actual use case. The user decides (must decide).<br></div><div>></div>> > Anna<br>> >><br>> >><br>> >> Second, axis order is no longer always easting, northing, e.g. for EPSG:4326 it is northing, easting, and an axisswap might need to be removed from operations provided by PROJ.<br>> >><br>> >> There are many more changes (see details in the PR), but these are the two most important ones.<br>> >><br>> >> Feedback welcome!<br>><br>> Thanks for your hard work.<br>><br>> I talked to Robert Bivand here at the GeoSTAT 2019 in Münster who gave<br>> a talk about the topic:<br>><br>> "Not just R-spatial: sustaining open source geospatial software stacks"<br>> <a href="https://github.com/rsbivand/geostat19_talk" target="_blank">https://github.com/rsbivand/geostat19_talk</a><br>><br>> You can quick-view the talk in rendered HTML like this (maybe there<br>> are better ways):<br>> <a href="http://htmlpreview.github.io/?https://raw.githubusercontent.com/rsbivand/geostat19_talk/master/docs/geostat_bivand_19.html" target="_blank">http://htmlpreview.github.io/?https://raw.githubusercontent.com/rsbivand/geostat19_talk/master/docs/geostat_bivand_19.html</a><br>><br>> It contains a section about PROJ (6).<br><div><br></div><div>Not so random citation:</div><div>"...we need to manipulate the CRS read in with the file to insert our choice of how to make the transformation..."</div><div><br></div><div>I will try to restrict the number of options based on the current region in a modified PR.<br></div><div><br></div><div>Markus M<br></div><br></div>
</blockquote></div>