Method for determining the direction of arrival ( θ), or the direction of arrival and of polarisation, of signals emanating from a plurality of radio sources by means of a single-polarisation or dual-polarisation array of spatially diverse sensors consisting of P receiving channels with P > 2 and it being possible for P to be less than the number of sources, when the degree of parsimony of the signals is such that locally, in the time/frequency representation thereof, the number of sources present is at most 2, and comprising the following steps: a) calculating the successive Discrete Fourier Transforms of the sampled received signal in order to obtain a P-vector time-frequency grid of the signal; each element of the grid containing a complex dimensional vector P, referred to as a measurement, the method being characterised in that it comprises the following additional steps: b) defining a division of the time-frequency grid into a set of windows; c) on each window, c1- calculating the covariance matrix of the measurements of said window, c2- counting the sources in order to identify a single-source or dual-source window, c3- determining the or the two largest eigenvalues, and the or the two associated eigenvectors of the covariance matrix; d) discretising the field of the directions of arrival ( θ) that are accessible to the array, e) generating directional vectors for the discretised directions of arrival ( θ) accessible to the array; f) associating the directional vectors with the windows, by comparing the magnitude of the projection of the directional vector with the vector subspace generated by the or the two eigenvectors of each considered window; g) for each set of windows associated with the same directional vector, obtaining the filtered directional vector from the corresponding source by calculating the eigenvector associated with the largest eigenvalue of a matrix constructed from the eigenvalues and eigenvectors of the covariance matrices of said windows. a) calculating the successive Discrete Fourier Transforms of the sampled received signal in order to obtain a P-vector time-frequency grid of the signal; each element of the grid containing a complex dimensional vector P, referred to as a measurement, the method being characterised in that it comprises the following additional steps: b) defining a division of the time-frequency grid into a set of windows; c) on each window, c1- calculating the covariance matrix of the measurements of said window, c2- counting the sources in order to identify a single-source or dual-source window, c3- determining the or the two largest eigenvalues, and the or the two associated eigenvectors of the covariance matrix; c1- calculating the covariance matrix of the measurements of said window, c2- counting the sources in order to identify a single-source or dual-source window, c3- determining the or the two largest eigenvalues, and the or the two associated eigenvectors of the covariance matrix; d) discretising the field of the directions of arrival ( θ) that are accessible to the array, e) generating directional vectors for the discretised directions of arrival ( θ) accessible to the array; f) associating the directional vectors with the windows, by comparing the magnitude of the projection of the directional vector with the vector subspace generated by the or the two eigenvectors of each considered window; g) for each set of windows associated with the same directional vector, obtaining the filtered directional vector from the corresponding source by calculating the eigenvector associated with the largest eigenvalue of a matrix constructed from the eigenvalues and eigenvectors of the covariance matrices of said windows.