Low-temperature scanning tunneling microscopy has been employed to analyze the formation of quantum well states (QWS) in two-dimensional gold islands, containing between 50 and 200 atoms, on MgO thin films. The energy position and symmetry of the eigenstates are revealed from conductance spectroscopy and imaging. The majority of the QWS originates from overlapping Au 6p orbitals in the individual atoms and is unoccupied. Their characteristic is already reproduced with simple particle-in-a-box models that account for the symmetry of the islands (rectangular, triangular, or linear). However, better agreement is achieved when considering the true atomic structure of the aggregates via a density functional tight-binding approach. Based on a statistically relevant number of single-island data, we have established a correlation between the island geometry and the gap between the highest-occupied and the lowest-unoccupied molecular orbital in the finite-sized islands. The linear eccentricity is identified as a suitable descriptor for this relationship, as it combines information on both island size and island shape. Finally, the depth of the confinement potential is determined from the spatial extension of QWS beyond the physical boundaries of the Au islands. Our paper demonstrates how electron quantization effects can be analyzed in detail in metal nanostructures. The results may help elucidating the interplay between electronic and chemical properties of oxide-supported clusters as used in heterogeneous catalysis.