[R-sig-Geo] subseting an spatial polygons data frame
Agustin Lobo
alobolistas at gmail.com
Thu Nov 12 15:51:24 CET 2009
I want to discard few polygons from an
sp. poly. DF. I'm trying things such as:
> str(x,max.level=2)
Formal class 'SpatialPolygonsDataFrame' [package "sp"] with 5 slots
..@ data :'data.frame': 42 obs. of 9 variables:
..@ polygons :List of 42
..@ plotOrder : int [1:42] 32 35 38 1 8 29 33 26 27 31 ...
..@ bbox : num [1:2, 1:2] 437951 4605986 459732 4633538
.. ..- attr(*, "dimnames")=List of 2
..@ proj4string:Formal class 'CRS' [package "sp"] with 1 slots
>
> x <- subset(x,row.names(x)!="1")
Error in `[.data.frame`(x at data, i, j, ..., drop = FALSE) :
undefined columns selected
but the output of row.names(x)!="1"
is correct, what am I doing wrong? Is subset() (from maptools) not
appropriate for objects of type SPDF ?
Thanks
Agus
More information about the R-sig-Geo
mailing list