A method is proposed to compute the response of highly heterogeneous structures by constructing homogenized models when no scale separation can be assumed. The technique is based on an extended homogenization scheme where the averaging operators are replaced by linear filters to cut off fine scale fluctuations, while maintaining locally varying mechanical fields with larger wavelength. As a result, the constitutive law at an intermediate scale, called mesoscale, arises to be naturally nonlocal by construction, where the kernel operator is fully constructed by computations on the unit cell associated with the microstructure.