我使用 R 作为 GIS 软件,感谢 mapview , gstat , sp 和其他包。
我用mapView()
函数
m <- vgm(psill=.49,model="Sph",range=600000,nugget=3.8)
idw <- krige(formula = temp~1, locations = data_test, newdata = grd, model=m)
idw.output = as.data.frame(idw)
names(idw.output)[1:3] <- c("long", "lat", "temp")
coordinates(idw.output) <- ~long+lat
morocco <- readOGR("/opt/lampp/htdocs/ardusky/public/data/TNG", "TNG")
proj4string(idw.output)<-proj4string(morocco)
tempData <- idw.output[morocco,]
proj4string(data_test)<-proj4string(morocco)
gridded(tempData) <- TRUE
m<-mapView(tempData, zcol = "temp") + data_test
m
结果:
我想控制着色,例如0->蓝色50->红色。
有什么方法可以做到吗?
答案 0 :(得分:3)
与spplot
类似,at
中的mapview
参数是您要使用的内容:
library(mapview)
library(sp)
data(meuse.grid)
coordinates(meuse.grid) <- ~x+y
proj4string(meuse.grid) <- CRS("+init=epsg:28992")
gridded(meuse.grid) <- TRUE
mapview(meuse.grid, zcol = "dist", at = seq(0, 1, 0.25))
注意:在CRAN版本中,这仅适用于Raster *对象,但我认为这就是你所追求的。