In this example, we'll combine the previous techniques in order to combine our terrain hillshade from Chapter 7, Python and Elevation Data, with the color classification that we used on the LIDAR. For this example, we'll need the ASCII Grid DEMs named dem.asc and relief.asc that we used in the previous chapter. We'll create a colorized DEM and a hillshade and then use PIL in order to blend them together for an enhanced elevation visualization. The code comments will guide you through the example as many of these steps are already familiar to you:
import gdal_array as gd
try:
import Image
except:
from PIL import Image
relief = "relief.asc"
dem = "dem.asc"
target = "hillshade.tif"
# Load the relief as the background image
bg = gd.numpy.loadtxt(relief, skiprows=6)
# Load the DEM into a numpy array as the foreground image
fg = gd.numpy.loadtxt(dem, skiprows=6)[:-2, :-2]
# Create a blank 3-band image to colorize the DEM
rgb = gd.numpy.zeros((3, len(fg), len...