<html><body><div style="color:#000; background-color:#fff; font-family:times new roman, new york, times, serif;font-size:12pt"><div style="font-family: 'times new roman', 'new york', times, serif; font-size: 12pt;"><span>Ludmila,</span></div><div style="font-family: 'times new roman', 'new york', times, serif; font-size: 16px; color: rgb(0, 0, 0); background-color: transparent; font-style: normal;"><span><br></span></div><div style="font-family: 'times new roman', 'new york', times, serif; font-size: 16px; color: rgb(0, 0, 0); background-color: transparent; font-style: normal;">Por que o seu código não funciona? Aqui no meu computador funcionou com um raster aleatório.</div><div style="font-family: 'times new roman', 'new york', times, serif; font-size: 16px; color: rgb(0, 0, 0); background-color: transparent; font-style: normal;"><br></div><div style="font-family: 'times new roman', 'new york', times, serif; font-size: 16px; color: rgb(0, 0, 0);
background-color: transparent; font-style: normal;">Você tentou usar a função clump? Ela pode ser mais próxima do que você precisa. Veja abaixo:</div><div style="font-family: 'times new roman', 'new york', times, serif; font-size: 16px; color: rgb(0, 0, 0); background-color: transparent; font-style: normal;"><br></div><div style="background-color: transparent;">library(raster)</div><div style="background-color: transparent;">library(igraph)</div><div style="background-color: transparent;"><br></div><div style="background-color: transparent;"># Create random raster</div><div style="background-color: transparent;">r <- raster(ncols=12, nrows=12)</div><div style="background-color: transparent;">set.seed(0)</div><div style="background-color: transparent;">r[] <- round(runif(ncell(r))*0.7 )</div><div style="background-color: transparent;"><br></div><div style="background-color: transparent;"># Detect clumps and plot result</div><div
style="background-color: transparent;">rc <- clump(r)</div><div style="background-color: transparent;">plot(rc)</div><div style="background-color: transparent;"><br></div><div style="background-color: transparent;"># Extract coordinates</div><div style="background-color: transparent;">xy <- coordinates(rc)</div><div style="background-color: transparent;">lon <- xFromCol(r,1:ncol(r))</div><div style="background-color: transparent;">lat <- yFromRow(r,1:nrow(r))</div><div style="font-family: 'times new roman', 'new york', times, serif; font-size: 12pt;"><br></div><div style="font-family: 'times new roman', 'new york', times, serif; font-size: 12pt;">Se o código acima não funcionar no seu computador, por favor retorne a mensagem de erro E o raster que você está usando para que a lista possa ajudar.</div><div style="font-family: 'times new roman', 'new york', times, serif; font-size: 12pt;"><span style="font-size:
12pt;"> </span><br></div><div style="font-family: 'times new roman', 'new york', times, serif; font-size: 12pt;">Saudações,<br>--<br>Thiago V. dos Santos<br>PhD student<br>Land and Atmospheric Science<br>University of Minnesota<br>http://www.laas.umn.edu/CurrentStudents/MeettheStudents/ThiagodosSantos/index.htm<br>Phone: (612) 323 9898</div><div style="font-family: 'times new roman', 'new york', times, serif; font-size: 12pt;"><br></div> <div style="font-family: 'times new roman', 'new york', times, serif; font-size: 12pt;"> <div style="font-family: 'times new roman', 'new york', times, serif; font-size: 12pt;"> <div dir="ltr"> <hr size="1"> <font size="2" face="Arial"> <b><span style="font-weight:bold;">From:</span></b> Ludmila Rattis <ludmilarattis@gmail.com><br> <b><span style="font-weight: bold;">To:</span></b> R-br@listas.c3sl.ufpr.br <br> <b><span style="font-weight: bold;">Sent:</span></b> Tuesday, August 20, 2013 1:54 PM<br>
<b><span style="font-weight: bold;">Subject:</span></b> [R-br] get coordinates of patches (clumps) in a raster map<br> </font> </div> <div class="y_msg_container"><br><div id="yiv1413123331"><div dir="ltr">Prezados,<div><br></div><div>Preciso extrair as coordenadas dos centroides de diferentes clumps de um raster. Para isso, usei o seguinte código: </div><div><br></div><div><div>p <- data.frame(rasterToPoints(meu_raster))</div>
<div>a<- p[p$layer > 0,]<br></div><div>c<-sapply(split(p[, c("x", "y")], p$layer), colMeans)<br></div><div><br></div><div>Contudo, esse código não está funcionando.</div><div><br></div><div>Alguém saberia me dizer como posso obter tais coordenadas?</div>
<div><br></div><div>[]'s</div><div><br></div><div><b>Ludmila Rattis</b><br></div>Programa de Pós-graduação em Ecologia/UNICAMP<br>
Conservation Biogeography Lab<br>
<a rel="nofollow" target="_blank" href="http://www.wix.com/rdloyola/lab">http://www.wix.com/rdloyola/lab</a><a rel="nofollow" target="_blank" href="http://sites.ffclrp.usp.br/ficus"></a>
</div></div></div><br>_______________________________________________<br>R-br mailing list<br><a ymailto="mailto:R-br@listas.c3sl.ufpr.br" href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><br><a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.<br><br></div> </div> </div> </div></body></html>