Veja este exemplo, fiz isto a algum tempo deve estar desatualizado<br><br>require(gpclib)<br>require(tripack)<br>require(sp)<br>require(spdep)<br>## Lendo o arquivo das bordas<br>bordas <- read.table("<a href="http://www.leg.ufpr.br/~wagner/ExemploBordaVoronoi/points_higi.csv">http://www.leg.ufpr.br/~wagner/ExemploBordaVoronoi/points_higi.csv</a>",header=TRUE,sep=";",dec=".")<br>
bordas = bordas[,-1]<br>bordas[44,] <- bordas[1,]<br><br>## Lendo o arquivo que contem as coordenadas<br>mosq = read.table("<a href="http://www.leg.ufpr.br/~wagner/ExemploBordaVoronoi/mosqhigi.csv">http://www.leg.ufpr.br/~wagner/ExemploBordaVoronoi/mosqhigi.csv</a>",header=TRUE,sep=";",dec=".")<br>
<br>## Funcoes extras<br>source("<a href="http://www.leg.ufpr.br/~wagner/ExemploBordaVoronoi/voronoi.R">http://www.leg.ufpr.br/~wagner/ExemploBordaVoronoi/voronoi.R</a>")<br>source("<a href="http://www.leg.ufpr.br/~wagner/ExemploBordaVoronoi/grp2sp.r">http://www.leg.ufpr.br/~wagner/ExemploBordaVoronoi/grp2sp.r</a>")<br>
<br>coords.x <- c()<br>coords.y <- c()<br>cod.mosq <- unique(mosq$mosquitrap)<br>for(i in 1:40){<br>coords.x[i] = unique(mosq[which(mosq$mosquitrap == cod.mosq[i]),]$x)<br>coords.y[i] = unique(mosq[which(mosq$mosquitrap == cod.mosq[i]),]$y)}<br>
<br>mapa.mosq <- voronoi(coords.x1=coords.x,coords.x2=coords.y,poligono = bordas)<br>mapa.mosq = grp2sp(mapa.mosq,ID = cod.mosq)<br>mapa.mosq.nb <- poly2nb(mapa.mosq)<br>plot(mapa.mosq)<br><br>-- <br>Wagner Hugo Bonat<br>
LEG - Laboratório de Estatística e Geoinformação<br>UFPR - Universidade Federal do Paraná<br>