[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