Multiband homogenization of metamaterials in real-space: Higher-order nonlocal models and scattering at external surfaces