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