Perspectiva 3D

GMT progs : grdview, grdcut

#!/bin/bash

paleta="/opt/gmt/share/cpt/GMT_relief.cpt"
grilla="/vox/201001/linux/topografia/ETOPO1_Bed_g_gdal.grd"
ps="cile.ps"
zona="-74/-71/-40/-30"
grdcut $grilla -R$zona -G$ps.grd
grdgradient $ps.grd -A50 -G$ps.grad -Nt
grdview $ps.grd -JM2.5i -C$paleta -Qs -JZ1i -I$ps.grad -E230/45 -Wc -B1a2 -K  > $ps
rm *.grad *.grd

Explicación
Generamos una perspectiva 3D del borde continental de Chile, donde se aprecia la fosa de subducción y la topografá continental. Formamos una sub-grilla de la zona a grafica con grdcut, luego añadimos iluminación con un acimut de 50 grados con el comando grdgradient para finalmente graficar la perspectiva con el comando grdview donde los parámetros funcionan del mismo modo como el ejemplo 2.



gonzalo 2011-09-16