传单中的单行

Question

我是 R 的新手，正在尝试在传单地图上显示区域供暖系统。从另一个系统我得到了 Gauss-Krueger-Coordinates，我将其转换为地理坐标。数据集基本上包含具有 4 列的数据帧（2 用于起点/2 用于端点）。我试图在传单中使用 addpolylines 绘制网格，但它只是一团糟...... 这是我的一部分数据（几千个管段）：

 PIPES <- data.frame(longitude_start = c(10.16598309, 10.1306068, 10.13709072, 10.11585874, 10.13122298, 
                                        10.13046081, 10.12480974, 10.09430826, 10.10316342),
                    latitude_start = c(53.59022654, 53.60358092, 53.60814442, 53.59990009, 53.60351254,
                                       53.6035185, 53.60134658, 53.59273843, 53.59777827),
                    longitude_end = c(10.16541999, 10.13051495, 10.13704627, 10.10316342, 10.13117327,
                                      10.12480932, 10.1231769, 10.09388444, 10.09789269),
                    latitude_end = c(53.59022261, 53.60351639, 53.60789653, 53.59777827, 53.6035107,
                                     53.60133096, 53.60104698, 53.59258596, 53.59664915))

这是我使用的代码：

mydf2 <- data.frame(group = c("max", "min"),
                    lat = c(PIPES$latitude_start, PIPES$latitude_end),
                    long = c(PIPES$longitude_start, PIPES$longitude_end))

leaflet() %>%
addTiles() %>%
addPolylines(data = mydf2, lng = ~long, lat = ~lat, group = ~group)

两行应解释为单行，但在地图上它们连在一起：

有人知道我能做什么吗？

（2020 年 4 月 20 日编辑）

问候丹尼尔

Answer 1

海吉·丹尼尔，

看起来小组没有像预期的那样工作。您可以改用循环：

library(magrittr)
library(leaflet)
#> Warning: Paket 'leaflet' wurde unter R Version 3.5.3 erstellt

PIPES <- data.frame(longitude_start = c(10.16598309, 10.1306068, 10.13709072, 10.11585874, 10.13122298, 
                                        10.13046081, 10.12480974, 10.09430826, 10.10316342),
                    latitude_start = c(53.59022654, 53.60358092, 53.60814442, 53.59990009, 53.60351254,
                                       53.6035185, 53.60134658, 53.59273843, 53.59777827),
                    longitude_end = c(10.16541999, 10.13051495, 10.13704627, 10.10316342, 10.13117327,
                                      10.12480932, 10.1231769, 10.09388444, 10.09789269),
                    latitude_end = c(53.59022261, 53.60351639, 53.60789653, 53.59777827, 53.6035107,
                                     53.60133096, 53.60104698, 53.59258596, 53.59664915))

m <- leaflet() %>%
  addTiles() 

for (i in 1:nrow(PIPES)) {
  m <- addPolylines(m, lng = c(PIPES$longitude_start[i], PIPES$longitude_end[i]), 
                    lat = c(PIPES$latitude_start[i], PIPES$latitude_end[i]))
}
m

^{由 reprex package (v0.3.0)}

于 2020-04-21 创建

传单中的单行

Single Lines in leaflet

r

line

leaflet