r.mapcalc not correctly determining null space around r.external data set
|Reported by:||dickeya||Owned by:|
Here's what I'm seeing:
I have a geotiff file that I link to using r.external (output=test_ras). I set the region to have 1px of null space around the input raster, and then run:
r.mapcalc "out_ras = if(isnull(test_ras),0,null())"
I expect out_rast to have all nulls set to 0, including the 1px buffer. The north and south borders are correctly detected as null and set to 0, but the east and west borders are still null in the output.
If I import the same geotiff using r.in.gdal and run the same routine, everything works as expected and I have a ring of 0s around my new data set.