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

[Yury Duarte]

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 em 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 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/b9d68457/attachment.html>