在R中使用传单绘制旅程路径

Drawing journey path using leaflet in R

我正在创建一个 Shiny 仪表板,其中 dataframe 开始 longitude/latitude 和结束 longitude/latitude cooridnated 我在 R 中使用 leaflet package:

`m=leaflet()%>%
      addTiles() %>%
      addMarkers(lng=(data$Start_long[i:j]), lat=(data$Start_lat[i:j]),popup="Start") %>%
      addCircleMarkers(lng=(data$End_long[i:j]), lat=(data$End_lat[i:j]),popup="End",clusterOptions=markerClusterOptions())`

我想知道是否有一种方法可以加入由 public 运输路线协调的起点和终点(可能通过 google 地图 API 或图书馆功能或失败, 用直线连接坐标?

您可以addPolylines()查看地图。

它需要两个向量作为参数,一个用于纬度,一个用于经度,其中每一行都是一个 'waypoint'。
在不了解数据结构的情况下很难为您提供帮助。 MRE:

library(leaflet)
cities <- read.csv(textConnection("
City,Lat,Long,Pop
Boston,42.3601,-71.0589,645966
Hartford,41.7627,-72.6743,125017
New York City,40.7127,-74.0059,8406000
Philadelphia,39.9500,-75.1667,1553000
Pittsburgh,40.4397,-79.9764,305841
Providence,41.8236,-71.4222,177994
"))

leaflet() %>% 
    addTiles() %>% 
    addPolylines(lat = cities$Lat, lng = cities$Long)

您可以使用我的 googleway 包来获取 directions/routes,并将其绘制在 Google 地图

要使用 Google 的 API,您需要为每个要使用的 API 提供一个有效密钥。在这种情况下,您需要 directions key, and for plotting the map you'll want a maps javascript key

(如果您愿意,您可以生成一个密钥并为两个 API 启用它)

要调用方向 API 并在 R 中绘制它,您可以这样做

library(googleway)

api_key <- "your_directions_api_key"
map_key <- "your_maps_api_key"

## set up a data.frame of locations
## can also use 'lat/lon' coordinates as the origin/destination
df_locations <- data.frame(
  origin = c("Melbourne, Australia", "Sydney, Australia")
  , destination = c("Sydney, Australia", "Brisbane, Australia")
  , stringsAsFactors = F
)

## loop over each pair of locations, and extract the polyline from the result
lst_directions <- apply(df_locations, 1, function(x){
  res <- google_directions(
    key = api_key
    , origin = x[['origin']]
    , destination = x[['destination']]
  )

  df_result <- data.frame(
    origin = x[['origin']]
    , destination = x[['destination']]
    , route = res$routes$overview_polyline$points
  )
  return(df_result)
})

## convert the results to a data.frame
df_directions <- do.call(rbind, lst_directions)

## plot the map
google_map(key = map_key ) %>%
  add_polylines(data = df_directions, polyline = "route")


在 Shiny 应用程序中也类似

library(shiny)
library(shinydashboard)
library(googleway)

ui <- dashboardPage(
  dashboardHeader(),
  dashboardSidebar(),
  dashboardBody(
    textInput(inputId = "origin", label = "Origin"),
    textInput(inputId = "destination", label = "Destination"),
    actionButton(inputId = "getRoute", label = "Get Rotue"),
    google_mapOutput("myMap")
  )
)

server <- function(input, output){

  api_key <- "your_directions_api_key"
  map_key <- "your_maps_api_key"

  df_route <- eventReactive(input$getRoute,{

    print("getting route")

    o <- input$origin
    d <- input$destination

    return(data.frame(origin = o, destination = d, stringsAsFactors = F))
  })


  output$myMap <- renderGoogle_map({

    df <- df_route()
    print(df)
    if(df$origin == "" | df$destination == "")
      return()

    res <- google_directions(
      key = api_key
      , origin = df$origin
      , destination = df$destination
    )

    df_route <- data.frame(route = res$routes$overview_polyline$points)

    google_map(key = map_key ) %>%
      add_polylines(data = df_route, polyline = "route")
  })
}

shinyApp(ui, server)

我用"for loop"来解决这个问题,一条一条的画多段线就行了。 (对不起我的中文表达^_^) 例如:

for(i in 1:nrow(sz)){
     if(i<=nrow(sz) ){
      a <- as.numeric(c(sz[i,c(8,10)])); 
      b <- as.numeric(c(sz[i,c(9,11)]));
      A <- A %>% addPolylines(a,b,group=NULL,weight = 1,color = "brown",
                              stroke = TRUE,fill = NULL,opacity = 0.8)}

或者更复杂的

for(j in 0:23){if(j<=23)
  #j--切每小时数据
  j1 <- as.character(paste(j,"点",sep='')) 
sz <- sz121[sz121$h==j,]
sz_4 <- sz121[sz121$bi_state==4 &sz121$h==j ,]
sz_8 <- sz121[sz121$bi_state==8&sz121$h==j,]
#还原A
A <- leaflet(sz121) %>% amap() %>% addLabelOnlyMarkers(~s_lon,~s_lat) %>% 
  addLegend(title=j1,colors=NULL,labels =NULL,position="topleft")
A <- A %>%addCircleMarkers(data=sz_8,~s_lon,~s_lat,color="orange",fill=TRUE,fillColor = "red", opacity = 1,fillOpacity=0.8,
                           weight =1,radius = 10) %>%addCircleMarkers(data=sz_4,~s_lon,~s_lat,color="black",fill=TRUE,fillColor = "red", 
                                                                      opacity = 1,fillOpacity=0.8,weight =5,radius = 10 ) %>% 
  addCircleMarkers(data=sz_8,~e_lon,~e_lat,color="orange",fill=TRUE,fillColor = "blue", opacity = 1,fillOpacity=0.8,weight=1,radius = 10) %>%
  addCircleMarkers(data=sz_4,~e_lon,~e_lat,color="black",fill=TRUE,fillColor = "blue", opacity = 1,fillOpacity=0.8,weight =5,radius = 10 ) 
for(i in 1:nrow(sz)){
  #i--画路径
  if(i<=nrow(sz) ){
    a <- as.numeric(c(sz[i,c(8,10)]));
    b <- as.numeric(c(sz[i,c(9,11)]));
    A <- A %>% addPolylines(a,b,group=NULL,weight = 1,color = "brown",stroke = TRUE,fill = NULL,opacity = 0.8)
  }
  if(i==nrow(sz)){print(A)}
}
Sys.sleep(3)
}