Neutron-deficient 177−185 Hg isotopes were studied using in-source laser resonance-ionization spectroscopy at the CERN-ISOLDE radioactive ion-beam facility, in an experiment combining different detection methods tailored to the studied isotopes. These include either α-decay tagging or Multireflection Time-of-Flight gating to identify the isotopes of interest. The endpoint of the odd-even nuclear shape staggering in mercury was observed directly by measuring for the first time the isotope shifts and hyperfine structures of 177−180 Hg. Changes in the mean-square charge radii for all mentioned isotopes, magnetic dipole and electric quadrupole moments of the odd-A isotopes and arguments in favor of I = 7/2 spin assignment for 177,179 Hg were deduced. Experimental results are compared with Density Functional Theory (DFT) and Monte-Carlo Shell Model (MCSM) calculations. DFT calculations with several Skyrme parameterizations predict a large jump in the charge radius around the neutron N = 104 mid shell, with an odd-even staggering pattern related to the coexistence of nearly-degenerate oblate and prolate minima. This near-degeneracy is highly sensitive to many aspects of the effective interaction, a fact that renders perfect agreement with experiment out of reach for current functionals. Despite this inherent difficulty, the SLy5s1 and a modified UNEDF1 SO parameterization predict a qualitatively correct staggering that is off by two neutron numbers. MCSM calculations of states with the experimental spins and parities show good agreement for both electromagnetic moments and the observed charge radii. A clear mechanism for the origin of shape staggering within this context is identified: a substantial change in occupancy of the proton πh 9/2 and neutron νi 13/2 orbitals.