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 <[email protected]> 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 < > [email protected]> 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 >> [email protected] >> 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 > ############################################################### >
_______________________________________________ R-br mailing list [email protected] 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.
