Hello all,<div><br></div><div>I&#39;m having an issue where r.mapcalc is incorrectly detecting nodata when using an r.external raster.  My basic work flow is something like this:</div><div><ul><li>import raster using r.external</li>
<li>set the region to include a 1 pixel buffer around the image</li><li>run r.mapcalc with isnull to output nodata as 1.  The 1px buffer ensures that I get a single, continuous nodata region around the raster.</li></ul><div>
With certain r.external rasters, r.mapcalc is only returning nodata for the first and last rows.  The first and last columns have the top and bottom pixels (first and last row), but nothing else.</div><div><br></div><div>
If I do the exact same routine with r.in.gdal, a continuous nodata region around the raster is identified by r.mapcalc.  </div><div><br></div><div>I reported this bug last year, but stopped seeing the problem and didn&#39;t follow up.  Looks like it&#39;s still there, I&#39;ll attach some files for testing.</div>
<div><br></div><div><a href="http://trac.osgeo.org/grass/ticket/369">http://trac.osgeo.org/grass/ticket/369</a></div><div><br></div><div>- Jamie</div><div><br></div><div><br></div></div>