TP Final - Introducción a la Ciencia de Datos y Ciencia de Datos Geográfica
Arq. Julieta Ercoli
Agosto 2020
Introducción
El presente trabajo pretende aplicar los conocimientos adquiridos en los módulos 2 y 3 del Posgrado: “Big Data e Inteligencia terririal”.
Para la aplicación de los mismos se abordará la temática de las obras iniciadas en la Ciudad Autónoma de Buenos Aires y se intentará explicar cómo se distribuyen las mismas y su relación con la densidad poblacional de las comunas en las que en mayor y menor medida se han iniciado obras.
Analizando los Datos
Los primeros datos que analizaremos son las obras civiles, públicas y privadas o remodelaciones registradas en el Ciudad Autónoma de Buenos Aires. (Fuente: data.buenosaires.gob.ar)
Obras iniciadas desde Mayo 2019 en la Ciudad Autónoma de Buenos Aires
library (tidyverse)## -- Attaching packages -------------------------- tidyverse 1.3.0 --## v ggplot2 3.3.0 v purrr 0.3.4
## v tibble 3.0.1 v dplyr 0.8.5
## v tidyr 1.1.0 v stringr 1.4.0
## v readr 1.3.1 v forcats 0.5.0## -- Conflicts ----------------------------- tidyverse_conflicts() --
## x dplyr::filter() masks stats::filter()
## x dplyr::lag() masks stats::lag()library (sf)## Linking to GEOS 3.8.0, GDAL 3.0.4, PROJ 6.3.1#llamo la base de datos de obras-inicidas
obras_iniciadas<-read.csv ("C:/Users/Luz/Desktop/BDeIT/TP Final/obras-iniciadas.csv")
#visualizo la base
names(obras_iniciadas)## [1] "periodo" "long"
## [3] "lat" "expediente"
## [5] "calle" "calle_nro"
## [7] "seccion" "manzana"
## [9] "parcela" "fecha_inicio_obra"
## [11] "estado" "ultimo_avo_solicitado"
## [13] "fecha_ultimo_avo" "estado_tramite"
## [15] "profesional" "metros_a_construir"
## [17] "fecha_estado_tramite" "tipo_obra"
## [19] "destino" "barrio"
## [21] "comuna" "codigo_postal"
## [23] "codigo_postal_argentino"head(obras_iniciadas)## periodo long lat expediente calle
## 1 202002 -58.41157 -34.60654 SARMIENTO
## 2 202002 NA NA 25871534/2016 SENDERO PEATONAL SIN NOMBRE OFICIAL
## 3 202002 -58.44741 -34.56938 13764166/2016 AGUILAR
## 4 202002 -58.48207 -34.56828 LUGONES
## 5 202002 -58.39544 -34.58798 11581348/2017 AZCUENAGA
## 6 202002 -58.49779 -34.63820 LEOPARDI
## calle_nro seccion manzana parcela fecha_inicio_obra estado
## 1 3208 13 057 001A 2019-10-12 Normal
## 2 4363 66 101 001 2017-03-31 Normal
## 3 2434 37 151 003a 2017-01-20 Normal
## 4 2748 53 052 026B Normal
## 5 1872 11 024 001H 2017-07-19 Clausurada
## 6 74 60 138 016c Normal
## ultimo_avo_solicitado fecha_ultimo_avo
## 1 Inicio de obra 2019-12-16
## 2 En AVO 4 2018-03-08
## 3 En AVO 4 2019-04-24
## 4 Demolición final 2019-11-29
## 5 Inicio de obra
## 6 Demolición
## estado_tramite profesional
## 1 Inicio de obra LALA NISIM RICARDO
## 2 Obra finalizada HATSERIAN FEDERICO
## 3 Obra finalizada PASQUARELLI PAOLA EMILIA
## 4 Pedido de excavación menor a 4 mts autorizado MACHIOTE IGNACIO LUIS
## 5 Inicio de obra MERLIN LEONARDO JOSE
## 6 Pedido de demolición autorizado CANO MONICA BEATRIZ
## metros_a_construir fecha_estado_tramite
## 1 973.44 2019-12-18
## 2 2400.00 2018-03-12
## 3 1479.60 2018-12-07
## 4 1301.00 2020-01-23
## 5 10.87 2018-11-05
## 6 86.00 2020-01-21
## tipo_obra destino
## 1 AMPLIACION Y MODIFICACION C/DEMOLICION PARCIAL OTROS
## 2 OBRA NUEVA VIVIENDA MULTIFAMILIAR
## 3 OBRA NUEVA VIVIENDA MULTIFAMILIAR
## 4 DEMOLICIÓN TOTAL Y OBRA NUEVA VIVIENDA MULTIFAMILIAR
## 5 MODIFICACION Y/O AMPLIACION OTROS
## 6 DEMOLICIÓN TOTAL VIVIENDA MULTIFAMILIAR
## barrio comuna codigo_postal codigo_postal_argentino
## 1 Comuna 3 1196 Balvanera C1196AAJ
## 2 NA
## 3 Comuna 13 1426 Colegiales C1426DSJ
## 4 Comuna 12 1430 Villa Urquiza C1430ECD
## 5 Comuna 2 1128 Recoleta C1128AAF
## 6 Comuna 10 1407 Villa Luro C1407IDB#ordeno la base para observar la cantidad de obras iniciadas por comuna
agrup_obras_comuna <- obras_iniciadas %>%
filter(!is.na(comuna)) %>%
group_by(barrio) %>%
summarise (cantidad = n())
head (agrup_obras_comuna)## # A tibble: 6 x 2
## barrio cantidad
## <chr> <int>
## 1 Comuna 1 282
## 2 Comuna 10 554
## 3 Comuna 11 821
## 4 Comuna 12 1035
## 5 Comuna 13 706
## 6 Comuna 14 650agrup_obras_comuna<-rename(agrup_obras_comuna, COMUNA = barrio)
head (agrup_obras_comuna)## # A tibble: 6 x 2
## COMUNA cantidad
## <chr> <int>
## 1 Comuna 1 282
## 2 Comuna 10 554
## 3 Comuna 11 821
## 4 Comuna 12 1035
## 5 Comuna 13 706
## 6 Comuna 14 650ggplot(agrup_obras_comuna) + geom_col(aes(x=cantidad, y=COMUNA), fill= "#FF6666") + labs (title = "Obras Iniciadas por Comuna", subtitle = "Ciudad Autónoma de Buenos Aires", x="Obras Iniciadas", y="Comunas", caption = "Fuente: data.buenosaires.gob.ar")
Del primer análisis realizado podemos observar que la mayor cantidad de obras iniciadas en la Ciudad Autónoma de Buenos Aires han sido en la comuna 12 y la menor en la comuna 8, ahora veremos qué tipo de obras son las que se iniciaron en la comuna 12 y qué tipo de obras son las que se iniciaron en la comuna 8.
Comuna_12 <- filter(obras_iniciadas, barrio == "Comuna 12")
Tipos_en_12<- Comuna_12 %>%
filter (!is.na(tipo_obra)) %>%
group_by(tipo_obra) %>%
summarise (cantidad = n())
head (Tipos_en_12)## # A tibble: 6 x 2
## tipo_obra cantidad
## <chr> <int>
## 1 "" 14
## 2 " " 1
## 3 "AMPLIACIÓN Y MODIFICACIÓN CON DEMOLICIÓN PARCIAL" 8
## 4 "AMPLIACION Y MODIFICACION C/DEMOLICION PARCIAL" 142
## 5 "DEMOLICIÓN TOTAL" 1
## 6 "DEMOLICIÓN TOTAL Y OBRA NUEVA" 21ggplot(Tipos_en_12) + geom_col(aes(x=cantidad, y=tipo_obra), fill= "#C77CFF") + labs (title = "Tipo de Obras Iniciadas en Comuna 12", subtitle = "Ciudad Autónoma de Buenos Aires", y="Tipo de Obras Iniciadas", x="Cantidad", caption = "Fuente: data.buenosaires.gob.ar")
Comuna_8 <- filter(obras_iniciadas, barrio == "Comuna 8")
Tipos_en_8<- Comuna_8 %>%
filter (!is.na(tipo_obra)) %>%
group_by(tipo_obra) %>%
summarise (cantidad = n())
head (Tipos_en_8)## # A tibble: 6 x 2
## tipo_obra cantidad
## <chr> <int>
## 1 "" 3
## 2 "AMPLIACIÓN Y MODIFICACIÓN CON DEMOLICIÓN PARCIAL" 4
## 3 "AMPLIACION Y MODIFICACION C/DEMOLICION PARCIAL" 31
## 4 "DEMOLICIÓN TOTAL" 1
## 5 "DEMOLICIÓN TOTAL Y OBRA NUEVA" 1
## 6 "DEMOLICION PARCIAL" 2ggplot(Tipos_en_8) + geom_col(aes(x=cantidad, y=tipo_obra), fill= "#00BA42") + labs (title = "Tipo de Obras Iniciadas en Comuna 8", subtitle = "Ciudad Autónoma de Buenos Aires", y="Tipo de Obras Iniciadas", x="Cantidad", caption = "Fuente: data.buenosaires.gob.ar")
Del segundo análisis observamos que, tanto en la comuna con mayor inicio de obras de la Ciudad (Comuna 12) y en la que posee menor inicio de obras (Comuna 8), el tipo de obra que más permisos solicitó es Demolición Total y Obra Nueva.
Visualizaciones territoriales
Observaremos a continuación la distribución de la totalidad de las obras inicias en la Ciudad de Buenos Aires por comuna:
comunas <- st_read('https://bitsandbricks.github.io/data/CABA_comunas.geojson')## Reading layer `CABA_comunas' from data source `https://bitsandbricks.github.io/data/CABA_comunas.geojson' using driver `GeoJSON'
## Simple feature collection with 15 features and 4 fields
## geometry type: MULTIPOLYGON
## dimension: XY
## bbox: xmin: -58.53152 ymin: -34.70529 xmax: -58.33514 ymax: -34.52754
## geographic CRS: WGS 84Obras_Inicio_Espacial <- obras_iniciadas %>%
filter(!is.na(lat), !is.na(long)) %>%
st_as_sf(coords = c("long", "lat"), crs = 4326)
ggplot() +
geom_sf(data = comunas) +
geom_sf(data = Obras_Inicio_Espacial, alpha = .2, color = "blue")
Obras_con_comunas <- st_join(Obras_Inicio_Espacial, comunas)## although coordinates are longitude/latitude, st_intersects assumes that they are planar
## although coordinates are longitude/latitude, st_intersects assumes that they are planarhead (Obras_con_comunas)## Simple feature collection with 6 features and 25 fields
## geometry type: POINT
## dimension: XY
## bbox: xmin: -58.49779 ymin: -34.6382 xmax: -58.39544 ymax: -34.56828
## geographic CRS: WGS 84
## periodo expediente calle calle_nro seccion manzana
## 1 202002 SARMIENTO 3208 13 057
## 2 202002 13764166/2016 AGUILAR 2434 37 151
## 3 202002 LUGONES 2748 53 052
## 4 202002 11581348/2017 AZCUENAGA 1872 11 024
## 5 202002 LEOPARDI 74 60 138
## 6 202002 22831554/2016 FRAY JUSTO SANTAMARIA DE ORO 2142 19 124
## parcela fecha_inicio_obra estado ultimo_avo_solicitado fecha_ultimo_avo
## 1 001A 2019-10-12 Normal Inicio de obra 2019-12-16
## 2 003a 2017-01-20 Normal En AVO 4 2019-04-24
## 3 026B Normal Demolición final 2019-11-29
## 4 001H 2017-07-19 Clausurada Inicio de obra
## 5 016c Normal Demolición
## 6 012a Normal En AVO 1 2020-08-01
## estado_tramite profesional
## 1 Inicio de obra LALA NISIM RICARDO
## 2 Obra finalizada PASQUARELLI PAOLA EMILIA
## 3 Pedido de excavación menor a 4 mts autorizado MACHIOTE IGNACIO LUIS
## 4 Inicio de obra MERLIN LEONARDO JOSE
## 5 Pedido de demolición autorizado CANO MONICA BEATRIZ
## 6 En AVO 1 CARPINACCI HORACIO EDUARDO
## metros_a_construir fecha_estado_tramite
## 1 973.44 2019-12-18
## 2 1479.60 2018-12-07
## 3 1301.00 2020-01-23
## 4 10.87 2018-11-05
## 5 86.00 2020-01-21
## 6 39372.85 2019-05-13
## tipo_obra destino
## 1 AMPLIACION Y MODIFICACION C/DEMOLICION PARCIAL OTROS
## 2 OBRA NUEVA VIVIENDA MULTIFAMILIAR
## 3 DEMOLICIÓN TOTAL Y OBRA NUEVA VIVIENDA MULTIFAMILIAR
## 4 MODIFICACION Y/O AMPLIACION OTROS
## 5 DEMOLICIÓN TOTAL VIVIENDA MULTIFAMILIAR
## 6 OBRA NUEVA VIVIENDA MULTIFAMILIAR
## barrio comuna codigo_postal codigo_postal_argentino
## 1 Comuna 3 1196 Balvanera C1196AAJ
## 2 Comuna 13 1426 Colegiales C1426DSJ
## 3 Comuna 12 1430 Villa Urquiza C1430ECD
## 4 Comuna 2 1128 Recoleta C1128AAF
## 5 Comuna 10 1407 Villa Luro C1407IDB
## 6 Comuna 14 1425 Palermo C1425FOD
## barrios
## 1 BALVANERA - SAN CRISTOBAL
## 2 BELGRANO - COLEGIALES - NUÑEZ
## 3 COGHLAN - SAAVEDRA - VILLA PUEYRREDON - VILLA URQUIZA
## 4 RECOLETA
## 5 FLORESTA - MONTE CASTRO - VELEZ SARSFIELD - VERSALLES - VILLA LURO - VILLA REAL
## 6 PALERMO
## perimetro area comunas geometry
## 1 10486.26 6385991 3 POINT (-58.41157 -34.60654)
## 2 26967.43 14571813 13 POINT (-58.44741 -34.56938)
## 3 17232.19 15570927 12 POINT (-58.48207 -34.56828)
## 4 21246.61 6140873 2 POINT (-58.39544 -34.58798)
## 5 18332.04 12656557 10 POINT (-58.49779 -34.6382)
## 6 21768.07 15772496 14 POINT (-58.42783 -34.58304)Obras_con_comunas %>%
group_by(comunas) %>%
summarise(cantidad = n())## Simple feature collection with 15 features and 2 fields
## geometry type: MULTIPOINT
## dimension: XY
## bbox: xmin: -58.52997 ymin: -34.69467 xmax: -58.34297 ymax: -34.53621
## geographic CRS: WGS 84
## # A tibble: 15 x 3
## comunas cantidad geometry
## <chr> <int> <MULTIPOINT [°]>
## 1 1 297 ((-58.39214 -34.60192), (-58.39173 -34.61106), (-58.39152 -~
## 2 10 556 ((-58.52997 -34.61767), (-58.52985 -34.62577), (-58.52978 -~
## 3 11 822 ((-58.52671 -34.61126), (-58.52604 -34.61096), (-58.52595 -~
## 4 12 1038 ((-58.51308 -34.58195), (-58.51286 -34.58351), (-58.5116 -3~
## 5 13 707 ((-58.47397 -34.54105), (-58.47348 -34.54028), (-58.47344 -~
## 6 14 650 ((-58.44834 -34.56809), (-58.44727 -34.56864), (-58.44598 -~
## 7 15 676 ((-58.50262 -34.59429), (-58.50235 -34.5947), (-58.50174 -3~
## 8 2 162 ((-58.41463 -34.59476), (-58.41463 -34.59565), (-58.41421 -~
## 9 3 249 ((-58.41393 -34.61077), (-58.41337 -34.60612), (-58.41273 -~
## 10 4 330 ((-58.43147 -34.65311), (-58.43118 -34.64183), (-58.43069 -~
## 11 5 379 ((-58.43267 -34.60342), (-58.43167 -34.60161), (-58.43122 -~
## 12 6 490 ((-58.46116 -34.60727), (-58.46106 -34.61066), (-58.46049 -~
## 13 7 450 ((-58.47748 -34.62298), (-58.47722 -34.62321), (-58.47678 -~
## 14 8 133 ((-58.49517 -34.67645), (-58.49508 -34.67257), (-58.49372 -~
## 15 9 464 ((-58.52874 -34.64413), (-58.52867 -34.63944), (-58.52859 -~ggplot() +
geom_sf(data = comunas) +
geom_sf(data = Obras_con_comunas, aes(color = comunas)) +
labs(title = "Cantidad de obras iniciadas por Comuna",
subtitle = "Ciudad Autónoma de Buenos Aires",
fill = "Cantidad",
caption= "Fuente:data.buenosaires.gob.ar",
y="",
x="")
comunas_obras <- Obras_con_comunas %>%
filter(!is.na(comunas)) %>%
group_by(comunas) %>%
summarise(cantidad=n())
comunas_obras <- comunas_obras %>%
st_set_geometry(NULL)
head(comunas_obras)## # A tibble: 6 x 2
## comunas cantidad
## <chr> <int>
## 1 1 297
## 2 10 556
## 3 11 822
## 4 12 1038
## 5 13 707
## 6 14 650comunas <- comunas %>%
left_join(comunas_obras, by="comunas")
head(comunas)## Simple feature collection with 6 features and 5 fields
## geometry type: MULTIPOLYGON
## dimension: XY
## bbox: xmin: -58.4627 ymin: -34.6625 xmax: -58.33514 ymax: -34.56935
## geographic CRS: WGS 84
## barrios
## 1 CONSTITUCION - MONSERRAT - PUERTO MADERO - RETIRO - SAN NICOLAS - SAN TELMO
## 2 RECOLETA
## 3 BALVANERA - SAN CRISTOBAL
## 4 BARRACAS - BOCA - NUEVA POMPEYA - PARQUE PATRICIOS
## 5 ALMAGRO - BOEDO
## 6 CABALLITO
## perimetro area comunas cantidad geometry
## 1 35572.65 17802807 1 297 MULTIPOLYGON (((-58.36854 -...
## 2 21246.61 6140873 2 162 MULTIPOLYGON (((-58.39521 -...
## 3 10486.26 6385991 3 249 MULTIPOLYGON (((-58.41192 -...
## 4 36277.44 21701236 4 330 MULTIPOLYGON (((-58.3552 -3...
## 5 12323.47 6660526 5 379 MULTIPOLYGON (((-58.41287 -...
## 6 10990.96 6851029 6 490 MULTIPOLYGON (((-58.43061 -...ggplot() +
geom_sf(data = comunas, aes(fill=cantidad)) +
geom_sf_text(data=comunas, aes(label = comunas), size=2.5, colour = "black") +
labs(title = "Cantidad de obras iniciadas por Comuna",
subtitle = "Ciudad Autónoma de Buenos Aires",
fill = "Cantidad",
caption= "Fuente:data.buenosaires.gob.ar",
y="",
x="") +
scale_fill_gradient(low="lightgoldenrod", high = "steelblue4")
En el plano de la Ciudad Autónoma de Buenos Aires dividido por comunas podemos visualizar rápidamente las comunas que han solicitado más cantidad de permisos de obra y las que menos lo han hecho.
Densidad Poblacional
Poblacion_radios <- st_read("http://cdn.buenosaires.gob.ar/datosabiertos/datasets/informacion-censal-por-radio/CABA_rc.geojson")## Reading layer `CABA_rc' from data source `http://cdn.buenosaires.gob.ar/datosabiertos/datasets/informacion-censal-por-radio/CABA_rc.geojson' using driver `GeoJSON'
## Simple feature collection with 3554 features and 8 fields
## geometry type: MULTIPOLYGON
## dimension: XY
## bbox: xmin: -58.53092 ymin: -34.70574 xmax: -58.33455 ymax: -34.528
## geographic CRS: WGS 84names (Poblacion_radios) ## [1] "RADIO_ID" "BARRIO" "COMUNA" "POBLACION" "VIVIENDAS"
## [6] "HOGARES" "HOGARES_NBI" "AREA_KM2" "geometry"# Densidad poblacional de cada radio censal
radios <- Poblacion_radios %>%
mutate(densidad=POBLACION/AREA_KM2)
#Densidad poblacional por radios censales
plot(ggplot() +
geom_sf(data = radios, aes(fill = densidad), color = NA) +
scale_fill_gradient(low="lightgoldenrod", high = "steelblue4") +
labs(title = "Densidad de población",
subtitle = "Ciudad Autónoma de Buenos Aires",
fill = "hab/km2")+
theme_void())
Volcando en el plano la combinación de densidad de población y obras iniciadas en Ciudad Autónoma de Buenos Aires
names (radios)## [1] "RADIO_ID" "BARRIO" "COMUNA" "POBLACION" "VIVIENDAS"
## [6] "HOGARES" "HOGARES_NBI" "AREA_KM2" "geometry" "densidad"comunas_geo <- radios %>%
group_by(COMUNA) %>%
summarise(POBLACION = sum(POBLACION),
VIVIENDAS = sum(VIVIENDAS),
HOGARES = sum(HOGARES),
HOGARES_NBI = sum(HOGARES_NBI),
AREA_KM2 = sum(AREA_KM2))
names (comunas_geo)## [1] "COMUNA" "POBLACION" "VIVIENDAS" "HOGARES" "HOGARES_NBI"
## [6] "AREA_KM2" "geometry"ggplot () + geom_sf (data=comunas_geo)
# Combinamos la base de radios censales con las de cantidad de obras por comuna
comu_geo <- comunas_geo %>% left_join(agrup_obras_comuna)## Joining, by = "COMUNA"head (comu_geo)## Simple feature collection with 6 features and 7 fields
## geometry type: POLYGON
## dimension: XY
## bbox: xmin: -58.53092 ymin: -34.64626 xmax: -58.33925 ymax: -34.528
## geographic CRS: WGS 84
## # A tibble: 6 x 8
## COMUNA POBLACION VIVIENDAS HOGARES HOGARES_NBI AREA_KM2
## <chr> <dbl> <dbl> <dbl> <dbl> <dbl>
## 1 1 205886 131382 84468 13429 17.8
## 2 10 166022 71664 61453 2149 12.6
## 3 11 189832 84734 71460 1444 14.1
## 4 12 200116 93502 78547 1335 15.6
## 5 13 230767 129482 100257 1879 14.4
## 6 14 226534 141710 103167 2425 16.0
## # ... with 2 more variables: geometry <POLYGON [°]>, cantidad <int>ggplot () + geom_sf (data = comu_geo, aes(fill= cantidad/POBLACION)) +
scale_fill_gradient(low="lightgoldenrod", high = "steelblue4") +
labs (title= "Cantidad de obras inicias per cápita",
subtitle= "Comunas de Ciudad Autónoma de Buenos Aires",
fill= "obras/habitantes")
De este plano podemos concluir que la densidad de población no es una variable significativamente relevante en la cantidad de obras iniciadas en la Ciudad Autónoma de Buenos Aires.
Conclusiones
Luego del análisis realizado podemos concluir que en la Ciudad Autónoma de Buenos Aires, la solicitud de permiso de obra más usual es: Demolición Total y Obra Nueva. Tanto en la comuna donde más obras se han iniciado (Comuna 12) como en la que menos se han iniciado (Comuna 8) las solicitudes siguen siendo en mayor medida de permisos para Demolición Total y Obra Nueva. Además podemos concluir que la densidad poblacional no es una variable lo suficientemente considerable como para definir si tiene alguna relación con la cantidad de obras iniciadas.
. .