Замените ненулевые разреженные матричные элементы на таковые
R = spones(S)
R = spones(S) генерирует матрицу R с такой же структурой разреженности, как и S, но с 1находится в ненулевых положениях.
c = sum(spones(S)) - количество ненулевых в каждом столбце.
r = sum(spones(S'))' - количество ненулей в каждой строке.
sum(c) и sum(r) равны и равны nnz(S).