rm(list=ls())

Quick tmap

A função qtm() pode ser utilizada para composição de mapas de maneira expressa (quick). O código a seguir permite carregar o dado georgia, recriar o dado georgia_sf e gerar um mapa expresso.

data("georgia")
ls()

## [1] "georgia"       "georgia.polys" "georgia2"

Para criar o objeto sf georgia_sf, faça:

georgia_sf <- st_as_sf(georgia)

Finalmente, crie o mapa expresso utilizando:

library(tmap)
qtm(georgia, fill="red", style = "natural")

O uso do parâmetro style permite configurações rápidas de mapas por meio do package tmap. Explore a função qtm para entender os demais temas possíveis.

O parâmetro fill pode ser utilizado para especificar uma cor ou a variável a ser mapeada.

qtm(georgia_sf, fill = "MedInc", text ="Name", text.size=0.5, format = "World_wide", style = "classic", text.root = 5, fill.title = "Median Income")

Mapa completo por meio do tmap

Para gerar mapas completos por meio do package tmap, é necessário utilizar a função tm_shapeacompanhada de uma outra função que especifique o precisa ser plotado. Antes de plotar, um objeto contendo a linha da divisa estadual da Georgia pode ser criada por meio da função st_union() do package sf. Uma alternativa para dados sp é função gUnaryUnion() do package rgeos, carregado juntamente com o package GISTools.

# união das feições para composição do limite estadual
g <- st_union(georgia_sf)

# plotar as camadas espaciais
tm_shape(georgia_sf) + 
  tm_fill("tomato") +
  tm_borders(lty="dashed", col = "gold")

Incluir alguns estilos:

# plotar as camadas espaciais
tm_shape(georgia_sf) + 
  tm_fill("tomato") +
  tm_borders(lty="dashed", col = "gold") +
  tm_style("natural", bg.color = "grey90")

Incluir limite estadual. Repare que a função tm_shape vai ser chamada novamente uma vez que outro objeto será incluído no mapa.

tm_shape(georgia_sf) + 
  tm_fill("tomato") +
  tm_borders(lty="dashed", col = "gold") +
  tm_style("natural", bg.color = "grey90") +
  tm_shape(g) + 
  tm_borders(lwd=2) # lwd - line width - espessura da linha

Incluindo todos os parâmetros juntos:

tm_shape(georgia_sf) + 
  tm_fill("tomato") +
  tm_borders(lty="dashed", col = "gold") +
  tm_style("natural", bg.color = "grey90") +
  tm_shape(g) + 
  tm_borders(lwd=2) +
  tm_layout(title = "O estado da Georgia",
            title.size = 1,
            title.position = c(0.55, "top"))

Experimente esses outros parâmetros da função qtm.

data(World, rivers, metro, land)

# just the map
qtm(World)

# choropleth
qtm(World, fill = "economy", format = "World", style = "col_blind")

# choropleth with more specifications
qtm(World, fill="HPI", fill.n = 9, fill.palette = "div",
    fill.title = "Happy Planet Index", fill.id = "name", 
    style = "gray", format = "World")

# Alternative to previous map - RECOMENDED
tm_shape(World) +
    tm_polygons("HPI", n = 9, palette = "div",
        title = "Happy Planet Index", id = "name") +
    tm_style("gray") +
    tm_format("World")

# Bubble map
qtm(World, borders = NULL) + 
qtm(metro, symbols.size = "pop2010", 
    symbols.title.size= "Metropolitan Areas", 
    symbols.id= "name",
    format = "World")

Para plotar dois mapas juntos com dados diferentes, experimente o seguinte código:

library(tmap)

t1 <- tm_shape(georgia_sf) +
              tm_fill("coral") +
              tm_borders() +
              tm_layout(bg.color = "grey85")

t2 <- tm_shape(georgia2) +
              tm_fill("orange") +
              tm_borders() +
              tm_layout(asp=0.86, bg.color = "grey95") #alinhar a segunda figura
t1

t2

Para especificar o layout da figura combinando os mapas:

library(grid)

# abrir uma nova página de plotagem
grid.newpage()

# definir o layout
pushViewport(viewport(layout = grid.layout(1,2)))

#  Plotar usando o comando print
print(t1, vp=viewport(layout.pos.col = 1, height = 5))
print(t2, vp=viewport(layout.pos.col = 2, height = 5))

Repare que t1 é um objeto sf e t2 é um objeto sp e apresentam sistemas de referência diferentes.

Para consultar os diferentes sistemas de projeção, verifique utilizando a função crs.

crs(georgia_sf)

## CRS arguments: +proj=longlat +ellps=WGS84 +no_defs

crs(georgia2)

## CRS arguments:
##  +proj=aea +lat_1=29.5 +lat_2=45.5 +lat_0=23 +lon_0=-96 +x_0=0.0
## +y_0=0.0 +ellps=GRS80 +units=m +datum=NAD83 +no_defs +towgs84=0,0,0

