"Todas" las terrazas de Madrid

Números

es un mapa en el que, en rojo, figuran todas (véase la coda) las terrazas de Madrid. Los datos están extraídos del censo de locales, sus actividades y terrazas de hostelería y restauración del ayuntamiento y están procesados con

terrazas <- fread("http://datos.madrid.es/egob/catalogo/200085-17-censo-locales.txt")
terrazas$coordenada_x_local <- as.numeric(gsub(",", ".", terrazas$coordenada_x_local))
terrazas$coordenada_y_local <- as.numeric(gsub(",", ".", terrazas$coordenada_y_local))
tmp <- terrazas[terrazas$coordenada_x_local > 1000, ]
tmp <- terrazas[terrazas$coordenada_y_local > 3e6,]

# UTM a siglo XXI
library(rgdal)
terrazas.utm     <- SpatialPoints(
    cbind(tmp$coordenada_x_local,
    tmp$coordenada_y_local),
    proj4string=CRS("+proj=utm +zone=30"))
terrazas.latlong <- spTransform(terrazas.utm,
    CRS("+proj=longlat"))

library(ggmap)
madrid <- get_map("madrid", zoom = 12)
tmp <- as.data.frame(terrazas.latlong)
ggmap(madrid) + geom_point(
    aes(x = coords.x1, y = coords.x2),
    data = tmp, size = .5,
    col = "red", alpha = 0.3)

Sobre las cursivas de todas:

  • Había coordenadas que apuntaban al golfo de Guinea. Ha habido que filtrarlas.
  • Según escribo, veo por la ventana una terraza que no figura en la lista. Habrá más, seguro.

En definitiva, ni están todas las que son, ni son todas las que están, ni son todas las que deberían serlo y las hay que no deberían ser pero que son y están. Hay que tener en cuenta de que para el ayuntamiento hay dos realidades:

  • La administrativa (licencias y autorizaciones)
  • La que se recoge a puros efectos estadísticos

Y luego, para la gente que tiene ojos en la cara, hay una tercera realidad:

  • La que es

Pero que no tiene necesariamente que coincidir con ninguna de las anteriores.