Approximating macroscopic observables in quantum spin systems with commuting matrices