[R-br] Obtenção de Valores em Pontos de Arquivo Raster x ASCII

[ASANTOS]

Boa noite Éder,

         Montei uma rotina abaixo para iniciarmos um novo tópico e é 
possível observar o que você havia comentado, sendo:

### <code r>
require(raster)
require(sp)
require(rgdal)

# RasterLayer inventado
r <- raster(nrows=10000, ncols=10000)
r <- setValues(r, 1:ncell(r))
r_SGDF<- as(r, 'SpatialGridDataFrame')
writeGDAL(r_SGDF,"r_test.tif")### Cria o Tiff
write.asciigrid(r_SGDF,"r_test.asc")### Cria o ASCII
#

# Lendo o tiff
r_tif<-raster(c("r_test.tif"))
s <- stack(r_tif)
POI <- SpatialPoints(cbind(lon=c(-47.5,-47.5,-44.5,-40), 
lat=c(-18.5,-19.5,-20.5,-22)))
proj4string(POI) <- proj4string(s)
e <- extract(s, POI, method= "bilinear", df=T) #ou method="simple"
proc.time()
#

# Lendo o ASCII
r_ascii<-read.asciigrid("r_test.asc")
POI <- SpatialPoints(cbind(lon=c(-47.5,-47.5,-44.5,-40), 
lat=c(-18.5,-19.5,-20.5,-22)))
proj4string(POI) <- proj4string(s)
e <- extract(s, POI, method= "bilinear", df=T) #ou method="simple"
proc.time()
### </code>

-- 
======================================================================
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
OrcID: orcid.org/0000-0001-8232-6722
Researchgate: https://www.researchgate.net/profile/Alexandre_Santos10
LinkedIn: https://br.linkedin.com/in/alexandre-dos-santos-87961635
======================================================================

Em 17/02/2016 16:07, Éder Comunello escreveu:
> Alexandre, boa tarde!
>
> Sobre o formato dos dados, fiz uns testes aqui e dá pra tirar umas 
> conclusões, mas acho melhor criar um tópico novo pra não embolar. 
> Resumidamente, ao acessar usando o pacote raster (antes que rgdal) 
> você mantém a informação no disco e economiza memória. Você poderia 
> usar um formato ascii, mas os binários otimizam a velocidade de 
> leitura e acesso.
>
> ​
> ================================================
> Éder Comunello
> Agronomist (UEM), MSc in Environ. Sciences (UEM)
> DSc in Agricultural Systems Engineering (USP/Esalq)
> Brazilian Agricultural Research Corporation (Embrapa)
> Dourados, MS, Brazil |<O>|
> ================================================
> GEO, -22.2752, -54.8182, 408m
> UTC-04:00 / DST: UTC-03:00
>
>

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