/usr/share/doc/gmt/examples/anim02/anim_02.sh is in gmt-examples 5.2.1+dfsg-3build1.
This file is owned by root:root, with mode 0o755.
The actual contents of the file can be viewed below.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 | #!/bin/bash
# GMT ANIMATION 02
# $Id: anim_02.sh 15178 2015-11-06 10:45:03Z fwobbe $
#
# Purpose: Make web page with simple animated GIF of a DEM grid
# GMT progs: gmt gmtset, gmt gmtmath, gmt grdgradient, gmt makecpt, gmt grdimage gmt psxy, gmt psconvert
# Unix progs: awk, mkdir, rm, mv, echo, convert, cat
# Note: Run with any argument to build movie; otherwise 1st frame is plotted only.
#
# 1. Initialization
# 1a) Assign movie parameters
. gmt_shell_functions.sh
width=3.5i
height=4.15i
dpi=72
n_frames=36
name=anim_02
ps=${name}.ps
# 1b) setup
del_angle=`gmt gmtmath -Q 360 ${n_frames} DIV =`
gmt makecpt -Crainbow -T500/4500/5000 -Z > $$.cpt
# 2. Main loop
mkdir -p $$
frame=0
while [ ${frame} -lt ${n_frames} ]; do
# Create file name using a name_##.tif format
file=`gmt_set_framename ${name} ${frame}`
angle=`gmt gmtmath -Q ${frame} ${del_angle} MUL =`
dir=`gmt gmtmath -Q ${angle} 180 ADD =`
gmt grdgradient us.nc -A${angle} -Nt2 -fg -G$$.us_int.nc
gmt grdimage us.nc -I$$.us_int.nc -JM3i -P -K -C$$.cpt -BWSne -B1 -X0.35i -Y0.3i \
--PS_MEDIA=${width}x${height} --FONT_ANNOT_PRIMARY=9p > $$.ps
gmt psxy -Rus.nc -J -O -K -Sc0.8i -Gwhite -Wthin >> $$.ps <<< "256.25 35.6"
gmt psxy -Rus.nc -J -O -Sv0.1i+e -Gred -Wthick >> $$.ps <<< "256.25 35.6 ${dir} 0.37"
[[ ${frame} -eq 0 ]] && cp $$.ps ${ps}
if [ $# -eq 0 ]; then
gmt_cleanup .gmt
gmt_abort "${0}: First frame plotted to ${name}.ps"
fi
# RIP to TIFF at specified dpi
gmt psconvert -E${dpi} -Tt $$.ps
mv -f $$.tif $$/${file}.tif
echo "Frame ${file} completed"
frame=`gmt_set_framenext ${frame}`
done
# 3. Create animated GIF file and HTML for web page
${GRAPHICSMAGICK-gm} convert -delay 10 -loop 0 $$/${name_}*.tif ${name}.gif
cat << END > ${name}.html
<HTML>
<TITLE>GMT shading: A tool for feature detection</TITLE>
<BODY bgcolor="#ffffff">
<CENTER>
<H1>GMT shading: A tool for feature detection</H1>
<IMG src="${name}.gif">
</CENTER>
<HR>
We make illuminated images of topography from a section of Colorado and
vary the azimuth of the illumination (see arrow). As the light-source sweeps around
the area over 360 degrees we notice that different features in the data
become hightlighted. This is because the illumination is based on data
gradients and such derivatives will high-light short-wavelength signal.
Again, our animation uses GraphicsMagick's convert tool to make an animated GIF file
with a 0.1 second pause between the 36 frames.
<HR>
<I>${name}.sh: Created by ${USER} on `date`</I>
</BODY>
</HTML>
END
# 4. Clean up temporary files
gmt_cleanup .gmt
|