It is not true that a measurement result of 0 distinguishes those two states. The probability of getting $a$ as a result of a measurement of the observable associated with the operator $\hat{A}$ is given by
$$
\Pr(m=a)=\operatorname{Tr}(\rho\hat{A})\,,
$$
and so if we compute these probabilities with the density matrices above, we have to get the same result, since the density matrices are mathematically the same.
To see this explicitly for this case, we first note that
$$
\lvert{\pm}\rangle=\frac{\lvert{0}\rangle\pm\lvert{1}\rangle}{\sqrt{2}}\,.
$$
Taking
$$
\hat{Z}=(1)\lvert{1}\rangle\langle1\rvert+(0)\lvert{0}\rangle\langle 0\rvert
=\lvert{1}\rangle\langle1\rvert\,
$$
as the operator we're measuring, we can then compute the probabilities. First,
$$
\Pr(m=0|\rho_1)=
\operatorname{Tr}\left(\rho\hat{A}\right)
=\operatorname{Tr}\left(
\left(
\frac{1}{2}|0\rangle\langle 0|+\frac{1}{2}| 1\rangle\langle 1|
\right)
\lvert{1}\rangle\langle1\rvert\right)
=\frac{1}{2}
\operatorname{Tr}\left(\lvert{1}\rangle\langle1\rvert\right)
=
\frac{1}{2}\,.
$$
A similar computation yields exactly the same result for the second state. (I leave it as an exercise. I didn't want to do it because as part of the computation, I would have just shown that the second density matrix can be written as the first.)
Any probability of an outcome for any measurement can be shown to yield the same result for these two density matrices (again, they are mathematically the same, so they have to yield the same results), and so these represent the same physical state.
I think perhaps the confusion comes in by interpreting
$$
\rho = \frac{1}{2}|0\rangle\langle0|+\frac{1}{2}|1\rangle\langle1|\,,
$$
as being similar in nature to
$$
\lvert\psi\rangle = \frac{1}{\sqrt{2}}\lvert0\rangle + \frac{1}{\sqrt{2}}\lvert1\rangle\,.
$$
That is, if we consider the two states
$$
\lvert\psi\rangle = \frac{1}{\sqrt{2}}\lvert0\rangle + \frac{1}{\sqrt{2}}\lvert1\rangle = \lvert+\rangle\,,
$$
and
$$
\lvert\psi\rangle = \frac{1}{\sqrt{2}}\lvert+\rangle + \frac{1}{\sqrt{2}}\lvert-\rangle = \lvert0\rangle\,,
$$
the probabilities of getting 0 as the result of a measurement of $\hat{Z}$ are clearly different! In one case the probability is 1/2 and in the other the probability is 1. Density matrices just don't work this way.
On the other hand, let's finally address the suggestion by the OP. Suppose we measured $\hat{Z}$ on only one of the systems in the ensemble. The OP suggests that upon getting the result 0, this would reveal that that particular system was one of the 50% of systems that was in state $|0\rangle$, thus revealing that the entire ensemble must consist of systems that are in either $|0\rangle$ or $|1\rangle$.
However, getting 0 as the result of the measurement doesn't tell us what the state before the measurement was. Rather, it only tells us that the system is in state $|0\rangle$ post-measurement*. Since there is a non-zero probability of getting 0 as the result of the measurement when half the systems are in state $|+\rangle$ and half are in state $|-\rangle$, this measurement cannot distinguish which statistical configuration the ensemble is in.
And, finally, in fact, in order to figure that out, you have to do a complete ensemble measurement, and the statistics of that sequence of measurements will be identical in both cases.