I've been reading about the von Neumann entropy of a state, as defined via $S(\rho)=-\operatorname{tr}(\rho\ln \rho)$. This equals the Shannon entropy of the probability distribution corresponding to the basis of eigenvectors of $\rho$.
Given a POVM measurement $\{\Pi_j\}$, this maps a state $\rho$ into a vector of probabilities $p_j^\Pi=\operatorname{tr}(\Pi_j \rho)$, and we can associate a Shannon entropy $H(p^\Pi)$ to any such probability distribution $p^\Pi$.
Is it always the case that $H(p^\Pi) \ge S(\rho)$? Or more precisely, can we say that, for any state $\rho$, the von Neumann entropy $S(\rho)$ equals the minimal Shannon entropy, minimized over the set of possible measurements performed on $\rho$?
This is clear for pure states, as for those we have $S(\rho)=0$, but I'm not sure how to see it in the general case. I'm not even sure this is actually true, as for a maximally mixed state I think this would mean that the Shannon entropy is equal regardless of the measurement basis. Maybe it holds if the minimization is restricted to measurements with a number of components equal to the dimension of the state?