El problema es que `ndvi` es un *ndarray* de 2 dimensiones y cuando asignas valores a ambos índices te devuelve el *item* para esa posición del *ndarray* y eso no es iterable y es un escalar (`np.uint8`) y, por tanto, el `for` se 'rompe' ya que no puede iterar sobre el escalar. 

Directamente, en los bucles deberías hacer:

    for i in range(rows):
        for j in range(cols):
            k = ndvi[i,j] ## QUITAMOS EL BUCLE PARA 'k'
            print (str([i,j,k]) + '\n')
            if k < 0.0 :
                y[i,j] = red
                break
            if 0.0 < k < 0.2 :
                y[i,j] = yellowred
                break
            if 0.2 < k < 0.4 :
                y[i,j] = yellow
                break
            if 0.4 < k < 0.6 :
                y[i,j] = greenyellow
                break
            if k > 0.6 :
                y[i,j] = green
                break
            else:
                y[i,j] = black
                break