### Plotting Simple Features

### 1 - Plot methods for 'sf' and 'sfc' objects

## 1.1 - Objetos somente com geometria: objeto 'sfc'

library(sf)
demo(nc,ask=F,echo=F) # lendo mapa 'nc' direto do pacote

# lendo o mapa 'nc'no diretório 'mapa_nc'

# nc <- st_read("./bancos_aula3/nc.shp")

plot(st_geometry(nc)) # 'st_geometry' - comando mostrar/criar um novo objeto 'sf' somente com geometria (sfc)

# inserção de cores, eixos, título e centroides
# 'axes': atrelado a CRS - para long/lat - graus e orientaçõo são adcioanados com 'axes = TRUE'

plot(st_geometry(nc),col=sf.colors(12,categorical=TRUE), border = 'grey', axes = TRUE,
     main="Map of North Caroline")
plot(st_geometry(st_centroid(nc)),pch=3, col='red', add=TRUE)

## 1.2 - Geometria com atributos: objeto 'sf'

# atributos de 'nc'

print(nc,n=3)
str(nc)

# plot default para um objeto 'sf': multi-plot com todos os atributos 
plot(nc) # plota 10, no máximo
plot(nc,max.plot=14) # plota todos os 14 atributos

# 'options' - controla o nmero de plotes

options(sf_max.plot=1)
plot(nc)

# escolha e plotagem de um único atributo - adição de legenda (color key)

plot(nc["AREA"])

## 1.3 Color key: posicionamento e tamanho

# posicionamento da color key: 1 = abaixo; 2 = esquerda; 3 = acima; 4 = direita

plot(nc["AREA"],key.pos = 4) # color key no lado direito

# transformando AREA em fatores e plotando

nc$f <- cut(nc$AREA, 10)
str(nc)

plot(nc["f"], axes = TRUE, key.pos = 4, pal = sf.colors(10))

## 1.3 Intervalos de classe

# color breaks (intervalos de classe): controlados pelos argumentos 'breaks' e 'nbreaks'
# 'breaks' - vetor com com valores das quebras
# 'nbreaks' - especifica o número de intervalos

plot(nc["AREA"],breaks=c(0,0.05,0.1,0.15,0.2,0.25))

# tipos de quebras

# default: 'pretty' - classe de quebras arredondadas - valores de quebras apromimados

plot(nc["AREA"], nbreaks=5)

# quantis

plot(nc["AREA"], nbreaks=5,breaks="quantile")

# intervalos iguais

plot(nc["AREA"], nbreaks=5,breaks="equal")