Elias,Pelo o que eu entendi, vc usou a funcao deldir(), no caso Voronoi para fazer os polígono da área. Se for isso acho que passei pelo mesmo processo.Vou passar algumas etapas e uma função que improvisei e no meu caso deu certo. Vou explicar os procedimentos:1) com a funcao deldir() vc faz seu polígonos na áreasss <- deldir( .....)entao vc tem o objeto sss2) Para gerar uma lista do objeto "sss" use a funcao tile.list()sss_list <- tile.list(sss)3) Agora desse objeto que é uma lista, vc tem que extrair o conjunto de polígonospara isso improvisei uma função, grave ela no seu ambiente:FUN_VOR_MAT <- function (A,TT) #TT é o arquivo, N o número do POL{AA <- 1BB <- 2NROW <- length (A[[TT]]$x) + 1 #add primeiro pontoVORONOI_MATRIX <- matrix(0,NROW,2)VORONOI_MATRIX[1:NROW,AA] <- A[[TT]]$x[1:NROW]VORONOI_MATRIX[1:NROW,BB] <- A[[TT]]$y[1:NROW]VORONOI_MATRIX[NROW,AA] <- A[[TT]]$x[1] #+ 1 add primeiro pontoVORONOI_MATRIX[NROW,BB] <- A[[TT]]$y[1] #+ 1 add primeiro pontoif(TT == NROW){print("17/10/2014, [hidden email] ")}return(VORONOI_MATRIX)}4) Agora vc vai usar a funcao:sss_areas <- lapply(1:length(sss_list), function (x) ( FUN_VOR_MAT(sss_list,x)) )5) #Formando Polygonssss_pols <- lapply(1:length(sss_areas), function(x) Polygons(list(Polygon(sss_areas[[x]])),ID = paste("reg", x, sep = "")))6) Para converter em SpatialPolyginssss_SPATIAL_POLS <- as.SpatialPolygons.PolygonsList(sss_pols, proj4string=CRS(as.character(NA)))Se tudo correu bem é só comecar a usar o poly2nb(sss_SPATIAL_POLS)Boa sorteRafael Tieppo
_______________________________________________
R-br mailing list
[hidden email]
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.
If you reply to this email, your message will be added to the discussion below:http://r-br.2285057.n4.nabble.com/R-br-skater-e-deldir-tp4663303p4663325.html