# NOT RUN {
require(sp)
grd <- GridTopology(c(1,1), c(1,1), c(5,5))
polys <- as(grd, "SpatialPolygons")
centroids <- coordinates(polys)
gdomain <- SpatialPolygonsDataFrame(polys, data=data.frame(x=centroids[,1],
y=centroids[,2], row.names=row.names(polys)))
feat <- gdomain$x*5+gdomain$y^2
plotMap(gdomain,feat)
# }
Run the code above in your browser using DataLab