quarta-feira, 11 de agosto de 2010

PSBASEMAP Parte 2

Projeção cartesiana de tempo –JXT/t

        Freqüentemente se faz necessária a produção de um gráfico de uma variável em função do tempo. Para essa finalidade existe a projeção cartesiana de tempo, onde podemos utilizar um eixo para a passagem de tempo, em geral “eixo x”, e outro eixo para uma variável genérica, em geral o “eixo y”. Nos exemplos abaixo utilizaremos somente o “eixo x” a fim de explicar a notação de tempo.
Primeiro exemplo:
gmtset PLOT_DATE_FORMAT -o ANNOT_FONT_SIZE_PRIMARY +9p TIME_LANGUAGE BR
psbasemap -R2010-6-1T/2010-8-1T/0/1 -JX6T/0.2 -Bpa7Rf1d -Bsa1OS -P > exemplotempo1.ps
Note que a formatação dos parâmetros –R e –B está diferente. Em –R foi utilizado 2010-6-1T/2010-8-1T/0/1, isto indica que o eixo começara no dia 1º de junho e termina no dia 1º de agosto, além disso, foi adicionado o”T” para indicar que se trata de tempo. Em –B foram utilizadas novas opções além das já conhecidas “a”, “f” e “g”. São as seguintes:


Segundo exemplo:
gmtset TIME_LANGUAGE BR
gmtset PLOT_CLOCK_FORMAT hh:mm ANNOT_FONT_SIZE_PRIMARY +9p
psbasemap -R2010-8-11T/2010-8-13T/0/1 -JX8t/0.2 -Bpa6Hf1h -Bsa1KS -P -K > exemplotempo2.ps
gmtset PLOT_DATE_FORMAT "o dd" ANNOT_FONT_SIZE_PRIMARY +9p
psbasemap -R -JX -Bpa6Hf1h -Bsa1DS -O -Y0.65i >> exemplotempo2.ps
Na fgura de baixo utilizamos como anotação, os dias da semana a1K e no exemplo de cima a data a1D. Note que na anotação de horario foram usados somente horas e minutos, PLOT_CLOCK_FORMAT hh:mm e que na anotação da data foram usados o nome do mês e os dígitos do dia do mês PLOT_DATE_FORMAT “o dd”.

Terceiro exemplo:
gmtset TIME_LANGUAGE BR
gmtset PLOT_DATE_FORMAT o TIME_FORMAT_PRIMARY Character
psbasemap -R2008T/2010T/0/1 -JX8T/0.2 -Bpa3Of1o -Bsa1YS -P 
> exemplotempo3.ps
Este exemplo apresenta um intervalo de dois anos ambos com anotações a cada três meses. O comando do gmtset PLOT_FORMAT_PRIMARY Character determinou que os nomes dos meses fossem abreviados com sua letra inicial.

Quarto exemplo:
gmtset PLOT_CLOCK_FORMAT -hham
psbasemap -R0t/0.25t/0/1 -JX-10/0.2 -Bpa15mf5m -Bsa1HS -P
> exemplotempo4.ps
Este exemplo apresenta o intervalo de algumas horas de um dia, perceba que o tempo corre da direita para a esquerda por causa do sinal negativo usado na dimensão do eixo x.

Quinto exemplo:
gmtset PLOT_DATE_FORMAT u TIME_FORMAT_SECONDARY full TIME_LANGUAGE BR
psbasemap -R2010-9-11T/2010-9-29T/0/1 -JX8t/0.2 -Bpa1K -Bsa1US -P -K > exemplotempo5.ps
gmtset PLOT_DATE_FORMAT o TIME_WEEK_START Sunday
gmtset TIME_FORMAT_PRIMARY Char TIME_LANGUAGE BR
psbasemap -R2010-9-8T/2010-9-29T/0/1 -JX8t/0.2 -Bpa3Kf1k -Bsa1rS -O -Y0.65i >> exemplotempo5.ps

            Este quinto exemplo mostra um intervalo de tempo de algumas semanas, sendo a figura de cima com um intervalo de semanas com a semana começando no domingo. Comando do gmtset: TIME_WEEK_START Sunday.
Já a próxima figura apresenta o mesmo intervalo de tempo, mas com anotações a cada dia da semana abreviados e as semanas marcadas com números.

Sexto exemplo:
gmtset PLOT_DATE_FORMAT "o yy" TIME_FORMAT_PRIMARY Abbreviated TIME_LANGUAGE BR
psbasemap -R1996T/1996-6T/0/1 -JX10T/0.2 -Ba1Of1dS -P > exemplotempo6.ps 


O sexto exemplo nos mostra o intervalo de tempo dos cinco primeiros meses de 2010, sendo que as anotações são de mês em mês com os nomes abreviados adicionados dos dois dígitos do ano.

Sétimo exemplo:
gmtset PLOT_DATE_FORMAT jjj TIME_INTERVAL_FRACTION 0.05
psbasemap -R2000-12-15T/2001-1-15T/0/1 -JX8T/0.2 -Bpa5Df1d -Bsa1YS -P > exemplotempo7.ps

           O sétimo e último exemplo nos mostra uma ilustração dos 15 últimos dias de 2000 e 15 primeiros dias de 2001 com anotação principal dos dias do ano e secundária do ano.

Agora que aprendemos como utilizar as variadas escalas de tempo. Criaremos um gráfico utilizando a escala de tempo no eixo x e outra dimensional no eixo y utilizando o script abaixo:
gmtset PLOT_DATE_FORMAT -o ANNOT_FONT_SIZE_PRIMARY +9p
gmtset CHAR_ENCODING ISOLatin1+ TIME_LANGUAGE BR
psbasemap -R2008-1-1T/2008-2-1T/-20/0 -JX15/10 -Bpa7Rf1d/a2f10S -Bsa1O/a2f10:" Temperatura M\351dia em \260C":SW -P > psbasemap2.ps

Por enquanto é isso, no próximo post entrarei no mérito das projeções de mapas!
Até mais!

0 comentários:

Postar um comentário

Translate

Seguidores

Total de visualizações de página

  ©Dicas de GMT - Generic Mapping Tools - Todos os direitos reservados.

Template by Dicas Blogger | Topo