久々に描画周りを使うとわすれていることが多くて悩ましい。
library(sp) #ポイントデータの場合 data(meuse) coordinates(meuse)=~x+y spplot(meuse, "cadmium", do.log = TRUE, key.space=list(x=0.2,y=0.9,corner=c(0,1)), cuts=c(1, 2, 3, 4, 10, 20),scales=list(draw=T)) #グリッドデータの場合 data(meuse.grid) coordinates(meuse.grid) = ~x+y proj4string(meuse.grid) <- CRS("+init=epsg:28992") gridded(meuse.grid) = TRUE spplot(meuse.grid["dist"], at = seq(-1,2,length=20), col.regions=bpy.colors()) #ポリゴンの場合 library(maptools) nc <- readShapePoly(system.file("shapes/sids.shp", package="maptools")[1], proj4string=CRS("+proj=longlat +datum=NAD27")) spplot(nc[,1],at=c(0,0.1,0.20)) #または ## plot of SpatialPolygonsDataFrame, using line densities names(nc) rrt <- nc$SID74/nc$BIR74 brks <- quantile(rrt, seq(0,1,1/7)) cols <- grey((length(brks):2)/length(brks)) dens <- (2:length(brks))*3 plot(nc, density=dens[findInterval(rrt, brks, all.inside=TRUE)]) #http://rspatial.r-forge.r-project.org/gallery/