[R-br] Diagrama de voronoi

[Wagner bonat]

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