elif np.isnan(rasterDEM[int(np.floor(Ly)), int(np.floor(Lx))]):