Learning Geospatial Analysis with Python - Second Edition

Creating a color hillshade

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 ...