Closed sutkss closed 4 years ago
missing value should be ignored when we calculate Y[x, :] - np.dot(W[x, :], H). To achieve this, after calculation of sparse code, masked position of Y have to be updated such that Y[masked, :] = np.dot(W[masked,:], H).
Y[x, :] - np.dot(W[x, :], H)
Y[masked, :] = np.dot(W[masked,:], H)
This operation should be added in line 71. https://github.com/hacarus/spm-image/blob/development/spmimage/decomposition/ksvd.py#L69-L71
Issue Description
missing value should be ignored when we calculate
Y[x, :] - np.dot(W[x, :], H). To achieve this, after calculation of sparse code, masked position of Y have to be updated such thatY[masked, :] = np.dot(W[masked,:], H).This operation should be added in line 71. https://github.com/hacarus/spm-image/blob/development/spmimage/decomposition/ksvd.py#L69-L71