extraindo dados de arquivo img
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
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... " ### 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.
Veja também <http://power.larc.nasa.gov/common/php/POWER_ParametersAgro.php> que disponibiliza muitos desses dados, com opções melhores de seleção e download em ASCII. ================================================ É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 14:16, Éder Comunello <comunello.eder@gmail.com> escreveu:
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... "
### 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.
participantes (2)
-
Hélio Gallo Rocha -
Éder Comunello