The microscopic valence electron density determines the optical, electronic, structural and thermal properties of materials. However, current techniques for measuring this electron charge density are limited: for example, scanning tunnelling microscopy is confined to investigations at the surface, and electron di raction requires very thin samples to avoid multiple scattering 1 . Therefore, an optical method is desirable for measuring the valence charge density of bulk materials. Since the discovery of high-harmonic generation (HHG) in solids 2 , there has been growing interest in using HHG to probe the electronic structure of solids 3-11 . Here, using single-crystal MgO, we demonstrate that high-harmonic generation in solids is sensitive to interatomic bonding. We find that harmonic e ciency is enhanced (diminished) for semi-classical electron trajectories that connect (avoid) neighbouring atomic sites in the crystal. These results indicate the possibility of using materials' own electrons for retrieving the interatomic potential and thus the valence electron density, and perhaps even wavefunctions, in an all-optical setting.High-harmonic generation (HHG) in bulk crystal is attributed to the sub-cycle electronic motion driven by an intense laser field [2][3][4][5][6][7][8][9][10][11] . There has been a growing interest in utilizing HHG to probe the electronic structure of solids 8,9,11 . Vampa et al. reconstructed the momentum-dependent bandgap of ZnO along the -M direction using HHG from a two-colour driving field 11 . Luu et al. retrieved the energy dispersion of the lowest conduction band of SiO 2 assuming that the harmonics are produced by the intraband currents 8 . The dependence of solid-state HHG on the coupling of multiple electronic bands has also been identified with the production of even harmonics in GaSe 9 and the emergence of a second plateau in rare-gas solids 12 . These findings show the possibility of using solidstate HHG to probe the electronic band structures in solids, but the analyses are so far limited to one dimension. For a complete electronic structure, it is desirable to exploit the microscopic process to measure the periodic potential in three dimensions (real space). This is analogous to tomographic imaging of a molecule, where the three-dimensional spatial information (that is, orbital wavefunction) of the target molecule is extracted [13][14][15] . Those measurement techniques are based critically on the dependence of HHG efficiency on molecular alignment with respect to the laser field 16 .In this letter, we demonstrate the strong sensitivity of HHG to the atomic-scale structure in the cubic wide-bandgap crystal MgO. First, using a linearly polarized field, we measure a highly anisotropic angular distribution in high-harmonic signal-despite the isotropic linear and weakly anisotropic nonlinear optical properties of the cubic crystal in the perturbative regime 17 . Second, we observe a strong ellipticity dependence of the HHG yield similar to the gas-phase HHG 18 for small elliptic...