单行中的单行

台风

我对R很陌生，想在传单地图上显示区域供暖系统。从另一个系统中，我获得了高斯-克鲁格坐标，并将其转换为地理坐标。数据集基本上包含带有4列的datatframe（起始点为2，端点为2）。我试图在传单中用addpolyline绘制网格，但这只是一团糟....这是我的数据的一部分（数千个管段）：

 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日编辑）

问候丹尼尔

比尔格

你好丹尼尔

看起来该小组无法正常工作。您可以改用循环：

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