Cargue inicial de la base de datos Comunas
comunas = shapefile("F:/MCD/Semestre 3/Análisis de información geográfica y espacial/Casos/cali/Comunas.shp")
comunas@data

Cargue inicial de la base de datos Encuesta origen y destino
require(readxl)
## Loading required package: readxl
Encuesta_ori_des= read_excel("F:/MCD/Semestre 3/Análisis de información geográfica y espacial/Casos/EncuestaOrigenDestino.xlsx", 
    sheet = "Sheet4")
## New names:
## • `MUNICIPIO` -> `MUNICIPIO...7`
## • `DEPARTAMENTO/LOCALIDAD/COMUNA/DISTRITO/BARRIO/VEREDA/HITO/DIRECCIÓN` ->
##   `DEPARTAMENTO/LOCALIDAD/COMUNA/DISTRITO/BARRIO/VEREDA/HITO/DIRECCIÓN...8`
## • `MUNICIPIO` -> `MUNICIPIO...11`
## • `DEPARTAMENTO/LOCALIDAD/COMUNA/DISTRITO/BARRIO/VEREDA/HITO/DIRECCIÓN` ->
##   `DEPARTAMENTO/LOCALIDAD/COMUNA/DISTRITO/BARRIO/VEREDA/HITO/DIRECCIÓN...12`
## • `OTRO¿CUÁL?` -> `OTRO¿CUÁL?...17`
## • `OTRO¿CUÁL?` -> `OTRO¿CUÁL?...24`
Encuesta_ori_des

Conversión de los datos de la columna a numericos
com=as.numeric(Encuesta_ori_des$comunaorigen)
## Warning: NAs introducidos por coerción

Datos Comuna Orígen


Filtro de la comuna 1 a la 22
posicion=which(com>=1&com<=22)

Conteo de cantidad por comuna origen y presentación en DataFrame
origenes=Encuesta_ori_des[posicion,]
head(table(origenes$comunaorigen))
## 
##   01   02   03   04   05   06 
##  794 2968 2121 1493  664  998
resultado=data.frame(table(origenes$comunaorigen))
head(resultado)

Se organiza las comunas de igual en las dos bases
as.character(resultado$Var1)
##  [1] "01" "02" "03" "04" "05" "06" "07" "08" "09" "10" "11" "12" "13" "14" "15"
## [16] "16" "17" "18" "19" "20" "21" "22"
orden=comunas$comuna
result_orden=resultado[orden,]
result_orden

Se crea la columna total origenes con sus respectivos valores
comunas@data$total_origenes=result_orden$Freq
comunas@data

Gráfico de los Datos origenes
lbls =paste(as.character(comunas$comuna),"-",comunas$total_origenes)
spl = list('sp.text', coordinates(comunas), lbls, cex=.5)
spl[[2]][7,]=spl[[2]][7,]+0.005

g1=spplot(comunas[,5], col.regions= terrain.colors(18), sp.layout = spl)
g1


Filtro de la comuna 1 a la 22 por Bicicleta
posicion=which(com>=1&com<=22&Encuesta_ori_des$TIPODEVEHÍCULO=="1")

Conteo de cantidad por comuna origen y presentación en DataFrame
origenes=Encuesta_ori_des[posicion,]
table(origenes$comunaorigen)
## 
##  01  02  03  04  05  06  07  08  09  10  11  12  13  14  15  16  17  18  19  20 
##  36 159  99  74  38  44  35  52  53  82  52  14  76  37  70  73 141 111 141  46 
##  21  22 
##  44  55
resultado=data.frame(table(origenes$comunaorigen))

Se organiza las comunas de igual en las dos bases
as.character(resultado$Var1)
##  [1] "01" "02" "03" "04" "05" "06" "07" "08" "09" "10" "11" "12" "13" "14" "15"
## [16] "16" "17" "18" "19" "20" "21" "22"
orden=comunas$comuna
result_orden=resultado[orden,]
result_orden

Se crea la columna Total origenes Bicicleta con sus respectivos valores
comunas@data$total_origen_bici=result_orden$Freq
comunas@data

Gráfico del Orígen de las Bicicletas
lbls =paste(as.character(comunas$comuna),"-",comunas$total_origen_bici)
spl = list('sp.text', coordinates(comunas), lbls, cex=.5)
spl[[2]][7,]=spl[[2]][7,]+0.005

g2=spplot(comunas[,6], col.regions= heat.colors(20,.95,.4), sp.layout = spl)
g2


