I have the following problem I am trying to solve. I have a shapefile with segments of varying lengths. For each segment, I have attributes including segment id, year, time period, and other variables like traffic volume. Right now a user is able to select a year, time period, and segment id from a drop down and it will populate a table with the associated values for that filter. I have that portion implemented successfully.
I also want the user to be able to click any of my custom segments (added via addPolylines() command) on the map and it should render the same table as above. I would prefer to not have the segment id as an item in my dropdown menu items list but rather as a selection from a map. However, my implementation does not do the following:
- Display my segments on the OpenStreets base map.
- Allow me to select the segment I want from the map and update the table accordingly for the selected segment.
Any insights to address this? My code, data, and screenshot of my output are below.
TIA,
Krishnan
Code
library(shiny)
library(leaflet)
library(dplyr)
library(sf)
library(scales)
library(bslib)
# Data
# sample_data
# segment_data
ui <- page_sidebar(
title = "Cap Analysis",
sidebar = sidebar(
selectInput("yr1", "Year", choices = sample_data$Year),
selectInput("per1", "Time Period", choices = sample_data$Period),
selectInput("segs","Segment ID",choices = unique(c(sample_data$SegID,segment_data$SegID))),
width = 300
),
card(
leafletOutput("map")
),
card(
tableOutput("attrib"),fill = TRUE
)
)
server <- function(input, output, session) {
selected <- reactive(sample_data %>% filter(Year == input$yr1 & Period==input$per1 & SegID==input$segs))
#selected <- reactive(sample_data %>% filter(Year == input$yr1 & Period==input$per1))
selectmap <- reactive(segment_data %>% filter(SegID == input$segs))
output$map <- renderLeaflet({
selectmap %>% leaflet() %>%
addTiles() %>% # Add OpenStreetMap tiles
setView(lng = -76.538192, lat = 39.310621, zoom = 13)
})
observe({
leafletProxy("map") %>%
clearShapes() %>%
addPolylines(
data = segment_data,
lng = ~c(start_lon,end_lon),
lat = ~c(start_lat,end_lat),
layerId = ~SegID,
color = "red",
weight = 3
)
})
observeEvent(input$map_shape_click, {
click <- input$map_shape_click
if (!is.null(click)) {
sample_data %>%
filter(SegID == click$SegID)
}
})
output$attrib <- renderTable(
selected() %>% count(`Segment ID`= SegID, Direction = Direction, `Segment Type` = SegType, Segment=Segment,`Number of Lanes`= NumLane, `Segment Length (ft)`= Segment_Length_ft, Type=Type,`FFS (mph)`= FFS_mph, `Volume (vph)` = scales::comma(round(Vol_vph,digits=0),big.mark = ","), `Truck Percent` = Truck_Percent, LOS = LOS, sort = TRUE, width="100%") %>% subset(select=c(1:11))
)
}
# Run the application
shinyApp(ui = ui, server = server)
Sample Data
sample_data <- structure(list(SegID = c("210", "211", "211A", "212", "213",
"213A", "214", "201", "202", "202A", "203", "210", "211", "211A",
"212", "213", "213A", "214", "201", "202", "202A", "203", "210",
"211", "211A", "212", "213", "213A", "214", "201", "202", "202A",
"203", "210", "211", "211A", "212", "213", "213A", "214", "201",
"202", "202A", "203", "210", "211", "211A", "212", "213", "213A",
"214", "201", "202", "202A", "203", "210", "211", "211A", "212",
"213", "213A", "214", "201", "202", "202A", "203", "210", "211",
"211A", "212", "213", "213A", "214", "201", "202", "202A", "203",
"210", "211", "211A", "212", "213", "213A", "214", "201", "202",
"202A", "203", "210", "211", "211A", "212", "213", "213A", "214",
"201", "202", "202A", "203"), Direction = c("NB", "NB", "NB",
"NB", "NB", "NB", "NB", "SB", "SB", "SB", "SB", "NB", "NB", "NB",
"NB", "NB", "NB", "NB", "SB", "SB", "SB", "SB", "NB", "NB", "NB",
"NB", "NB", "NB", "NB", "SB", "SB", "SB", "SB", "NB", "NB", "NB",
"NB", "NB", "NB", "NB", "SB", "SB", "SB", "SB", "NB", "NB", "NB",
"NB", "NB", "NB", "NB", "SB", "SB", "SB", "SB", "NB", "NB", "NB",
"NB", "NB", "NB", "NB", "SB", "SB", "SB", "SB", "NB", "NB", "NB",
"NB", "NB", "NB", "NB", "SB", "SB", "SB", "SB", "NB", "NB", "NB",
"NB", "NB", "NB", "NB", "SB", "SB", "SB", "SB", "NB", "NB", "NB",
"NB", "NB", "NB", "NB", "SB", "SB", "SB", "SB"), SegType = c("Freeway",
"Freeway", "Ramp Diverge", "Freeway", "Freeway", "Ramp Diverge",
"Freeway", "Freeway", "Freeway", "Ramp Merge", "Freeway", "Freeway",
"Freeway", "Ramp Diverge", "Freeway", "Freeway", "Ramp Diverge",
"Freeway", "Freeway", "Freeway", "Ramp Merge", "Freeway", "Freeway",
"Freeway", "Ramp Diverge", "Freeway", "Freeway", "Ramp Diverge",
"Freeway", "Freeway", "Freeway", "Ramp Merge", "Freeway", "Freeway",
"Freeway", "Ramp Diverge", "Freeway", "Freeway", "Ramp Diverge",
"Freeway", "Freeway", "Freeway", "Ramp Merge", "Freeway", "Freeway",
"Freeway", "Ramp Diverge", "Freeway", "Freeway", "Ramp Diverge",
"Freeway", "Freeway", "Freeway", "Ramp Merge", "Freeway", "Freeway",
"Freeway", "Ramp Diverge", "Freeway", "Freeway", "Ramp Diverge",
"Freeway", "Freeway", "Freeway", "Ramp Merge", "Freeway", "Freeway",
"Freeway", "Ramp Diverge", "Freeway", "Freeway", "Ramp Diverge",
"Freeway", "Freeway", "Freeway", "Ramp Merge", "Freeway", "Freeway",
"Freeway", "Ramp Diverge", "Freeway", "Freeway", "Ramp Diverge",
"Freeway", "Freeway", "Freeway", "Ramp Merge", "Freeway", "Freeway",
"Freeway", "Ramp Diverge", "Freeway", "Freeway", "Ramp Diverge",
"Freeway", "Freeway", "Freeway", "Ramp Merge", "Freeway"), Segment = c("I-95 NB North of Ramp from Eastern Ave WB",
"I-95 NB Ramp to Moravia Rd NB", "I-95 NB Ramp to Moravia Rd NB",
"I-95 NB North of Ramp to Moravia Rd NB", "I-95 NB Ramp to Pulaski Hwy EB",
"I-95 NB Ramp to Pulaski Hwy EB", "I-95 NB North of Ramp to Pulaski Hwy EB",
"I-95 SB South of Ramp to I-895 SB", "I-95 SB Ramp from Pulaski Hwy WB",
"I-95 SB Ramp from Pulaski Hwy WB", "I-95 SB South of Ramp from Pulaski Hwy WB",
"I-95 NB North of Ramp from Eastern Ave WB", "I-95 NB Ramp to Moravia Rd NB",
"I-95 NB Ramp to Moravia Rd NB", "I-95 NB North of Ramp to Moravia Rd NB",
"I-95 NB Ramp to Pulaski Hwy EB", "I-95 NB Ramp to Pulaski Hwy EB",
"I-95 NB North of Ramp to Pulaski Hwy EB", "I-95 SB South of Ramp to I-895 SB",
"I-95 SB Ramp from Pulaski Hwy WB", "I-95 SB Ramp from Pulaski Hwy WB",
"I-95 SB South of Ramp from Pulaski Hwy WB", "I-95 NB North of Ramp from Eastern Ave WB",
"I-95 NB Ramp to Moravia Rd NB", "I-95 NB Ramp to Moravia Rd NB",
"I-95 NB North of Ramp to Moravia Rd NB", "I-95 NB Ramp to Pulaski Hwy EB",
"I-95 NB Ramp to Pulaski Hwy EB", "I-95 NB North of Ramp to Pulaski Hwy EB",
"I-95 SB South of Ramp to I-895 SB", "I-95 SB Ramp from Pulaski Hwy WB",
"I-95 SB Ramp from Pulaski Hwy WB", "I-95 SB South of Ramp from Pulaski Hwy WB",
"I-95 NB North of Ramp from Eastern Ave WB", "I-95 NB Ramp to Moravia Rd NB",
"I-95 NB Ramp to Moravia Rd NB", "I-95 NB North of Ramp to Moravia Rd NB",
"I-95 NB Ramp to Pulaski Hwy EB", "I-95 NB Ramp to Pulaski Hwy EB",
"I-95 NB North of Ramp to Pulaski Hwy EB", "I-95 SB South of Ramp to I-895 SB",
"I-95 SB Ramp from Pulaski Hwy WB", "I-95 SB Ramp from Pulaski Hwy WB",
"I-95 SB South of Ramp from Pulaski Hwy WB", "I-95 NB North of Ramp from Eastern Ave WB",
"I-95 NB Ramp to Moravia Rd NB", "I-95 NB Ramp to Moravia Rd NB",
"I-95 NB North of Ramp to Moravia Rd NB", "I-95 NB Ramp to Pulaski Hwy EB",
"I-95 NB Ramp to Pulaski Hwy EB", "I-95 NB North of Ramp to Pulaski Hwy EB",
"I-95 SB South of Ramp to I-895 SB", "I-95 SB Ramp from Pulaski Hwy WB",
"I-95 SB Ramp from Pulaski Hwy WB", "I-95 SB South of Ramp from Pulaski Hwy WB",
"I-95 NB North of Ramp from Eastern Ave WB", "I-95 NB Ramp to Moravia Rd NB",
"I-95 NB Ramp to Moravia Rd NB", "I-95 NB North of Ramp to Moravia Rd NB",
"I-95 NB Ramp to Pulaski Hwy EB", "I-95 NB Ramp to Pulaski Hwy EB",
"I-95 NB North of Ramp to Pulaski Hwy EB", "I-95 SB South of Ramp to I-895 SB",
"I-95 SB Ramp from Pulaski Hwy WB", "I-95 SB Ramp from Pulaski Hwy WB",
"I-95 SB South of Ramp from Pulaski Hwy WB", "I-95 NB North of Ramp from Eastern Ave WB",
"I-95 NB Ramp to Moravia Rd NB", "I-95 NB Ramp to Moravia Rd NB",
"I-95 NB North of Ramp to Moravia Rd NB", "I-95 NB Ramp to Pulaski Hwy EB",
"I-95 NB Ramp to Pulaski Hwy EB", "I-95 NB North of Ramp to Pulaski Hwy EB",
"I-95 SB South of Ramp to I-895 SB", "I-95 SB Ramp from Pulaski Hwy WB",
"I-95 SB Ramp from Pulaski Hwy WB", "I-95 SB South of Ramp from Pulaski Hwy WB",
"I-95 NB North of Ramp from Eastern Ave WB", "I-95 NB Ramp to Moravia Rd NB",
"I-95 NB Ramp to Moravia Rd NB", "I-95 NB North of Ramp to Moravia Rd NB",
"I-95 NB Ramp to Pulaski Hwy EB", "I-95 NB Ramp to Pulaski Hwy EB",
"I-95 NB North of Ramp to Pulaski Hwy EB", "I-95 SB South of Ramp to I-895 SB",
"I-95 SB Ramp from Pulaski Hwy WB", "I-95 SB Ramp from Pulaski Hwy WB",
"I-95 SB South of Ramp from Pulaski Hwy WB", "I-95 NB North of Ramp from Eastern Ave WB",
"I-95 NB Ramp to Moravia Rd NB", "I-95 NB Ramp to Moravia Rd NB",
"I-95 NB North of Ramp to Moravia Rd NB", "I-95 NB Ramp to Pulaski Hwy EB",
"I-95 NB Ramp to Pulaski Hwy EB", "I-95 NB North of Ramp to Pulaski Hwy EB",
"I-95 SB South of Ramp to I-895 SB", "I-95 SB Ramp from Pulaski Hwy WB",
"I-95 SB Ramp from Pulaski Hwy WB", "I-95 SB South of Ramp from Pulaski Hwy WB"
), NumLane = c(4, 4, 1, 4, 4, 1, 4, 4, 4, 1, 4, 4, 4, 1, 4, 4,
1, 4, 4, 4, 1, 4, 4, 4, 1, 4, 4, 1, 4, 4, 4, 1, 4, 4, 4, 1, 4,
4, 1, 4, 4, 4, 1, 4, 4, 4, 1, 4, 4, 1, 4, 4, 4, 1, 4, 4, 4, 1,
4, 4, 1, 4, 4, 4, 1, 4, 4, 4, 1, 4, 4, 1, 4, 4, 4, 1, 4, 4, 4,
1, 4, 4, 1, 4, 4, 4, 1, 4, 4, 4, 1, 4, 4, 1, 4, 4, 4, 1, 4),
Segment_Length_ft = c(5280, 1500, 1500, 5280, 1500, 420,
5280, 5280, 1500, 1070, 5280, 5280, 1500, 1500, 5280, 1500,
420, 5280, 5280, 1500, 1070, 5280, 5280, 1500, 1500, 5280,
1500, 420, 5280, 5280, 1500, 1070, 5280, 5280, 1500, 1500,
5280, 1500, 420, 5280, 5280, 1500, 1070, 5280, 5280, 1500,
1500, 5280, 1500, 420, 5280, 5280, 1500, 1070, 5280, 5280,
1500, 1500, 5280, 1500, 420, 5280, 5280, 1500, 1070, 5280,
5280, 1500, 1500, 5280, 1500, 420, 5280, 5280, 1500, 1070,
5280, 5280, 1500, 1500, 5280, 1500, 420, 5280, 5280, 1500,
1070, 5280, 5280, 1500, 1500, 5280, 1500, 420, 5280, 5280,
1500, 1070, 5280), Type = c("Level", "Level", "Level", "Level",
"Level", "Level", "Level", "Level", "Level", "Level", "Level",
"Level", "Level", "Level", "Level", "Level", "Level", "Level",
"Level", "Level", "Level", "Level", "Level", "Level", "Level",
"Level", "Level", "Level", "Level", "Level", "Level", "Level",
"Level", "Level", "Level", "Level", "Level", "Level", "Level",
"Level", "Level", "Level", "Level", "Level", "Level", "Level",
"Level", "Level", "Level", "Level", "Level", "Level", "Level",
"Level", "Level", "Level", "Level", "Level", "Level", "Level",
"Level", "Level", "Level", "Level", "Level", "Level", "Level",
"Level", "Level", "Level", "Level", "Level", "Level", "Level",
"Level", "Level", "Level", "Level", "Level", "Level", "Level",
"Level", "Level", "Level", "Level", "Level", "Level", "Level",
"Level", "Level", "Level", "Level", "Level", "Level", "Level",
"Level", "Level", "Level", "Level"), FFS_mph = c(60, 60,
50, 60, 60, 40, 60, 60, 60, 30, 60, 60, 60, 50, 60, 60, 40,
60, 60, 60, 30, 60, 60, 60, 50, 60, 60, 40, 60, 60, 60, 30,
60, 60, 60, 50, 60, 60, 40, 60, 60, 60, 30, 60, 60, 60, 50,
60, 60, 40, 60, 60, 60, 30, 60, 60, 60, 50, 60, 60, 40, 60,
60, 60, 30, 60, 60, 60, 50, 60, 60, 40, 60, 60, 60, 30, 60,
60, 60, 50, 60, 60, 40, 60, 60, 60, 30, 60, 60, 60, 50, 60,
60, 40, 60, 60, 60, 30, 60), Year = c(2015, 2015, 2015, 2015,
2015, 2015, 2015, 2015, 2015, 2015, 2015, 2018, 2018, 2018,
2018, 2018, 2018, 2018, 2018, 2018, 2018, 2018, 2021, 2021,
2021, 2021, 2021, 2021, 2021, 2021, 2021, 2021, 2021, 2015,
2015, 2015, 2015, 2015, 2015, 2015, 2015, 2015, 2015, 2015,
2018, 2018, 2018, 2018, 2018, 2018, 2018, 2018, 2018, 2018,
2018, 2021, 2021, 2021, 2021, 2021, 2021, 2021, 2021, 2021,
2021, 2021, 2015, 2015, 2015, 2015, 2015, 2015, 2015, 2015,
2015, 2015, 2015, 2018, 2018, 2018, 2018, 2018, 2018, 2018,
2018, 2018, 2018, 2018, 2021, 2021, 2021, 2021, 2021, 2021,
2021, 2021, 2021, 2021, 2021), Period = c("AM Peak", "AM Peak",
"AM Peak", "AM Peak", "AM Peak", "AM Peak", "AM Peak", "AM Peak",
"AM Peak", "AM Peak", "AM Peak", "AM Peak", "AM Peak", "AM Peak",
"AM Peak", "AM Peak", "AM Peak", "AM Peak", "AM Peak", "AM Peak",
"AM Peak", "AM Peak", "AM Peak", "AM Peak", "AM Peak", "AM Peak",
"AM Peak", "AM Peak", "AM Peak", "AM Peak", "AM Peak", "AM Peak",
"AM Peak", "PM Peak", "PM Peak", "PM Peak", "PM Peak", "PM Peak",
"PM Peak", "PM Peak", "PM Peak", "PM Peak", "PM Peak", "PM Peak",
"PM Peak", "PM Peak", "PM Peak", "PM Peak", "PM Peak", "PM Peak",
"PM Peak", "PM Peak", "PM Peak", "PM Peak", "PM Peak", "PM Peak",
"PM Peak", "PM Peak", "PM Peak", "PM Peak", "PM Peak", "PM Peak",
"PM Peak", "PM Peak", "PM Peak", "PM Peak", "WE Peak", "WE Peak",
"WE Peak", "WE Peak", "WE Peak", "WE Peak", "WE Peak", "WE Peak",
"WE Peak", "WE Peak", "WE Peak", "WE Peak", "WE Peak", "WE Peak",
"WE Peak", "WE Peak", "WE Peak", "WE Peak", "WE Peak", "WE Peak",
"WE Peak", "WE Peak", "WE Peak", "WE Peak", "WE Peak", "WE Peak",
"WE Peak", "WE Peak", "WE Peak", "WE Peak", "WE Peak", "WE Peak",
"WE Peak"), Vol_vph = c(2800, 2800, 425, 2375, 2375, 375,
2000, 5550, 5550, 600, 6150, NA, NA, NA, NA, NA, NA, NA,
NA, NA, NA, NA, 3225, 3225, 450, 2775, 2775, 450, 2325, 4725,
4725, 500, 5225, 6150, 6150, 600, 5550, 5500, 550, 5000,
2700, 2700, 325, 3025, NA, NA, NA, NA, NA, NA, NA, NA, NA,
NA, NA, 5975, 5975, 650, 5325, 5325, 650, 4675, 3300, 3300,
325, 3625, 4125, 4125, 350, 3775, 3775, 275, 3500, 3450,
3450, 250, 3700, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA,
NA, 3950, 3950, 300, 3650, 3650, 350, 3300, 3725, 3725, 275,
4000), Truck_Percent = c(9.7, 9.7, 2, 10.7, 10.7, 6.1, 7.3,
7.9, 7.9, 6.1, 7.7, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA,
NA, 12.7, 12.7, 1.7, 14, 14, 5.3, 15.3, 14.7, 14.7, 8, 14.1,
9.7, 9.7, 2, 10.7, 10.7, 6.1, 7.3, 7.9, 7.9, 6.1, 7.7, NA,
NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 12.7, 12.7, 1.7,
14, 14, 5.3, 15.3, 14.7, 14.7, 8, 14.1, 4.4, 4.4, 1, 4.7,
4.7, 5.8, 3.6, 3.2, 3.2, 2.3, 3.1, NA, NA, NA, NA, NA, NA,
NA, NA, NA, NA, NA, 6.7, 6.7, 0.8, 7.3, 7.3, 2.2, 7.9, 6.5,
6.5, 2.7, 6.2), LOS = c("B", "B", "B", "B", "B", "B", "B",
"D", "C", "C", "D", "B", "A", "A", "B", "B", "B", "A", "C",
"C", "C", "C", "B", "A", "A", "B", "B", "B", "B", "C", "C",
"C", "D", "D", "F", "F", "E", "D", "D", "D", "B", "B", "B",
"B", "C", "C", "C", "D", "D", "D", "C", "B", "B", "B", "B",
"D", "C", "C", "D", "D", "D", "C", "B", "B", "B", "C", "C",
"B", "B", "C", "C", "C", "C", "B", "B", "B", "B", "B", "A",
"A", "B", "B", "B", "B", "B", "B", "B", "C", "C", "A", "A",
"B", "B", "B", "B", "B", "B", "B", "C")), class = c("tbl_df",
"tbl", "data.frame"), row.names = c(NA, -99L))
Segment Data
segment_data <- structure(list(Shape_Leng = c(1188.14747869, 604.391081074, 845.943032786,
924.419150361, 1170.61433301, 1284.04786287, 592.277680455, 713.523651468,
95.6789913537, 110.510041468, 184.887359789), ID = c("214", "213",
"212", "211", "210", "203", "202", "201", "211A", "213A", "202A"
), SegID = c("214", "213", "212", "211", "210", "203", "202",
"201", "211A", "213A", "202A"), start_lat = c(39.306108, 39.306108,
39.305346, 39.299475, 39.298813, 39.299187, 39.309301, 39.309301,
39.298813, 39.306108, 39.309301), start_lon = c(-76.530686, -76.530686,
-76.530567, -76.530471, -76.530557, -76.530727, -76.531848, -76.531848,
-76.530557, -76.530686, -76.531848), end_lat = c(39.31416, 39.310007,
39.299475, 39.30554, 39.290808, 39.308053, 39.309029, 39.314113,
39.299475, 39.305346, 39.308053), end_lon = c(-76.530686, -76.530686,
-76.530567, -76.530471, -76.530557, -76.530727, -76.531848, -76.531848,
-76.530557, -76.530686, -76.531848), geometry = structure(list(
structure(c(-8519356.9876, -8519390.6853, -8519438.3123,
-8519512.6105, -8519617.4122, 4765614.3042, 4765830.5091,
4766061.9766, 4766350.5966, 4766772.7688), dim = c(5L, 2L
), class = c("XY", "LINESTRING", "sfg")), structure(c(-8519356.9876,
-8519347.8763, -8519357.7182, -8519372.9345, -8519376.8361,
-8519374.1049, -8519364.7411, -8519349.1346, -8519325.3348,
-8519267.7883, 4765614.3042, 4765620.2143, 4765731.774, 4765929.5856,
4765997.0834, 4766027.9061, 4766058.3386, 4766087.2106, 4766115.3021,
4766175.2652), dim = c(10L, 2L), class = c("XY", "LINESTRING",
"sfg")), structure(c(-8519343.8012, -8519339.2718, -8519321.1497,
-8519317.3396, -8519321.1497, -8519333.0071, 4765504.5837,
4765458.0392, 4765214.2149, 4764986.5575, 4764837.0086, 4764660.0102
), dim = c(6L, 2L), class = c("XY", "LINESTRING", "sfg")),
structure(c(-8519333.0071, -8519318.855, -8519305.4834, -8519287.0054,
-8519266.2801, -8519229.7454, -8519220.7717, -8519219.2111,
-8519223.893, -8519234.0372, -8519248.4731, -8519281.6368,
-8519324.1643, -8519346.5252, 4764660.0102, 4764686.1622,
4764829.5659, 4764937.4376, 4765026.8313, 4765186.3291, 4765235.4893,
4765271.7743, 4765317.8133, 4765365.4129, 4765405.9896, 4765462.1728,
4765515.2347, 4765532.5401), dim = c(14L, 2L), class = c("XY",
"LINESTRING", "sfg")), structure(c(-8519342.6305, -8519357.7182,
-8519375.7875, -8519549.531, 4764564.8164, 4764415.7613,
4764286.3254, 4763413.3404), dim = c(4L, 2L), class = c("XY",
"LINESTRING", "sfg")), structure(c(-8519361.5967, -8519350.3051,
-8519344.0869, -8519345.3062, -8519362.109, -8519396.4005,
-8519411.5147, -8519441.7311, 4764618.5569, 4764759.8832,
4764926.4331, 4765122.7328, 4765444.7299, 4765692.3906, 4765770.9106,
4765894.1418), dim = c(8L, 2L), class = c("XY", "LINESTRING",
"sfg")), structure(c(-8519486.3647, -8519495.4003, -8519513.7825,
-8519541.3392, -8519555.0558, -8519570.144, -8519589.6902,
-8519605.8072, -8519621.2383, -8519639.7557, -8519651.4148,
-8519659.6448, -8519665.4743, -8519668.5606, -8519664.7885,
-8519650.3861, -8519622.9529, -8519590.7189, -8519571.5157,
-8519550.9408, -8519528.3084, -8519505.676, 4766073.5608,
4766076.6252, 4766147.8562, 4766216.6312, 4766238.9207, 4766253.666,
4766262.2389, 4766263.9534, 4766262.9247, 4766254.3518, 4766244.0644,
4766231.0336, 4766215.6024, 4766199.1425, 4766178.9106, 4766148.0482,
4766105.8697, 4766065.7487, 4766049.2887, 4766039.6871, 4766034.8863,
4766034.5434), dim = c(22L, 2L), class = c("XY", "LINESTRING",
"sfg")), structure(c(-8519486.3647, -8519658.6161, 4766073.5608,
4766765.9808), dim = c(2L, 2L), class = c("XY", "LINESTRING",
"sfg")), structure(c(-8519342.6305, -8519333.0071, 4764564.8164,
4764660.0102), dim = c(2L, 2L), class = c("XY", "LINESTRING",
"sfg")), structure(c(-8519356.9876, -8519343.8012, 4765614.3042,
4765504.5837), dim = c(2L, 2L), class = c("XY", "LINESTRING",
"sfg")), structure(c(-8519486.3647, -8519441.7311, 4766073.5608,
4765894.1418), dim = c(2L, 2L), class = c("XY", "LINESTRING",
"sfg"))), n_empty = 0L, crs = structure(list(input = "WGS 84 / Pseudo-Mercator",
wkt = "PROJCRS[\"WGS 84 / Pseudo-Mercator\",\n BASEGEOGCRS[\"WGS 84\",\n ENSEMBLE[\"World Geodetic System 1984 ensemble\",\n MEMBER[\"World Geodetic System 1984 (Transit)\"],\n MEMBER[\"World Geodetic System 1984 (G730)\"],\n MEMBER[\"World Geodetic System 1984 (G873)\"],\n MEMBER[\"World Geodetic System 1984 (G1150)\"],\n MEMBER[\"World Geodetic System 1984 (G1674)\"],\n MEMBER[\"World Geodetic System 1984 (G1762)\"],\n MEMBER[\"World Geodetic System 1984 (G2139)\"],\n ELLIPSOID[\"WGS 84\",6378137,298.257223563,\n LENGTHUNIT[\"metre\",1]],\n ENSEMBLEACCURACY[2.0]],\n PRIMEM[\"Greenwich\",0,\n ANGLEUNIT[\"degree\",0.0174532925199433]],\n ID[\"EPSG\",4326]],\n CONVERSION[\"Popular Visualisation Pseudo-Mercator\",\n METHOD[\"Popular Visualisation Pseudo Mercator\",\n ID[\"EPSG\",1024]],\n PARAMETER[\"Latitude of natural origin\",0,\n ANGLEUNIT[\"degree\",0.0174532925199433],\n ID[\"EPSG\",8801]],\n PARAMETER[\"Longitude of natural origin\",0,\n ANGLEUNIT[\"degree\",0.0174532925199433],\n ID[\"EPSG\",8802]],\n PARAMETER[\"False easting\",0,\n LENGTHUNIT[\"metre\",1],\n ID[\"EPSG\",8806]],\n PARAMETER[\"False northing\",0,\n LENGTHUNIT[\"metre\",1],\n ID[\"EPSG\",8807]]],\n CS[Cartesian,2],\n AXIS[\"easting (X)\",east,\n ORDER[1],\n LENGTHUNIT[\"metre\",1]],\n AXIS[\"northing (Y)\",north,\n ORDER[2],\n LENGTHUNIT[\"metre\",1]],\n USAGE[\n SCOPE[\"Web mapping and visualisation.\"],\n AREA[\"World between 85.06°S and 85.06°N.\"],\n BBOX[-85.06,-180,85.06,180]],\n ID[\"EPSG\",3857]]"), class = "crs"), class = c("sfc_LINESTRING",
"sfc"), precision = 0, bbox = structure(c(xmin = -8519668.5606,
ymin = 4763413.3404, xmax = -8519219.2111, ymax = 4766772.7688
), class = "bbox"))), row.names = c(NA, -11L), class = c("sf",
"data.frame"), sf_column = "geometry", agr = structure(c(Shape_Leng = NA_integer_,
ID = NA_integer_, SegID = NA_integer_, start_lat = NA_integer_,
start_lon = NA_integer_, end_lat = NA_integer_, end_lon = NA_integer_
), class = "factor", levels = c("constant", "aggregate", "identity"
)))
Output
