Mapas Geológicos: 13.5. Bloque 3D + Imagen Satelital

jueves, 17 de mayo de 2018

13.5. Bloque 3D + Imagen Satelital

Superponer una imagen satelital sobre un bloque 3D.

Introducción: Se explica como superponer una imagen satelital sobre un bloque 3D. Para poder superponer la imagen satelital es necesario descomponer en grillas RGB.

Script

0. Datos:

DEM: Grilla GMRT3.1.grd. Ver anexo 4.
Imagen Satelital: SanJuan_Geo.tif. Ver anexo 2.

1. Imagen a RGB: con el software Mirone (ver 2.2) se descompone la imagen en grillas 3 grillas:

2. Crear bloque 3D: grdview se crea el bloque 3D con la superposición de la imagen satelital.

-Qi300: indica la resolución de 300 dpi de la imagen superpuesta.
-G: Grillas RGB.

gmt grdview %CUT% -R -J -JZ -O -K -p -Qi300 >> %OUT% -G"SanJuan_SAT_r.grd","SanJuan_SAT_g.grd","SanJuan_SAT_b.grd" -N-10000+glightgray -Wf0.5 -BnSwEZ+b -Baf -Bzaf+l"Altura (m)"

3. Complementar bloque 3D: Se pueden agregar datos de océanos, línea de costa, límites, rutas, pueblos y ciudades.

gmt pscoast -R -J -JZ -O -K -p160/30/0 -Df >> %OUT% -Sdodgerblue2 -A0/0/1
gmt pscoast -R -J -JZ -O -K -p160/30/0 -Df >> %OUT% -W1/0.3,black
gmt pscoast -R%REGION% -Df -M -N1/ | grdtrack -G%CUT% -sa | psxyz -R%REGION3D% -J -JZ -O -K -p >> %OUT% -W0.5,black
gmt pscoast -R%REGION% -Df -M -N2/ | grdtrack -G%CUT% -sa | psxyz -R%REGION3D% -J -JZ -O -K -p >> %OUT% -W0.2,black,-
gmt grdtrack -R%REGION% "E:\Facultad\Datos_Geofisicos\IGN\1_GMT\003_Red_Ferroviaria.gmt" -G%CUT% -sa | psxyz -R%REGION3D% -J -JZ -O -K -p >> %OUT% -Wthin,blue
gmt grdtrack -R%REGION% "E:\Facultad\Datos_Geofisicos\IGN\1_GMT\005_Centros_Poblados.gmt" -G%CUT% -sa | psxyz -R%REGION3D% -J -JZ -O -K -p >> %OUT% -W0.2,black,- -Sc0.1 -Gred
gmt grdtrack -R%REGION% "E:\Facultad\Datos_Geofisicos\IGN\1_GMT\010_Ejidos_Urbanos.gmt" -G%CUT% -sa | psxyz -R%REGION3D% -J -JZ -O -K -p >> %OUT% -Wfaint -Ggreen