[R-br] loops no R

[Humberto Hazin]

Bom dia Benilton,

Obrigado por responder!

É 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

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


hdf4to5("S20040012004031.hdf");m0401 <- 
satinGet("S20040012004031.hdf5",lats=c(10, -30), lons=c(-20, -50), 
itype="oceancolor")
  m0401<-(simplifySat(m0401, 
extent=1));m0401<-extract(m0401);m0401<-data.frame(m0401);m0401$mes<-1; 
m0401$year<-2004

  hdf4to5("S20040322004060.hdf");m0402 <- 
satinGet("S20040322004060.hdf5",lats=c(10, -30), lons=c(-20, -50), 
itype="oceancolor")
  m0402<-(simplifySat(m0402, 
extent=1));m0402<-extract(m0402);m0402<-data.frame(m0402);m0402$mes<-2; 
m0402$year<-2004

Eu não sei se ficou mais claro agora!

Humberto


Em 3/29/2012 9:25 PM, Benilton Carvalho escreveu:
> 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)
>
> b
>
> 2012/3/29 Humberto Hazin <hghazin em hotmail.com 
> <mailto:hghazin em hotmail.com>>
>
>     Prezados,
>
>     Venho mais uma vez pedir ajuda a voces para resolver um problema!
>
>     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:
>
>     hdf4to5("S20040012004031.hdf");m0401 <-
>     satinGet("S20040012004031.hdf5",lats=c(10, -30), lons=c(-20, -50),
>     itype="oceancolor")
>      m0401<-(simplifySat(m0401,
>     extent=1));m0401<-extract(m0401);m0401<-data.frame(m0401);m0401$mes<-1;
>     m0401$year<-2004
>
>      hdf4to5("S20040322004060.hdf");m0402 <-
>     satinGet("S20040322004060.hdf5",lats=c(10, -30), lons=c(-20, -50),
>     itype="oceancolor")
>      m0402<-(simplifySat(m0402,
>     extent=1));m0402<-extract(m0402);m0402<-data.frame(m0402);m0402$mes<-2;
>     m0402$year<-2004
>
>     Alguem poderia me ajudar!
>
>     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
>
>     Humberto
>
>     _______________________________________________
>     R-br mailing list
>     R-br em listas.c3sl.ufpr.br <mailto: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.
>
>
>
>
> _______________________________________________
> 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/20120330/0655f9b7/attachment.html>