我制作了一个闪亮的应用程序,它使用ggmap返回静态地图。 但是,当我想覆盖邮政编码边界时,我遇到一个错误,ggplot无法找到数据集。
数据集poa是邮政编码边界的数据帧,即具有多边形ID的拉特和隆起。
我已经尝试添加environment = environment()
,但这并没有解决我的问题。
我知道数据存在,因为我调用print(str(poa))
打印到R控制台。
有人可以为我建议一个解决方案,以便ggplot可以访问poa数据帧吗? 我很抱歉这不是一个非常可重复的例子。
更新: 当我使用这段代码时,ggplot能够访问poa数据帧:
print(ggmap(map)) + geom_polygon(data = poa, aes(x = LON, y = LAT, group = order), alpha = .5, colour = "black", fill = NA))
但我需要对地图的基础层进行嵌套调用ggplot,当我这样做时ggplot无法找到数据
这是我的server.R代码
我正在使用isolate,因为我的ui.R中有一个动作按钮,我只希望在点击它时更新绘图。
library(shiny)
library(ggmap)
library(RODBC)
# Define server logic required to summarize and view the selected dataset
shinyServer(function(input, output) {
output$searchString <- renderText({
if (input$searchButton == 0)
return()
isolate({input$searchString})
})
mapSourceInput <- reactive({
switch(input$mapSource
, "google" = "google"
, "stamen" = "stamen")
})
mapTypeInput <- reactive({
switch(input$mapType
, "terrain" = "terrain"
, "satellite" = "satellite"
, "roadmap" = "roadmap"
, "hybrid" = "hybrid"
, "toner" = "toner"
, "watercolor" = "watercolor")
})
overlayInput <- reactive({
switch(input$overlay
, "postcodes" = "postcodes"
, "states" = "states"
, "nothing" = "nothing")
})
output$map <- renderPlot({
if (input$searchButton == 0)
return()
isolate({
if (overlayInput() == "nothing"){
map <- get_map(location = input$searchString, zoom = input$zoom, source = mapSourceInput(), maptype = mapTypeInput())
mapPlot <- ggmap(map)
print(mapPlot)
#return()
} else {
if (overlayInput() == "postcodes"){
#postcode boundaries
map <- get_map(location = input$searchString, zoom = input$zoom, source = mapSourceInput(), maptype = mapTypeInput())
poa <- structure(list(POAOBS = c(1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L)
, COORD_REF = 1:10
, COORD_POL = 1:10
, POLYGON = c(1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L)
, LON = c(144.951431274414, 144.956451416016, 144.95539855957, 144.955993652344, 144.958465576172, 144.956634521484, 144.956817626953, 144.954727172852, 144.957550048828, 144.958831787109)
, LAT = c(-37.8131675720215, -37.8117561340332, -37.8094863891602, -37.8058776855469, -37.8061485290527, -37.8021659851074, -37.8010902404785, -37.7994079589844, -37.7997169494629, -37.799861907959)
, POA = c("3000", "3000", "3000", "3000", "3000", "3000", "3000", "3000", "3000", "3000"))
, .Names = c("POAOBS", "COORD_REF", "COORD_POL", "POLYGON", "LON", "LAT", "POA")
, row.names = c(NA, 10L)
, class = "data.frame")
print(str(poa))
print(ggmap(map, base_layer = ggplot(data = poa, aes(x = LON, y = LAT), environment = environment()), extent = "normal", maprange = FALSE, environment = .GlobalEnv) +
geom_polygon(data = poa, aes(x = LON, y = LAT, group = order), alpha = .5, colour = "black", fill = NA) +
coord_map(projection = "mercator",
xlim = c(attr(map, "bb")$ll.lon, attr(map, "bb")$ur.lon),
ylim = c(attr(map, "bb")$ll.lat, attr(map, "bb")$ur.lat)))
} else {
if(overlayInput() == "states"){
return()
}}}
})
})
})
编辑:添加了ui.R
library(shiny)
# Define UI for dataset viewer application
shinyUI(pageWithSidebar(
# Application title.
headerPanel("The New Map App"),
# Sidebar with controls
sidebarPanel(
textInput("searchString", "Get Map Of", value = "melbourne, australia")
, selectInput("mapSource", "Choose a Map Source", choices = c("google", "stamen"))
, numericInput("zoom", "Zoom Level", 10)
, helpText("Note: An integer from 3 (continent) to 21 (building), default value 10 (city)")
, selectInput("mapType", "Choose a Map Type", choices = c("terrain", "satellite", "roadmap", "hybrid", "toner", "watercolor"))
, helpText("Note: Options available are 'terrain', 'satellite', 'roadmap', and 'hybrid' (google maps), 'watercolor', and 'toner' (stamen maps)")
, radioButtons("overlay", "Overlay Polygon",
list("Postcodes" = "postcodes"
, "States" = "states"
, "Nothing" = "nothing"))
, actionButton("searchButton", "Get Map")
, tags$style(type='text/css', "button#searchButton { margin-bottom: 9px; }")
),
#output panel
mainPanel(
h3(textOutput("searchString"))
, plotOutput("map")
)
))
sessionInfo的输出:
R version 3.0.1 (2013-05-16)
Platform: x86_64-pc-linux-gnu (64-bit)
locale:
[1] LC_CTYPE=en_AU.UTF-8 LC_NUMERIC=C LC_TIME=en_AU.UTF-8 LC_COLLATE=en_AU.UTF-8 LC_MONETARY=en_AU.UTF-8
[6] LC_MESSAGES=en_AU.UTF-8 LC_PAPER=C LC_NAME=C LC_ADDRESS=C LC_TELEPHONE=C
[11] LC_MEASUREMENT=en_AU.UTF-8 LC_IDENTIFICATION=C
attached base packages:
[1] stats graphics grDevices utils datasets methods base
other attached packages:
[1] RODBC_1.3-7 ggmap_2.3 ggplot2_0.9.3.1 shiny_0.6.0
loaded via a namespace (and not attached):
[1] bitops_1.0-5 caTools_1.14 colorspace_1.2-2 dichromat_2.0-0 digest_0.6.3 grid_3.0.1
[7] gtable_0.1.2 httpuv_1.0.6.3 labeling_0.1 mapproj_1.2-1 maps_2.3-2 MASS_7.3-26
[13] munsell_0.4 plyr_1.8 png_0.1-5 proto_0.3-10 RColorBrewer_1.0-5 Rcpp_0.10.4
[19] reshape2_1.2.2 RgoogleMaps_1.2.0.3 rjson_0.2.12 RJSONIO_1.0-3 scales_0.2.3 stringr_0.6.2
[25] tools_3.0.1 xtable_1.7-1
答案 0 :(得分:5)
所以我终于想出了这个。
要让ggplot能够在嵌套的ggplot查询中找到数据集,需要使用&lt;&lt; -
分配有问题的数据集例如,在函数内部调用时,无法在此嵌套查询中找到数据集poa。
print(ggmap(map)) + geom_polygon(data = poa, aes(x = LON, y = LAT, group = order), alpha = .5, colour = "black", fill = NA))
因此,在您需要在函数中使用poa之前,请使用此行
poa <<- poa
来自帮助: “运算符&lt;&lt; - 和 - &gt;&gt;通常仅在函数中使用,并导致通过父环境进行搜索,以查找正在分配的变量的现有定义。”
答案 1 :(得分:2)
这是一个棘手的问题。我可能会编辑这个答案,因为我对ggplot的环境嵌套有了更好的了解。但是这里的修改后的块似乎让它对我起作用了:
env <- environment()
print(ggmap(map, base_layer = ggplot(data = poa, aes(x = LON, y = LAT), environment=env), extent = "normal", maprange = FALSE, environment=environment()) +
geom_polygon(data = poa, aes(x = LON, y = LAT), alpha = .5, colour = "black", fill = NA, environment=env) +
coord_map(projection = "mercator",
xlim = c(attr(map, "bb")$ll.lon, attr(map, "bb")$ur.lon),
ylim = c(attr(map, "bb")$ll.lat, attr(map, "bb")$ur.lat)))
这里有几点需要注意。
group = order
指令,因为这会导致“行数不同”错误。不确定你在那之后是什么,但看起来不对。poa
存在的环境存储在名为env
的变量中。这是您要在引用poa
。geom_polygon
参数,指示对ggplot
和environment=env
的调用使用我刚刚创建的环境。ggmap
作为该命令中的最后一个参数,在调用环境中执行外部environment=environment()
调用。没有它,它将无法找到我们创建的env
变量。通过这些更改,事情似乎正常。
答案 2 :(得分:0)
当我将group = order
更改为例如group = POAOBS
时,应用就可以运行。