setValues for a RasterStack
* edit -- the title of this post should have been "setValues with a RasterBrick".
On Wed, Jun 13, 2012 at 12:23 PM, Jonathan Greenberg <jgrn at illinois.edu> wrote:
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
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