<div dir="ltr"><br><br>On Thu, Mar 14, 2019 at 3:16 PM Nikos Alexandris <<a href="mailto:nik@nikosalexandris.net">nik@nikosalexandris.net</a>> wrote:<br>><br>> Following up, why are there differences between GDAL and GRASS GIS in<br>> the following example?<br>><br>> This <a href="ftp://ftp.soilgrids.org/data/aggregated/5km/OCDENS_M_sl1_5km_ll.tif">ftp://ftp.soilgrids.org/data/aggregated/5km/OCDENS_M_sl1_5km_ll.tif</a><br>> raster map, subject to `gdalinfo`:<br>> ```<br>> gdalinfo OCDENS_M_sl1_5km_ll.tif -nogcp -nomd -norat -noct -nofl<br>> Driver: GTiff/GeoTIFF<br>> Files: OCDENS_M_sl1_5km_ll.tif<br>> Size is 7200, 2987<br>> Coordinate System is:<br>> GEOGCS["WGS 84",<br>>     DATUM["WGS_1984",<br>>         SPHEROID["WGS 84",6378137,298.257223563,<br>>             AUTHORITY["EPSG","7030"]],<br>>         AUTHORITY["EPSG","6326"]],<br>>     PRIMEM["Greenwich",0],<br>>     UNIT["degree",0.0174532925199433],<br>>     AUTHORITY["EPSG","4326"]]<br>> Origin = (-180.000000000000000,87.370000000000005)<br>> Pixel Size = (0.050000000000000,-0.050000000000000)<br>> Corner Coordinates:<br>> Upper Left  (-180.0000000,  87.3700000) (180d 0' 0.00"W, 87d22'12.00"N)<br>> Lower Left  (-180.0000000, -61.9800000) (180d 0' 0.00"W, 61d58'48.00"S)<br>> Upper Right ( 180.0000000,  87.3700000) (180d 0' 0.00"E, 87d22'12.00"N)<br>> Lower Right ( 180.0000000, -61.9800000) (180d 0' 0.00"E, 61d58'48.00"S)<br>> Center      (   0.0000000,  12.6950000) (  0d 0' 0.00"E, 12d41'42.00"N)<br>> Band 1 Block=7200x1 Type=Int16, ColorInterp=Gray<br>>   NoData Value=-32768<br>> ```<br>><br>> and GRASS GIS<br>> ```<br>> # create a new Location<br>> grass -c OCDENS_M_sl1_5km_ll.tif /geoyeux/grassdb/global/soil_grids/<br>><br>> # projection info<br>> g.proj -g<br>><br>> name=WGS 84<br>> datum=wgs84<br>> ellps=wgs84<br>> proj=ll<br>> no_defs=defined<br>> epsg=4326<br>> unit=degree<br>> units=degrees<br>> meters=1.0<br>><br>> # raster info<br>> <a href="http://r.info">r.info</a> -g OCDENS_M_sl1_5km_ll<br>><br>> north=87.37<br>> south=-61.98<br>> east=180<br>> west=-180<br>> nsres=0.05<br>> ewres=0.05<br>> rows=2987<br>> cols=7200<br>> cells=21506400<br>> datatype=CELL<br>> ncats=0<br>><br>> # report non-NULL cells and their x, y grid location<br>> r.stats OCDENS_M_sl1_5km_ll -n -x > stats_x<br>> ```<br>><br>> Comparing a few single pixels via:<br>> ```<br>> while read LINE;do<br>>     set -- $LINE<br>>     echo " ($1,$2)<br>>     echo "GDAL:  $(gdallocationinfo -valonly OCDENS_M_sl1_5km_ll.tif $1 $2)"<br>>     echo "GRASS: $3"<br>>     echo<br>> done < stats_x_head<br>> ```<br>><br>> gives<br>> ```<br>>  (2930,77)<br>> GDAL:  2090<br>> GRASS: 2096<br>><br>>  (2931,77)<br>> GDAL:  2055<br>> GRASS: 2090<br>><br>>  (2932,77)<br>> GDAL:  2063<br>> GRASS: 2055<br>><br>>  (2933,77)<br>> GDAL:  2093<br>> GRASS: 2063<br>><br>>  (2934,77)<br>> GDAL:  2240<br>> GRASS: 2093<br>><br>>  (2935,77)<br>> GDAL:  2332<br>> GRASS: 2240<br>><br>>  (2936,77)<br>> GDAL:  2296<br>> GRASS: 2332<br>><br>>  (2937,77)<br>> GDAL:  2253<br>> GRASS: 2296<br>><br>>  (2938,77)<br>> GDAL:  2179<br>> GRASS: 2253<br>><br>>  (2939,77)<br>> GDAL:  2115<br>> GRASS: 2179<br>> ```<br>><br><div>> Why these differences?</div><div><br></div><div>With r.stats -x, indexing starts with 1 (first row is 1).</div><div>With gdallocationinfo, indexing starts with 0 (first row is 0).</div><div><br></div><div>Markus M<br></div>><br>> Nikos<br>><br>><br>> # meta<br>><br>> GRASS 7.7.svn (2019)<br>> libgis Revision: 74118<br>> libgis Date: 2019-02-21 10:38:28 +0100 (Thu, 21 Feb 2019)<br>><br>> GDAL 2.3.1, released 2018/06/22<br>> _______________________________________________<br>> grass-user mailing list<br>> <a href="mailto:grass-user@lists.osgeo.org">grass-user@lists.osgeo.org</a><br>> <a href="https://lists.osgeo.org/mailman/listinfo/grass-user">https://lists.osgeo.org/mailman/listinfo/grass-user</a></div>