RStudio_SII_acquire_data.R

library(readr)
library(dplyr)
library(stringr)
library(sp)
library(rgdal)
library(rgeos)
library(spatialEco)

##%######################################################%##
#                                                          #
####                    NUM DEPTOS                      ####
#                                                          #
##%######################################################%##

sii.data <- read_delim("BRORGA2441N_00000_2014.txt", delim = "|", col_names = c("CodigoSIIdelaComuna",
"NumerodeManzana",
"NumerodePredial",
"Direcciononombredelpredio",
"Avaluofiscaltotal",
"Contribucionsemestral",
"Codigodedestinoprincipal",
"Avalúoexentodelapropiedad",
"CodigoSIIdelaComunaRolBienComun1",
"NumerodeManzanaRolBienComun1",
"NumerodepredioRolBienComun1",
"CodigoSIIdelaComunaRolBienComun2",
"NumerodeManzanaRolBienComun2",
"NumerodepredioRolBienComun2",
"Superficietotaldelterrenodelapropiedad"))

sii <- select(sii.data,"CodigoSIIdelaComuna","NumerodeManzana","NumerodePredial","Direcciononombredelpredio","Codigodedestinoprincipal","Superficietotaldelterrenodelapropiedad")
sii$CodigoSIIdelaComuna <- as.integer(sii$CodigoSIIdelaComuna)
sii$Superficietotaldelterrenodelapropiedad <- as.integer(sii$Superficietotaldelterrenodelapropiedad)
sii$CMN_MZ <- paste0(as.integer(sii$CodigoSIIdelaComuna),"-",as.integer(sii$NumerodeManzana))

deptos <- sii[grep("DEPTO|DPTO|DP", sii$Direcciononombredelpredio), ]

deptos.unique <- deptos
deptos.unique$buildings <- str_extract(deptos$Direcciononombredelpredio, "^(\\d+[A-Z\\s\\.]+|[A-Z\\s\\.]+)\\d+")
deptos.unique <- select(deptos.unique, CodigoSIIdelaComuna, NumerodeManzana, CMN_MZ, buildings) %>%
group_by(CodigoSIIdelaComuna,CMN_MZ,buildings) %>%
summarize(ndeptos=n())

##%######################################################%##
#                                                          #
####                    INTERSECTION                    ####
#                                                          #
##%######################################################%##


mz_valdivia.sii <- readOGR(dsn= "MANZANAS_SII/MZ_DR_VALDIVIA/", layer = "MZ_DR_VALDIVIA")
mz_valdivia.sii <- spTransform(mz_valdivia.sii,
CRS("+proj=utm +zone=19 +south +datum=WGS84 +units=m +no_defs +ellps=WGS84 +towgs84=0,0,0"))

mz_valdivia.censo <- readOGR(dsn="/Users/robsalasco/Downloads/Mapa Manzanas sin CODSII (geometría buena) - Cartografia Precenso 2011 - Público/Carto_Region_14.gdb/", layer="MANZANA")
mz_valdivia.censo <- spTransform(mz_valdivia.censo,
CRS("+proj=utm +zone=19 +south +datum=WGS84 +units=m +no_defs +ellps=WGS84 +towgs84=0,0,0"))

mz_valdivia.sii.centroid <- gCentroid(mz_valdivia.sii,byid = T)
mz_valdivia.sii.centroid <- SpatialPointsDataFrame(mz_valdivia.sii.centroid,
mz_valdivia.sii@data, match.ID=FALSE)

t <- point.in.poly(mz_valdivia.sii.centroid, mz_valdivia.censo)
t$MANZENT <-  as.character(t$MANZENT)

##%######################################################%##
#                                                          #
####                        JOIN                        ####
#                                                          #
##%######################################################%##


joined <- left_join(deptos.unique,t@data,by=c("CMN_MZ"="CMN_MZ"))
joined.final <-joined[joined$CodigoSIIdelaComuna==10101,]