Preciso extrair as coordenadas dos centroides de diferentes clumps de um raster. Para isso, usei o seguinte código:
p <- data.frame(rasterToPoints(meu_raster))
a<- p[p$layer > 0,]
c<-sapply(split(p[, c("x", "y")], p$layer), colMeans)
Contudo, esse código não está funcionando.
Alguém saberia me dizer como posso obter tais coordenadas?
[]'s
Ludmila Rattis
Programa de Pós-graduação em Ecologia/UNICAMP
Conservation Biogeography Lab
http://www.wix.com/rdloyola/lab