gpt4 book ai didi

r - Shiny 的 map View

转载 作者:行者123 更新时间:2023-12-03 20:16:25 24 4
gpt4 key购买 nike

我使用 mapView 创建了一个交互式 map
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) <- TRUE

idw <- idw(formula = temp ~ 1, locations = data_test, newdata = grd)
idw.output = as.data.frame(idw) # output is defined as a data table

names(idw.output)[1:3] <- c("long", "lat", "temp")
idw.output <- idw.output[,1:3]

coordinates(idw.output) <- ~long+lat
morocco <- readOGR("Data/morocco/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

结果:

working part

现在

我想转向 Shiny ,问题是 map View 没有渲染功能。
我尝试使用 fpView() bView() 但没有结果。

我对 Shiny 的尝试

代码:

用户界面
library(shiny)
library(shinydashboard)
library(mapview)

header <- dashboardHeader(title="Ardusky")

sidebar <- dashboardSidebar(
)

body <- dashboardBody(
# Define UI for application
fluidPage(
mainPanel(
mapview:::fpViewOutput("mapplot"),
mapview:::plainViewOutput("test")
))
)

ui <- dashboardPage(header, sidebar, body, skin="black")

服务器.R
library(shiny)
library(mapview)
library(ggplot2)
library(sp)
library(gstat)
library(rgdal)
library(RMySQL)

shinyServer(function(input, output, session) {

repInput <- reactive({
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) <- TRUE
idw <- idw(formula = temp ~ 1, locations = data_test, newdata = grd)
idw.output = as.data.frame(idw) # output is defined as a data table
names(idw.output)[1:3] <- c("long", "lat", "temp")
idw.output <- idw.output[,1:3]
coordinates(idw.output) <- ~long+lat
morocco <- readOGR("/home/bloodesu/Data/morocco/TNG", "TNG")
proj4string(idw.output)<-proj4string(morocco)
tempData <- idw.output[morocco,]
proj4string(data_test)<-proj4string(morocco)
gridded(tempData) <- TRUE
tempData
})

output$mapplot <- mapview:::renderfpView({
mapview:::fpView(repInput(), zcol = "temp")
})

output$test <- mapview:::renderPlainView({
mapview:::plainview(repInput(), zcol = "temp")
})


})

结果

enter image description here

结论

如您所见,只有 plainView 给出了一些可接受的结果,但没有传单支持

最佳答案

map 查看 Shiny 不是天生为彼此而生的。但是, map 查看 基于 传单所以我们可以利用 Shiny 来自 的支持传单 .诀窍是使用 设置 map 对象。 map 查看 然后调用@map renderLeaflet() 内的插槽( 传单 部分)

ui.R

library(shiny)
library(shinydashboard)
library(mapview)

header <- dashboardHeader(title="Ardusky")

sidebar <- dashboardSidebar(
)

body <- dashboardBody(
# Define UI for application
fluidPage(
mainPanel(
leafletOutput("mapplot"),
mapview:::plainViewOutput("test")
))
)

ui <- dashboardPage(header, sidebar, body, skin="black")

server.ui
library(shiny)
library(mapview)
library(sp)

shinyServer(function(input, output, session) {

data(meuse)
coordinates(meuse) <- ~x+y
proj4string(meuse) <- CRS("+init=epsg:28992")

data(meuse.grid)
coordinates(meuse.grid) <- ~x+y
proj4string(meuse.grid) <- CRS("+init=epsg:28992")
gridded(meuse.grid) <- TRUE

m <- mapview(meuse.grid, zcol = "dist") + meuse

output$mapplot <- renderLeaflet({
m@map
})

})

这能解决你的问题吗?

更新:
我刚刚添加了 mapviewOutputrenderMapview到github上的开发版本。这意味着我们现在可以跳过对 @map 的显式调用。 mapview 对象的插槽。所以像 output$mapplot <- renderMapview(m)现在应该可以工作了。

的开发版 map 查看 可以安装 devtools::install_github("environmentalinformatics-marburg/mapview", ref = "develop")
更新 (2018/04/01): renderMapviewmapviewOutput目前不工作!因此,调用 renderLeaflet({ m@map })目前是使用 Shiny 的 map View 的方式。

关于r - Shiny 的 map View ,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/36679944/

24 4 0
Copyright 2021 - 2024 cfsdn All Rights Reserved 蜀ICP备2022000587号
广告合作:1813099741@qq.com 6ren.com