Obrigado Jônatan,
Perfeita essa função usando getValuesBlocke consegui fazer a extração de 10 valores na horizontal à partir de uma coordenada arbitrária como xy(-150,80), porém não consigo agora capturar além dos valores dos pixel a coordenada central desses pixel pelo output dado pela função, poderia me ajudar?
No meu caso tentei usar a função xyFromCell():
require(raster); require(sp)
set.seed(1)
## Criando um raster
r <- raster(nc=10, nr=10)
r <- setValues(r, round(runif(ncell(r))* 255))
#Consigo extrair os valores dos pixels de interesse dadas as
coordenadas centrais dos mesmos
x <- c(-150)
y <- c(-80)
p <- data.frame(x,y)
pontos <- SpatialPoints(p)
p$cel <- cellFromXY(r, pontos)
p$col <- colFromCell(r, p$cel)
p$row <- rowFromCell(r, p$cel)
p
plot(r)
text(r)
points(pontos, pch = 4, col = 2)
library(plyr)
vals2cols <- ldply(1:nrow(p),
function(ir){
getValuesBlock(r,
col = p$col[ir],
ncols = 10,
row = p$row[ir],
nrows = 1)
}# end fun
)
df <- data.frame(p, vals2cols)
df
p.cen <- xyFromCell(r, vals2cols); p.cen
#END-------------------------------------------------------------------------------------------------------
-- ====================================================================== 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 99686-6970 (VIVO) (+55) 65 3221-2674 (FIXO) e-mails:alexandresantosbr@yahoo.com.br alexandre.santos@cas.ifmt.edu.br Lattes: http://lattes.cnpq.br/1360403201088680 OrcID: orcid.org/0000-0001-8232-6722 - ResearcherID: A-5790-2016 Researchgate: www.researchgate.net/profile/Alexandre_Santos10 LinkedIn: br.linkedin.com/in/alexandre-dos-santos-87961635 Mendeley:www.mendeley.com/profiles/alexandre-dos-santos6/ ======================================================================
use a função getValuesBlock do raster.
require(raster); require(sp)set.seed(1)## Criando um rasterr <- raster(nc=10, nr=10)r <- setValues(r, round(runif(ncell(r))* 255))#Consigo extrair os valores dos pixels de interesse dadas as coordenadas centrais dos mesmosx <- c(54,18,54,54)y <- c(27,81,9,63)p <- data.frame(x,y)pontos <- SpatialPoints(p)
p$cel <- cellFromXY(r, pontos)p$col <- colFromCell(r, p$cel)p$row <- rowFromCell(r, p$cel)pplot(r)text(r)points(pontos, pch = 4, col = 2)
library(plyr)vals2cols <- ldply(1:nrow(p),function(ir){getValuesBlock(r,col = p$col[ir],ncols = 2,row = p$row[ir],nrows = 1)}# end fun)df <- data.frame(p, vals2cols)df
2017-01-26 17:43 GMT-02:00 ASANTOS via R-br <r-br@listas.c3sl.ufpr.br>:
Prezados Membros,
Tenho um raster hipotético r:
require(raster); require(sp)
## Criando um raster
r <- raster(nc=10, nr=10)
r <- setValues(r, round(runif(ncell(r))* 255))
Consigo extrair os valores dos pixels de interesse dadas as coordenadas centrais dos mesmos:
x <- c(54,18,54,54)
y <- c(27,81,9,63)
pontos <- SpatialPoints(cbind(x,y))
cells <- cellFromXY(r, pontos); cells
p.cen <- xyFromCell(r, cells); p.cen
values <- extract(r,p.cen)
dados<-cbind(p.cen,cells,values)
No entanto, eu gostaria de encontrar um modo para extrair o valor dos pixels, no qual, eu selecionaria uma coordenada e a função realizaria a extração dos pixels de linhas verticais inteiras do raster a cada dois pixels à partir da coordenada dada, por exemplo?
Alguém poderia me dar uma ideia de como chegar a esse resultado?
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 99686-6970 (VIVO) (+55) 65 3221-2674 (FIXO)
e-mails:alexandresantosbr@yahoo.com.br
alexandre.santos@cas.ifmt.edu.br
Lattes: http://lattes.cnpq.br/1360403201088680
OrcID: orcid.org/0000-0001-8232-6722 - ResearcherID: A-5790-2016
Researchgate: www.researchgate.net/profile/Alexandre_Santos10
LinkedIn: br.linkedin.com/in/alexandre-dos-santos-87961635
Mendeley:www.mendeley.com/profiles/alexandre-dos-santos6/
============================================================ ==========
_______________________________________________
R-br mailing list
R-br@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 fornea cdigo mnimo reproduzvel.
--
###############################################################
## Jônatan Dupont Tatsch## Professor do Departamento de Física## Centro de Ciências Exatas e Naturais (CCNE)
## Universidade Federal de Santa Maria - UFSM## Faixa de Camobi, Prédio 13 - Campus UFSM - Santa Maria, RS, Brasil - 97105-900
## Telefone: +55(55)33012083
###############################################################