Several methods have been proposed in the past years to extract the Rayleigh wave ellipticity from horizontal-to-vertical spectral ratios of single station ambient noise recordings. The disadvantage of this set of techniques is the difficulty in clearly identifying and separating the contribution of higher modes. In most cases, only the fundamental mode of ellipticity can be identified. Moreover, it is generally difficult to correct for the energy of SH and Love waves present in the horizontal components of the ambient vibration wavefield. We introduce a new methodology to retrieve Rayleigh wave ellipticity using high-resolution frequency-wavenumber array analysis. The technique is applied to the three components of motion and is based on the assumption that an amplitude maximum in the f-k cross-spectrum must represent the true power amplitude of the corresponding signal. In the case of Rayleigh waves, therefore, the ratio between maxima obtained from the horizontal (radial-polarized) and vertical components of motion will also represent the frequency-dependent ellipticity function. Consequently, if we can identify the Rayleigh dispersion curves of several modes on the f-k plane, then the corresponding modal ellipticity patterns can also be separated and extracted. To test the approach, synthetic and real data sets were processed. In all tested cases, a reliable estimation of segments of the fundamental mode ellipticity was obtained. The identification of higher modes is possible in most cases. The quality of results depends on the selected array geometry and the signal-to-noise ratio, with a major improvement achieved by increasing the number of receivers employed during the survey. An experiment conducted in the town of Visp (Switzerland) allowed the retrieval of portions of ellipticity curves up to the second Rayleigh higher mode, using two concentric circular array configurations of 14 and 11 receivers eac