I have a vector of data p_reconstructed that is zero everywhere except a few places. The nonzero indices are obtained with interpolant_points = find(p_reconstructed) and the corresponding values interpolant_values = p_reconstructed(interpolant_points).