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!