I have CDL state files for 2010 at 30' meter resolution. I'm using
functions from the *raster* package:

After using '*reclass*' function to extract only the ~30 crops I care to
locate (approximately) (step 1),
I want to somehow (*step 2*) create a raster brick file with 30 layers, one
identifying the presence/absence (i.e. 1 vs. NA) for each crop, and then *
aggregate* (step 3) to a higher resolution (say 900 meter resolution), so
that location information for each crop type is retained in the lower
resolution raster brick. Then I will *merge* (step 4) the raster bricks
together to have a raster (mostly consisting of NA's - uncropped land),
with rough locational information for those 30 major crops.

The point is, I can implement all the steps *except for step 2*.
I read through the raster vignettes, and have found no elegant way to
aggregate to a raster with categorical data without (arbitrarily) losing
lots of information. *

The point is just to identify whether or not each crop (of the ~30 key
crops that I care about) is grown in each 900 meter raster cell.*

I hope the question has been clearly posed.
Thanks very much in advance!

Jacob Teter
Graduate Student Researcher
University of California, Davis

	[[alternative HTML version deleted]]

