Bug or wishful thinking on my part? from scipy import sparse import numpy as np X = np.random.rand(50,3) spX = sparse.csr_matrix(X) # any sparse matrix really newX = np.asanyarray(spX) spX is newX spX.ndim == newX.ndim Cheers, Skipper