
Olá Obrigado a todos pela ajuda!! Abraços 2014-04-15 17:22 GMT-03:00 Éder Comunello <comunello.eder@gmail.com>:
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.
-- Sérgio Henrique Almeida da Silva Junior Doutorando em Epidemiologia em Saúde Pública Escola Nacional de Saúde Pública Sérgio Arouca - ENSP/FIOCRUZ http://lattes.cnpq.br/1611345552843383 Tel: (21) 968463637 http://www.linkedin.com/profile/view?id=250437145&trk=tab_pro