Convergence of square tilings to the Riemann map

A well-known theorem of Rodin \& Sullivan, previously conjectured by Thurston, states that the circle packing of the intersection of a lattice with a simply connected planar domain $\Omega$ into the unit disc $\mathbb{D}$ converges to a Riemann map from $\Omega$ to $\mathbb{D}$ when the mesh size converges to 0. We prove the analogous statement when circle packings are replaced by the square tilings of Brooks et al.