r - 使用 gcIntermediate 在 R 传单/Shiny 中映射跨越日期变更线的最短飞行路径

Question

This question already has answers here:

Plotting routes that cross the international dateline using leaflet library in R

(2 个回答)


2年前关闭。




我正在使用 R-Leaflet 创建澳大利亚机场及其国际目的地的 map 。

这是我的示例数据:

df<-data.frame("Australian_Airport" = "Brisbane", 
           "International" =  c("Auckland", "Bandar Seri Begawan","Bangkok","Christchurch","Denpasar","Dunedin","Hamilton","Hong Kong","Honiara","Kuala Lumpur"),
           "Australian_lon" = c(153.117, 153.117,153.117,153.117,153.117,153.117, 153.117, 153.117, 153.117, 153.117),
           "Australian_lat" = c(-27.3842,-27.3842,-27.3842,-27.3842,-27.3842,-27.3842, -27.3842, -27.3842, -27.3842, -27.3842),
           "International_lon" = c(174.7633, 114.9398, 100.5018, 172.6362, 115.2126,-82.77177, -84.56134, 114.10950, 159.97290, 101.68685),
           "International_lat" = c(-36.848460, 4.903052, 13.756331, -43.532054,-8.670458,28.019740, 39.399501, 22.396428, -9.445638,  3.139003)
           )

我认为使用 gcIntermediate 使用弯曲的飞行路径会很酷，所以我创建了一个 SpatialLines 对象:

library(rgeos)
library(geosphere)

p1<-as.matrix(df[,c(3,4)])

p2<-as.matrix(df[,c(5,6)])

df2 <-gcIntermediate(p1, p2, breakAtDateLine=F, 
                    n=100, 
                    addStartEnd=TRUE,
                    sp=T) 

然后我使用传单和 Shiny 绘制它:

server <-function(input, output) {

airportmap<- leaflet() %>% addTiles() %>% 
    addCircleMarkers(df, lng = df$Australian_lon, lat = df$Australian_lat, 
    radius = 2, label = paste(df$Australian_Airport, "Airport"))%>% 
    addPolylines(data = df2, weight = 1)

output$mymap <- renderLeaflet({airportmap}) # render the base map
  }


ui<-  navbarPage("International flight path statistics - top routes",
      tabPanel("Interactive map",

      leafletOutput('mymap',  width="100%", height=900)

         )
         )

# Run the application 
shinyApp(ui = ui, server = server)

它看起来像这样:



因此，如果它们越过日期变更线，则路径是不正确的。将 breakAtDateLine 更改为 FALSE 并不能修复它(行消失但路径仍然断开)。在这个阶段，我怀疑我可能需要使用不同的映射系统或其他东西，但如果有人有一些建议，我将不胜感激。

提前致谢。

最佳答案

如果您对另一个 map 库感兴趣，那么 googleway 使用 Google Maps，根据我的经验，它更擅长处理跨越日期变更线的线。

笔记

要使用 Google map ，您需要一个 API key

目前只支持 sf 对象，不支持 sp

这也适用于 Shiny 的；我只是在这里向您展示基本 map

我编写了 googleway

library(sf)
library(googleway)

## convert the sp object to sf
sf <- sf::st_as_sf(df2)

set_key("your_api_key")

google_map() %>%
    add_polylines(data = sf)

关于r - 使用 gcIntermediate 在 R 传单/Shiny 中映射跨越日期变更线的最短飞行路径，我们在Stack Overflow上找到一个类似的问题： https://stackoverflow.com/questions/48857748/

上一篇：regex - 导入名称不同月份的文件

下一篇：perl - Perl $$(双美元)符号-它取消引用什么？

Answer 1

如果您对另一个 map 库感兴趣，那么 googleway 使用 Google Maps，根据我的经验，它更擅长处理跨越日期变更线的线。

笔记

要使用 Google map ，您需要一个 API key

目前只支持 sf 对象，不支持 sp

这也适用于 Shiny 的；我只是在这里向您展示基本 map

我编写了 googleway

library(sf)
library(googleway)

## convert the sp object to sf
sf <- sf::st_as_sf(df2)

set_key("your_api_key")

google_map() %>%
    add_polylines(data = sf)