Olá Mauro Sznelwar, 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. abs, ____________________________________________________________________________ MATEUS DIAS NUNES MESTRANDO DO PROGRAMA DE PÓS-GRADUAÇÃO EM METEOROLOGIA - PPGMET UNIVERSIDADE FEDERAL DE PELOTAS - UFPEL TELEFONE: +55 (53) 81125154 ____________________________________________________________________________ Em 24 de fevereiro de 2016 12:00, <r-br-request@listas.c3sl.ufpr.br> escreveu:

Enviar submissões para a lista de discussão R-br para r-br@listas.c3sl.ufpr.br

Para se cadastrar ou descadastrar via WWW, visite o endereço https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br ou, via email, envie uma mensagem com a palavra 'help' no assunto ou corpo da mensagem para r-br-request@listas.c3sl.ufpr.br

Você poderá entrar em contato com a pessoa que gerencia a lista pelo endereço r-br-owner@listas.c3sl.ufpr.br

Quando responder, por favor edite sua linha Assunto assim ela será mais específica que "Re: Contents of R-br digest..."

Tópicos de Hoje:

1. RES: Cores em mapas (Mauro Sznelwar) 2. Artigo que não tem no periodicos capes (Pedro Emmanuel Alvarenga Americano do Brasil) 3. Re: Artigo que não tem no periodicos capes (Marcelo Laia)

----------------------------------------------------------------------

Message: 1 Date: Wed, 24 Feb 2016 05:45:58 -0300 From: "Mauro Sznelwar" <sznelwar@uol.com.br> To: <r-br@listas.c3sl.ufpr.br>, "'Thiago V. dos Santos'" <thi_veloso@yahoo.com.br> Subject: [R-br] RES: Cores em mapas Message-ID: <004f01d16edf$c8d41680$5a7c4380$@com.br> Content-Type: text/plain; charset="UTF-8"

Teria o arquivo teste.nc para rodar?

filled.contour( sort(lon), lat, media_ColumnAmountO3[ order(lon), ]) map( xlim=c(-100,-10), ylim=c(-60,10), add=T)

Ou então eu deixo aqui algumas sugestões usando o levelplot (vai precisar baixar o shapefile da AS em anexo):

---------------------- library(raster) library(rasterVis)

# Abre Arquivo r <- brick("~/Downloads/teste.nc")

# Tira média dos dois mapas r1 <- calc(r, fun=mean)

# Uma das vantagens de usar raster/rasterVis: depois de duas linhas de código, já estamos plotando os dados levelplot(r1, margin=F, contour=T)

# Descompacta e abre contorno da América do Sul para adicionar ao plot #(precisa baixar arquivo em anexo ou do link http://thematicmapping.org/downloads/TM_WORLD_BORDERS_SIMPL-0.3.zip)

unzip("~/Downloads/TM_WORLD_BORDERS_SIMPL-0.3.zip") mapaSHP <- shapefile('TM_WORLD_BORDERS_SIMPL-0.3.shp')

# Plota com contorno da America do Sul levelplot(r1, margin=F, contour=T, main="Campo médio de Ozônio Novembros") + layer(sp.lines(mapaSHP, lwd=0.7, col='black'))

# Quer mudar a cor? Basta mudar as configurações de 'par.settings' 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'))

# Também é possível criar a sua própria paleta de cores, usando o pacote Rcolorbrewer # Veja aconselhamento para cores de mapa aqui: http://colorbrewer2.org/ # Vou criar uma escala sequencial, com 5 tons de azul # Você pode ver todas as paletas disponíveis no pacote digitando display.brewer.all() my.Theme=rasterTheme(region=brewer.pal('Blues', n=5))

# Aqui eu personalizo a escala usando os valores do raster e um numero maior de quebras intervalos <- seq(trunc(min(values(r1))), trunc(max(values(r1))), length.out=length(my.Theme$regions$col)-1)

# Plota p <- levelplot(r1, margin=F, contour=T, par.settings=my.Theme, at=intervalos, main="Campo médio de Ozônio Novembros") p + layer(sp.lines(mapaSHP, lwd=0.7, col='black'))

# Para salvar, eu recomendo copiar a imagem direto do display para um arquivo pdf dev.copy2pdf(file="~/Desktop/ozone.pdf", width = 11, height = 8.5)

# Depois fica fácil convertê-la para png em altíssima resolução (600dpi) usando o ImageMagick (que precisa estar instalado): # Assim a sua imagem está de acordo com as normas de publicação das principais revistas científicas system("convert -density 600 ~/Desktop/ozone.pdf ~/Desktop/ozone.png")

Veja o resultado final: https://dl.dropboxusercontent.com/u/27700634/ozone.png

Greetings, -- Thiago V. dos Santos

PhD student Land and Atmospheric Science University of Minnesota

On Friday, February 19, 2016 1:18 PM, Mateus Dias Nunes < nunes.mateusdias@gmail.com> wrote:

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"; abaixo o link da figura que consegui gerar através dessa função.

https://www.dropbox.com/s/xi622mbpqg7tgch/campo_medio_O3.png?dl=0

abaixo o script que gerou o grafico com os contornos.

