Analysing a large neural population concurrently in space and time through non-negative matrix factorization