Extrair coordenadas do polígono do Brasil - maptools
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/ <http://augustoribas.heliohost.org> Github: https://github.com/Squiercg Lattes: http://lattes.cnpq.br/7355685961127056
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/ <http://augustoribas.heliohost.org> Github: https://github.com/Squiercg Lattes: http://lattes.cnpq.br/7355685961127056
-- Grato Augusto C. A. Ribas Site Pessoal: http://recologia.com.br/ <http://augustoribas.heliohost.org> Github: https://github.com/Squiercg Lattes: http://lattes.cnpq.br/7355685961127056
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
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 <comunello.eder@gmail.com>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.
Augusto, nesse mapa o Brasil esta definido por 62 poligonos plot(wrld_simpl[21, ]) length(wrld_simpl@polygons[[21]]@Polygons) o primeiro sendo o principal (maior area) which.max(sapply(wrld_simpl@polygons[[21]]@Polygons, function(x) x@area)) voce deseja as coordenadas de todos os poligonos? Para pegar as coordenadas do primeiro use coo1 <- wrld_simpl@polygons[[21]]@Polygons[[1]]@coords Elias.
participantes (4)
-
Augusto Ribas -
Elias T. Krainski -
Paulo Justiniano -
Éder Comunello