Monday, December 1, 2014


In Chapter 8 of my book, the shapefile doesn't have a .prj file which isn't a big deal for the example

but can be annoying if you try to open it in QGIS or ArcGIS.  There was also an unused numpy array in the example so I updated the whole package.  You can download the updated script and data here:

This example is the very first example in the chapter.  The second example shows you how to classify the output data to create the image shown here.


  1. Joel, great book!
    Just in case, in Chapter 8 you say:
    ndvi = 1.0 * (irClip - rClip) / irClip + rClip + 1.0

    when I guess you meant to say:
    ndvi = 1.0 * (irClip - rClip) / (irClip + rClip + 1.0)


  2. Thank you so much for sharing the information. I hope to hear more updates from you.
    outlook email login