We propose a model reduction method for positive systems that ensures the positivity of the reduced model. For both, continuous-time and discrete-time systems, our approach is based on constructing diagonal solutions of Lyapunov inequalities. These are linear matrix inequalities (LMIs), which are shown to be feasible. Stability is preserved and an error bound in the $\mathcal{H}_\infty$-norm is provided.