[R-br] Ler .MAP no R - Shapefile

[Éder Comunello]

Senhores, boa tarde!

Supondo que o mapa de regiões não estivesse disponível, ainda seria
possível 'derivá-lo' do mapa de munícipios fornecido pelo IBGE.

A rotina que segue une polígonos com uma característica (atributo) comum.
Nesse caso específico utilizei a informação referente às mesoregiões do RJ,
disponíveis no shape 1:500.000 com projeção Policônica.

### <code r>
# install.packages('gpclib', dep=T)
# library(gpclib)     # General Polygon Clipping library
library(maptools)   # for geospatial services; also loads foreign and sp
library(rgdal)      # for map projection work; also loads sp

setwd('c:/Lab/Datasus'); getwd() ### ALTERAR!!!

URL1    <- 'ftp://geoftp.ibge.gov.br/malhas_digitais'
URL2    <-
'municipio_2005/escala_500mil/proj_policonica/arcview_shp/uf/rj/33mu500pc.zip'
fName   <-  basename(URL2)
dirName <-  gsub("\\.zip", "", fName)

if(!file.exists(fName)) download.file(file.path(URL1, URL2), fName,
mode='wb')
unzip(fName, exdir=dirName)

RJMun <- readShapePoly(file.path(dirName, gsub("\\.zip", "\\.shp", fName)))
proj4string(RJMun) <- CRS("+init=epsg:29101") # policônica
proj4string(RJMun)
# print(CRSargs(CRS("+init=epsg:29101")))
plot(RJMun, col=terrain.colors(8))

names(RJMun em data)
IDNew   <- unlist(RJMun em data["MESOREGIAO"]) ### Novos IDs
RJMeso  <- unionSpatialPolygons(RJMun, IDNew)

plot(RJMeso, col=heat.colors(6))
names(RJMeso)
### </code>


================================================
Éder Comunello
Agronomist, MSc in Environ. Sciences (UEM)

PhD Student in Agricultural Systems Engineering (USP/ESALQ)
Piracicaba, SP, Brazil [22 42.7'S, 47 37.8'W]

Researcher at Embrapa Western Region Agriculture
Dourados, MS, Brazil [22 16.5'S, 54 49.0'W]
================================================
UTC-04:00 / DST: UTC-03:00
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20140415/9727bfb7/attachment.html>