Diagrama de voronoi
Boa tarde pssoal, Estou usando o pacote tripack para a construção de diagrama de voronoi.
data(tritest) vm <- voronoi.mosaic(tritest) plot.voronoi(vm, do.points=T) points(tritest)
Nesse exemplo, a "borda" do diagrama é um retângulo. É possível utilizar alguma função já existente para criar o diagrama de forma que a "borda" seja o polígono que eu especificar? (Algo como o que é feito no software MapInfo, para os que o conhecem) Fabrício. http://vizzual.net
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á
Wagner, funcionou perfeitamente! Muito obrigado! Fabrício. http://vizzual.net 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
participantes (2)
-
Fabrício Tavares -
Wagner bonat