extrair bordas de arquivos .shp

Bom dia pessoal, Alguem poderia ajudar com a seguinta dúvida? Estou usando o pacote sp pra visualização e manipulação de dados espaciais no R. Um dos arquivos que estou trabalhando é do tipo .shp Qual procedimento para extrair as bordas dos polígonos deste arquivo que não seja tão lusitânica quanto o seguinte código: dados <- readShapePoly('dados.shp') length(dados@polygons) #supondo que o resultado deste comando seja == 2 length(dados@polygons[[1]]@Polygons) #e supondo que o resultado aqui seja == 10 length(dados@polygons[[2]]@Polygons) #e supondo que o resultado aqui seja == 10 Entao meu codigo para extrair as coordenadas ditas como bordas ou limites seria dados@polygons[[1]]@Polygons[[1]])@coords dados@polygons[[1]]@Polygons[[2]])@coords . . . dados@polygons[[1]]@Polygons[[10]])@coords e estes mesmos comandos para o poligono 2 o que totaliza 20 linhas de comando Talvez existisse algo como dados@borders mas não existe pois é possível verificar o que se pode extrair do objeto dados que é um objeto da classe "SpatialPolygonsDataFrame" pelo comando getSlots('SpatialPolygonsDataFrame') Obrigado Samuel ==================================== Samuel P. C. Carvalho Mestre em Ciências Florestais [UFLA] Doutorando em Recursos Florestais [ESALQ/USP] =============================================

Estou seguro de que a família apply() poderia ser usada para automatizar. Na falta de um CMR não tive como testar. À disposição. Walmes. ========================================================================== Walmes Marques Zeviani LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) Departamento de Estatística - Universidade Federal do Paraná fone: (+55) 41 3361 3573 VoIP: (3361 3600) 1053 1173 e-mail: walmes@ufpr.br skype: walmeszeviani twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ==========================================================================

Não sei muito bem o que você quer, mas vê se isso ajuda! plot(mapa) w <- owin(poly=locator()) #Desenhar o polígono da borda desejada, clicando no sentindo anti-horário coord.ppp <- as.ppp(coord.dt, W=w) Em 4 de novembro de 2012 13:53, Walmes Zeviani <walmeszeviani@gmail.com>escreveu:
Estou seguro de que a família apply() poderia ser usada para automatizar. Na falta de um CMR não tive como testar.
À disposição. Walmes.
========================================================================== Walmes Marques Zeviani LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) Departamento de Estatística - Universidade Federal do Paraná fone: (+55) 41 3361 3573 VoIP: (3361 3600) 1053 1173 e-mail: walmes@ufpr.br skype: walmeszeviani twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ==========================================================================
_______________________________________________ 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) 94429486/78101651 id: 123*20942

O pacote rgdal também é uma opção extrair os limites. ?readGDAL 2012/11/4 Samuel Carvalho <samukajm@yahoo.com.br>
Bom dia pessoal, Alguem poderia ajudar com a seguinta dúvida? Estou usando o pacote sp pra visualização e manipulação de dados espaciais no R. Um dos arquivos que estou trabalhando é do tipo .shp Qual procedimento para extrair as bordas dos polígonos deste arquivo que não seja tão lusitânica quanto o seguinte código:
dados <- readShapePoly('dados.shp') length(dados@polygons) #supondo que o resultado deste comando seja == 2 length(dados@polygons[[1]]@Polygons) #e supondo que o resultado aqui seja == 10 length(dados@polygons[[2]]@Polygons) #e supondo que o resultado aqui seja == 10
Entao meu codigo para extrair as coordenadas ditas como bordas ou limites seria dados@polygons[[1]]@Polygons[[1]])@coords dados@polygons[[1]]@Polygons[[2]])@coords . . . dados@polygons[[1]]@Polygons[[10]])@coords
e estes mesmos comandos para o poligono 2 o que totaliza 20 linhas de comando Talvez existisse algo como dados@borders mas não existe pois é possível verificar o que se pode extrair do objeto dados que é um objeto da classe "SpatialPolygonsDataFrame" pelo comando getSlots('SpatialPolygonsDataFrame')
Obrigado Samuel
*====================================* *Samuel P. C. Carvalho *Mestre em Ciências Florestais [UFLA] Doutorando em Recursos Florestais [ESALQ/USP] =============================================
_______________________________________________ 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.
-- "Small steps toward a much better world" \begin{signature} Daniel Marcelino Land Phone 1+514 343 6111 #3799 3200 Jean Brillant, Office C5071 Montreal, QC; H3T 1N8 Canada \end{signature}
participantes (4)
-
Daniel Marcelino
-
Samuel Carvalho
-
Sérgio Henrique almeida da silva ju
-
Walmes Zeviani