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ígonoplot(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!AttYury DuarteEngenheiro 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###############################################################