Senhores, bom dia!

Aproveitando a discussão, encaminho uma sugestão de outro pacote para fins similares. A principal diferença é poder optar por resoluções diferentes e a possibilidade de excluir a representação por múltiplos polígonos utilizando ("less islands" ou "li")

### <code r>

require(rworldmap)

require(rworldxtra) # necessário para "high resolution"

world <- getMap(resolution = "li")

### resolution: "coarse","low","less islands","li","high" (require rworldxtra)

grep("Bra.il", world@data$ADMIN) #23

map <- world@polygons[[23]]@Polygons[[1]]@coords

map.sp <- SpatialPolygons(list(world@polygons[[23]]))

map.spdf <- world[23,]

parOri <- par(no.readonly=T)

par(mfrow=c(1,3))

plot(map, asp=1, type="l", col=2)

plot(map.sp, axes=T, bor=2)

plot(map.spdf, axes=T, bor=2)

par(parOri)

sapply(ls(patt="^map"), object.size)

### </ code>

Éder Comunello <c omunello.eder@gmail.com>
Dourados, MS - [22 16.5'S, 54 49'W]

Em 6 de março de 2015 22:12, Paulo Justiniano <paulojus@leg.ufpr.br> escreveu:

isto é pq certos municipios e estados podem ser definidos por mais de um poligono

Por isto p polygon que defini a unidade é formada por uma lista de Polygons, que é o desesho de cada (sub) área

Pense por exemplo em um estado com ilhas que pertençam a ele.

On Fri, 6 Mar 2015, Augusto Ribas wrote:

Consegui.

O que eu queria era o seguinte

brasil<-wrld_simpl[21,]
brasil@polygons[[1]]@Polygons[[1]]@coords
plot(brasil@polygons[[1]]@Polygons[[1]]@coords,type="l")

Acho que não estava entendendo o tal do Polygons dentro de polygons. Estrutura complicada.
Mas de qualquer forma agradeço a atenção. Bom fim de semana a todos.

Em 6 de março de 2015 17:11, Augusto Ribas <ribas.aca@gmail.com> escreveu:
Ola pessoal, estou com uma duvida aqui, como eu retiro as coordenadas que formam o polígono
do Brasil de um spatialpolygondataframe.

No pacote maptools tem o wrld_simpl que é um mapa do mundo, um dos polígonos é o mapa do Brasil,
como eu consigo somente os pontos, coordenadas na ordem para fazer o polígono, não estou
conseguindo entender como extrai essa informação

######
library(maptools)
data(wrld_simpl)
which(wrld_simpl$NAME=="Brazil")

plot(wrld_simpl[21,])
wrld_simpl[21,]
#######

Eu gostaria ter os pontos para escrever polygon(pontos) ou plot(pontos,type="l") e desenhar o
mapa em um plot igual acima, tem como fazer isso?

--
Grato
Augusto C. A. Ribas

Site Pessoal: http://recologia.com.br/
Github: https://github.com/Squiercg
Lattes: http://lattes.cnpq.br/7355685961127056

--
Grato
Augusto C. A. Ribas

Site Pessoal: http://recologia.com.br/
Github: https://github.com/Squiercg
Lattes: http://lattes.cnpq.br/7355685961127056

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