[R-br] Selecionar coordenadas e valor de pixels contidos no interior de parcelas

ASANTOS alexandresantosbr em yahoo.com.br
Terça Dezembro 24 16:29:09 BRST 2013


Boa tarde Pessoal,

         Estou quebrando a cabeça em um desafio natalino e gostaria de 
compartilhar o desafio com vocês, tenho uma imagem simulada stackIm1, 
sobre a qual lancei pontos aleatórios position.parc que são o centro de 
parcelas circulares CIR de 5 de raio, sendo meu CRM:

require(raster)
require(spatstat)
require(sp)

##Criando um raster 
-----------------------------------------------------------
r <- raster(nc=100, nr=100)
set.seed(20130622)
stackIm1 <- stack(lapply(1, function(x) setValues(r, 
round(runif(ncell(r))* 255))))
plot(stackIm1)
#


## Criando 15 pontos amostrais 
-------------------------------------------------
att.Im1<-as(stackIm1, "SpatialGridDataFrame")
df.Im1<-as.data.frame(att.Im1)
colnames(df.Im1)<-c("px","x","y")
gridded(df.Im1) = ~x+y
position.parc<-spsample(df.Im1,15, type="random")
points(position.parc, col="red",pch=20)
#

#Criar as parcelas circulares em torno do ponto amostral com 5 de raio
#
CIR <- vector("list",length(position.parc$x))
for(i in 1:length(position.parc$x)) CIR[[i]] <- disc(radius=5, 
centre=c(position.parc$x[i],position.parc$y[i]))
#

Gostaria de alguma luz para selecionar as informações("px","x","y") para 
os pixels contidos no interior de cada parcela
circular CIR na imagem stackIm1.

Obrigado,

-- 
======================================================================
Alexandre dos Santos
Proteção Florestal
IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso
Campus Cáceres
Caixa Postal 244
Avenida dos Ramires, s/n
Bairro: Distrito Industrial
Cáceres - MT                      CEP: 78.200-000
Fone: (+55) 65 8132-8112 (TIM)   (+55) 65 9686-6970 (VIVO)
e-mails:alexandresantosbr em yahoo.com.br
         alexandre.santos em cas.ifmt.edu.br
Lattes: http://lattes.cnpq.br/1360403201088680
======================================================================



Mais detalhes sobre a lista de discussão R-br