[R-sig-Geo] error using filename argument with focalFilter() in package raster

[William McCoy]

Robert,

Thanks very much for the very quick bug fix.  And thanks for the raster 
package -- I find it very useful.

Bill


On 09/29/2010 04:45 PM, Robert J. Hijmans wrote:
> Bill, you are the not the one doing anything wrong. I fixed this bug
> in 'raster' version 1.5-12. Thanks for reporting it, Robert
>
> On Wed, Sep 29, 2010 at 12:17 PM, William McCoy<wdmccoy at geo.umass.edu>  wrote:
>> I have a raster DEM loaded in R and I am using the raster package.  My DEM
>> is called lakeside:
>>
>>> library(raster)
>> Loading required package: sp
>> raster version 1.5-8 (20-September-2010)
>>> lakeside
>> class       : RasterLayer
>> filename    :
>> nrow        : 341
>> ncol        : 321
>> ncell       : 109461
>> min value   : 1288.399
>> max value   : 1704.749
>> projection  : +proj=utm +zone=12 +ellps=WGS84 +datum=WGS84 +units=m +no_defs
>> +towgs84=0,0,0
>> xmin        : 344459.7
>> xmax        : 347669.7
>> ymin        : 4532817
>> ymax        : 4536227
>> xres        : 10
>> yres        : 10
>>
>>
>> When I use the focalFilter function with a filename argument, I get the
>> following error:
>>
>>> xslopeLS<- focalFilter(lakeside, matrix(filtHorn, nrow = 3)
>> + /(8*xres(lakeside)), filename = "xslopeLS")
>> Error in .local(object, values) :
>>   setValues no longer supports setting rows of values
>>
>>
>> However, if I call the same function without the filename argument, it seems
>> to work fine:
>>
>>> xslopeLS<- focalFilter(lakeside, matrix(filtHorn, nrow = 3)
>> + /(8*xres(lakeside)))
>>> xslopeLS
>> class       : RasterLayer
>> filename    :
>> nrow        : 341
>> ncol        : 321
>> ncell       : 109461
>> min value   : -1.468210
>> max value   : 1.196593
>> projection  : +proj=utm +zone=12 +ellps=WGS84 +datum=WGS84 +units=m +no_defs
>> +towgs84=0,0,0
>> xmin        : 344459.7
>> xmax        : 347669.7
>> ymin        : 4532817
>> ymax        : 4536227
>> xres        : 10
>> yres        : 10
>>
>>
>> Here is my filter:
>>
>>> filtHorn
>> [1] -1 -2 -1  0  0  0  1  2  1
>>
>>
>> And here is my sessionInfo():
>>
>>> sessionInfo()
>> R version 2.11.1 (2010-05-31)
>> x86_64-redhat-linux-gnu
>>
>> locale:
>>   [1] LC_CTYPE=en_US.UTF-8       LC_NUMERIC=C
>>   [3] LC_TIME=en_US.UTF-8        LC_COLLATE=en_US.UTF-8
>>   [5] LC_MONETARY=C              LC_MESSAGES=en_US.UTF-8
>>   [7] LC_PAPER=en_US.UTF-8       LC_NAME=C
>>   [9] LC_ADDRESS=C               LC_TELEPHONE=C
>> [11] LC_MEASUREMENT=en_US.UTF-8 LC_IDENTIFICATION=C
>>
>> attached base packages:
>> [1] stats     graphics  grDevices utils     datasets  methods   base
>>
>> other attached packages:
>> [1] raster_1.5-8 sp_0.9-69
>>
>> loaded via a namespace (and not attached):
>> [1] grid_2.11.1     lattice_0.19-11
>>
>>
>> Any idea what I'm doing wrong?
>>
>> Thanks for your help,
>> Bill
>>
>> --
>> William D. McCoy
>> Geosciences
>> University of Massachusetts
>> Amherst, MA  01003
>>
>> _______________________________________________
>> R-sig-Geo mailing list
>> R-sig-Geo at stat.math.ethz.ch
>> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
>>
>
> _______________________________________________
> R-sig-Geo mailing list
> R-sig-Geo at stat.math.ethz.ch
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo


-- 
William D. McCoy
Geosciences
University of Massachusetts
Amherst, MA 01003