我正在使用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
Run Code Online (Sandbox Code Playgroud)
结果:
我想控制颜色,例如0->蓝色50->红色。
有什么办法可以做到这一点吗?
与您想要使用的参数spplot类似:atmapview
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))
Run Code Online (Sandbox Code Playgroud)
注意:在 CRAN 版本中,这只适用于 Raster* 对象,但我认为这就是您所追求的。