Как управлять цветами карты в mapview (пакете) в R?

Я работаю с R как с программным обеспечением ГИС, благодаря 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-> красный.

есть ли способ сделать это?


r r-mapview
person Serhan    schedule 03.05.2016    source источник
comment
Почтовый индекс .... не картинки.   -  person IRTFM    schedule 04.05.2016
comment
нет кода, чтобы показать человека, только одна функция-ведьма: mapView (tempData, zcol = temp)   -  person Serhan    schedule 04.05.2016
comment
просто SpatialPixelsDataFrame   -  person Serhan    schedule 04.05.2016
comment
›Summary (tempData) Объект класса SpatialPixelsDataFrame Координаты: min max long -5,991 -5,741 lat 35,575 35,811 Прогнозируется: FALSE proj4string: [+ proj = longlat + datum = WGS84 + no_defs + ellps = WGS84 + towgs84 = 0,0,0 ] Количество точек: 8833 Атрибуты сетки: cellcentre.offset размер ячейки Cell.dim long -5.990 0.002 125 lat 35.576 0.002 118 Атрибуты данных: temp Мин. : 22,91 1-я четверть: 22,91 Медиана: 22,91 Среднее значение: 22,91 3-я четверть: 22,91 Макс. : 22.91   -  person Serhan    schedule 04.05.2016

Ответы (1)


arrow_upward
4
arrow_downward

Подобно 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 *, но я предполагаю, что это то, что вам нужно.

person TimSalabim    schedule 18.05.2016
comment
Я получаю эту ошибку! : Ошибка в листовке :: addCircleMarkers (m, lng = координаты (lst [[i]]) [, 1],: неиспользуемый аргумент (at = c (0, 0,25, 0,5, 0,75, 1)) - person Serhan; 20.05.2016
comment
Не могли бы вы предоставить минимальный воспроизводимый пример для воспроизведения ошибки, а также сообщить мне, какие версии mapview и брошюры вы используете? - person TimSalabim; 20.05.2016