[R-br] extraindo dados de arquivo img

Éder Comunello comunello.eder em gmail.com
Quinta Janeiro 21 15:16:04 BRST 2016


Boa tarde, Hélio!

Segue uma ideia inicial...

### <code r>
setwd("D:/Temp") ### alterar aqui!

url <- "
http://spirits.jrc.ec.europa.eu/files/ecmwf/ope/south-america/tmin/ope_south-america_tmin_20151121.zip
"

### Faz download, caso ainda não tenha sido feito....
if (!file.exists(basename(url))) download.file(url, basename(url),
mode="wb")

unzip(basename(url), list=T) -> fnames
# 1 ope_south-america_tmin_20151121.img 341796 2015-12-04 11:22:00
# 2 ope_south-america_tmin_20151121.hdr    547 2015-12-04 11:22:00

### descompacta
for (fn in fnames$Name) unzip(basename(url), fn)
file.info(fn)

require(rgdal)

### .hdr -> header (cabeçalho) em formato texto
GDALinfo(fnames$Name[1])                 ### informações a partir do
arquivo .img
cat(readLines(fnames$Name[2]), sep="\n") ### informações do .hdr
(basicamente, as mesmas)

img <- readGDAL(fnames$Name[1]) ### armazena em um SGDF (pacote sp)
spplot(img)
head(coordinates(img))
head(img em data) ### dados - nesse caso temp. mínima
### </code>

​
================================================
Éder Comunello
PhD Student in Agricultural Systems Engineering (USP/Esalq)
Brazilian Agricultural Research Corporation (Embrapa)
Dourados, MS, Brazil [22 16.5'S, 54 49.0'W]




Em 21 de janeiro de 2016 12:12, Hélio Gallo Rocha <heliogallorocha em gmail.com
> escreveu:

> Boa tarde Listeiros,
>
> Alguém poderia me dar uma mão nesta tarefa:
>
> O site : http://spirits.jrc.ec.europa.eu/
>
> dispõem  arquivos climatológicos com as extenções : img e hdr
>
> Quem duipõem um CRM ou indicaria um pacote que extraia os dados e suas
> coordenadas?
>
> Grato
>
>
> --
> Hélio Gallo Rocha
> IFSULDEMINAS - Câmpus Muzambinho
>
> _______________________________________________
> 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.
>
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20160121/97936247/attachment.html>


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