[R-br] Diagrama de voronoi

[Fabrício Tavares]

Wagner, funcionou perfeitamente! Muito obrigado!

Fabrício.
http://vizzual.net



Em 30 de maio de 2011 12:18, Wagner bonat <wbonat em gmail.com> escreveu:

> Veja este exemplo, fiz isto a algum tempo deve estar desatualizado
>
> require(gpclib)
> require(tripack)
> require(sp)
> require(spdep)
> ## Lendo o arquivo das bordas
> bordas <- read.table("
> http://www.leg.ufpr.br/~wagner/ExemploBordaVoronoi/points_higi.csv
> ",header=TRUE,sep=";",dec=".")
> bordas = bordas[,-1]
> bordas[44,] <- bordas[1,]
>
> ## Lendo o arquivo que contem as coordenadas
> mosq = read.table("
> http://www.leg.ufpr.br/~wagner/ExemploBordaVoronoi/mosqhigi.csv
> ",header=TRUE,sep=";",dec=".")
>
> ## Funcoes extras
> source("http://www.leg.ufpr.br/~wagner/ExemploBordaVoronoi/voronoi.R")
> source("http://www.leg.ufpr.br/~wagner/ExemploBordaVoronoi/grp2sp.r")
>
> coords.x <- c()
> coords.y <- c()
> cod.mosq <- unique(mosq$mosquitrap)
> for(i in 1:40){
> coords.x[i] = unique(mosq[which(mosq$mosquitrap == cod.mosq[i]),]$x)
> coords.y[i] = unique(mosq[which(mosq$mosquitrap == cod.mosq[i]),]$y)}
>
> mapa.mosq <- voronoi(coords.x1=coords.x,coords.x2=coords.y,poligono =
> bordas)
> mapa.mosq = grp2sp(mapa.mosq,ID = cod.mosq)
> mapa.mosq.nb <- poly2nb(mapa.mosq)
> plot(mapa.mosq)
>
> --
> Wagner Hugo Bonat
> LEG - Laboratório de Estatística e Geoinformação
> UFPR - Universidade Federal do Paraná
>
> _______________________________________________
> R-br mailing list
> R-br em listas.c3sl.ufpr.br
> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
>
>
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110530/644eef83/attachment.html>