script2.sh

#!/usr/bin/env bash

# Script 2: vamos pegar o grid do script 1 e fazer uma imagem bonita
# Imagem de área de estudo, grid com sombra e globinho

# modulos: makecpt, grdgradient, grdhisteq, grdimage, grdcontour, pscoast, psscale, psxy

# entrada: topo_srtm15_btaub.grd
# saida: myetopo1.cpt, hillshade-grad.nc, hillshade-hist.nc, topo_srtm15_btaub2.ps

# hoje a ideia é mexer com os parâmetros e tirar as dúvidas até aqui!

# primeiro sugiro nomear as coisas mais repetitivas no script, nome do ps, dimensões (-R) e formato do mapa (-J), depois é só chamar $nomevariavel

r="-R-48/-40/-25/-20"         # xmin/xmax/ymin/ymax
j="-JM20"                     # tem trocentas projeções, vamos ver algumas em sala, essa M é a mercartor, boa para regiões pequenas
g="topo_srtm15_btaub.grd"     # aceita varios formatos de grid
m="topo_srtm15_btaub2.ps"     # imprime em ps.. por isso precisa instalar o ghost!

#tabela de cores
gmt makecpt -Cglobe -T-3500/3500/10 > myetopo1.cpt                         

# fazer a sombra do grid
gmt grdgradient $g -Ghillshade-grad.nc -A335                 # -A é o azimute
gmt grdhisteq hillshade-grad.nc -Ghillshade-hist.nc -N.7     # -N é normalização, tentem varia-lo

# mapa
gmt grdimage $g $r $j -Cmyetopo1.cpt -Ihillshade-hist.nc -Ba2f1g1wEnS -Xc -Yc -K > $m   # com o -I colocamos o grid de sombra, -Ba<texto>f<borda>g<linhas>, -Xc centraliza
gmt grdcontour topo_btaub_30s.grd -R -J -C500m -W.001,50 -K -O >> $m                    # estou usando o grid com menos resolução feito no script1, tentem usar o $g para ver a diferença
gmt pscoast -R -J -B+t"Topografia" -Cdarkblue -Da -K -O >> $m                           # -B+t coloca o titulo, -C colore as regiões alagadas, -Smar, -Wlinhadecosta, -Gareaseca, -Da resolução auto
gmt pscoast -R -J -Na/0.3p -Df -Lf-41/-24.5/10/100+l -T-40/-26/1 -W -K -O >> $m         # -Lf escala +l unidade, -T norte, se coloca mais /1 rosa dos ventos
gmt psscale -D10/-1/18/0.5h -Ba1000:"m": -Cmyetopo1.cpt -E -K -O >> $m                  # -B intervalo dos numeros na paleta :textonapaleta:, -E coloca triangulos nas extremidades, ele reclama pq é do GMT4...
gmt pscoast -Rg -JG-50/5/4 -Bg30 -Swhite -Da -X18 -Y12 -A10000 -G50 -K -O >> $m         # -Rg projeção global, -JG é projeção global:x/y centrais e tamanho do globinho em cm (default)
                                                                                        # -Bg indica a distancia das linhas de latlong, -X e -Y sempre posição em cm (default)
                                                                                        # -A area minima em km2 para desenhar a coast, se n colocar ele desenha detalhes q ficam enormes em um globinho
                                                                                        # como vimos, -G pinta area seca, a cor pode ser o nome da cor, o numero RGB (3 num) ou escala de cinza (1 num) (0-255)
                                                                                        # no globinho não dá pra ver, mas no -GP/p tem 90 texturas

                                                             # EOF significa end of file, podemos criar um file daqui também com cat << EOF > nomefile.txt
gmt psxy -W0.5,red -Gred -R -J -O << EOF >> $m               # podemos entrar com um arquivo de texto, da no mesmo, só colocar o que está antes do EOF dentro do arquivo
-48 -25
-48 -20
-40 -20
-40 -25
-48 -25
EOF