If you want to see more info of the dataset, you can also view it with
str(map@data) or head(map@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 ~~~~~~~~~~~~~~~~~~~~~~~~~~ _______________________________________________ R-sig-Geo mailing list R-sig-Geo@r-project.org https://stat.ethz.ch/mailman/listinfo/r-sig-geo