Bom dia Jônatan e colegas listeiros, como vão? Muito obrigado pela sua ajuda, Jônatan, acho que vão ser bastante úteis esses comandos que você mostrou. Delimitar o perímetro da minha área também é um objetivo desse trabalho, mas a intenção principal dele é de formar polígonos unindo cada ponto com seus pontos adjacentes (vizinhos), de forma que eu fique com N polígonos dentro da minha área e que cada polígono gerado carregue consigo as informações dos pontos utilizados para gerar o polígono. Não sei se utilizando a função buffer eu teria um resultado satisfatório. Também me foi indicado utilizar a função 'voronoi' para atingir esse objetivo, mas não consigo aplica-la de forma assertiva. Mais uma vez, agradeço pela ajuda e atenção de todos! Att Yury Duarte Engenheiro Agrônomo - ESALQ/USP Em qui, 23 de ago de 2018 às 00:53, Jônatan <jdtatsch@gmail.com> escreveu:
Supondo que o polígono de interesse seja a envoltória convexa dos pontos.
library(sp) library(rgdal) library(raster)
pts_shp <- shapefile("X/x.shp") plot(pts_shp) p <- coordinates(pts_shp)
ch <- chull(p) coords <- p[c(ch, ch[1]), ] # fechando o polígono
plot(p, pch = 19) lines(coords, col = "red", lwd = 3)
sp_pol <- SpatialPolygons( list(Polygons( list(Polygon(coords)), ID = 1 )) ) plot(sp_pol, axes = TRUE) points(pts_shp, pch = 20, cex = 0.2)
sp_pol_df <- SpatialPolygonsDataFrame( sp_pol, data = data.frame(ID = 1) ) #shapefile(sp_pol_df, "pol-from-pts.shp", overwrite = TRUE) #x <- shapefile("pol-from-pts.shp") #plot(x, axes = TRUE)
On Wed, Aug 22, 2018 at 2:23 PM Yury Duarte via R-br < r-br@listas.c3sl.ufpr.br> wrote:
Boa tarde colegas listeiros, como vão?
Estou começando a trabalhar com dados georreferenciados e geoestatistica no R e me deparei com uma dificuldade que não estou conseguindo contornar. Tenho um shapefile de pontos que eu gostaria de transformar em um shape de polígonos e não estou sabendo como fazer essa conversão. Depois de fazer essa transformação de pontos para polígonos, queria saber se também é possível eu unir os polígonos adjacentes quando seus valores forem iguais, afim de reduzir o numero total de polígonos do meu shapefile.
Desde já, agradeço pela ajuda e atenção de todos!
Att
Yury Duarte Engenheiro Agrônomo - 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.
-- ############################################################### ## Jônatan Dupont Tatsch ## Professor do Departamento de Física ## Centro de Ciências Exatas e Naturais (CCNE) ## Universidade Federal de Santa Maria - UFSM ## Faixa de Camobi, Prédio 13 - Campus UFSM - Santa Maria, RS, Brasil - 97105-900 ## Telefone: +55(55)33012083 ## www.ufsm.br/meteorologia ###############################################################