2

由于mapviewgstatsp和其他软件包,我正在使用R作为GIS软件。

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->blue 50->red。

有什么办法吗?

4

1 回答 1

4

spplot与 in 中的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* 对象,但我认为这就是您所追求的。

于 2016-05-18T11:14:32.413 回答