Learning nonnegative matrix factorizations from compressed data