We present a new robust method for identifying three dynamically distinct regions in a stratified turbulent flow, which we characterise as quiescent flow, intermittent layers, and turbulent patches. The method uses the cumulative filtered distribution function of the local density gradient to identify each region. We apply it to data from direct numerical simulations of homogeneous stratified turbulence, with unity Prandtl number, resolved on up to 8192×8192×4096 grid points. In addition to classifying regions consistently with contour plots of potential enstrophy, our method identifies quiescent regions as regions where /νN 2 ∼ O(1), layers as regions where /νN 2 ∼ O(10), and patches as regions where /νN 2 ∼ O(100). Here is the dissipation rate of turbulence kinetic energy, ν is the kinematic viscosity, and N is the (overall) buoyancy frequency. By far the highest local dissipation and mixing rates, and the majority of dissipation and mixing, occur in patch regions, even when patch regions occupy only 5% of the flow volume. We conjecture that treating stratified turbulence as an instantaneous assemblage of these different regions in varying proportions may explain some of the apparently highly scattered flow dynamics and statistics previously reported in the literature.