[R-sig-Geo] ASCII grid file - extracting data after import with readGDAL()

Matthew Landis landis at isciences.com
Fri Aug 26 16:26:50 CEST 2011


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
~~~~~~~~~~~~~~~~~~~~~~~~~~



More information about the R-sig-Geo mailing list