Ciclovías de la Ciudad: cuánto se tarda, qué verde queda cerca y hasta dónde se llega

Cada cifra trae la línea de R que la produjo, copiable, tal como está en el script. Tocá cualquier función y se abre su ficha completa acá abajo: para qué sirve, cómo se llama, qué argumentos toma, qué versión y qué licencia. Las fichas se le piden al diccionario de funciones de R al abrir esta página, una consulta por función, y quedan guardadas para cuando las abrís. Lo único escrito por nosotros es la glosa en castellano: va rotulada como tal dentro de cada ficha.

== 1. construir_red (toda la ciudad) ==
Red construida: 2505 nodos, 2613 tramos, 286.8 km de ciclovia (componente conexa mayor, EPSG:32721).
== 2. matriz_od Obelisco / Pque Centenario / Pque Lezama ==
Matriz OD en metros por ciclovia (origenes x destinos):
     [,1] [,2] [,3]
[1,]    0 6009 3895
[2,] 6009    0 9024
[3,] 3895 9024    0
== 3. instalacion_mas_cercana desde el Obelisco ==
Paquita Bernardo: 0 m por ciclovia
La Cumparsita: 0 m por ciclovia
Plaza De La Republica: 91 m por ciclovia
== 4. area_servicio 2000 m desde el Obelisco ==
Area de servicio a 2000 m por red: 197 nodos alcanzables, poligono de 573.1 ha.
== tiempo total: 29.8 s ==

La obra no es anónima. Tiene una sección responsable, con nombre propio y dirección permanente en el registro público de ANS.

Lo que hay que construir no lo inventamos: lo escribió Esri. La extensión Network Analyst de ArcGIS es la especificación. Acá está medido cuánto de esa especificación cubre esta obra — con las 35 herramientas listadas una por una, incluidas las 20 que no cubre.

Todo lo que se usó es software libre y datos abiertos, y cada uno viene con su licencia y su versión exacta. Nada acá tiene costo de licencia ni ata a un proveedor.

Es la pregunta que decide si un municipio necesita un servidor o no. La respuesta la da el diccionario, función por función, consultada en vivo al abrir esta página.

Un recibo que sólo dice lo que salió bien no es un recibo. Esto es lo que quedó flojo en esta corrida, dicho por quien la hizo.

• Dos resultados dan 0 m y no es que estén encima del Obelisco. «Paquita Bernardo» y «La Cumparsita» salen a 0 metros porque el cálculo pega cada espacio verde al nodo de ciclovía más cercano, y esos dos cayeron en el mismo nodo que el origen. El número honesto de la pregunta «¿qué verde tengo cerca?» es el tercero: Plaza de la República, 91 m — que además cualquiera que conozca el Obelisco puede verificar caminando.
• El área de servicio es una envolvente, no la mancha real. Las 573,1 ha salen del polígono convexo que envuelve los 197 nodos alcanzables. ArcGIS dibuja un polígono ajustado a las calles; esta versión todavía no. Por eso esa herramienta figura como parcial y no como cubierta en la matriz de las 35.
• Las fichas del diccionario pueden ir una versión atrás de lo que corrió. Cuando pasa, la ficha lo dice arriba de todo, al lado del nombre de la función, y la cuenta exacta está en Los papeles. Es una diferencia menor, pero la página la muestra en vez de taparla.
• Siete de las 19 funciones del recibo no tienen ficha. El diccionario cubre 19 paquetes a nivel función, e igraph, tidygraph, dplyr y ellmer todavía no están adentro. Esas siete se usaron igual y están listadas con nombre y apellido, pero su firma no la verificamos nosotros. La página lo dice donde pasa, en lugar de inventar la ficha.
• El mapa del encabezado es una silueta, no un mapa. Está dibujado con los mismos 2.764 tramos que leyó la corrida, pero con una proyección simple para poder mirarlo. Las distancias nunca se midieron sobre ese dibujo: se midieron en EPSG:32721.
• Al abrir esta página, la consola del navegador muestra dos errores 404. Son las dos consultas al diccionario por paquetes que no cubre (igraph y tidygraph): el 404 es la respuesta correcta del servicio, no una falla. Lo decimos acá para que nadie confunda ese ruido con un problema real.
• Nadie la aprobó todavía. Esta obra está corrida y verificada, pero no tiene la firma de una persona de la comuna. Hasta que la tenga, es un borrador con recibo.
• Si la aprobás desde acá, la firma queda en este navegador y en ningún lado más. Los botones de arriba hacen lo que dicen —cambian el estado de la obra, y se puede deshacer—, pero esta página guarda ese estado en el almacenamiento local de tu equipo. En la plataforma, la misma acción se asienta en el tablero de la Comuna 1 con la fecha y con quién la puso. Lo decimos porque un recibo que se miente sobre su propia firma no sirve para nada.

