Abstract: We propose an improved image simulation procedure for atomic-resolution annular dark-field scanning transmission electron microscopy (STEM) based on the multislice formulation, which takes thermal diffuse scattering fully into account. The improvement with regard to the classical frozen phonon approach is realized by separating the lattice configuration statistics from the dynamical scattering so as to avoid repetitive calculations. As an example, the influence of phonon scattering on the image contrast is calculated and investigated. STEM image simulation of crystals can be applied with reasonable computing times to problems involving a large number of atoms and thick or large supercells.