我使用创建交互式地图 mapView
该mapview()
功能适用于我的网格数据(SpatialPixelsDataFrame
):
library(sp)
library(ggplot2)
library(gstat)
library(rgdal)
library(mapview)
library(RMySQL)
con <- dbConnect(MySQL(),
user="root",
password="",
host="127.0.0.1",
dbname="rstudio")
data<-dbReadTable(con,"data")
on.exit(dbDisconnect(con))
data_test <- data
data_test$x <- data$long
data_test$y <- data$lat
coordinates(data_test) = ~x + y
x.range <- as.numeric(c(-5.99, -5.74))
y.range <- as.numeric(c(35.57, 35.81))
grd <- expand.grid(x = seq(from = x.range[1], to = x.range[2], by = 0.002),
y = seq(from = y.range[1], to = y.range[2], by = 0.002)) # expand points to grid
coordinates(grd) <- ~x + y
gridded(grd) …
Run Code Online (Sandbox Code Playgroud) 我正在使用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->红色。
有什么办法可以做到这一点吗?
我想将我的数据框的一些列保存为特定的格式(JavaScript格式).我试图从rjson包中使用toJSON(),但这样做不起作用.
我的结果应该是这样的:http://leaflet.github.io/Leaflet.markercluster/example/realworld.388.js
请帮助,并提前感谢.