
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}