Below we calculate the GICs of the Rice-Mele model, which is a lattice model (in other words, the Wannier functions are completely localized on particular sites), and approximately reconstruct the probability distribution of the polarization. Our reconstruction is based [@Bandyopadhyay05; @Collins77; @Mead84] on maximizing the information entropy under the constraints provided by the moments calculated. The expression for the entropy we use is $$\label{eqn:SPx}
S[P(x)] = - \int d x P(x) \ln P(x),$$ minimized as a functional of $P(x)$ under the constraints $$\mu_P^{(k)} = \int d x P(x) x^k,$$ as well as the constraint that $P(x)$ is normalized. The functional minimization of Eq. (\[eqn:SPx\]) under the constraints results in the functional differential equation $$\label{eqn:Sopt}
\frac{\delta}{\delta P(x)} \left[ S[P(x)]
- \sum_k A_k (\mu_P^{(k)} - \mu_C^{(k)})
 \right] = 0,$$ where $\mu_C^{(k)}$ are the moments obtained from the cumulants of the Berry phase (see Eq. (\[eqn:muC\])), and $A_k$ are Lagrange multipliers. The solution of Eq. (\[eqn:Sopt\]) is $$P(x) = C\exp\left( - \sum_k A_k x^k\right),$$ where $C$ is the normalization constant. We determine the constants $A_k$ by numerically minimizing the quantity $$\label{eqn:chisquared}
\chi^2 = \sum_k (\mu_P^{(k)}-\mu_C^{(k)})^2,$$ as a function of $A_k$. As our initial guess in all cases studied below, we take the Gaussian 