2016-05-03 44 views
0

我与[R工作作为一个GIS软件,感谢MapView的gstatSP包等。如何在R中控制mapview(package)中的地图颜色?

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 

结果绘制结果:

enter image description here

我想控制着色,0->蓝色50->红色例子。

有什么办法可以做到这一点?

+1

邮政编码....不是图片。 –

+0

没有代码显示man,只有一个函数是:mapView(tempData,zcol =“temp”) – Serhan

+0

什么是tempData? – MLavoie

回答

2

spplotat论点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版本只针对光栅作品*的对象,但我认为这是你所追求的。

+0

我得到这个错误! : 错误小叶:: addCircleMarkers(米,LNG =坐标(LST [[I]])[1],: 未使用的参数(在= C(0,0.25,0.5,0.75,1)) – Serhan

+0

莫非你提供一个最小可重现的例子来重现错误,并让我知道你正在使用哪个版本的mapview和传单? – TimSalabim