Soil moisture percentage 2002-08-25 (vertical exaggeration 50)
#!/bin/bashfor ANGLE in `seq 0 5 360` doprintf -v NAME "out_%03d" $ANGLEecho $NAMEX=`python -c "import math;print 0.5*math.cos(${ANGLE}./180.*math.pi)+0.5"`Y=`python -c "import math;print -0.5*math.sin(${ANGLE}./180.*math.pi)+0.5"`printf -v XX "%0.2f" ${X}printf -v YY "%0.2f" ${Y}echo $XX, $YY m.nviz.image -b elevation_map=dem_small@soils,dem_small@soils -a mode=fine resolution_fine=1 color_map=dem_small@soils,dem_small@soils
volume=interp_2002_08_25@soils volume_shading=gouraud volume_resolution=1 volume_position=0,0,19.5 slice=1:z,1:z,1:z,1:z
slice_position=0,1,1,1,0,1,0,1,0,0,0,1,1,1,0,1,0,1,0,0,1,0,0,1 slice_transparency=127,127,127,127 position=$XX,$YY
height=30 perspective=19 twist=0 zexag=50. focus=226,158,20 light_position=-0.47,0.61,0.80 light_brightness=80
light_ambient=20 light_color=255:255:255 output=${NAME} format=ppm size=798,533
vpoint=observation3D_2002_08_25@soils,observation3D_2002_08_25@soils vpoint_width=2,2 vpoint_color=0:0:255,255:0:0
vpoint_size=2,6 vpoint_marker=sphere,x vpoint_position=0,0,0,0,0,0 vpoint_layer=1,1 vpoint_mode=3D,surface convert ${NAME}.ppm legend.png -gravity SouthWest -geometry +10+10 -compose over -composite ${NAME}.png done convert -delay 1x5 `seq -f out_%03g.png 0 5 360` -coalesce -layers OptimizeTransparency animation_rotate.gif
Soil moisture percentage 2002-2005 (vertical exaggeration 50)
#!/bin/bashfor VOXEL in `g.mlist --quiet type=rast3d pattern=interp* exclude=*nomask* separator=" "`doDATE=`python -c "print '${VOXEL}'.split('_',1)[1]"`DATE2=`python -c "print '${DATE}'.replace('_', '-')"`echo ${DATE}
m.nviz.image -b elevation_map=dem_small@soils -a mode=fine resolution_fine=1 color_map=cross_dem_interp_nomask_${DATE}
volume=${VOXEL} volume_shading=gouraud volume_resolution=1 volume_position=0,0,19.5 slice=1:z,1:z,1:z,1:z
slice_position=0,1,1,1,0,1,0,1,0,0,0,1,1,1,0,1,0,1,0,0,1,0,0,1 slice_transparency=127,127,127,127 position=0,1 height=35
perspective=19 twist=0 zexag=50 focus=226,158,20 light_position=-0.47,0.61,0.80 light_brightness=80 light_ambient=20
light_color=255:255:255 output=out_${DATE} format=ppm size=798,533 vpoint=observation3D_${DATE}@soils,observation3D_${DATE}@soils
vpoint_width=2,2 vpoint_color=0:0:255,0:0:0 vpoint_size=2,6 vpoint_marker=sphere,x vpoint_position=0,0,0,0,0,0 vpoint_layer=1,1 vpoint_mode=3D,surface
convert out_${DATE}.ppm legend.png -gravity SouthWest -geometry +10+10 -compose over -composite out_${DATE}.png mogrify -fill black -pointsize 16 -annotate +20+50 ${DATE2} out_${DATE}.png done convert +dither -delay 1x2 `ls out*.png` -coalesce -layers OptimizeTransparency animation_time-series.gif
Soil moisture percentage 2003 (without taking into account landscape)
done in Animation Tool
isosurface of value 28 combined with slices, interpolation tension 20, done in animation tool