Top panels: Ionization maps and galaxies colored by their Lyα transmissions, at $\overline{x}=0.5$, 0.7, and 0.9 for our z = 8 simulation. Lower panels: Recovered bubble edges (red contours) using the approach described in Section 3 and Appendix A.1 on top of the input ionization maps (white – ionized, gray – neutral). For our mock observations we use a number density of n_gal = 0.002/cMpc³ (corresponding to a detection limit of M_UV<−17.8 at z = 8) and a Lyα 5σ equivalent width limit of ≥30 Å for the faintest galaxies. We see the bubble edges are recovered well at all neutral fractions.