Below is a chunk of code that plots a rasterStack with 6 layers; I am trying to
put labels a-f in the bottom right corner of each plot. The code below only
shows the letter “a”, how do I get different letters to show in the bottom left
corner of each of the 6 plots as I use levelplot?
I have a rasterBrick (temp) with global temperatures, and a rasterLayer
(regions) of global regions. I would like to create a rasterBrick of regional
average temperatures. So far I have the following:
region1-mask(temp,regions,maskvalue=1,inverse=TRUE)
mean1-cellStats(region1,mean)
for an alternative function/package.
Thanks so much,
Tiffany Smith
___
R-sig-Geo mailing list
R-sig-Geo@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
map
Robert
On Mon, Mar 25, 2013 at 9:32 AM, Tiffany Smith smithtiffa...@gmail.com
wrote:
Robert,
Thanks. I guess my question now is, how can I map the output from using
'extract'? I'd like to be able to map is as I would a Raster* layer.
Thanks,
Tiffany
On Mar 25, 2013, at 12:07
polygon. See
the examples. Alternatively, you can use rasterize and then use crosstab.
Robert
On Fri, Mar 22, 2013 at 10:13 AM, Tiffany Smith smithtiffa...@gmail.com
wrote:
Hi all,
I'm looking to use a shapefile to aggregate a raster layer. Here's what I
have so far:
indexAL-raster
Hi all,
I'm looking to use a shapefile to aggregate a raster layer. Here's what I have
so far:
indexAL-raster(myfile.nc)
ALshp-readShapePoly(ALzipcodes.shp)
ALrst-rasterize(ALshp,indexAL,fun=mean)
indexAL has a finer resolution than ALrst, and I'd like to aggregate
indexAL up to the scale of
Hello,
I'm trying to write out the values of a RasterBrick to a .txt file. Here's what
I have:
b-brick(file1.nc)
vals-getValues(b)
write.table(vals,file1.txt)
Where 'b' is a RasterBrick with dimensions [224,464,33] which corresponds to
[lat,lon,time]. The output, 'vals' has dimensions
and read in the file, unless it's
convenient to you.
Hph,
Thiago.
- Original Message -
From: Tiffany Smith smithtiffa...@gmail.com
To: r-sig-geo@r-project.org
Cc:
Sent: Thursday, August 2, 2012 3:26 PM
Subject: Re: [R-sig-Geo] Mann-Kendall test on raster data
An update, I
work:
myraster.raster-as.raster(myraster.table)
Error in UseMethod(as.raster) :
no applicable method for 'as.raster' applied to an object of class
data.frame
Thanks,
Tiffany
On Aug 1, 2012, at 6:24 PM, Tiffany Smith wrote:
Hello,
I'm trying to do a Mann-Kendall test on my RasterBrick data
Hello,
I'm trying to do a Mann-Kendall test on my RasterBrick data. So far, I've been
able to complete the Mann-Kendall test (via the Kendall package) by
extracting a series like so:
onelocation-rev(extract(myraster,SpatialPoints(cbind(-86,33
onelocation.mk-MannKendall(onelocation)
I need
10 matches
Mail list logo