<div dir="ltr"><div>## usando rotate</div><div>require(raster)</div><div>b <- brick("~/Downloads/SST_81_14.nc", <b>stopIfNotEqualSpaced=FALSE</b>)</div><div>b</div><div>extent(b)</div><div>## fix long</div><div>extent(b) <- c(0, 360, -51, 11)</div><div>b</div><div>plot(b)</div><div>## rotate</div><div>br <- rotate(b)</div><div>extent(br)</div><div>br <- setValues(br, getValues(b))</div><div>plot(br)</div></div><div class="gmail_extra"><br><div class="gmail_quote">2014-12-05 0:55 GMT-02:00 Humberto <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 lang="PT-BR" link="#0563C1" vlink="#954F72"><div><p class="MsoNormal">Boa noite pessoal,<u></u><u></u></p><p class="MsoNormal"><u></u> <u></u></p><p class="MsoNormal">Estou extraindo de um arquivo ncdf os valores de latitude (lat) longitude( lon)  e temperatura da superfície da agua do mar( SST), entretanto a longitude (lon) ela está na projeção 0 a 360 graus. Minha pergunta é como faço para converter ela para -180 a 180 graus? Segue o script que estou utilizando juntamente com o link dos dados.<u></u><u></u></p><p class="MsoNormal">Desde já agradeço a ajuda <u></u><u></u></p><p class="MsoNormal">Humberto <u></u><u></u></p><p class="MsoNormal"><u></u> <u></u></p><table border="0" cellspacing="0" cellpadding="0" width="450" style="width:337.5pt;border-collapse:collapse"><tbody><tr style="height:15.0pt"><td width="64" nowrap valign="bottom" style="width:48.0pt;padding:0in 3.5pt 0in 3.5pt;height:15.0pt"><p class="MsoNormal"><span style="color:black">          lon<u></u><u></u></span></p></td><td width="64" nowrap valign="bottom" style="width:48.0pt;padding:0in 3.5pt 0in 3.5pt;height:15.0pt"><p class="MsoNormal"><span style="color:black">             lat<u></u><u></u></span></p></td><td width="77" nowrap valign="bottom" style="width:58.0pt;padding:0in 3.5pt 0in 3.5pt;height:15.0pt"><p class="MsoNormal"><span style="color:black">     data<u></u><u></u></span></p></td><td width="64" nowrap valign="bottom" style="width:48.0pt;padding:0in 3.5pt 0in 3.5pt;height:15.0pt"><p class="MsoNormal"><span style="color:black">          sst<u></u><u></u></span></p></td><td width="64" nowrap valign="bottom" style="width:48.0pt;padding:0in 3.5pt 0in 3.5pt;height:15.0pt"></td><td width="117" nowrap valign="bottom" style="width:87.5pt;padding:0in 3.5pt 0in 3.5pt;height:15.0pt"><p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Times New Roman",serif">               Lon como deveria ser<u></u><u></u></span></p></td></tr><tr style="height:15.0pt"><td width="64" nowrap valign="bottom" style="width:48.0pt;padding:0in 3.5pt 0in 3.5pt;height:15.0pt"><p class="MsoNormal" align="right" style="text-align:right"><span style="color:black">299.5<u></u><u></u></span></p></td><td width="64" nowrap valign="bottom" style="width:48.0pt;padding:0in 3.5pt 0in 3.5pt;height:15.0pt"><p class="MsoNormal" align="right" style="text-align:right"><span style="color:black">10.5<u></u><u></u></span></p></td><td width="77" nowrap valign="bottom" style="width:58.0pt;padding:0in 3.5pt 0in 3.5pt;height:15.0pt"><p class="MsoNormal" align="right" style="text-align:right"><span style="color:black">11/30/1981<u></u><u></u></span></p></td><td width="64" nowrap valign="bottom" style="width:48.0pt;padding:0in 3.5pt 0in 3.5pt;height:15.0pt"><p class="MsoNormal" align="right" style="text-align:right"><span style="color:black">27.85<u></u><u></u></span></p></td><td width="64" nowrap valign="bottom" style="width:48.0pt;padding:0in 3.5pt 0in 3.5pt;height:15.0pt"><p class="MsoNormal" align="right" style="text-align:right"><span style="color:black"><u></u> <u></u></span></p></td><td width="117" nowrap valign="bottom" style="width:87.5pt;padding:0in 3.5pt 0in 3.5pt;height:15.0pt"><p class="MsoNormal" align="right" style="text-align:right"><span style="color:black">-60.000<u></u><u></u></span></p></td></tr><tr style="height:15.0pt"><td width="64" nowrap valign="bottom" style="width:48.0pt;padding:0in 3.5pt 0in 3.5pt;height:15.0pt"><p class="MsoNormal" align="right" style="text-align:right"><span style="color:black">300.5<u></u><u></u></span></p></td><td width="64" nowrap valign="bottom" style="width:48.0pt;padding:0in 3.5pt 0in 3.5pt;height:15.0pt"><p class="MsoNormal" align="right" style="text-align:right"><span style="color:black">10.5<u></u><u></u></span></p></td><td width="77" nowrap valign="bottom" style="width:58.0pt;padding:0in 3.5pt 0in 3.5pt;height:15.0pt"><p class="MsoNormal" align="right" style="text-align:right"><span style="color:black">11/30/1981<u></u><u></u></span></p></td><td width="64" nowrap valign="bottom" style="width:48.0pt;padding:0in 3.5pt 0in 3.5pt;height:15.0pt"><p class="MsoNormal" align="right" style="text-align:right"><span style="color:black">27.96<u></u><u></u></span></p></td><td width="64" nowrap valign="bottom" style="width:48.0pt;padding:0in 3.5pt 0in 3.5pt;height:15.0pt"><p class="MsoNormal" align="right" style="text-align:right"><span style="color:black"><u></u> <u></u></span></p></td><td width="117" nowrap valign="bottom" style="width:87.5pt;padding:0in 3.5pt 0in 3.5pt;height:15.0pt"><p class="MsoNormal" align="right" style="text-align:right"><span style="color:black">-59.500<u></u><u></u></span></p></td></tr><tr style="height:15.0pt"><td width="64" nowrap valign="bottom" style="width:48.0pt;padding:0in 3.5pt 0in 3.5pt;height:15.0pt"><p class="MsoNormal" align="right" style="text-align:right"><span style="color:black">301.5<u></u><u></u></span></p></td><td width="64" nowrap valign="bottom" style="width:48.0pt;padding:0in 3.5pt 0in 3.5pt;height:15.0pt"><p class="MsoNormal" align="right" style="text-align:right"><span style="color:black">10.5<u></u><u></u></span></p></td><td width="77" nowrap valign="bottom" style="width:58.0pt;padding:0in 3.5pt 0in 3.5pt;height:15.0pt"><p class="MsoNormal" align="right" style="text-align:right"><span style="color:black">11/30/1981<u></u><u></u></span></p></td><td width="64" nowrap valign="bottom" style="width:48.0pt;padding:0in 3.5pt 0in 3.5pt;height:15.0pt"><p class="MsoNormal" align="right" style="text-align:right"><span style="color:black">28.02<u></u><u></u></span></p></td><td width="64" nowrap valign="bottom" style="width:48.0pt;padding:0in 3.5pt 0in 3.5pt;height:15.0pt"><p class="MsoNormal" align="right" style="text-align:right"><span style="color:black"><u></u> <u></u></span></p></td><td width="117" nowrap valign="bottom" style="width:87.5pt;padding:0in 3.5pt 0in 3.5pt;height:15.0pt"><p class="MsoNormal" align="right" style="text-align:right"><span style="color:black">-58.500<u></u><u></u></span></p></td></tr><tr style="height:15.0pt"><td width="64" nowrap valign="bottom" style="width:48.0pt;padding:0in 3.5pt 0in 3.5pt;height:15.0pt"><p class="MsoNormal" align="right" style="text-align:right"><span style="color:black">302.5<u></u><u></u></span></p></td><td width="64" nowrap valign="bottom" style="width:48.0pt;padding:0in 3.5pt 0in 3.5pt;height:15.0pt"><p class="MsoNormal" align="right" style="text-align:right"><span style="color:black">10.5<u></u><u></u></span></p></td><td width="77" nowrap valign="bottom" style="width:58.0pt;padding:0in 3.5pt 0in 3.5pt;height:15.0pt"><p class="MsoNormal" align="right" style="text-align:right"><span style="color:black">11/30/1981<u></u><u></u></span></p></td><td width="64" nowrap valign="bottom" style="width:48.0pt;padding:0in 3.5pt 0in 3.5pt;height:15.0pt"><p class="MsoNormal" align="right" style="text-align:right"><span style="color:black">28.02<u></u><u></u></span></p></td><td width="64" nowrap valign="bottom" style="width:48.0pt;padding:0in 3.5pt 0in 3.5pt;height:15.0pt"><p class="MsoNormal" align="right" style="text-align:right"><span style="color:black"><u></u> <u></u></span></p></td><td width="117" nowrap valign="bottom" style="width:87.5pt;padding:0in 3.5pt 0in 3.5pt;height:15.0pt"><p class="MsoNormal" align="right" style="text-align:right"><span style="color:black">-57.500<u></u><u></u></span></p></td></tr><tr style="height:15.0pt"><td width="64" nowrap valign="bottom" style="width:48.0pt;padding:0in 3.5pt 0in 3.5pt;height:15.0pt"><p class="MsoNormal" align="right" style="text-align:right"><span style="color:black">303.5<u></u><u></u></span></p></td><td width="64" nowrap valign="bottom" style="width:48.0pt;padding:0in 3.5pt 0in 3.5pt;height:15.0pt"><p class="MsoNormal" align="right" style="text-align:right"><span style="color:black">10.5<u></u><u></u></span></p></td><td width="77" nowrap valign="bottom" style="width:58.0pt;padding:0in 3.5pt 0in 3.5pt;height:15.0pt"><p class="MsoNormal" align="right" style="text-align:right"><span style="color:black">11/30/1981<u></u><u></u></span></p></td><td width="64" nowrap valign="bottom" style="width:48.0pt;padding:0in 3.5pt 0in 3.5pt;height:15.0pt"><p class="MsoNormal" align="right" style="text-align:right"><span style="color:black">27.97<u></u><u></u></span></p></td><td width="64" nowrap valign="bottom" style="width:48.0pt;padding:0in 3.5pt 0in 3.5pt;height:15.0pt"><p class="MsoNormal" align="right" style="text-align:right"><span style="color:black"><u></u> <u></u></span></p></td><td width="117" nowrap valign="bottom" style="width:87.5pt;padding:0in 3.5pt 0in 3.5pt;height:15.0pt"><p class="MsoNormal" align="right" style="text-align:right"><span style="color:black">-56.500<u></u><u></u></span></p></td></tr><tr style="height:15.0pt"><td width="64" nowrap valign="bottom" style="width:48.0pt;padding:0in 3.5pt 0in 3.5pt;height:15.0pt"><p class="MsoNormal" align="right" style="text-align:right"><span style="color:black">304.5<u></u><u></u></span></p></td><td width="64" nowrap valign="bottom" style="width:48.0pt;padding:0in 3.5pt 0in 3.5pt;height:15.0pt"><p class="MsoNormal" align="right" style="text-align:right"><span style="color:black">10.5<u></u><u></u></span></p></td><td width="77" nowrap valign="bottom" style="width:58.0pt;padding:0in 3.5pt 0in 3.5pt;height:15.0pt"><p class="MsoNormal" align="right" style="text-align:right"><span style="color:black">11/30/1981<u></u><u></u></span></p></td><td width="64" nowrap valign="bottom" style="width:48.0pt;padding:0in 3.5pt 0in 3.5pt;height:15.0pt"><p class="MsoNormal" align="right" style="text-align:right"><span style="color:black">27.84<u></u><u></u></span></p></td><td width="64" nowrap valign="bottom" style="width:48.0pt;padding:0in 3.5pt 0in 3.5pt;height:15.0pt"><p class="MsoNormal" align="right" style="text-align:right"><span style="color:black"><u></u> <u></u></span></p></td><td width="117" nowrap valign="bottom" style="width:87.5pt;padding:0in 3.5pt 0in 3.5pt;height:15.0pt"><p class="MsoNormal" align="right" style="text-align:right"><span style="color:black">-55.500<u></u><u></u></span></p></td></tr><tr style="height:15.0pt"><td width="64" nowrap valign="bottom" style="width:48.0pt;padding:0in 3.5pt 0in 3.5pt;height:15.0pt"><p class="MsoNormal" align="right" style="text-align:right"><span style="color:black">305.5<u></u><u></u></span></p></td><td width="64" nowrap valign="bottom" style="width:48.0pt;padding:0in 3.5pt 0in 3.5pt;height:15.0pt"><p class="MsoNormal" align="right" style="text-align:right"><span style="color:black">10.5<u></u><u></u></span></p></td><td width="77" nowrap valign="bottom" style="width:58.0pt;padding:0in 3.5pt 0in 3.5pt;height:15.0pt"><p class="MsoNormal" align="right" style="text-align:right"><span style="color:black">11/30/1981<u></u><u></u></span></p></td><td width="64" nowrap valign="bottom" style="width:48.0pt;padding:0in 3.5pt 0in 3.5pt;height:15.0pt"><p class="MsoNormal" align="right" style="text-align:right"><span style="color:black">27.81<u></u><u></u></span></p></td><td width="64" nowrap valign="bottom" style="width:48.0pt;padding:0in 3.5pt 0in 3.5pt;height:15.0pt"><p class="MsoNormal" align="right" style="text-align:right"><span style="color:black"><u></u> <u></u></span></p></td><td width="117" nowrap valign="bottom" style="width:87.5pt;padding:0in 3.5pt 0in 3.5pt;height:15.0pt"><p class="MsoNormal" align="right" style="text-align:right"><span style="color:black">-54.500<u></u><u></u></span></p></td></tr><tr style="height:15.0pt"><td width="64" nowrap valign="bottom" style="width:48.0pt;padding:0in 3.5pt 0in 3.5pt;height:15.0pt"><p class="MsoNormal" align="right" style="text-align:right"><span style="color:black">306.5<u></u><u></u></span></p></td><td width="64" nowrap valign="bottom" style="width:48.0pt;padding:0in 3.5pt 0in 3.5pt;height:15.0pt"><p class="MsoNormal" align="right" style="text-align:right"><span style="color:black">10.5<u></u><u></u></span></p></td><td width="77" nowrap valign="bottom" style="width:58.0pt;padding:0in 3.5pt 0in 3.5pt;height:15.0pt"><p class="MsoNormal" align="right" style="text-align:right"><span style="color:black">11/30/1981<u></u><u></u></span></p></td><td width="64" nowrap valign="bottom" style="width:48.0pt;padding:0in 3.5pt 0in 3.5pt;height:15.0pt"><p class="MsoNormal" align="right" style="text-align:right"><span style="color:black">27.78<u></u><u></u></span></p></td><td width="64" nowrap valign="bottom" style="width:48.0pt;padding:0in 3.5pt 0in 3.5pt;height:15.0pt"><p class="MsoNormal" align="right" style="text-align:right"><span style="color:black"><u></u> <u></u></span></p></td><td width="117" nowrap valign="bottom" style="width:87.5pt;padding:0in 3.5pt 0in 3.5pt;height:15.0pt"><p class="MsoNormal" align="right" style="text-align:right"><span style="color:black">-53.500<u></u><u></u></span></p></td></tr><tr style="height:15.0pt"><td width="64" nowrap valign="bottom" style="width:48.0pt;padding:0in 3.5pt 0in 3.5pt;height:15.0pt"><p class="MsoNormal" align="right" style="text-align:right"><span style="color:black">307.5<u></u><u></u></span></p></td><td width="64" nowrap valign="bottom" style="width:48.0pt;padding:0in 3.5pt 0in 3.5pt;height:15.0pt"><p class="MsoNormal" align="right" style="text-align:right"><span style="color:black">10.5<u></u><u></u></span></p></td><td width="77" nowrap valign="bottom" style="width:58.0pt;padding:0in 3.5pt 0in 3.5pt;height:15.0pt"><p class="MsoNormal" align="right" style="text-align:right"><span style="color:black">11/30/1981<u></u><u></u></span></p></td><td width="64" nowrap valign="bottom" style="width:48.0pt;padding:0in 3.5pt 0in 3.5pt;height:15.0pt"><p class="MsoNormal" align="right" style="text-align:right"><span style="color:black">27.75<u></u><u></u></span></p></td><td width="64" nowrap valign="bottom" style="width:48.0pt;padding:0in 3.5pt 0in 3.5pt;height:15.0pt"><p class="MsoNormal" align="right" style="text-align:right"><span style="color:black"><u></u> <u></u></span></p></td><td width="117" nowrap valign="bottom" style="width:87.5pt;padding:0in 3.5pt 0in 3.5pt;height:15.0pt"><p class="MsoNormal" align="right" style="text-align:right"><span style="color:black">-52.500<u></u><u></u></span></p></td></tr></tbody></table><p class="MsoNormal"><u></u> <u></u></p><p class="MsoNormal"><u></u> <u></u></p><p class="MsoNormal">Link para acesso dos dados <a href="https://www.dropbox.com/s/xnab2g365mduykx/SST_81_14.nc?dl=0" target="_blank">https://www.dropbox.com/s/xnab2g365mduykx/SST_81_14.nc?dl=0</a><u></u><u></u></p><p class="MsoNormal"><u></u> <u></u></p><p class="MsoNormal">Script:<u></u><u></u></p><p class="MsoNormal"><span lang="EN-US">require(ncdf) #'ncdf' versao 1.6.6<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US">require(reshape)<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US">require(chron)<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p><p class="MsoNormal"><span lang="EN-US">###sst<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US">nc <- open.ncdf("SST_81_14.nc")<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US">nc<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p><p class="MsoNormal"><span lang="EN-US">lat <- get.var.ncdf(nc, "lat")<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US">lon <- get.var.ncdf(nc, "lon")<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US">tunits <- att.get.ncdf(nc,"time","units")<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p><p class="MsoNormal"><span lang="EN-US"># print tunits<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US">tunits$value<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p><p class="MsoNormal"><span lang="EN-US"># criar tempo "real" Times<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US">time<- get.var.ncdf(nc,"time")<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US">Time<-chron(time,origin=c(month=1, day=1, year=1900))<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p><p class="MsoNormal">## array da variavel sst<u></u><u></u></p><p class="MsoNormal">wm <- get.var.ncdf(nc, "sst")<u></u><u></u></p><p class="MsoNormal">str(wm)<u></u><u></u></p><p class="MsoNormal"><u></u> <u></u></p><p class="MsoNormal">## atribuindo dimnames conforme lon, lat e tempo<u></u><u></u></p><p class="MsoNormal"><span lang="EN-US">dimnames(wm)[[1]] <- as.character(lon)<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US">dimnames(wm)[[2]] <- as.character(lat)<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US">dimnames(wm)[[3]] <- as.character(Time)<u></u><u></u></span></p><p class="MsoNormal">str(wm)<u></u><u></u></p><p class="MsoNormal"><u></u> <u></u></p><p class="MsoNormal">## melt para reestruturação dos dados<u></u><u></u></p><p class="MsoNormal">mwm <- melt(wm)<u></u><u></u></p><p class="MsoNormal">head(mwm)<u></u><u></u></p></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><div class="gmail_signature"><div dir="ltr"><div><font style="font-family:arial">###############################################################</font><br><span style="font-family:arial;font-size:small">##  </span><span style="font-size:11px;text-align:center;background-color:rgb(239,239,239);color:rgb(150,150,150);font-family:Verdana,Arial,Helvetica,sans-serif;line-height:13px">Jônatan Dupont Tatsch</span></div><div><span style="font-family:arial;font-size:small">##  </span><span style="color:rgb(150,150,150);font-family:Verdana,Arial,Helvetica,sans-serif;font-size:10.909090995788574px;line-height:12.997159004211426px;text-align:center;background-color:rgb(239,239,239)">Professor do Departamento de Física</span><font face="arial" style="font-size:small"><br></font></div><div><span style="font-size:small;font-family:arial">##  </span><span style="color:rgb(150,150,150);font-family:Verdana,Arial,Helvetica,sans-serif;font-size:10.909090995788574px;line-height:12.997159004211426px;text-align:center;background-color:rgb(239,239,239)">Centro de Ciências Exatas e Naturais (CCNE)</span><br></div><div><span style="font-size:small;font-family:arial">##  </span><span style="color:rgb(150,150,150);font-family:Verdana,Arial,Helvetica,sans-serif;font-size:10.909090995788574px;line-height:12.997159004211426px;text-align:center;background-color:rgb(239,239,239)">Universidade Federal de Santa Maria</span></div><div><span style="font-family:arial;font-size:small">##  </span><span style="font-size:11px;text-align:center;background-color:rgb(239,239,239);color:rgb(150,150,150);font-family:Verdana,Arial,Helvetica,sans-serif;line-height:13px">Faixa de Camobi, Prédio 13 - Campus UFSM - </span><span style="font-size:11px;text-align:center;background-color:rgb(239,239,239);color:rgb(150,150,150);font-family:Verdana,Arial,Helvetica,sans-serif;line-height:13px">Santa Maria, RS, Brasil - 97105-900</span><br></div><div style="font-family:arial;font-size:small">##  <span style="font-size:11px;text-align:center;background-color:rgb(239,239,239);color:rgb(150,150,150);font-family:Verdana,Arial,Helvetica,sans-serif;line-height:13px">Telefone: +55(55)33012083</span><br></div><div style="font-family:arial;font-size:small">##  <span style="font-size:11px;text-align:center;background-color:rgb(239,239,239);color:rgb(150,150,150);font-family:Verdana,Arial,Helvetica,sans-serif;line-height:13px"><a href="http://www.ufsm.br/meteorologia" target="_blank">www.ufsm.br/meteorologia</a></span></div><div style="font-family:arial;font-size:small">###############################################################</div></div></div>
</div>