Hi Edzer, Roger and list, is it intended that subset.SpatialPointsDataFrame() returns two additional columns with the coordinates in the data-slot of it's return value?
data(meuse.grid) names(meuse.grid) coordinates(meuse.grid) <- ~x+y names(meuse.grid) meuse.grid <- subset(meuse.grid, meuse.grid$dist > 0.01) names(meuse.grid) Regards, Tom _______________________________________________ R-sig-Geo mailing list [email protected] https://stat.ethz.ch/mailman/listinfo/r-sig-geo
