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!

Leia mais...

PSBASEMAP Parte 1

O PSBASEMAP  é o responsável por criar um código postscript que produz a borda das ilustrações geradas com o GMT. Apesar de seu uso não ser obrigatório, uma vez que é possível produzir a borda das ilustrações usando os demais comandos, é recomendável sempre utilizá-lo para tal fim de modo a permitir possíveis alterações nas configurações do GMTDEFAULT.
Além de gerar a base da ilustração que você deseja fazer, com o PSBASEMAP podemos adicionar escalas -L, definir a cor do fundo -G, definir a área a ser utilizada no mapa -R, o tipo de projeção utilizada -J, adicionar eixos e anotações -B ou mesmo para desenhar uma rosa dos ventos -T, definir a posição da ilustração gerada -X/-Y entre outras opções. 
Sintaxe básica:  
psbasemap -R,-J mais pelo menos UMA dessas opções obrigatoriamente -B,-G,-L,-T, 
onde: 


-R: Define a área a ser focada na ilustração ou do mapa.
Sintaxe básica:    -Rxmin/xmax/ymin/ymax/zmin/zmax
Podendo ser em números decimais (e.g., -147.456710) ou em graus minutos e segundos [graus:minutos:segundos][W|E|S|N]} (e.g., 123:45:35W ou -123:45:35).Pode-se também usar os atalhos –Rg para –R0/360/-90/90 e –Rd para –R-180/180/-90/90.  


-B: Define as marcações, anotações e nomes de cada eixo, além de definir o 
título e quais eixos serão mostrados na ilustração. 
Sintaxe básica:  
  -B[p/s][eixox/eixoy/eixoz][:.”Título”:][W|w][S|s][E|e][N|n]ou[Z|z], 
onde [p/s] correspondem a anotação principal e secundária respectivamente. 

Para cada eixo teremos as seguintes regras:  
a[x] – define o intervalo em que serão anotadas as marcações do eixo. 
f [x] – define o intervalo das marcações do eixo. 
g[x] – define o intervalo em que será desenhada a grade. 
Nas anotações [W|w][S|s][E|e][N|n]ou[Z|z] a letra maiúscula serve para indicar que o eixo será marcado e anotado, enquanto que a letra minúscula indica apenas marcação e caso seja suprimido de algum o eixo, o mesmo não terá anotação alguma. 


-E: Define o ponto de vista, azimute e elevação da figura em visões em perspectiva. 
Sintaxe básica:  
-E[grau do azimute]/[grau da elevação] 


 -J: Define o tipo de projeção a ser utilizada:  


-Jz: Define um eixo em z.
Sintaxe básica: 
-Jz[relação do eixo z com os demais] 




Projeção Cartesiana Linear


Sintaxe básica:

-JX[tamanho do eixo x]/[tamanho do eixo y]

Abaixo  segue um exemplo de imagem gerada utilizando-se apenas o psbasemap e seus respectivos parâmetros comentados a projeção utilizada foi a Cartesiana Linear.

psbasemap -R0/10/0/10/0/10 -JX10 -Jz0.5 -Glightblue -Bpa1f0.5g2:"eixo-x":/a1f0.5g2:"eixo-y":/a1f0.5g2:"eixo-z":WeSnZ -E210/30 > psbasemap.ps

Detalhes de cada parâmetro usado:


-R0/10/0/10/0/10 – Área representada no gráfico 


-JX10/10 – Projeção adotada: linear “X”


-Jz0.5 – Indica que foi adicionado um eixo em “z” com a metade dos demais eixos. 


-Glightblue – Cor do fundo da ilustração.


-Bpa1f0.5g2:"eixo-x":/a1f0.5g2:"eixo-y":/a1f0.5g2:"eixo-z":WeSnZ – Este talvez seja o parâmetro mais complicado e por que não, CHATO de se explicar. Para todos os eixos foi usado a1f0.5g2 que indica marcações a cada meia unidade, anotações a cada uma e grade a cada duas unidades de medida. WeSnZ indica que todos os eixos foram marcados mas somente os eixos W, S e Z foram anotados.


-E210/30 – Inclinação de 210º e elevação de 30º. 


Aproveitem para exercitar a criatividade criando figuras utilizando o script dado acima como ponto de partida alterando os valores dos parâmetros...


Boa Sorte!

Leia mais...