<div dir="ltr"><div><div><div><div><br><br>On Fri, Nov 10, 2017 at 6:08 PM, Vaclav Petras <<a href="mailto:wenzeslaus@gmail.com">wenzeslaus@gmail.com</a>> wrote:<br>><br>> Hi all,<br>><br>> I'm trying to understand a change of behavior between 7.2 and trunk. With 7.2, I can do the following:<br>><br>> curl -SL <a href="http://fatra.cnr.ncsu.edu/foss4g2017/nc_tile_0793_016_spm.zip">http://fatra.cnr.ncsu.edu/foss4g2017/nc_tile_0793_016_spm.zip</a> \<br>>     > nc_tile_0793_016_spm.zip  && unzip nc_tile_0793_016_spm.zip<br>><br>> grass72 -c EPSG:3358 ~/grassdata/crs_test -e<br>> grass72 ~/grassdata/crs_test2/PERMANENT/ --exec r.in.lidar input=nc_tile_0793_016_spm.las output=count_10 method=n -e -n resolution=10<br>><br>> The last step (r.in.lidar) fails in current trunk because the projection don't match. Here is the difference (name, datum and swapped lat_1 and lat_2):<br>><br>>        GRASS LOCATION PROJ_INFO is:<br>>        name: NAD83(HARN) / North Carolina<br>>        datum: nad83harn<br>>        lat_1: 36.16666666666666<br>>        lat_2: 34.33333333333334<br>><br>>        Import dataset PROJ_INFO is:<br>>        name: NAD_1983_StatePlane_North_Carolina_FIPS_3200<br>>        datum: nad83<br>>        lat_1: 34.33333333333334<br>>        lat_2: 36.16666666666666<br>><br>> 7.2 considers this OK while trunk considers this different. I'm not sure which one is correct.<br><br></div>In this case, 7.2 is correct, but 7.2 does not compare datums, i.e. projections would be regarded as matching even with e.g. nad27 and nad83.<br></div>The comparison of datums is new in 7.3 and needs some refinement.<br><br></div>I will fix the comparison of swapped lat_1 and lat_2.<br><br></div>Markus M<br><div><div><div><br><div><br>><br>> Here is the lasinfo reference:<br>><br>> PROJCS["NAD_1983_StatePlane_North_Carolina_FIPS_3200",<br>>     GEOGCS["GCS_North_American_1983",<br>>         DATUM["North_American_Datum_1983",<br>>             SPHEROID["GRS_1980",6378137,298.257222101]],<br>>         PRIMEM["Greenwich",0],<br>>         UNIT["Degree",0.0174532925199432955]],<br>>     PROJECTION["Lambert_Conformal_Conic_2SP"],<br>>     PARAMETER["False_Easting",609601.22],<br>>     PARAMETER["False_Northing",0],<br>>     PARAMETER["Central_Meridian",-79],<br>>     PARAMETER["Standard_Parallel_1",34.33333333333334],<br>>     PARAMETER["Standard_Parallel_2",36.16666666666666],<br>>     PARAMETER["Latitude_Of_Origin",33.75],<br>>     UNIT["Meter",1]]<br>><br>> `r.in.lidar -p` in both 7.2 and trunk gives:<br>><br>> +proj=lcc +lat_1=34.33333333333334 +lat_2=36.16666666666666 +lat_0=33.75 +lon_0=-79 +x_0=609601.22 +y_0=0 +datum=NAD83 +units=m +no_defs<br>><br>> `g.proj -p` gives the following in 7.2 and trunk in locations created by 7.2 and trunk:<br>><br>> -PROJ_INFO-------------------------------------------------<br>> name       : NAD83(HARN) / North Carolina<br>> datum      : nad83harn<br>> ellps      : grs80<br>> proj       : lcc<br>> lat_1      : 36.16666666666666<br>> lat_2      : 34.33333333333334<br>> lat_0      : 33.75<br>> lon_0      : -79<br>> x_0        : 609601.22<br>> y_0        : 0<br>> no_defs    : defined<br>> -PROJ_EPSG-------------------------------------------------<br>> epsg       : 3358<br>> -PROJ_UNITS------------------------------------------------<br>> unit       : meter<br>> units      : meters<br>> meters     : 1<br>><br>> I don't know how to resolve it besides using -o but I don't know how to justify it and why there is a difference (I don't see it in the recent proj-related commits).<br>><br>> Thank you,<br>> Vaclav<br>><br>> _______________________________________________<br>> grass-dev mailing list<br>> <a href="mailto:grass-dev@lists.osgeo.org">grass-dev@lists.osgeo.org</a><br>> <a href="https://lists.osgeo.org/mailman/listinfo/grass-dev">https://lists.osgeo.org/mailman/listinfo/grass-dev</a><br><br></div></div></div></div></div>