Éder <br><br>seu arquivo binario (que está no formato do software Grads - <a href="http://www.iges.org/grads/">http://www.iges.org/grads/</a>) tem somente 49 variáveis (da "temp" até a "tsfc"), embora o arquivo descritor tenha 58 variáveis. Eu abri o arquivo no grads e verifiquei-o.<br>
<br>Portanto, a leitura dos dados no R está ok (195*200*49 = 1911000).<br><br>[ ]'s<br><br><br><div class="gmail_quote">2011/4/14 <span dir="ltr"><<a href="mailto:eder@leg.ufpr.br">eder@leg.ufpr.br</a>></span><br>
<blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">Isso que é minha maior duvida, pois no readbin o argumento size aceita no<br>
meu pc no maximo 4 byte. não sou especialista nisso pois creio que o<br>
numero 1E20 não poder ser representado com 4 bytes.<br>
Éder<br>
<div><div></div><div class="h5"><br>
<br>
> Arquivos GRIB padrao tem cabecalho obrigatorio, mas pode nao ser esse<br>
> o seu caso.<br>
><br>
> De qualquer forma, vc precisa dos specs do seu .dat... pode ser que o<br>
> tamanho das variaveis sejam diferentes do padrao (4 bytes para<br>
> inteiros e 8 para doubles)... pode ser que a representacao seja<br>
> diferente (o que acho ser o caso de GRIB para medidas em double).<br>
><br>
> b<br>
><br>
> 2011/4/14 <<a href="mailto:eder@leg.ufpr.br">eder@leg.ufpr.br</a>>:<br>
>> Benilton,<br>
>> O formato dos números creio ser double, visto que quanto abri os<br>
>> arquivos<br>
>> que deram certo, os valores são coerentes com a variável, e pelo o que<br>
>> entendi este bin não tem cabeçalho, visto que quem diz como é o bin é o<br>
>> arquivo .ctl que faz o papel do cabeçalho.<br>
>> Desta forma ainda não entendi o que pode estar acontecendo.<br>
>> Att.<br>
>><br>
>>> muito pouco provavel que o seu readBin() esteja correto.<br>
>>><br>
>>> arquivos binarios contem cabecalhos e campos que nao sao<br>
>>> necessariamente 'double'... entao, antes de tudo, vc precisa ter uma<br>
>>> descricao do formato para poder usar o readBin() apropriadamente....<br>
>>> e, como vc deve ter notado, o uso incorreto vai te trazer resultados<br>
>>> inesperados.<br>
>>><br>
>>> b<br>
>>><br>
>>> 2011/4/14 <<a href="mailto:eder@leg.ufpr.br">eder@leg.ufpr.br</a>>:<br>
>>>> Pessoal,<br>
>>>><br>
>>>> Minha duvida agora e com arquivos binarios, contextualizando:<br>
>>>> Quero abrir arquivos de saida de modelos climáticos (Eta), estes<br>
>>>> modelos<br>
>>>> soltam arquivo no formato GRIB, em em palavras R-anas uma array onde<br>
>>>> cada<br>
>>>> variável meteorológica seria um z dento do array (dados[x,y,z]), no<br>
>>>> arquivo ctl (indica como o binario é) indica qual o código falha (NA),<br>
>>>> em<br>
>>>> arquivo que o código era -9999 a função readBin funcionou, em outro<br>
>>>> arquivo que o código era 1.0E20 ela não funcionou ou seja estou<br>
>>>> suspeitando que ela não esta entendo este código falha, exemplo:<br>
>>>><br>
>>>> formato <-<br>
>>>> readLines('<a href="http://www.leg.ufpr.br/%7Eeder/ll20gfsa2006120100+2006120100.ctl" target="_blank">http://www.leg.ufpr.br/~eder/ll20gfsa2006120100+2006120100.ctl</a>')<br>
>>>> formato<br>
>>>> cone <-<br>
>>>> file('<a href="http://www.leg.ufpr.br/%7Eeder/ll20gfsa2006120100+2006120100.dat" target="_blank">http://www.leg.ufpr.br/~eder/ll20gfsa2006120100+2006120100.dat</a>','rb')<br>
>>>> dados <- readBin(cone,numeric(),195*200*58,4)<br>
>>>> length(dados)<br>
>>>> head(dados)<br>
>>>> close(cone)<br>
>>>><br>
>>>> pelas informações em formato era para ter no arquivo binário<br>
>>>> 195*200*58=<br>
>>>> 2262000 números, porem ela carrega apenas 1911000, faltando dados,<br>
>>>> desta<br>
>>>> forma. Alguma dica?<br>
>>>> Atenciosamente<br>
>>>><br>
>>>><br>
>>>><br>
>>>><br>
>>>> _______________________________________________<br>
>>>> R-br mailing list<br>
>>>> <a href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><br>
>>>> <a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>
>>>><br>
>>><br>
>><br>
>><br>
>><br>
>><br>
>><br>
><br>
<br>
<br>
_______________________________________________<br>
R-br mailing list<br>
<a href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><br>
<a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>
</div></div></blockquote></div><br><br clear="all"><br>-- <br>#------------------------------------------------------------------------------#<br># Jônatan Dupont Tatsch #<br>
# Climate and Biosphere Laboratory #<br># Department of Atmospheric Sciences #<br># Institute of Astronomy, Geophysics and Atmospheric Sciences #<br>
# University of São Paulo #<br># Rua do Matão, 1226 #<br># Cid. Universitária, São Paulo, SP, Brazil, CEP: 05508-090 #<br>
# Phone:+55 11 3091-4772, Fax:+55 11 3091-4714 #<br># <a href="http://jonatandupont.weebly.com/index.html" target="_blank">http://jonatandupont.weebly.com/index.html</a> #<br>#------------------------------------------------------------------------------#<br>