Para nomear os atributos, verifique o data.frame do objeto georgia_sf. Utilize a função names() para ter acesso a uma lista de nomes de todos os atributos.

names(georgia_sf)

##  [1] "Latitude" "Longitud" "TotPop90" "PctRural" "PctBach"  "PctEld"  
##  [7] "PctFB"    "PctPov"   "PctBlack" "X"        "Y"        "ID"      
## [13] "Name"     "MedInc"   "geometry"

# Construir um data frame com os nomes. 
data.frame(georgia_sf)[,13]

##   [1] "Appling"       "Atkinson"      "Bacon"         "Baker"        
##   [5] "Baldwin"       "Banks"         "Barrow"        "Bartow"       
##   [9] "Ben Hill"      "Berrien"       "Bibb"          "Bleckley"     
##  [13] "Brantley"      "Brooks"        "Bryan"         "Bulloch"      
##  [17] "Burke"         "Butts"         "Calhoun"       "Camden"       
##  [21] "Candler"       "Carroll"       "Catoosa"       "Charlton"     
##  [25] "Chatham"       "Chattahoochee" "Chattooga"     "Cherokee"     
##  [29] "Clarke"        "Clay"          "Clayton"       "Clinch"       
##  [33] "Cobb"          "Coffee"        "Colquitt"      "Columbia"     
##  [37] "Cook"          "Coweta"        "Crawford"      "Crisp"        
##  [41] "Dade"          "Dawson"        "Decatur"       "DeKalb"       
##  [45] "Dodge"         "Dooly"         "Dougherty"     "Douglas"      
##  [49] "Early"         "Echols"        "Effingham"     "Elbert"       
##  [53] "Emanuel"       "Evans"         "Fannin"        "Fayette"      
##  [57] "Floyd"         "Forsyth"       "Franklin"      "Fulton"       
##  [61] "Gilmer"        "Glascock"      "Glynn"         "Gordon"       
##  [65] "Grady"         "Greene"        "Gwinnett"      "Habersham"    
##  [69] "Hall"          "Hancock"       "Haralson"      "Harris"       
##  [73] "Hart"          "Heard"         "Henry"         "Houston"      
##  [77] "Irwin"         "Jackson"       "Jasper"        "Jeff Davis"   
##  [81] "Jefferson"     "Jenkins"       "Johnson"       "Jones"        
##  [85] "Lamar"         "Lanier"        "Laurens"       "Lee"          
##  [89] "Liberty"       "Lincoln"       "Long"          "Lowndes"      
##  [93] "Lumpkin"       "McDuffie"      "McIntosh"      "Macon"        
##  [97] "Madison"       "Marion"        "Meriwether"    "Miller"       
## [101] "Mitchell"      "Monroe"        "Montgomery"    "Morgan"       
## [105] "Murray"        "Muscogee"      "Newton"        "Oconee"       
## [109] "Oglethorpe"    "Paulding"      "Peach"         "Pickens"      
## [113] "Pierce"        "Pike"          "Polk"          "Pulaski"      
## [117] "Putnam"        "Quitman"       "Rabun"         "Randolph"     
## [121] "Richmond"      "Rockdale"      "Schley"        "Screven"      
## [125] "Seminole"      "Spalding"      "Stephens"      "Stewart"      
## [129] "Sumter"        "Talbot"        "Taliaferro"    "Tattnall"     
## [133] "Taylor"        "Telfair"       "Terrell"       "Thomas"       
## [137] "Tift"          "Toombs"        "Towns"         "Treutlen"     
## [141] "Troup"         "Turner"        "Twiggs"        "Union"        
## [145] "Upson"         "Walker"        "Walton"        "Ware"         
## [149] "Warren"        "Washington"    "Wayne"         "Webster"      
## [153] "Wheeler"       "White"         "Whitfield"     "Wilcox"       
## [157] "Wilkes"        "Wilkinson"     "Worth"

Para incluir os nomes, faça:

tm_shape(georgia_sf) +
              tm_fill("white") +
              tm_borders() +
              tm_text("Name", size=0.3) +
              tm_layout(frame = FALSE)

Para criar um subconjunto dos dados, é possível fazer uma seleção criando uma camada com o subconjunto por meio dos índices:

index <- c(81, 82, 83, 150, 62, 53, 21, 16, 124, 121, 17)
georgia_sub <- georgia_sf[index, ]

Para sobrepor os mapas com as duas camadas, faça:

tm_shape(georgia_sf) +
              tm_fill("white") +
              tm_borders("grey", lwd=0.5) +
tm_shape(g) +
              tm_borders(lwd=2) +
tm_shape(georgia_sub) +
              tm_fill("lightblue") +
              tm_borders() +
              tm_layout(frame = TRUE, title = "Georgia com subconjunto", title.size=1, title.position = c(0.02, "bottom"))

---

Mapas coropléticos

---

Mapas de pontos

---

Representação de fluxos

---

Mapas de calor

---

Cartogramas

---

Símbolos proporcionais

---