Introducción al GMT (Generic Mapping Mapping Tools)
Por: Cristobal Condori Lima-2016
Generic Mapping Tools (GMT)
Página oficial: http://gmt.soest.hawaii.ed/ !toriales: http://gmt.soest.hawaii.ed/doc/".1.0/#aller$.html
¿Qué es GMT?
% El GMT es un programa gratuito y de código abierto para generar imágenes a partir de un conjunto de datos geográficos y cartesianas.
%
¿Qué es GMT?
% Funciona en plataformas de: Windows Mac y !inu". % !enguaje de programación: GMT#fortran GMT#$yt%on y GMT#Matlab
Instalación ' ersiones recientes de 3bnt/4ebian $a 5ienen con repositorios de #*! gmt
&' https://www.soest.hawaii.ed/gmt/gmt/gmt(install(form.html gedit .bashrc )#*! e+port #*!,*/+++/+++/++++ e+port P!,P!,:#*!,*/bin sorce .bashrc
Parámetros de formato (GMT defaults) %)/bin/bash %) C*7!89 %)89&L %gmtset.......
#*!!: define el aspecto general del mapa; por e
neas de comando gmtset 7!(?7! 2 7!(?7!(9@ 16 gmtset PP8(*49 A gmtset P#(897!!97 portrait
Script GMT psbasemap psbasemap -8-B"/-6"/-20/0 -*1"c -&a"g2 -Dc -Ec - F map.ps
-8-B"/-6"/-20/0: -8min.lon./ma+.lon/min.lat./ma+.lat -*1"c: Pro$ecciGn $ tama=o del mapa. -* es *ercator $ 1"c es el ancho total del mapa en cm. -&a"g2: tiHetado en los e
Proyección !
Script GMT pscoast pscoast -8-B"/-6"/-20/0 -*1"c -&a"g2 -4h -71 -J1;brown -#2A"/2A"/200 -1A0/2K"/2"" -Dc -Ec - -P F map.ps -4h: 4efine la resolciGn. -7: *estra los l>mites de frontera de pa>ses; departamentos $ mar>timas. -J1;brown: ncho $ color de las l>neas de costa. -#2A"/2A"/200: 4efine el color Itierra' de los continentes. -1A0/2K"/2"": 4efine el color Iaga' de los ocanos; r>os $ lagos. -P: rientaciGn de página de la figra
Prácticas %*apa >smico 1a(Plot(ismos(Per.gmt 1b(Plot(ismos(Per(Color(&ar.gmt
%Perfil >smico 2(Perfil(ismico.gmt
%*apa !opográfico K(*apa(!opografico.gmt
%*apa de *ecanismos ?ocales A(*ecanismos(?ocales.gmt
Script "# Mapa S$smico psxy
)/bin/bash )crip para generar mapa con distribciGn de sismos )tor: C. Condori )Ciencias de la !ierra Glida C! gmtset *38(379! cm gmtset 7!(?7!(9@ 10 gmtset L&L(?7!(9@ 10 gmtset ,48(?7!(9@ 10 gmtset !9CM(L7#!, 0.1 gmtset ?8*(P7 K.2" gmtset &*P(!EP ?7CE gmtset J7!(38(?7! !83 gmtset PL!(4#8(?8*! -4 gmtset 4!(P8(97C, K00 gmtset P#(897!!97 landscape name*apa(ismico(del(Per2 titloN*P O*9C 4L P8N figname.ps 8plot-8/-BA/-6B/-20/00 psbasemap 8plot -*1"." -&a2f2Jen:N.titloN: -P -Dc -Ec - -M F fig pscoast -8 - -4h -J1p;blacQ -71/1p;blacQ -lightble -#lightbrown -20 -!f-R0/-2/1/2 -Lf-R".0/-1S/-20.00/K00 - - -M -P FF fig
Script "# Mapa S$smico psxy
)4epartamentos PerT awQ UVprint 1; 2WU deptos.dat X ps+$ -8 - -* -J0.1p;blacQ;( - -P -M FF fig )?osa; 4orsal de 7aYca $ ?ractra de *enda=a ps+$ ?osa.gmt -8 - -* -f0.Bc/0.1"clt:0.Bc -J2/2""/2""/2"" -#2""/2""/2"" - - -P -M FF fig ps+$ *endana.gmt -8 - -* -J1p;white;( - -M FF fig ps+$ 7aYca.gmt -8 - -* -J1p;white;( - -M FF fig )plotear sismos cat Catalogo(79C-1S06(2016.t+t X awQ UVif I" Z 60' print A;K;6[0.0K;" WU F sperficiales.gmt ps+$ sperficiales.gmt -8 - -c -J0.01 -#2""/0/0 -J0.01 - - -M FF fig cat Catalogo(79C-1S06(2016.t+t X awQ UVif I"F60.1 \\ "ZK00' print A;K;6[0.0K;"WU F intermedios.gmt ps+$ intermedios.gmt -8 - -c -J0.01 -#0/2""/0 -J0.01 - - -M FF fig cat Catalogo(79C-1S06(2016.t+t X awQ UVif I" F K00.1' print A;K;6[0.0K;" WU F profndos.gmt ps+$ profndos.gmt -8 - -c -J0.01 -#0/0/2"" -J0.01 - - -M FF fig
Script "# Mapa S$smico psxy )Le$enda ) # is 5ertical gap; is 5ertical line; 7 sets ) of colmns; 4 draws horiYontal line. ) , is header; L is label; is s$mbol; ! is paragraph te+t; * is map scale. pslegend ZZ74 -4+2.R/0/2."/A.0/8& - -8 -? -#light$ellow - FF fig # 0.2c , 10 1 *agnitd # 0.2 71 0.K"c c 0.12c red 0.Ap 0.Ki A.0 # 0.1 0.K"c c 0.1"c red 0.Ap 0.Ki ".0 # 0.1 0.K"c c 0.1Bc red 0.Ap 0.Ki 6.0 # 0.1 0.K"c c 0.21c red 0.Ap 0.Ki R.0 # 0.1 0.K"c c 0.2Ac red 0.Ap 0.Ki B.0 # 0.1 74 ) *estra la ?igra
Script "# Mapa S$smico psxy
)plotear sismos cat Catalogo(79C-1S06(2016.t+t X awQ UV print A;K;"[-1;6 WU F tmp.t+t maQecpt -!-R00/0/" -@ -Crainbow F depth.cpt ps+$ tmp.t+t -8 - -c0.2c -JK -Cdepth.cpt - - -M FF fig psscale -Cdepth.cpt -46."i/K.Ki/Ai/0.2"i5 -&a"0:NProfndidad IMm'N: - - FF fig
Script %# Perfil S$smico
pscoupe
Script %# Perfil S$smico )Pntos iniciales del perfil pro
Script &# Mapa Topo'ráfico grdimage CPT color para topografía y batimetría
http://soliton.5m.b$temarQ.co.Q/pb/cpt-cit$/ grdgradient !P1(&ed(g(gmtA.grd -K0" -#mshadow.grd -7t - grdimage !P1(&ed(g(gmtA.grd 8plot -*1Rc -&a2f2Jen -Cele5ation.cpt -M -P F fig grd5iew -8 - !P1(&ed(g(gmtA.grd -&02hg0Ah -Y0.000"c -1B0/0S0 -]i2"0 -Cele5ation.cpt -9mshadow.grd -P - -M FF fig psscale -Cele5ation.cpt -4B/-1/10/0."h -&2000/:*!8: - -P -M FF fig pscoast -8 - -4h -J1/0/0/2"" -71/0.0Kc -J0.0Kc/0/0/0 -20 - -!f-R0/-2/1/2 -Lf-R0.0/-1S."/-20.00/K00 -P - FF fig ) *estra la ?igra g5 fig
Script # Mecanismos ocales psmeca http://www.globalcmt.org/C*!search.html
)plotear mecanismos focales psmeca mecaper.t+t - - -8 -* -J1 -L -#red -4000/60 -d.60/-1/0 -M FF fig psmeca mecaper.t+t - - -8 -* -J1 -L -#green -460.1/K00 -d.60/-1/0 -M FF fig psmeca mecaper.t+t - - -8 -* -J1 -L -#ble -4K00.1/R00 -d.60/-1/0 FF fig ) *estra la ?igra g5 fig ps2pdf *apa(*eca(Per.ps
