Aggregate using Shapefile - R-SIG-Geo

Fri, Mar 22, 2013 10:13 AM #

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 "ALrst" in a way that I'm able to create a map in the end. I've used "extract", but am not sure then how to map that output. "aggregate" doesn't seem flexible enough to use a shapefile (or raster layer). Any suggestions?

Thanks,
Tiffany

Robert J. Hijmans

Mon, Mar 25, 2013 9:07 AM #

An embedded and charset-unspecified text was scrubbed...
Name: not available
URL: <https://stat.ethz.ch/pipermail/r-sig-geo/attachments/20130325/48422516/attachment.pl>

Tiffany Smith

Mon, Mar 25, 2013 9:32 AM #

An embedded and charset-unspecified text was scrubbed...
Name: not available
URL: <https://stat.ethz.ch/pipermail/r-sig-geo/attachments/20130325/f1fe2fc9/attachment.pl>

Robert J. Hijmans

Mon, Mar 25, 2013 12:33 PM #

An embedded and charset-unspecified text was scrubbed...
Name: not available
URL: <https://stat.ethz.ch/pipermail/r-sig-geo/attachments/20130325/cced0e00/attachment.pl>

Tiffany Smith

Tue, Mar 26, 2013 8:49 AM #

An embedded and charset-unspecified text was scrubbed...
Name: not available
URL: <https://stat.ethz.ch/pipermail/r-sig-geo/attachments/20130326/c18a3573/attachment.pl>