<div dir="ltr"><div class="gmail_default" style=""><div class="gmail_default" style="font-family:verdana,sans-serif;font-size:12.8px">Olá eu gostaria de colocar cores (com a barra de cores ao lado da figura) no meu mapa ao invés de usar os contornos com a função "contour";<br>abaixo o link da figura que consegui gerar através dessa função.</div><div class="gmail_default" style="font-family:arial,sans-serif;font-size:12.8px"><br></div><div class="gmail_default" style="font-family:arial,sans-serif;font-size:12.8px"><font face="verdana, sans-serif"><a href="https://www.dropbox.com/s/xi622mbpqg7tgch/campo_medio_O3.png?dl=0" target="_blank">https://www.dropbox.com/s/xi622mbpqg7tgch/campo_medio_O3.png?dl=0</a></font></div><div class="gmail_default" style="font-family:arial,sans-serif;font-size:12.8px"><br></div><div class="gmail_default" style=""><span style="font-size:12.8px">abaixo o script que gerou o grafico com os contornos.</span><br><br><div class="gmail_default" style=""><br></div><div class="gmail_default" style=""><span style="font-size:12.8px"><br></span></div><div class="gmail_default" style=""><span style="font-size:12.8px"># Carregando biblioteca para manipular arquivos netCDF</span></div><div class="gmail_default" style=""><span style="font-size:12.8px"><br></span></div><div class="gmail_default" style=""><span style="font-size:12.8px">library(maps)</span></div><div class="gmail_default" style=""><span style="font-size:12.8px">library(ncdf4) </span></div><div class="gmail_default" style=""><span style="font-size:12.8px">  </span></div><div class="gmail_default" style=""><span style="font-size:12.8px">#==========================================================================</span></div><div class="gmail_default" style=""><br></div><div class="gmail_default" style=""><span style="font-size:12.8px"><br></span></div><div class="gmail_default" style=""><span style="font-size:12.8px">#BIBLIOTECA "ncdf4".</span></div><div class="gmail_default" style=""><span style="font-size:12.8px"># PARA ESTE EXEMPLO SÃO USADOS RECURSOS PARA ABRIR A BIBLIOTECA )</span></div><div class="gmail_default" style=""><span style="font-size:12.8px">#POIS RNetCDF, ncdf e ncdf4 APRESENTAM COMANDOS DIFERENTES PARA ABRIRMOS AS VARIÁVEIS</span></div><div class="gmail_default" style=""><span style="font-size:12.8px"><br></span></div><div class="gmail_default" style=""><span style="font-size:12.8px"><br></span></div><div class="gmail_default" style=""><span style="font-size:12.8px"> dados <- nc_open('<a href="http://teste.nc">teste.nc</a>')  </span></div><div class="gmail_default" style=""><span style="font-size:12.8px"># lendo coordenadas espaço-temporal</span></div><div class="gmail_default" style=""><span style="font-size:12.8px">lat <- ncvar_get( dados, 'lat' )</span></div><div class="gmail_default" style=""><span style="font-size:12.8px">lon <- ncvar_get( dados, 'lon' )</span></div><div class="gmail_default" style=""><span style="font-size:12.8px">time <- ncvar_get( dados, 'time' )</span></div><div class="gmail_default" style=""><span style="font-size:12.8px"><br></span></div><div class="gmail_default" style=""><span style="font-size:12.8px">#=======================================================================================================</span></div><div class="gmail_default" style=""><span style="font-size:12.8px"><br></span></div><div class="gmail_default" style=""><span style="font-size:12.8px"># lendo dados coluna total de Ozônio</span></div><div class="gmail_default" style=""><span style="font-size:12.8px">ColumnAmountO3 <- ncvar_get( dados, 'ColumnAmountO3' )</span></div><div class="gmail_default" style=""><span style="font-size:12.8px"><br></span></div><div class="gmail_default" style=""><span style="font-size:12.8px"># dimensoes da variavel ColumnAmountO3</span></div><div class="gmail_default" style=""><span style="font-size:12.8px">dims_ColumnAmountO3 <- dim(ColumnAmountO3)</span></div><div class="gmail_default" style=""><span style="font-size:12.8px"><br></span></div><div class="gmail_default" style=""><span style="font-size:12.8px"># tornando o arranjo 3D (ColumnAmountO3) em um 2D, organizado em ptos de grade X tempo</span></div><div class="gmail_default" style=""><span style="font-size:12.8px"><br></span></div><div class="gmail_default" style=""><span style="font-size:12.8px">dim(ColumnAmountO3) <- c( dims_ColumnAmountO3[1]*dims_ColumnAmountO3[2], dims_ColumnAmountO3[3] )</span></div><div class="gmail_default" style=""><span style="font-size:12.8px"><br></span></div><div class="gmail_default" style=""><span style="font-size:12.8px"># calculando a média e retornado-a em 2D</span></div><div class="gmail_default" style=""><span style="font-size:12.8px"><br></span></div><div class="gmail_default" style=""><span style="font-size:12.8px">media_ColumnAmountO3 <- rowMeans( ColumnAmountO3)</span></div><div class="gmail_default" style=""><span style="font-size:12.8px">dim(media_ColumnAmountO3) <- c( dims_ColumnAmountO3[1], dims_ColumnAmountO3[2] )</span></div><div class="gmail_default" style=""><span style="font-size:12.8px"><br></span></div><div class="gmail_default" style=""><span style="font-size:12.8px">#==========================================================================================================</span></div><div class="gmail_default" style=""><span style="font-size:12.8px"># longitude varia de 0 a 360, convertendo para -180 a 180, essa conversão é feita para plotagem sobre o mapa</span></div><div class="gmail_default" style=""><span style="font-size:12.8px">for (i in 1:dim(lon)) { if (lon[i]>180) { lon[i] <- lon[i]-360 } }</span></div><div class="gmail_default" style=""><span style="font-size:12.8px"><br></span></div><div class="gmail_default" style=""><span style="font-size:12.8px"><br></span></div><div class="gmail_default" style=""><span style="font-size:12.8px"># criando arquivo PNG que receberá o campo com o mapa</span></div><div class="gmail_default" style=""><span style="font-size:12.8px"><br></span></div><div class="gmail_default" style=""><span style="font-size:12.8px">#png( filename="campo_medio_O3_jan2005.png",width=600,height=800 )</span></div><div class="gmail_default" style=""><span style="font-size:12.8px"><br></span></div><div class="gmail_default" style=""><span style="font-size:12.8px"># plotando mapa da America do Sul</span></div><div class="gmail_default" style=""><span style="font-size:12.8px">map( xlim=c(-100,-10), ylim=c(-60,10) )</span></div><div class="gmail_default" style=""><span style="font-size:12.8px">map.axes()              # plotando eixos</span></div><div class="gmail_default" style=""><span style="font-size:12.8px">title( main="Campo médio de ozonio janeiros" )   # título do gráfico</span></div><div class="gmail_default" style=""><span style="font-size:12.8px"><br></span></div><div class="gmail_default" style=""><span style="font-size:12.8px"># definindo intervalo de 5 Dobson Units (DU)</span></div><div class="gmail_default" style=""><span style="font-size:12.8px">intervalos = seq( trunc(min(ColumnAmountO3)), trunc(max(ColumnAmountO3)), 5 )</span></div><div class="gmail_default" style=""><span style="font-size:12.8px"><br></span></div><div class="gmail_default" style=""><span style="font-size:12.8px"># adicionando campo de coluna de ozonio</span></div><div class="gmail_default" style=""><span style="font-size:12.8px"><br></span></div><div class="gmail_default" style=""><span style="font-size:12.8px"><br></span></div><div class="gmail_default" style=""><span style="font-size:12.8px">contour( sort(lon), lat, media_ColumnAmountO3[ order(lon), ], add=T, levels=intervalos, lwd=2, labcex=1.3, col="black" )</span></div><div class="gmail_default" style=""><span style="font-size:12.8px"><br></span></div><div class="gmail_default" style=""><span style="font-size:12.8px"><br></span></div><div class="gmail_default" style=""><span style="font-size:12.8px"># fechando arquivo PNG</span></div><div class="gmail_default" style=""><span style="font-size:12.8px">#dev.off() </span></div><div class="gmail_default" style=""><span style="font-size:12.8px"><br></span></div><div class="gmail_default" style=""><span style="font-size:12.8px"><br></span></div><div class="gmail_default" style=""><span style="font-size:12.8px">obrigado</span></div></div></div><div><div class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><font face="tahoma, sans-serif">____________________________________________________________________________</font></div><font face="tahoma, sans-serif"><div dir="ltr"><font face="tahoma, sans-serif"><br></font></div>MATEUS DIAS NUNES</font><div><div><font face="tahoma, sans-serif">MESTRANDO DO PROGRAMA DE PÓS-GRADUAÇÃO EM METEOROLOGIA - PPGMET</font></div><div><span style="font-family:tahoma,sans-serif">UNIVERSIDADE FEDERAL DE PELOTAS - UFPEL</span><font face="tahoma, sans-serif"><br></font></div><div><font face="tahoma, sans-serif"><span style="background-color:rgb(255,255,255)"><span style="color:rgb(0,0,0)">TELEFONE: +55 (53) 81125154           </span>  </span><br>____________________________________________<span style="background-color:rgb(255,255,255)"><span></span></span>________________________________</font></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div>
</div>