In this work, we explore the relevant methodology for the investigation of interacting systems with contact interactions, and we introduce a class of zonal estimators for path-integral Monte Carlo methods, designed to provide physical information about limited regions of inhomogeneous systems. We demonstrate the usefulness of zonal estimators by their application to a system of trapped bosons in a quasiperiodic potential in two dimensions, focusing on finite temperature properties across a wide range of values of the potential. Finally, we comment on the generalization of such estimators to local fluctuations of the particle numbers and to magnetic ordering in multi-component systems, spin systems, and systems with nonlocal interactions.