You don't miss anything! It's just a wrong explanation for the "phase-space-cell factor" h^3=(2 \pi \hbar)^3 in statistical mechanics.
The correct explanation is as follows. Suppose we have a gas of particles in a container of volume V. It's important to have a container of finite size. The shape doesn't matter. To keep the argument simple, suppose it's a cube with length L. Take this container as a typical box in a very large box containing the gas. Then it makes sense to assume that the container is as good as any other box within the big box, and we can impose periodic boundary conditions. Then the momentum operator for particles in the container is well defined, and the eigenvalues of the three momentum components are discrete. The normalized eigenvectors are
u_{\vec{p}}(\vec{x})=\frac{1}{L^{3/2}} \exp \left (\frac{\mathrm{i} \vec{p} \cdot \vec{x}}{\hbar} \right )
with the boundary conditions enforcing that
\vec{p} \in \frac{2 \pi \hbar}{L} \mathbb{Z}^3.
Now to get the thermodynamic potentials you have to sum over the states or count states. In this case you have to sum over all these momenta. Calculating such sums is much more tedious than doing integrals, and often you are not interested on finite-size effects at all, i.e., you make the container very large and take the "thermodynamic" limit. If you make the container very large, i.e., L becomes very large, the grid of momenta becomes very fine, and it is appropriate to replace the sum by an integral. To get this limit right, we need a density of states in a given momentum volume \mathrm{d}^3 \vec{p}. In the discretized version in such a momentum volume you have
\mathrm{d} \Omega = \mathrm{d}^3 \vec{p} \frac{L^3}{(2 \pi \hbar)^3}=\mathrm{d}^3 \vec{p} \frac{V}{(2 \pi \hbar)^3}
momentum states for your single particle.
So in the thermodynamic limit you can substitute the sum over momentum states by an integral with the just calculated integration measure:
\sum_{\vec{p} \in \frac{2 \pi \hbar}{L} \mathbb{Z}^3} \rightarrow \int_{\mathbb{R}^3} \mathrm{d}^3 \vec{p} \frac{V}{(2 \pi \hbar)^3}.
Note that here also the dimensions are correct, because both the sum symbol and the integral are formally of dimension 1.
That's the true reason for the phase-space volume (2 \pi \hbar)^3 in statistical mechanics. It has nothing to do with the uncertainty relation in a direct way, but with the number of states in a given momentum range for a particle in a container of volume V in the thermodynamical limit!