Speed up computation of nodata mask.

Proposal:

For the computation of the nodata mask it is not needed to run np.all/np.any on all pixels in a multispectral image. Just run it on the first that contains some data to get a first guess and then only run it on those pixels that are labelled as nodata due to the first band.