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 ###############################################################