Filtro de la comuna 1 a la 22 por Moto
posicion=which(com>=1&com<=22&Encuesta_ori_des$TIPODEVEHÍCULO=="2")

Conteo de cantidad por comuna origen y presentación en DataFrame
origenes=Encuesta_ori_des[posicion,]
table(origenes$comunaorigen)
## 
##   01   02   03   04   05   06   07   08   09   10   11   12   13   14   15   16 
##  375 1414  970  720  319  473  350  552  413  675  431  176  554  310  536  450 
##   17   18   19   20   21   22 
## 1136  746 1280  339  378  515
resultado=data.frame(table(origenes$comunaorigen))      

Se organiza las comunas de igual en las dos bases
as.character(resultado$Var1)
##  [1] "01" "02" "03" "04" "05" "06" "07" "08" "09" "10" "11" "12" "13" "14" "15"
## [16] "16" "17" "18" "19" "20" "21" "22"
orden=comunas$comuna
result_orden=resultado[orden,]
result_orden

Se crea la columna Total origenes Moto con sus respectivos valores
comunas@data$total_origen_moto=result_orden$Freq
comunas@data

Gráfico del Orígen de las Motos
lbls =paste(as.character(comunas$comuna),"-",comunas$total_origen_moto)
spl = list('sp.text', coordinates(comunas), lbls, cex=.5)
spl[[2]][7,]=spl[[2]][7,]+0.005

g3=spplot(comunas[,7], col.regions= terrain.colors(18), sp.layout = spl)
g3


Filtro de la comuna 1 a la 22 por Carro
posicion=which(com>=1&com<=22&Encuesta_ori_des$TIPODEVEHÍCULO=="3")

Conteo de cantidad por comuna origen y presentación en DataFrame
origenes=Encuesta_ori_des[posicion,]
table(origenes$comunaorigen)
## 
##   01   02   03   04   05   06   07   08   09   10   11   12   13   14   15   16 
##  301 1121  854  551  241  389  283  420  382  529  361  123  486  272  456  385 
##   17   18   19   20   21   22 
##  899  559 1165  272  328  489
resultado=data.frame(table(origenes$comunaorigen))      

Se organiza las comunas de igual en las dos bases
as.character(resultado$Var1)
##  [1] "01" "02" "03" "04" "05" "06" "07" "08" "09" "10" "11" "12" "13" "14" "15"
## [16] "16" "17" "18" "19" "20" "21" "22"
orden=comunas$comuna
result_orden=resultado[orden,]
result_orden

Se crea la columna Total origenes Carro con sus respectivos valores
comunas@data$total_origen_carro=result_orden$Freq
comunas@data

Gráfico del Orígen de los Carros
lbls =paste(as.character(comunas$comuna),"-",comunas$total_origen_carro)
spl = list('sp.text', coordinates(comunas), lbls, cex=.5)
spl[[2]][7,]=spl[[2]][7,]+0.005

g4=spplot(comunas[,8], col.regions= heat.colors(20,.95,.4), sp.layout = spl)
g4


Datos Comuna Destino


Conversión de los datos de la columna a numericos
com=as.numeric(Encuesta_ori_des$comunadestino)
## Warning: NAs introducidos por coerción

Filtro de la comuna 1 a la 22
posicion=which(com>=1&com<=22)

Conteo de cantidad por comuna destino y presentación en DataFrame
destino=Encuesta_ori_des[posicion,]
table(destino$comunadestino)
## 
##   01   02   03   04   05   06   07   08   09   10   11   12   13   14   15   16 
##  230 4810 3859 1856  522  716  748 1025 1279  927  564  246  673  379  579  743 
##   17   18   19   20   21   22 
## 2164  662 3158  351  515 1860
resultado=data.frame(table(destino$comunadestino))

Se organiza las comunas de igual en las dos bases
as.character(resultado$Var1)
##  [1] "01" "02" "03" "04" "05" "06" "07" "08" "09" "10" "11" "12" "13" "14" "15"
## [16] "16" "17" "18" "19" "20" "21" "22"
orden=comunas$comuna
result_orden=resultado[orden,]
result_orden

Se crea la columna Total Destino con sus respectivos valores
comunas@data$total_destino=result_orden$Freq
comunas@data

Gráfico de los Datos destino
lbls =paste(as.character(comunas$comuna),"-",comunas$total_destino)
spl = list('sp.text', coordinates(comunas), lbls, cex=.5)
spl[[2]][7,]=spl[[2]][7,]+0.005

g5=spplot(comunas[,9], col.regions= terrain.colors(18), sp.layout = spl)
g5


