<div dir="ltr"><div class="gmail_default" style="font-family:verdana,sans-serif"><span style="font-family:arial,sans-serif;font-size:12.8px">Olá Mauro Sznelwar,<br></span><br>obrigado pelas sujestões. No assunto 27 eu respondi que obtive êxito com as Sujestões propostas pelo Thiago Santos. Mas mesmo assim muito obrigado.</div><div class="gmail_default" style="font-family:verdana,sans-serif"><br></div><div class="gmail_default" style="font-family:verdana,sans-serif">abs,</div></div><div class="gmail_extra"><br clear="all"><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>
<br><div class="gmail_quote">Em 24 de fevereiro de 2016 12:00, <span dir="ltr"><<a href="mailto:r-br-request@listas.c3sl.ufpr.br" target="_blank">r-br-request@listas.c3sl.ufpr.br</a>></span> escreveu:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Enviar submissões para a lista de discussão R-br para<br>
<a href="mailto:r-br@listas.c3sl.ufpr.br">r-br@listas.c3sl.ufpr.br</a><br>
<br>
Para se cadastrar ou descadastrar via WWW, visite o endereço<br>
<a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" rel="noreferrer" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>
ou, via email, envie uma mensagem com a palavra 'help' no assunto ou<br>
corpo da mensagem para<br>
<a href="mailto:r-br-request@listas.c3sl.ufpr.br">r-br-request@listas.c3sl.ufpr.br</a><br>
<br>
Você poderá entrar em contato com a pessoa que gerencia a lista pelo<br>
endereço<br>
<a href="mailto:r-br-owner@listas.c3sl.ufpr.br">r-br-owner@listas.c3sl.ufpr.br</a><br>
<br>
Quando responder, por favor edite sua linha Assunto assim ela será<br>
mais específica que "Re: Contents of R-br digest..."<br>
<br>
<br>
Tópicos de Hoje:<br>
<br>
1. RES: Cores em mapas (Mauro Sznelwar)<br>
2. Artigo que não tem no periodicos capes<br>
(Pedro Emmanuel Alvarenga Americano do Brasil)<br>
3. Re: Artigo que não tem no periodicos capes (Marcelo Laia)<br>
<br>
<br>
----------------------------------------------------------------------<br>
<br>
Message: 1<br>
Date: Wed, 24 Feb 2016 05:45:58 -0300<br>
From: "Mauro Sznelwar" <<a href="mailto:sznelwar@uol.com.br">sznelwar@uol.com.br</a>><br>
To: <<a href="mailto:r-br@listas.c3sl.ufpr.br">r-br@listas.c3sl.ufpr.br</a>>, "'Thiago V. dos Santos'"<br>
<<a href="mailto:thi_veloso@yahoo.com.br">thi_veloso@yahoo.com.br</a>><br>
Subject: [R-br] RES: Cores em mapas<br>
Message-ID: <004f01d16edf$c8d41680$5a7c4380$@<a href="http://com.br" rel="noreferrer" target="_blank">com.br</a>><br>
Content-Type: text/plain; charset="UTF-8"<br>
<br>
Teria o arquivo <a href="http://teste.nc" rel="noreferrer" target="_blank">teste.nc</a> para rodar?<br>
<br>
<br>
<br>
filled.contour( sort(lon), lat, media_ColumnAmountO3[ order(lon), ])<br>
map( xlim=c(-100,-10), ylim=c(-60,10), add=T)<br>
<br>
<br>
Ou então eu deixo aqui algumas sugestões usando o levelplot (vai precisar baixar o shapefile da AS em anexo):<br>
<br>
----------------------<br>
library(raster)<br>
library(rasterVis)<br>
<br>
# Abre Arquivo<br>
r <- brick("~/Downloads/<a href="http://teste.nc" rel="noreferrer" target="_blank">teste.nc</a>")<br>
<br>
# Tira média dos dois mapas<br>
r1 <- calc(r, fun=mean)<br>
<br>
# Uma das vantagens de usar raster/rasterVis: depois de duas linhas de código, já estamos plotando os dados<br>
levelplot(r1, margin=F, contour=T)<br>
<br>
# Descompacta e abre contorno da América do Sul para adicionar ao plot<br>
#(precisa baixar arquivo em anexo ou do link <a href="http://thematicmapping.org/downloads/TM_WORLD_BORDERS_SIMPL-0.3.zip" rel="noreferrer" target="_blank">http://thematicmapping.org/downloads/TM_WORLD_BORDERS_SIMPL-0.3.zip</a>)<br>
<br>
unzip("~/Downloads/TM_WORLD_BORDERS_SIMPL-0.3.zip")<br>
mapaSHP <- shapefile('TM_WORLD_BORDERS_SIMPL-0.3.shp')<br>
<br>
# Plota com contorno da America do Sul<br>
levelplot(r1, margin=F, contour=T, main="Campo médio de Ozônio Novembros") + layer(sp.lines(mapaSHP, lwd=0.7, col='black'))<br>
<br>
# Quer mudar a cor? Basta mudar as configurações de 'par.settings'<br>
levelplot(r1, margin=F, contour=T, main="Campo médio de Ozônio Novembros", par.settings=RdBuTheme) + layer(sp.lines(mapaSHP, lwd=0.7, col='black'))<br>
<br>
# Também é possível criar a sua própria paleta de cores, usando o pacote Rcolorbrewer<br>
# Veja aconselhamento para cores de mapa aqui: <a href="http://colorbrewer2.org/" rel="noreferrer" target="_blank">http://colorbrewer2.org/</a><br>
# Vou criar uma escala sequencial, com 5 tons de azul<br>
# Você pode ver todas as paletas disponíveis no pacote digitando display.brewer.all()<br>
my.Theme=rasterTheme(region=brewer.pal('Blues', n=5))<br>
<br>
# Aqui eu personalizo a escala usando os valores do raster e um numero maior de quebras<br>
intervalos <- seq(trunc(min(values(r1))), trunc(max(values(r1))), length.out=length(my.Theme$regions$col)-1)<br>
<br>
# Plota<br>
p <- levelplot(r1, margin=F, contour=T, par.settings=my.Theme, at=intervalos, main="Campo médio de Ozônio Novembros")<br>
p + layer(sp.lines(mapaSHP, lwd=0.7, col='black'))<br>
<br>
# Para salvar, eu recomendo copiar a imagem direto do display para um arquivo pdf<br>
dev.copy2pdf(file="~/Desktop/ozone.pdf", width = 11, height = 8.5)<br>
<br>
# Depois fica fácil convertê-la para png em altíssima resolução (600dpi) usando o ImageMagick (que precisa estar instalado):<br>
# Assim a sua imagem está de acordo com as normas de publicação das principais revistas científicas<br>
system("convert -density 600 ~/Desktop/ozone.pdf ~/Desktop/ozone.png")<br>
<br>
<br>
Veja o resultado final:<br>
<a href="https://dl.dropboxusercontent.com/u/27700634/ozone.png" rel="noreferrer" target="_blank">https://dl.dropboxusercontent.com/u/27700634/ozone.png</a><br>
<br>
Greetings,<br>
-- Thiago V. dos Santos<br>
<br>
PhD student<br>
Land and Atmospheric Science<br>
University of Minnesota<br>
<br>
<br>
On Friday, February 19, 2016 1:18 PM, Mateus Dias Nunes <<a href="mailto:nunes.mateusdias@gmail.com">nunes.mateusdias@gmail.com</a>> wrote:<br>
<br>
<br>
<br>
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.<br>
<br>
<a href="https://www.dropbox.com/s/xi622mbpqg7tgch/campo_medio_O3.png?dl=0" rel="noreferrer" target="_blank">https://www.dropbox.com/s/xi622mbpqg7tgch/campo_medio_O3.png?dl=0</a><br>
<br>
abaixo o script que gerou o grafico com os contornos.<br>
<br>
<br>
<br>
<br>
# Carregando biblioteca para manipular arquivos netCDF<br>
<br>
library(maps)<br>
library(ncdf4)<br>
<br>
#==========================================================================<br>
<br>
<br>
#BIBLIOTECA "ncdf4".<br>
# PARA ESTE EXEMPLO SÃO USADOS RECURSOS PARA ABRIR A BIBLIOTECA )<br>
#POIS RNetCDF, ncdf e ncdf4 APRESENTAM COMANDOS DIFERENTES PARA ABRIRMOS AS VARIÁVEIS<br>
<br>
<br>
dados <- nc_open('<a href="http://teste.nc" rel="noreferrer" target="_blank">teste.nc</a>')<br>
# lendo coordenadas espaço-temporal<br>
lat <- ncvar_get( dados, 'lat' )<br>
lon <- ncvar_get( dados, 'lon' )<br>
time <- ncvar_get( dados, 'time' )<br>
<br>
#=======================================================================================================<br>
<br>
# lendo dados coluna total de Ozônio<br>
ColumnAmountO3 <- ncvar_get( dados, 'ColumnAmountO3' )<br>
<br>
# dimensoes da variavel ColumnAmountO3<br>
dims_ColumnAmountO3 <- dim(ColumnAmountO3)<br>
<br>
# tornando o arranjo 3D (ColumnAmountO3) em um 2D, organizado em ptos de grade X tempo<br>
<br>
dim(ColumnAmountO3) <- c( dims_ColumnAmountO3[1]*dims_ColumnAmountO3[2], dims_ColumnAmountO3[3] )<br>
<br>
# calculando a média e retornado-a em 2D<br>
<br>
media_ColumnAmountO3 <- rowMeans( ColumnAmountO3)<br>
dim(media_ColumnAmountO3) <- c( dims_ColumnAmountO3[1], dims_ColumnAmountO3[2] )<br>
<br>
#==========================================================================================================<br>
# longitude varia de 0 a 360, convertendo para -180 a 180, essa conversão é feita para plotagem sobre o mapa<br>
for (i in 1:dim(lon)) { if (lon[i]>180) { lon[i] <- lon[i]-360 } }<br>
<br>
<br>
# criando arquivo PNG que receberá o campo com o mapa<br>
<br>
#png( filename="campo_medio_O3_jan2005.png",width=600,height=800 )<br>
<br>
# plotando mapa da America do Sul<br>
map( xlim=c(-100,-10), ylim=c(-60,10) )<br>
map.axes() # plotando eixos<br>
title( main="Campo médio de ozonio janeiros" ) # título do gráfico<br>
<br>
# definindo intervalo de 5 Dobson Units (DU)<br>
intervalos = seq( trunc(min(ColumnAmountO3)), trunc(max(ColumnAmountO3)), 5 )<br>
<br>
# adicionando campo de coluna de ozonio<br>
<br>
<br>
contour( sort(lon), lat, media_ColumnAmountO3[ order(lon), ], add=T, levels=intervalos, lwd=2, labcex=1.3, col="black" )<br>
<br>
<br>
# fechando arquivo PNG<br>
#dev.off()<br>
<br>
<br>
obrigado<br>
____________________________________________________________________________<br>
MATEUS DIAS NUNES<br>
MESTRANDO DO PROGRAMA DE PÓS-GRADUAÇÃO EM METEOROLOGIA - PPGMET<br>
UNIVERSIDADE FEDERAL DE PELOTAS - UFPEL<br>
<br>
TELEFONE: +55 (53) 81125154<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" rel="noreferrer" 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" rel="noreferrer" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forne� c�igo m�imo reproduz�el.<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" rel="noreferrer" 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" rel="noreferrer" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e fornea cdigo mnimo reproduzvel.<br>
<br>
<br>
---<br>
Este email foi escaneado pelo Avast antivírus.<br>
<a href="https://www.avast.com/antivirus" rel="noreferrer" target="_blank">https://www.avast.com/antivirus</a><br>
<br>
<br>
<br>
------------------------------<br>
<br>
Message: 2<br>
Date: Wed, 24 Feb 2016 08:31:03 -0300<br>
From: Pedro Emmanuel Alvarenga Americano do Brasil<br>
<<a href="mailto:emmanuel.brasil@gmail.com">emmanuel.brasil@gmail.com</a>><br>
To: R BR lista <<a href="mailto:r-br@listas.c3sl.ufpr.br">r-br@listas.c3sl.ufpr.br</a>><br>
Subject: [R-br] Artigo que não tem no periodicos capes<br>
Message-ID:<br>
<<a href="mailto:CAFfGvy%2BdQnZQ1XKJX3zUR56oH3ouV118FwW1Pq%2B6KQ73%2BcK84A@mail.gmail.com">CAFfGvy+dQnZQ1XKJX3zUR56oH3ouV118FwW1Pq+6KQ73+cK84A@mail.gmail.com</a>><br>
Content-Type: text/plain; charset="utf-8"<br>
<br>
Amigos de R,<br>
<br>
Será que alguem tem na mão que possa me enviar o pdf desse artigo aqui:<br>
<br>
<a href="http://onlinelibrary.wiley.com/doi/10.1002/sim.4780141909/full" rel="noreferrer" target="_blank">http://onlinelibrary.wiley.com/doi/10.1002/sim.4780141909/full</a><br>
<br>
Grato,<br>
<br>
Pedro Brasil<br>
-------------- Próxima Parte ----------<br>
Um anexo em HTML foi limpo...<br>
URL: <<a href="http://listas.inf.ufpr.br/pipermail/r-br/attachments/20160224/c6a4909e/attachment-0001.html" rel="noreferrer" target="_blank">http://listas.inf.ufpr.br/pipermail/r-br/attachments/20160224/c6a4909e/attachment-0001.html</a>><br>
<br>
------------------------------<br>
<br>
Message: 3<br>
Date: Wed, 24 Feb 2016 08:33:57 -0300<br>
From: Marcelo Laia <<a href="mailto:marcelolaia@gmail.com">marcelolaia@gmail.com</a>><br>
To: r-br <<a href="mailto:r-br@listas.c3sl.ufpr.br">r-br@listas.c3sl.ufpr.br</a>><br>
Subject: Re: [R-br] Artigo que não tem no periodicos capes<br>
Message-ID:<br>
<<a href="mailto:CAEEYVUD8wuMn%2BN9ymLOd8_5wiP7uLFbxhx8WBUkExhRaPHc99g@mail.gmail.com">CAEEYVUD8wuMn+N9ymLOd8_5wiP7uLFbxhx8WBUkExhRaPHc99g@mail.gmail.com</a>><br>
Content-Type: text/plain; charset="utf-8"<br>
<br>
<a href="http://sci-hub.io/" rel="noreferrer" target="_blank">http://sci-hub.io/</a><br>
Em 24 de fev de 2016 08:31, "Pedro Emmanuel Alvarenga Americano do Brasil" <<br>
<a href="mailto:emmanuel.brasil@gmail.com">emmanuel.brasil@gmail.com</a>> escreveu:<br>
<br>
> Amigos de R,<br>
><br>
> Será que alguem tem na mão que possa me enviar o pdf desse artigo aqui:<br>
><br>
> <a href="http://onlinelibrary.wiley.com/doi/10.1002/sim.4780141909/full" rel="noreferrer" target="_blank">http://onlinelibrary.wiley.com/doi/10.1002/sim.4780141909/full</a><br>
><br>
> Grato,<br>
><br>
> Pedro Brasil<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" rel="noreferrer" 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" rel="noreferrer" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça<br>
> código mínimo reproduzível.<br>
><br>
-------------- Próxima Parte ----------<br>
Um anexo em HTML foi limpo...<br>
URL: <<a href="http://listas.inf.ufpr.br/pipermail/r-br/attachments/20160224/11e2033d/attachment-0001.html" rel="noreferrer" target="_blank">http://listas.inf.ufpr.br/pipermail/r-br/attachments/20160224/11e2033d/attachment-0001.html</a>><br>
<br>
------------------------------<br>
<br>
Subject: Legenda do Digest<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" rel="noreferrer" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>
<br>
<br>
------------------------------<br>
<br>
Fim da Digest R-br, volume 62, assunto 29<br>
*****************************************<br>
</blockquote></div><br></div>