We analyze seven year and nine year WMAP temperature maps for signatures of three finite flat topologies M_0=T^3, M_1=T^2 x R^1, and M_2=S^1 x R^2. We use Monte-Carlo simulations with the Feldman-Cousins method to obtain confidence intervals for the size of the topologies considered. We analyze the V, W, and Q frequency bands along with the ILC map and find no significant difference in the results. The 95.5% confidence level lower bound on the size of the topology is 1.5L_0 for M_0, 1.4L_0 for M_1, and 1.1L_0 for M_2, where L_0 is the radius of the last scattering surface. Our results agree very well with the recently released results from the Planck temperature data. We show that the likelihood function is not Gaussian in the size, and therefore simulations are important for obtaining accurate bounds on the size. We then introduce the formalism for including polarization data in the analysis. The improvement that we find from WMAP polarization maps is small because of the high level of instrumental noise, but our forecast for Planck maps shows a much better improvement on the lower bound for L. For the M_0 topology we expect an improvement on the lower bound of L from 1.7L_0 to 1.9L_0 at 95.5% confidence level. Using both polarization and temperature data is important because it tests the hypothesis that deviations in the TT spectrum at small l originate in the primordial perturbation spectrum.