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@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@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@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.