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

ASANTOS alexandresantosbr em yahoo.com.br
Quinta Dezembro 26 14:46:15 BRST 2013


Fantástico Éder,

      Muito obrigado,

Resolveu meu problema,


Em 26/12/2013 11:20, Éder Comunello escreveu:
> Alexandre, bom dia!
>
> Não tenho experiência com a {spatstat}, mas tentei trabalhar seu 
> código e acredito que possa lhe ser útil.
>
> ### <BEGIN>
> require(raster); require(spatstat); require(sp)
>
> ## Raster
> -----------------------------------------------------------
> set.seed(20130622)
> r <- raster(nc=100, nr=100)
> r <- setValues(r, round(runif(ncell(r))* 255))
> plot(r)
>
> r.spdf <- data.frame(coordinates(r), px=getValues(r)); names(r.spdf)
> gridded(r.spdf) <- ~x+y
> spplot(r.spdf)
>
> ## Amostras
> -------------------------------------------------
> parc.p <- spsample(r.spdf, 15, type="random")
> parc.c <- vector("list",length(parc.p))
> for(i in 1:length(parc.p)) {
> parc.c[[i]] <- disc(radius=5, centre=c(parc.p$x[i], parc.p$y[i]))}
>
> W    <- do.call(union.owin,parc.c)
> W.df <- as.data.frame(as.im <http://as.im>(W))
>
> W.sp <- W.df; coordinates(W.sp) <- ~x+y
> plot(W.sp)
>
> W.df$px <- over(W.sp, r.spdf)
> head(W.df)
> #
> #          x         y value  px
> #1 -179.9856 -39.83543     1  11
> #2 -179.9856 -39.29907     1  76
> #3 -179.9856 -38.76271     1  76
> #4 -179.9856 -38.22635     1  76
> #5 -179.9856 -37.68998     1  34
> #6 -178.7263 -41.98087     1 157
> #
> ### <END>
>
> Éder Comunello <c 
> <mailto:comunello.eder em gmail.com>omunello.eder em gmail.com 
> <mailto:omunello.eder em gmail.com>>
> Dourados, MS - [22 16.5'S, 54 49'W]
>
>
> Em 26 de dezembro de 2013 08:37, ASANTOS 
> <alexandresantosbr em yahoo.com.br 
> <mailto:alexandresantosbr em yahoo.com.br>> escreveu:
>
>     Bom dia Éder,
>
>       Conforme sugerido, adicionei  ao meu CRM:
>
>     W <- do.call(union.owin,CIR)
>     #
>     imW<-as.im <http://as.im>(W)
>     grids<-as(imW, "SpatialGridDataFrame")
>     fullgrid(grids)=FALSE
>     plot(grids)
>     points(df.Im1)
>     over(grids,df.Im1))##
>
>         O over fez a listagem do valor dos pixels na imagem (stackIm1)
>     contidos no interior das parcelas circulares CIR, mas meu objetivo
>     é listar qual a coordenada e o valor de cada pixel em cada parcela
>     circular CIR, com isso terei um resultado em colunas contendo:
>     CIR, x, y e px
>
>     Alexandre
>
>
>     Em 26/12/2013 05:53, Éder Comunello escreveu:
>>     Bom dia, Alexandre!
>>
>>     Já deu uma olhada em sp::over?
>>
>>     -- 
>>     Éder Comunello <c
>>     <mailto:comunello.eder em gmail.com>omunello.eder em gmail.com
>>     <mailto:omunello.eder em gmail.com>>
>>     Dourados, MS - [22 16.5'S, 54 49'W]
>>
>>
>>
>>     _______________________________________________
>>     R-br mailing list
>>     R-br em listas.c3sl.ufpr.br  <mailto:R-br em listas.c3sl.ufpr.br>
>>     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.
>
>     -- 
>     ======================================================================
>     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  <tel:%28%2B55%29%2065%208132-8112>  (TIM)(+55) 65 9686-6970  <tel:%28%2B55%29%2065%209686-6970>  (VIVO)
>     e-mails:alexandresantosbr em yahoo.com.br  <mailto:e-mails:alexandresantosbr em yahoo.com.br>  
>              alexandre.santos em cas.ifmt.edu.br  <mailto:alexandre.santos em cas.ifmt.edu.br>  
>     Lattes:http://lattes.cnpq.br/1360403201088680                    
>     ======================================================================
>
>
>     _______________________________________________
>     R-br mailing list
>     R-br em listas.c3sl.ufpr.br <mailto:R-br em listas.c3sl.ufpr.br>
>     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.
>
>
>
>
> _______________________________________________
> R-br mailing list
> R-br em listas.c3sl.ufpr.br
> 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.

-- 
======================================================================
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
======================================================================

-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20131226/e773a8cf/attachment.html>


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