[GRASS-stats] Error rgrass7 using readRAST()

[Francesco Tonini]

Roger,

Apparently the error was due to a mismatch with the default g.region. 
After checking in GRASS for g.region -p I found that it was set to 100 x 
101 columns. If any of you ever encounter the same problem, just make 
sure to set the g.region using the same input raster of interest, i.e. 
using g.region raster=/your_raster_name/.

Are there any plans to make use of the 'raster' package so that 
readRAST() reads a 'raster' R obj instead of a sp griddataframe? Same 
for writeRAST().

Thank you,
Francesco




On 3/18/2015 11:32 AM, Roger Bivand wrote:
> On Wed, 18 Mar 2015, Francesco Tonini wrote:
>
>> Dear all,
>>
>> I am testing rgrass7 and trying to read a single raster file calling 
>> R from GRASS. I am using the function readRAST(). Currently, this 
>> return a SpatialGridDataFrame R object but it seems to have a bug 
>> when reading the correct number of observations. In my case the input 
>> raster is 100 x 100 so it should have 10,000 observations (ncell). 
>> Instead, readRAST() read 100 x 101 with 10,100 observations. Did 
>> anybody encounter this issue or know what it could be caused by?
>> I can send over my input raster if that helps to test on your machines.
>
> Please make an archive file of (parts of) your location, including all 
> the upper-level files, but maybe excluding other data than the raster 
> file in question. In the packaged location, the same problem should 
> occur. Put the archive file somewhere from which it can be downloaded. 
> I'll take a look.
>
> Roger
>
>>
>> Any help is much appreciated.
>>
>> Thank you.
>> Francesco
>> _______________________________________________
>> grass-stats mailing list
>> grass-stats at lists.osgeo.org
>> http://lists.osgeo.org/mailman/listinfo/grass-stats
>>
>

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/grass-stats/attachments/20150318/880991cd/attachment.html>