[R-br] Converter Shape de Pontos para Shape de Polígonos

Jônatan jdtatsch em gmail.com
Quinta Agosto 23 00:53:20 -03 2018


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 em 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 em 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
###############################################################
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20180823/343904e5/attachment.html>


Mais detalhes sobre a lista de discussão R-br