Filtro de la comuna 1 a la 22 por Bicicletas
posicion=which(com>=1&com<=22&Encuesta_ori_des$TIPODEVEHÍCULO=="1")

Conteo de cantidad por comuna destino y presentación en DataFrame
destino=Encuesta_ori_des[posicion,]
table(destino$comunadestino)
## 
##  01  02  03  04  05  06  07  08  09  10  11  12  13  14  15  16  17  18  19  20 
##  12 218 189  85  41  34  41  54  65  43  29  15  33  22  20  35 115  31 147  25 
##  21  22 
##  26 100
resultado=data.frame(table(destino$comunadestino))

Se organiza las comunas de igual en las dos bases
as.character(resultado$Var1)
##  [1] "01" "02" "03" "04" "05" "06" "07" "08" "09" "10" "11" "12" "13" "14" "15"
## [16] "16" "17" "18" "19" "20" "21" "22"
orden=comunas$comuna
result_orden=resultado[orden,]
result_orden

Se crea la columna Total destino Bicicleta con sus respectivos valores
comunas@data$total_destino_bici=result_orden$Freq
comunas@data

Gráfico del Destino de las Bicicletas
lbls =paste(as.character(comunas$comuna),"-",comunas$total_destino_bici)
spl = list('sp.text', coordinates(comunas), lbls, cex=.5)
spl[[2]][7,]=spl[[2]][7,]+0.005

g6=spplot(comunas[,10], col.regions= heat.colors(20,.95,.4), sp.layout = spl)
g6


Filtro de la comuna 1 a la 22 por Motos
posicion=which(com>=1&com<=22&Encuesta_ori_des$TIPODEVEHÍCULO=="2")

Conteo de cantidad por comuna destino y presentación en DataFrame
destino=Encuesta_ori_des[posicion,]
table(destino$comunadestino)
## 
##   01   02   03   04   05   06   07   08   09   10   11   12   13   14   15   16 
##  112 2177 1734  846  215  323  353  470  608  419  254  128  346  188  271  327 
##   17   18   19   20   21   22 
##  927  302 1421  151  243  911
resultado=data.frame(table(destino$comunadestino))

Se organiza las comunas de igual en las dos bases
as.character(resultado$Var1)
##  [1] "01" "02" "03" "04" "05" "06" "07" "08" "09" "10" "11" "12" "13" "14" "15"
## [16] "16" "17" "18" "19" "20" "21" "22"
orden=comunas$comuna
result_orden=resultado[orden,]
result_orden

Se crea la columna Total destino Moto con sus respectivos valores
comunas@data$total_destino_moto=result_orden$Freq
comunas@data

Gráfico del Destino de las Motos
lbls =paste(as.character(comunas$comuna),"-",comunas$total_destino_moto)
spl = list('sp.text', coordinates(comunas), lbls, cex=.5)
spl[[2]][7,]=spl[[2]][7,]+0.005

g7=spplot(comunas[,11], col.regions= terrain.colors(18), sp.layout = spl)
g7


Filtro de la comuna 1 a la 22 por Carros
posicion=which(com>=1&com<=22&Encuesta_ori_des$TIPODEVEHÍCULO=="3")

Conteo de cantidad por comuna destino y presentación en DataFrame
destino=Encuesta_ori_des[posicion,]
table(destino$comunadestino)
## 
##   01   02   03   04   05   06   07   08   09   10   11   12   13   14   15   16 
##   86 2006 1557  757  225  275  286  403  489  384  239   85  226  138  245  321 
##   17   18   19   20   21   22 
##  959  276 1302  144  198  694
resultado=data.frame(table(destino$comunadestino))

Se organiza las comunas de igual en las dos bases
as.character(resultado$Var1)
##  [1] "01" "02" "03" "04" "05" "06" "07" "08" "09" "10" "11" "12" "13" "14" "15"
## [16] "16" "17" "18" "19" "20" "21" "22"
orden=comunas$comuna
result_orden=resultado[orden,]
result_orden

Se crea la columna Total destino Carros con sus respectivos valores
comunas@data$total_destino_carros=result_orden$Freq
comunas@data

Gráfico del Destino de los Carros
lbls =paste(as.character(comunas$comuna),"-",comunas$total_destino_carros)
spl = list('sp.text', coordinates(comunas), lbls, cex=.5)
spl[[2]][7,]=spl[[2]][7,]+0.005

g8=spplot(comunas[,12], col.regions= heat.colors(20,.95,.4), sp.layout = spl)
g8