<div dir="ltr"><div><br><br>On Mon, Apr 16, 2018 at 12:28 PM, Nikos Alexandris <<a href="mailto:nik@nikosalexandris.net">nik@nikosalexandris.net</a>> wrote:<br>><br>> * Markus Neteler <<a href="mailto:neteler@osgeo.org">neteler@osgeo.org</a>> [2018-04-15 22:49:09 +0200]:<br>><br>>> On Sun, Apr 15, 2018 at 9:46 PM, Markus Neteler <<a href="mailto:neteler@osgeo.org">neteler@osgeo.org</a>> wrote:<br>>>><br>>>> On Sun, Apr 15, 2018 at 8:48 PM, Veronica Andreo <<a href="mailto:veroandreo@gmail.com">veroandreo@gmail.com</a>> wrote:<br>>>>><br>>>>> Hi devs,<br>>>>><br>>>>> I just found something odd. I have a CHIRPS tif for Africa [0], but I would<br>>>>> only need to import Ghana, so I used r.in.gdal -r. However when I visualize<br>>>>> the data, most of it is nodata, while the full map of Africa has data all<br>>>>> over the continent. Then, I tried also with r.import extent=region. Same<br>>>>> result. The only thing providing what I expected to be the right result (a<br>>>>> subset from the full map) is importing the whole map and then using<br>>>>> r.mapcalc to clip (i need the space in my hard disk)<br>>><br>>><br>>> trying harder I can reproduce the issue (of course not NULL related).<br>>><br>>> So, I have no clue right now why it behaves differently.<br>><br>><br>> Hi,<br>><br>> no intention to add noise. There is something weird, here too.<br>><br>> # first things first<br>> grass -c /geo/geodata/examples/chirps-v2.0.2018.01.01.tif /geo/grassdb/onsight/chirps<br>><br>> # importing the whole map looks visually ok<br>> r.in.gdal input=/geo/geodata/examples/chirps-v2.0.2018.01.01.tif output=chirps_v2<br>><br>> # then, using the following box coordinates for Ghana from<br>> # <a href="https://en.wikipedia.org/wiki/User%3aThe_Anome/country_bounding_boxes">https://en.wikipedia.org/wiki/User%3aThe_Anome/country_bounding_boxes</a><br>> g.region w=-4.000 e=4.733 n=20.000 s=11.167 -p<br><br></div>This bbox is north of Ghana, not Ghana.<br><div><br>><br>> # to draw the box on-display<br>> v.in.region output=ghana_box<br>><br>> # and doing a partial import<br>> r.in.gdal -r input=/geo/geodata/examples/chirps-v2.0.2018.01.01.tif output=chirps_v2_subset<br>><br>> Regardless whether the box figures are not correct, while the imported<br>> portion is, indeed, a box fraction of the whole map, it isn't, visually<br>> the "right" part. It appears as a part from the coast has been imported.<br><br></div><div>Apparently you did not update your local copy of trunk or relbr74, see my previous post.<br></div><div><br></div><div>Markus M<br></div><div><br></div></div>