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á