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@data)

IDNew   <- unlist(RJMun@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

_______________________________________________
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.