If you want to see more info of the dataset, you can also view it with
str(map at data)
or
head(map at data)
Matt
On 8/26/2011 9:57 AM, Forrest Stevens wrote:
As you said, readGDAL will give you a SpatialGridDataFrame object, which
behaves much like a normal data.framet. ?Your best bet is to run it
through
names to see what data are available to work with:
names(map)
Like a regular data.frame this will give you the bands or data read in
from
your ASCII grid, on which you can pull subsets or index like you would a
normal vector (e.g. pretend your map object has "band1" as its only data,
you can pull a subset of the first 10 grid values like this:
map$band1[1:10]
I hope that gets you closer to solving the rest of your question,
--
~~~~~~~~~~~~~~~~~~~~~~~~~~
Matthew Landis, Ph.D.
Research Scientist
ISciences, LLC
61 Main St. Suite 200
Burlington VT 05401
802.864.2999
~~~~~~~~~~~~~~~~~~~~~~~~~~