<div dir="ltr"><div><div><br><br>On Fri, Nov 10, 2017 at 6:43 PM, Markus Metz <<a href="mailto:markus.metz.giswork@gmail.com">markus.metz.giswork@gmail.com</a>> wrote:<br>><br>><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>> 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>> The comparison of datums is new in 7.3 and needs some refinement.<br>><br>> I will fix the comparison of swapped lat_1 and lat_2.<br><br></div>Fixed in trunk r71656.<br><br></div>The test for different datum names has been disabled again in trunk r71657. There are several different datum names in lib/gis/datum.table that apparently mean the same: same ellipsoid and same transformation parameters. Apparently, GRASS does not provide a datum name when converting projection information from GRASS to proj4 for reprojecting data.<br><div><div><br></div><div>Markus M</div><div><br></div><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>