Humberto sugiro que dê uma olhada nesse material do Walmes.<div><br></div><div><a href="http://www.leg.ufpr.br/~walmes/ensino/ce223-2011-01/ce223-2011-01-aula22.R">http://www.leg.ufpr.br/~walmes/ensino/ce223-2011-01/ce223-2011-01-aula22.R</a></div>
<div><br></div><div><br></div><div>[ ]'s</div><div><br></div><div><br><br><div class="gmail_quote">2012/10/1 Humberto Hazin <span dir="ltr"><<a href="mailto:hghazin@hotmail.com" target="_blank">hghazin@hotmail.com</a>></span><br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div><div dir="ltr">Jonatan,<br><br>Muito obrigado! Problema resolvido!!! Porém tenho uma pergunta em relação o "time" que na verdade sao datas mm/dd/yy, contudo aparece como numeros (ex. 731962) como posso fazer para converter-las para datas mesmo ou colocar em colunas dia, mes, ano?<div class="im">
<br><br>Humberto<br><br><font style="font-weight:bold;font-size:10pt"><span style="color:rgb(0,0,255)"><span style="color:rgb(0,0,0)">__________________________________________________</span><br></span></font><pre>Prof. Dr. Humberto Gomes Hazin<br>
Universidade Federal Rural de Pernambuco - UFRPE<br>Unidade Acadêmica de Serra Talhada - UAST<br>Coordenador do Curso de Engenharia de Pesca<br>Fazenda Saco, S/N - Caixa Postal 063<br>CEP 56900-000<br>Serra Talhada - PE<br>
Tel.(87) 3831-2206/3831-1927<br>Cel.(81) 9158 0564</pre><br><br></div><div><div></div><hr>Date: Sun, 30 Sep 2012 18:58:13 -0300<br>From: <a href="mailto:jdtatsch@gmail.com" target="_blank">jdtatsch@gmail.com</a><br>To: <a href="mailto:r-br@listas.c3sl.ufpr.br" target="_blank">r-br@listas.c3sl.ufpr.br</a><br>
Subject: Re: [R-br] Organizacao de dados Netcdf<div><div class="h5"><br><br>Humberto, voce pode fazer isso utilizando a pacote ncdf (versao 1.6.6) e a funcao melt do pacote reshape, conforme CMR abaixo.<div><br></div><div>
<br></div><div><div>require(ncdf) #‘ncdf’ version 1.6.6</div><div>require(reshape)</div>
<div> </div><div>nc <- open.ncdf("DML05_12.nc")</div><div>nc</div><div># [1] "file DML05_12.nc has 4 dimensions:"</div><div># [1] "LONN55_12 Size: 68"</div><div># [1] "LAT46_103 Size: 58"</div>
<div># [1] "TIME Size: 96"</div><div># [1] "bnds Size: 2"</div><div># [1] "------------------------"</div><div># [1] "file DML05_12.nc has 2 variables:"</div><div># [1] "double TIME_bnds[bnds,TIME] Longname:TIME_bnds Missval:1e+30"</div>
<div># [1] "float MLD[LONN55_12,LAT46_103,TIME] Longname:argo mixed layer depth (m) Missval:-999"</div><div> </div><div>lat <- get.var.ncdf(nc, "LAT46_103")</div><div>lon <- get.var.ncdf(nc, "LONN55_12")</div>
<div>time <- get.var.ncdf(nc, "TIME")</div><div><br></div><div> ## array da variavel MLD</div><div> dml <- get.var.ncdf(nc, "MLD")</div><div> str(dml)</div><div> # num [1:68, 1:58, 1:96] 25.4 25.9 26.4 26.9 27.3 ...</div>
<div> </div><div> ## atribuindo dimnames conforme lon, lat e tempo</div><div> dimnames(dml)[[1]] <- as.character(lon)</div><div> dimnames(dml)[[2]] <- as.character(lat)</div><div> dimnames(dml)[[3]] <- as.character(time)</div>
<div> str(dml)</div><div> #num [1:68, 1:58, 1:96] 25.4 25.9 26.4 26.9 27.3 ...</div><div> #- attr(*, "dimnames")=List of 3</div><div> #..$ : chr [1:68] "-55.5" "-54.5" "-53.5" "-52.5" ...</div>
<div> #..$ : chr [1:58] "-44.5" "-43.5" "-42.5" "-41.5" ...</div><div> #..$ : chr [1:96] "731962" "731993" "732021" "732052" ... </div><div>
</div><div> ## melt para reestruturacao dos dados</div><div> mdml <- melt(dml)</div><div> </div><div> head(mdml)</div><div> # X1 X2 X3 value</div><div> #1 -55.5 -44.5 731962 25.41</div><div> #2 -54.5 -44.5 731962 25.91</div>
<div> #3 -53.5 -44.5 731962 26.42</div><div> #4 -52.5 -44.5 731962 26.89</div><div> #5 -51.5 -44.5 731962 27.30</div><div> #6 -50.5 -44.5 731962 27.63 </div><div> </div><div> tail(mdml)</div><div> # X1 X2 X3 value</div>
<div> #378619 6.5 12.5 734853 NA</div><div> #378620 7.5 12.5 734853 NA</div><div> #378621 8.5 12.5 734853 NA</div><div> #378622 9.5 12.5 734853 NA</div><div> #378623 10.5 12.5 734853 NA</div><div>
#378624 <a href="tel:11.5%2012.5%20734853" value="+15125734853" target="_blank">11.5 12.5 734853</a> NA</div>
<div><br></div><div>require(lattice)</div>levelplot(value~X1+X2|factor(X3),data=mdml)</div><div><br></div><div><br><div><br></div><br><div>2012/9/30 Humberto Hazin <span dir="ltr"><<a href="mailto:hghazin@hotmail.com" target="_blank">hghazin@hotmail.com</a>></span><br>
<blockquote style="border-left:1px #ccc solid;padding-left:1ex">
<div><div dir="ltr">
<br>Bom dia pessoal,<br><br>Estou tentando extrair variáveis ambientais de arquivos formato netcdf. Até um certo momento consigo separtar em arquivos as Latitudes, Longitudes, time (mes/ano) e a dml (exemplo que mando aqui). <br>
<br>Link para baixar os dados do exemplo<br><a href="https://www.dropbox.com/s/66dp4rgei791r52/DML05_12.nc" target="_blank">https://www.dropbox.com/s/66dp4rgei791r52/DML05_12.nc</a><br><br>nc <- <a href="http://open.nc" target="_blank">open.nc</a>("DML05_12.nc", write=TRUE)<br>
<a href="http://print.nc" target="_blank">print.nc</a>(nc)<br>lat<-<a href="http://var.get.nc" target="_blank">var.get.nc</a>(nc, "LAT46_103")<br>lon<-<a href="http://var.get.nc" target="_blank">var.get.nc</a>(nc, "LONN55_12")<br>
time<-<a href="http://var.get.nc" target="_blank">var.get.nc</a>(nc, "TIME")<br>dml<-<a href="http://var.get.nc" target="_blank">var.get.nc</a>(nc, "MLD")<br><br><br>Meu problema é que nao estou conseguindo organizar eles em colunas com
seus respectivos valores para cada latitude, longitude, time e dml. <br><br>o output final deveria ficar assim:<br><span style="text-indent:0px;letter-spacing:normal;font-variant:normal;text-align:-webkit-left;font-style:normal;background-color:rgb(225,226,229);font-weight:normal;line-height:15px;border-collapse:separate;text-transform:none;font-size:13px;white-space:pre-wrap;font-family:'Lucida Console';word-spacing:0px"><pre style="font-family:'Lucida Console';font-size:10pt!important;outline-style:none;outline-width:initial;outline-color:initial;border-top-style:none;border-right-style:none;border-bottom-style:none;border-left-style:none;border-width:initial;border-color:initial;white-space:pre-wrap!important;margin-right:0px;margin-bottom:0px;margin-left:0px;line-height:1.2">
lat lon time dml
1 -44.5 -55.5
2 -43.5 -55.5
3 -42.5 -55.5
4 -41.5 -55.5
5 -40.5 -55.5
6 -39.5 -55.5</pre></span><br>Agradeço a ajuda desde já<br><br>Humberto<br><br><br><br><br><br><br><font style="font-weight:bold"><span style="color:rgb(0,0,255)"><span style="color:rgb(0,0,0)">__________________________________________________</span><br>
</span></font><pre>Prof. Dr. Humberto Gomes Hazin<br>Universidade Federal Rural de Pernambuco - UFRPE<br>Unidade Acadêmica de Serra Talhada - UAST<br>Coordenador do Curso de Engenharia de Pesca<br>Fazenda Saco, S/N - Caixa Postal 063<br>
CEP 56900-000<br>Serra Talhada - PE<br>Tel.(87) 3831-2206/3831-1927<br>Cel.(81) 9158 0564</pre> </div></div>
<br>_______________________________________________<br>
R-br mailing list<br>
<a href="mailto:R-br@listas.c3sl.ufpr.br" target="_blank">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>
Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.<br></blockquote></div><br><br clear="all"><div><br></div>-- <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:<a href="tel:%2B55%2011%203091-4772" value="+551130914772" target="_blank">+55 11 3091-4772</a>, Fax:<a href="tel:%2B55%2011%203091-4714" value="+551130914714" target="_blank">+55 11 3091-4714</a> #<br>
# <a href="http://jonatandupont.weebly.com/index.html" target="_blank">http://jonatandupont.weebly.com/index.html</a> #<br>#------------------------------------------------------------------------------#<br>
</div>
<br></div></div>_______________________________________________
R-br mailing list
<a href="mailto:R-br@listas.c3sl.ufpr.br" target="_blank">R-br@listas.c3sl.ufpr.br</a>
<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>
Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forne�a c�digo m�nimo reproduz�vel.</div> </div></div>
<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>
Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.<br></blockquote></div><br><br clear="all"><div><br></div>-- <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>
</div>