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

Sérgio Henrique almeida da silva ju sergio.edfisica em gmail.com
Terça Abril 15 17:35:03 BRT 2014


Olá

Obrigado a todos pela ajuda!!

Abraços


2014-04-15 17:22 GMT-03:00 Éder Comunello <comunello.eder em 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 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
>
> _______________________________________________
> R-br mailing list
> R-br em 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
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20140415/50d782d8/attachment.html>


Mais detalhes sobre a lista de discussão R-br