7

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コードです

ui.Rにアクションボタンがあり、クリックされたときにのみプロットを更新したいので、isolateを使用しています。

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   
4

3 に答える 3

5

だから私は最終的にこれを理解しました。

ネストされた ggplot クエリで ggplot がデータ セットを見つけられるようにするには、問題のデータ セットを <<- を使用して割り当てる必要があります。

たとえば、データセット poa は、関数内で呼び出されたときに、このネストされたクエリで見つけることができません。

print(ggmap(map)) + geom_polygon(data = poa, aes(x = LON, y = LAT, group = order), alpha = .5, colour = "black", fill = NA))

したがって、関数で poa を使用する前に、次の行を使用してください

poa <<- poa

ヘルプから: 「演算子 <<- および ->> は通常、関数でのみ使用され、割り当てられている変数の既存の定義を親環境で検索します。」

于 2013-09-09T00:23:06.590 に答える
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)))

ここで注意すべき点がいくつかあります。

  1. group = order「行数が異なります」というエラーが発生するため、ディレクティブを削除しました。あなたが何を求めているのかわかりませんが、正しく見えません。
  2. poaという名前の変数に存在する環境を保存しますenv。これは、 を参照するすべての関数で使用する環境ですpoa
  3. パラメータを指定して、先ほど作成した環境を呼び出して使用するように指示しgeom_polygonます。ggplotenvironment=env
  4. そのコマンドの最後のパラメーターとして使用して、呼び出し環境で外部ggmap呼び出しを実行するように指示します。environment=environment()これがないと、env作成した変数を見つけることができません。

これらの変更により、物事は適切に機能しているように見えます。

于 2013-07-08T05:32:11.047 に答える
0

group = orderたとえば、に変更するとgroup = POAOBS、アプリは機能します。

于 2013-07-11T14:20:01.213 に答える