Prezados e Prezadas,

 

Estou trabalhando com o IVS (Índice de Vulnerabilidade Social), mais precisamente com as UDHs (Unidades de Habitação) para o Município de Jaboatão dos Guararapes - PE. A base  Completa + Shapefiles estão no site do Ipea (http://ivs.ipea.gov.br/index.php/pt/biblioteca). Gostaria de fazer um mapa do Índice de Vulnerabilidade Social do Município para saber quais as localidades mais vulneráveis.

Tentei unir pelo código do município e pela variável  UDH_Atlas no shapefile e udh e travei.

Ao que parece:

1. UDHs são colocados como "strings" pelo R, ao colocar a função as.numeric(), há uma mudança nas linhas.

2. UDHs são coordenadas únicas (não se separam em Long. e Lat.), ou temos que separá-las em Long e Lat?

 

Agradeço desde já a ajuda,

 

Marcelo Nery Magalhães

 

Script:

#Bibliotecas

library(ggplot2)

library(rgdal)

 

#Shape

Jab<-readOGR("C:/Users/usuario/Documents/IVS", "UDHs_RM_Todas")

head(Jab@data)

 

#Transformando string em numérico

Jab$UDH_ATLAS<-as.numeric(Jab$UDH_ATLAS)

 

#Base de dados/tabela

ivs<-read.csv("C:/Users/usuario/Documents/IVS/jab.csv", sep = ";", head=T, dec = ",")

head(ivs)

 

ivs<-na.omit(ivs)

 

names(Jab@data)

names(ivs)

 

#Transformando os nomes das variáveis em um só para cada objeto

names(ivs)[12]<-c("Cod_mun_6")

names(Jab)[1]<-c("udh")  # tentei criar uma chave com essa variável, já que com o código do município não estava funcionando

names(Jab)

names(ivs)

 

#Ordenando

ivs<-ivs[order(ivs$Cod_mun_6),]

Jab<-Jab@data[order(Jab$Cod_mun_6),]

 

head(ivs)

head(Jab)

 

#Combinando através da Chave do Código do Município

j<-merge(Jab,ivs,by.y="Cod_mun_6") # Houve problema ao criar “j”

head(j)

 

dim(Jab)

dim(ivs)

dim(j)

 

fortify(Jab)

j.f<-fortify(Jab,region="udh")

 

j$ivs

#Categorizando

j$ivscat<-cut(j$ivs, breaks = c(0,0.200,0.300,0.400,0.500,1), labels=c("0 a 0,200 - Muito Baixa", "0,200 a 0,300 - Baixa", "0,300 a 0,400 - Média", "0,400 a 0,500 - Alta", "0,500 a 1 - Muito Alta"), include.lowest=T)

names(j)

names(j.f)

#Problema: udh não possui latitude e lobgitude (como fazer, então?)

jabo<-merge(j, j.f, by.y="udh")

 

 

_______________________________________________
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código 
mínimo reproduzível.

Responder a