Un recibo sirve si otro puede repetir la cuenta y obtener lo mismo. Acá está el código que corrió, no una descripción del código que corrió.

https://cdn.buenosaires.gob.ar/datosabiertos/datasets/transporte-y-obras-publicas/ciclovias/ciclovias.csv
https://cdn.buenosaires.gob.ar/datosabiertos/datasets/secretaria-de-desarrollo-urbano/espacios-verdes/espacio_verde_publico.geojson

# R 4.3.3 o superior
install.packages(c("sf", "sfnetworks", "tidygraph", "dplyr", "igraph", "ellmer"))
# el paquete propio, desde nuestro r-universe
install.packages("ansgis.geostat", repos = "https://emri-setup.r-universe.dev")

library(sf); library(sfnetworks); library(tidygraph); library(dplyr)

URL_CICLOVIAS <- "https://cdn.buenosaires.gob.ar/datosabiertos/datasets/transporte-y-obras-publicas/ciclovias/ciclovias.csv"
URL_VERDES    <- "https://cdn.buenosaires.gob.ar/datosabiertos/datasets/secretaria-de-desarrollo-urbano/espacios-verdes/espacio_verde_publico.geojson"
CRS_M <- 32721   # UTM 21S: metros para CABA

puntos_de_texto <- function(txt) {
  m <- do.call(rbind, lapply(strsplit(strsplit(txt, ";")[[1]], ","), as.numeric))
  st_transform(st_as_sf(data.frame(lon = m[, 1], lat = m[, 2]),
                        coords = c("lon", "lat"), crs = 4326), CRS_M)
}

# 1 · Build Network / Create Network Dataset  ->  2.505 nodos y 286,8 km
f_construir_red <- function(comuna = "") {
  df     <- utils::read.csv(URL_CICLOVIAS, encoding = "UTF-8")
  lineas <- st_as_sf(df, wkt = "geometry", crs = 4326)
  if (nzchar(comuna)) lineas <- lineas[lineas$comuna == as.integer(comuna), ]
  lineas <- st_transform(lineas, CRS_M)
  red <- as_sfnetwork(lineas, directed = FALSE) |>
    activate("edges") |> mutate(peso = edge_length()) |>
    activate("nodes") |> filter(group_components() == 1L)   # componente conexa mayor
  .ans$red <- red
  km <- sum(as.numeric(st_length(st_as_sf(red, "edges")))) / 1000
  sprintf("Red construida: %d nodos, %d tramos, %.1f km", igraph::vcount(red), igraph::ecount(red), km)
}

# 2 · Make OD Cost Matrix + Solve  ->  6.009 m y 3.895 m
f_matriz_od <- function(puntos) {
  p <- puntos_de_texto(puntos)
  m <- st_network_cost(.ans$red, from = p, to = p, weights = "peso")
  round(matrix(as.numeric(m), nrow = nrow(m)))
}

# 3 · Make Closest Facility + Solve  ->  91 m
f_instalacion_mas_cercana <- function(lon, lat, n = "3") {
  v      <- st_transform(st_make_valid(st_read(URL_VERDES, quiet = TRUE)), CRS_M)
  origen <- puntos_de_texto(paste0(lon, ",", lat))
  costos <- as.numeric(st_network_cost(.ans$red, from = origen,
                                       to = st_centroid(v), weights = "peso"))
  ord <- order(costos)[seq_len(min(as.integer(n), length(costos)))]
  sprintf("%.0f m por ciclovia", costos[ord])
}

# 4 · Make Service Area + Solve  ->  197 nodos y 573,1 ha
f_area_servicio <- function(lon, lat, corte_m = "1000") {
  origen <- puntos_de_texto(paste0(lon, ",", lat))
  costos <- as.numeric(st_network_cost(.ans$red, from = origen, weights = "peso"))
  alcanz <- which(costos <= as.numeric(corte_m))
  nodos  <- st_as_sf(.ans$red, "nodes")[alcanz, ]
  hull   <- st_convex_hull(st_union(nodos))
  sprintf("%d nodos alcanzables, poligono de %.1f ha",
          length(alcanz), as.numeric(st_area(hull)) / 10000)
}

cat(f_construir_red(""))
cat(f_matriz_od("-58.3816,-34.6037;-58.4358,-34.6064;-58.3697,-34.6280"))
cat(f_instalacion_mas_cercana("-58.3816", "-34.6037", "3"))
cat(f_area_servicio("-58.3816", "-34.6037", "2000"))