<div dir="ltr"><div><div><div>Hi all,<br><br></div>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></div>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></div>7.2 considers this OK while trunk considers this different. I'm not sure which one is correct.<br><br>Here is the lasinfo reference:<br><div><div><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></div><div><br></div><div>`r.in.lidar -p` in both 7.2 and trunk gives:</div><div><br></div><div>+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</div><div><br></div><div>`g.proj -p` gives the following in 7.2 and trunk in locations created by 7.2 and trunk:<br></div><div><br></div><div>-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></div><div><br></div><div>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></div><div><br></div><div>Thank you,</div><div>Vaclav<br></div></div></div>