<html>
<head>
<meta content="text/html; charset=ISO-8859-1"
http-equiv="Content-Type">
</head>
<body bgcolor="#FFFFFF" text="#000000">
como demorei demais, imagino que já deve ter solucionado, mas em
caso contrario<br>
fica alguns comandos para inspirar a solução do problema<br>
<br>
<br>
########################################################<br>
# tente pegar a lista de arquivos com o comando dir<br>
#arqs <- dir(pattern='*.hdf')<br>
arqs <- c("S20040012004031.hdf", "S20040322004060.hdf")<br>
<br>
# d1 <- as.numeric(gsub( 'S2004(...)2004(...)\\.hdf','\\1',
arqs[2] ))<br>
d2 <- as.numeric(gsub( 'S2004(...)2004(...)\\.hdf','\\2',
arqs[2] ))<br>
<br>
# descobre qual é o mês correspondente (em numero de 1 a 12)<br>
mes_number <- table((d2 - 30*(1:12)+5>0))[2], sep='')<br>
<br>
nome <- paste('mes', '04', sprintf("%02.0f", mes_number),
sep='')<br>
assign( nome, 'alguma coisa a ser salva' )<br>
########################################################<br>
<br>
qualquer duvido, escreva com mais detalhes<br>
<br>
t+<br>
cleber<br>
<br>
<br>
<br>
Em 30/03/2012 10:51, Humberto Hazin escreveu:
<blockquote cite="mid:BLU0-SMTP58BD3CD4E2D61680A7BE87B4490@phx.gbl"
type="cite">
<meta content="text/html; charset=ISO-8859-1"
http-equiv="Content-Type">
Bom dia Benilton,<br>
<br>
Obrigado por responder!<br>
<br>
É o seguinte quando baixo a imagem ela já vem com esse número todo
(S20040322004060) onde S=seawifs, em seguida vem o ano 2004032
2004060 significa o periodo da coleta do satélite ou seja 032-60
dias de 2004, ou seja mes 2. Assim para o outro arquivo
S20040012004031 001-31 dias de 2004 ou seja mes 1<br>
<br>
Esse arquivos estão em hdf4 porém preciso converter para hdf5 para
poder trabalhar! Então criei uma função no R (hdf4to5) para para
ler um executavel de um programa
(H4H5TOOLS-2.2.1-win764-vs10static) em segundo plano para esse
fim! Uso o pacote Satin para ler o arquivo convertido em hdf5 e
fazer as extrações da área que quero ter esses dados
"satinGet("S20040012004031.hdf5",lats=c(10, -30), lons=c(-20,
-50), itype="oceancolor")" como preciso de quadrantes de 1x1 grau
uso a função "simplifySat(m0401, extent=1)" em seguida criei uma
função no R "extract(m0401)" para extrair tudo em coluna Lat,
Long, Parametro...<br>
<br>
<br>
hdf4to5("S20040012004031.hdf");m0401 <-
satinGet("S20040012004031.hdf5",lats=c(10, -30), lons=c(-20, -50),
itype="oceancolor")<br>
m0401<-(simplifySat(m0401,
extent=1));m0401<-extract(m0401);m0401<-data.frame(m0401);m0401$mes<-1;
m0401$year<-2004<br>
<br>
hdf4to5("S20040322004060.hdf");m0402 <-
satinGet("S20040322004060.hdf5",lats=c(10, -30), lons=c(-20, -50),
itype="oceancolor")<br>
m0402<-(simplifySat(m0402,
extent=1));m0402<-extract(m0402);m0402<-data.frame(m0402);m0402$mes<-2;
m0402$year<-2004<br>
<br>
Eu não sei se ficou mais claro agora!<br>
<br>
Humberto<br>
<br>
<br>
Em 3/29/2012 9:25 PM, Benilton Carvalho escreveu:
<blockquote
cite="mid:CAO-arWOXyWSntOwFq3VtjhKEU2aovd8bVRC7nu94BMgv48Rv7A@mail.gmail.com"
type="cite">Sim, alguem pode te ajudar... mas depende muito de
vc ajudar o alguem primeiro.... Qual o truque para sair
de S20040322004060 para "mes 2"? (minha bola de cristal diz que
o ano corresponde aos 4 primeiros digitos apos S, mas isso eh
soh uma hipotese, jah que vc nao disse explicitamente)
<div> <br>
</div>
<div>b<br>
<br>
<div class="gmail_quote">2012/3/29 Humberto Hazin <span
dir="ltr"><<a moz-do-not-send="true"
href="mailto:hghazin@hotmail.com">hghazin@hotmail.com</a>></span><br>
<blockquote class="gmail_quote" style="margin:0 0 0
.8ex;border-left:1px #ccc solid;padding-left:1ex">
Prezados,<br>
<br>
Venho mais uma vez pedir ajuda a voces para resolver um
problema!<br>
<br>
baixei uma série de images de satelite, por ano e mes
(2004 a 2010) totalizando 163 arquivos ate ai tudo bem!
Por termos de comodidade gostaria de criar um loop para
nao ter que digitar essas 163 linhas ou mais para extrair
latitude, longitude e SST. conforme o codigo abaixo:<br>
<br>
hdf4to5("S20040012004031.hdf");m0401 <-
satinGet("S20040012004031.hdf5",lats=c(10, -30),
lons=c(-20, -50), itype="oceancolor")<br>
m0401<-(simplifySat(m0401,
extent=1));m0401<-extract(m0401);m0401<-data.frame(m0401);m0401$mes<-1;
m0401$year<-2004<br>
<br>
hdf4to5("S20040322004060.hdf");m0402 <-
satinGet("S20040322004060.hdf5",lats=c(10, -30),
lons=c(-20, -50), itype="oceancolor")<br>
m0402<-(simplifySat(m0402,
extent=1));m0402<-extract(m0402);m0402<-data.frame(m0402);m0402$mes<-2;
m0402$year<-2004<br>
<br>
Alguem poderia me ajudar!<br>
<br>
Nao coloquei o arquivo pois é muito grande quase 36gigas
porém caso precisem posso colocar o parte do codigo para
baixarem um ou mais arquivos e testarem<br>
<br>
Humberto<br>
<br>
</blockquote>
</div>
</div>
</blockquote>
</blockquote>
</body>
</html>