# Carregando biblioteca para manipular arquivos netCDF

library(maps) library(ncdf4)

#==========================================================================

#BIBLIOTECA "ncdf4". # PARA ESTE EXEMPLO SÃO USADOS RECURSOS PARA ABRIR A BIBLIOTECA ) #POIS RNetCDF, ncdf e ncdf4 APRESENTAM COMANDOS DIFERENTES PARA ABRIRMOS AS VARIÁVEIS

dados <- nc_open('teste.nc') # lendo coordenadas espaço-temporal lat <- ncvar_get( dados, 'lat' ) lon <- ncvar_get( dados, 'lon' ) time <- ncvar_get( dados, 'time' )

#=======================================================================================================

# lendo dados coluna total de Ozônio ColumnAmountO3 <- ncvar_get( dados, 'ColumnAmountO3' )

# dimensoes da variavel ColumnAmountO3 dims_ColumnAmountO3 <- dim(ColumnAmountO3)

# tornando o arranjo 3D (ColumnAmountO3) em um 2D, organizado em ptos de grade X tempo

dim(ColumnAmountO3) <- c( dims_ColumnAmountO3[1]*dims_ColumnAmountO3[2], dims_ColumnAmountO3[3] )

# calculando a média e retornado-a em 2D

media_ColumnAmountO3 <- rowMeans( ColumnAmountO3) dim(media_ColumnAmountO3) <- c( dims_ColumnAmountO3[1], dims_ColumnAmountO3[2] )

#========================================================================================================== # longitude varia de 0 a 360, convertendo para -180 a 180, essa conversão é feita para plotagem sobre o mapa for (i in 1:dim(lon)) { if (lon[i]>180) { lon[i] <- lon[i]-360 } }

# criando arquivo PNG que receberá o campo com o mapa

#png( filename="campo_medio_O3_jan2005.png",width=600,height=800 )

# plotando mapa da America do Sul map( xlim=c(-100,-10), ylim=c(-60,10) ) map.axes() # plotando eixos title( main="Campo médio de ozonio janeiros" ) # título do gráfico

# definindo intervalo de 5 Dobson Units (DU) intervalos = seq( trunc(min(ColumnAmountO3)), trunc(max(ColumnAmountO3)), 5 )

# adicionando campo de coluna de ozonio

contour( sort(lon), lat, media_ColumnAmountO3[ order(lon), ], add=T, levels=intervalos, lwd=2, labcex=1.3, col="black" )

# fechando arquivo PNG #dev.off()

obrigado

____________________________________________________________________________ MATEUS DIAS NUNES MESTRANDO DO PROGRAMA DE PÓS-GRADUAÇÃO EM METEOROLOGIA - PPGMET UNIVERSIDADE FEDERAL DE PELOTAS - UFPEL

TELEFONE: +55 (53) 81125154

____________________________________________________________________________ _______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forne� c�igo m�imo reproduz�el. _______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e fornea cdigo mnimo reproduzvel.

--- Este email foi escaneado pelo Avast antivírus. https://www.avast.com/antivirus

------------------------------

Message: 2 Date: Wed, 24 Feb 2016 08:31:03 -0300 From: Pedro Emmanuel Alvarenga Americano do Brasil <emmanuel.brasil@gmail.com> To: R BR lista <r-br@listas.c3sl.ufpr.br> Subject: [R-br] Artigo que não tem no periodicos capes Message-ID: < CAFfGvy+dQnZQ1XKJX3zUR56oH3ouV118FwW1Pq+6KQ73+cK84A@mail.gmail.com> Content-Type: text/plain; charset="utf-8"

Amigos de R,

Será que alguem tem na mão que possa me enviar o pdf desse artigo aqui:

http://onlinelibrary.wiley.com/doi/10.1002/sim.4780141909/full

Grato,

Pedro Brasil -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: < http://listas.inf.ufpr.br/pipermail/r-br/attachments/20160224/c6a4909e/attac...

...

------------------------------

Message: 3 Date: Wed, 24 Feb 2016 08:33:57 -0300 From: Marcelo Laia <marcelolaia@gmail.com> To: r-br <r-br@listas.c3sl.ufpr.br> Subject: Re: [R-br] Artigo que não tem no periodicos capes Message-ID: < CAEEYVUD8wuMn+N9ymLOd8_5wiP7uLFbxhx8WBUkExhRaPHc99g@mail.gmail.com> Content-Type: text/plain; charset="utf-8"

http://sci-hub.io/ Em 24 de fev de 2016 08:31, "Pedro Emmanuel Alvarenga Americano do Brasil" < emmanuel.brasil@gmail.com> escreveu:

...

Amigos de R,

Será que alguem tem na mão que possa me enviar o pdf desse artigo aqui:

http://onlinelibrary.wiley.com/doi/10.1002/sim.4780141909/full

Grato,

Pedro Brasil

_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.

-------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: < http://listas.inf.ufpr.br/pipermail/r-br/attachments/20160224/11e2033d/attac...

...

------------------------------

Subject: Legenda do Digest

_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br

------------------------------

Fim da Digest R-br, volume 62, assunto 29 *****************************************