[R-sig-Geo] setValues for a RasterStack

[Jonathan Greenberg]

R-sig-geo'ers:

I'm trying to use setValues to replace values of a layer in a
RasterBrick, but I'm getting the following error:

Error in setValues(x, values[startidx[i]:endidx[i]], layer = i) :
  you can only setValues for a single layer if all values are in
memory. But values could not be loaded

The "empty" raster brick (note the "values" is none) is on disk:

class       : RasterBrick
dimensions  : 319, 10853, 3462107, 6  (nrow, ncol, ncell, nlayers)
resolution  : 2, 2  (x, y)
extent      : 748005.1, 769711.1, 4307056, 4307694  (xmin, xmax, ymin, ymax)
coord. ref. : +proj=utm +zone=10 +datum=WGS84 +units=m +no_defs
+ellps=WGS84 +towgs84=0,0,0
values      : none
layer names : layer.1, layer.2, layer.3, layer.4, layer.5, layer.6

> inMemory(x)
FALSE

"values" is a "squashed" vector (BSQ) of all values I wish to write to
the RasterBrick of length ncell(x) * nlayers(x), i is 1:nlayers(x),
and startidx and endidx are the start/stop points for this vector for
each band i:

startidx:
[1]        1  3462108  6924215 10386322 13848429 17310536
endidx:
[1]  3462107  6924214 10386321 13848428 17310535 20772642

I'm trying to setValues(x) one band at a time:

for(i in 1:nlayers(x))
{
	out <- setValues(x,values[startidx[i]:endidx[i]],layer=i)
}

What am I doing wrong?  Does setValues ONLY work if inMemory(x)==TRUE?  Thanks!

--j

-- 
Jonathan A. Greenberg, PhD
Assistant Professor
Department of Geography and Geographic Information Science
University of Illinois at Urbana-Champaign
607 South Mathews Avenue, MC 150
Urbana, IL 61801
Phone: 415-763-5476
AIM: jgrn307, MSN: jgrn307 at hotmail.com, Gchat: jgrn307, Skype: jgrn3007
http://www.geog.illinois.edu/people/JonathanGreenberg.html