Fabrício.
Em 30 de maio de 2011 12:18, Wagner bonat
<wbonat@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@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br