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 comment:

  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)