shape ## Be sure that we can do 'floating operation' imf = im.astype(sp.float64) ndvi = sp.empty((nl, nc)) if nb < 2: print "Two bands are needed to compute the NDVI" return None else: den = imf[:, :, ir - 1] + imf[:, :, r - 1] # Pre compute the denominator t...