The mammalian liver is composed of repeating hexagonal units termed lobules. Spatially resolved single-cell transcriptomics revealed that about half of hepatocyte genes are differentially expressed across the lobule, yet technical limitations impeded reconstructing similar global spatial maps of other hepatocyte features. Here, we show how zonated surface markers can be used to sort hepatocytes from defined lobule zones with high spatial resolution. We apply transcriptomics, miRNA array measurements and mass spectrometry proteomics to reconstruct spatial atlases of multiple zonated features. We demonstrate that protein zonation largely overlaps with mRNA zonation, with the periportal HNF4α as an exception. We identify zonation of miRNAs such as miR-122, and inverse zonation of miRNAs and their hepatocyte target genes, highlighting potential regulation of protein levels through zonated mRNA degradation. Among the targets we find the pericentral Wnt receptors Fzd7 and Fzd8 and the periportal Wnt inhibitors Tcf7l1 and Ctnnbip1. Our approach facilitates reconstructing spatial atlases of multiple cellular features in the liver and other structured tissues.