Shiny,R - 错误:无法从非数字矩阵派生坐标

时间:2017-09-09 17:52:04

标签: r shiny sp

我试图使用Riny Shiny包来绘制地图上的点,因为它们在一天中的不同点出现。

当Shiny服务器被告知采用子集化数据帧的坐标(lon和lat)时,会出现意外问题:

coordinates(selected.tweet.points.spdf) <- ~ lon + lat

有趣的是,这在Shiny之外起作用,但在Shiny中并不适合我。

Shiny窗口中显示的错误消息是:

Error: cannot derive coordinates from non-numeric matrix

如果你在一个不是空间对象的对象上调用coordinates()而没有告诉它要考虑什么坐标,这就是 出现的错误信息,如下所示:< / p>

coordinates(selected.tweet.points.spdf)

这是一个最低限度可重复的例子(无实际绘图):

library(shiny)
library(tigris)
library(dplyr)
library(leaflet)
library(sp)
library(ggmap)
library(maptools)
library(broom)
library(httr)
library(rgdal)
library(tidyr)

nyc.neighborhoods.url <- GET('http://data.beta.nyc//dataset/0ff93d2d-90ba-457c-9f7e-39e47bf2ac5f/resource/35dd04fb-81b3-479b-a074-a27a37888ce7/download/d085e2f8d0b54d4590b1e7d1f35594c1pediacitiesnycneighborhoods.geojson')
nyc.neighborhoods.spdf <- readOGR(content(nyc.neighborhoods.url, 'text'), 'OGRGeoJSON', verbose = F)
nyc.neighborhoods.df <- tidy(nyc.neighborhoods.spdf)

tweet.points.df <- data.frame(class = c("Not Private", "Not Private", "Private", "Private", "Private", "Private"),
                        lat = c(40.65514, 40.65514, 42.74662, 42.74662, 40.65514, 40.57238),
                        lon = c(-73.94878, -73.94878, -75.77004, -75.77004, -73.94878, -74.15395),
                        time = c("14:00", "14:00", "14:30", "14:30", "14:30", "14:30"))

shinyApp(
  # User interface
  ui = fluidPage(
    titlePanel("Social Media Post Privacy Classifier"),
    sidebarLayout(
      sidebarPanel(
        sliderInput("time", "Time of day (by half hour)",   
                    min = as.POSIXlt("2017-01-01 00:00:00", tz = "GMT"),   
                    max = as.POSIXlt("2017-01-01 23:30:00", tz = "GMT"),   
                    value = as.POSIXct("2017-01-01 00:00:00", tz = "GMT"),   
                    timeFormat="%H:%M", timezone = "+0000", step = 60 * 30, animate = T)
      ),
      mainPanel(
        textOutput("out")
      )
    )
  ),
  # Server
  server = function(input, output) {
    choropleth.selection <- reactive({
      selected.time <- format(input$time, '%H:%M')
      selected.tweet.points.df <- tweet.points.df[tweet.points.df$time == selected.time, ]
      selected.tweet.points.spdf <- selected.tweet.points.df
      coordinates(selected.tweet.points.spdf) <- ~ lon + lat
      proj4string(selected.tweet.points.spdf) <- proj4string(nyc.neighborhoods.spdf)
      selected.tweet.points.df.matches <- over(selected.tweet.points.spdf, nyc.neighborhoods.spdf)
      selected.tweet.points.df <- cbind(selected.tweet.points.df, selected.tweet.points.df.matches)
      selected.tweet.points.df <- selected.tweet.points.df %>% drop_na()
      print(selected.tweet.points.df$class)
    })
    output$out <- renderPrint({ choropleth.selection() })
  }
)

这是一个证明相同技术在Shiny之外工作的例子:

selected.time <- "14:30"
points.df <- data.frame(class = c("Not Private", "Not Private", "Private", "Private", "Private", "Private"),
                        lat = c(40.65514, 40.65514, 42.74662, 42.74662, 40.65514, 40.57238),
                        lon = c(-73.94878, -73.94878, -75.77004, -75.77004, -73.94878, -74.15395),
                        time = c("14:00", "14:00", "14:30", "14:30", "14:30", "14:30"))
points.df <- data.frame(points.df[points.df$time == selected.time, ])
coordinates(points.df) <- ~ lon + lat

1 个答案:

答案 0 :(得分:0)

selected.tweet.points.spdf为空时,即没有匹配时间时,会出现错误。在selected.time00:00时,初始运行就是这种情况 - 没有匹配的tweet.points.df$time

一种可能的解决方案是改变您的服务器功能:

server = function(input, output) {
    choropleth.selection <- reactive({
      selected.time <- format(input$time, '%H:%M')
      selected.tweet.points.df <- tweet.points.df[tweet.points.df$time == selected.time, ]
      selected.tweet.points.spdf <- selected.tweet.points.df
      validate(need(nrow(selected.tweet.points.df)>0, "No match")) # <=======
      coordinates(selected.tweet.points.spdf) <- ~ lon + lat
      proj4string(selected.tweet.points.spdf) <- proj4string(nyc.neighborhoods.spdf)
      selected.tweet.points.df.matches <- over(selected.tweet.points.spdf, nyc.neighborhoods.spdf)
      selected.tweet.points.df <- cbind(selected.tweet.points.df, selected.tweet.points.df.matches)
      selected.tweet.points.df <- selected.tweet.points.df %>% drop_na()
      print(selected.tweet.points.df$class)
    })
    output$out <- renderPrint({ choropleth.selection() })
  }