### AULA 10 - GEOESTATISTICA

## AULA 10.1 - CRIACAO DE UM OBJETO 'geoda' A PARTIR DE UM SHAPE

#getwd()
#rm(list=ls())

### 1. INSTALACAO E CARREGAMENTO DO  PACOTE 'geoR'

# Instalacao do pacote geoR

install.packages("geoR",dep=T) # 'dep=T' - instalacao dos demais pctes necessarios

library(geoR)

# No 'geoR" - objetos da classe 'geodata'

### 2. INSTALA??O DE OUTROS PACOTES NECESS?RIOS

# instalacao de 'maptools', 'sf' e 'spdep'
# 'mapttols' - funcoes para imp/exp e manipulacao
# 'sf' - representacao de dados espaciais (exp/imp de shape, por exemplo)
# 'spdep' - analise de dados de area

# carregamento dos pacotes instalados

library(maptools) # pode tb usar a funcao 'require'
library(sf)
require(spdep)

### 3. ALGUMAS OBSERVA??ES SOBRE O 'geoR'

# bancos de dados (bd) disponiveis em cada pacote

data()

# listagem dos bd do pacote 'geoR'

data(package="geoR")

# carregar bd e listar seus objetos

data(s100)

# s100 - pontos com coord
class(sic.100) # objeto da classe 'geodata'

# funcao 'plot' - funcao generica - metodos espec para cd classe
# para objeto da classe 'geodata' - 'plot' roda funcao especifica

points(sic.100)

### 4. IMPORTA??O DE SHAPEs E CRIACAO DE OBJETO 'geodata'

### 4.1 IMPORTACAO DO SHAPE DOS SET CENS COM AS VARIAVEIS

### usando pacote 'sf'

setores <- st_read("./bd_aula_12/inc_den_var_pol.shp")
plot(st_geometry(setores))

# olhando o banco de setores

class(setores)
print(setores)

### 4.2 CRIACAO DE UM ARQ DE CENTROIDES DOS SET CENS

set.centr <- st_coordinates(st_centroid(setores))
class(set.centr)

### 4.3 IMPORT COORD DA BORDA DA AREA DE ESTUDO

borda.rp <- st_read("./bd_aula_12/perimetro_urbano_pol.shp")
class(borda.rp)
plot(st_geometry(borda.rp))

# extraindo coordenadas dos pontos da borda

head(borda.rp)
str(borda.rp)

borda.rp.geom <- st_geometry(borda.rp)
str(borda.rp.geom)
borda.coord <- st_coordinates(borda.rp.geom)
class(borda.coord)
head(borda.coord)
borda.coord[,c(1,2)]
plot(borda.coord[,1],borda.coord[,2])

### 4.4 CRIACAO DE UM OBJETO 'geodata'

deng <- list() # criacao de uma lista vazia (como abrir uma caixa!)

# colocacao dos centroides no obj 'deng'

deng$coords <- set.centr
head(deng$coords)

# colocacao dos valores da variavel depend (incid) no objeto 'deng'

head(setores)

# transformado data frame espacial em data frame n?o espacial

set.nao.esp <- as.data.frame(setores)
class(set.nao.esp)
head(set.nao.esp)
str(set.nao.esp)

deng$data <- set.nao.esp[,4] # informando q a var da col 4 eh a depend (incidencia)
head(deng$data) # objeto do tipo 'list' com os valores da incidencia

# colocacao dos valores das covariaveis no obj 'deng'

deng$covariates <- set.nao.esp[,-c(4,26)] # informando quem sao as covariaveis (estamos excluindo a var da col 4 - q eh a var depend e a coluna 26 q eh geometria)
head(deng) # todo o bd
head(deng$covariates) # so mostra as covariaveis
head(deng$data) # so a var depend (incid)
head(deng$coords) # so mostra as coordenadas
str(deng)

# colocando a borda no objeto 'deng'

deng$borders <- borda.coord[,c(1,2)] # coord do poligono de borda
deng$borders # mostra todos as coord de todos os pontos da borda
head(deng$borders) # mostra as coord dos seis primeiros pontos

# transf 'deng' em um objeto do tipo 'geodata'

class(deng) # eh objeto do tipo lista
class(deng) <- 'geodata' # para dizer q eh da classe 'geodata'
class(deng) # agora eh objeto 'geodata'

### 'olhando o obj 'deng'  com a funcao 'plot'

deng$data

plot(deng) # sempre que usar so o nome do bd - R faz operacao so para a var depend
names(deng) # componentes do bd

### 5. EXPORTACAO DO 'geodata' - FUNCAO 'dump'
# para leitura pelo proprio R em, por exemplo, outra workspace

dump("deng","./resol_aulas_12_1a4/deng.R") # exporta o geodata 'deng'