A method implemented by computer means for estimating interference on a radiofrequency system using a set of channels, CHi, said interference being caused by interferers of an interfering system using a set of transmission bands, TBi, each of said transmission bands extending on a plurality of contiguous channels of said set of channels, Wherein the method comprises: - Determining a set Ω of all possible configurations of occupation or non-occupation of said set of transmission bands, defined as a set of possible vectors Z k = [ Z 1,k,..., Z i,k,... Z I,k] satisfying at a time instant k a non-overlapping condition of said radiofrequency system, said non-overlapping condition corresponding to the fact that only one interferer i, among a set of I possible interferers, can be active at a same time k on each channel of said set of channels and forming, with contiguous channels, a transmission band i, - Obtaining measurements X 1,..., X k,..., X K of occupation of at least a part of said set of channels, at respective time instants k: 0 < k ≤ K, where K defines a given observation time window, - Performing probabilities calculations so as to determine, for each transmission band, an estimated activation rate τ, on the basis of said measurements X 1,..., X k,..., X K, said estimated activation rate τ corresponding to an occupation rate of a transmission band i by an interferer within said given observation time window, wherein said probabilities calculations follow an expectation-maximum approach which iteratively approximates a maximum likelihood solution, according to two steps at each iteration t: a) For a fixed estimate τ ( t), define the expectation Q ( τ|τ ( t)) such that Q τ | τ t = E Z ∈ Ω | X, τ t log P X, Z | τ b) Find the estimate τ ( t +1) that maximizes Q ( τ|τ ( t)): τ t + 1 = argmax Q τ | τ t, and wherein successive t iterations of steps a) and b) lead to a converging solution giving said estimated activation rate τ ( t +1) for each transmission band. - Determining a set Ω of all possible configurations of occupation or non-occupation of said set of transmission bands, defined as a set of possible vectors Z k = [ Z 1,k,..., Z i,k,... Z I,k] satisfying at a time instant k a non-overlapping condition of said radiofrequency system, said non-overlapping condition corresponding to the fact that only one interferer i, among a set of I possible interferers, can be active at a same time k on each channel of said set of channels and forming, with contiguous channels, a transmission band i, - Obtaining measurements X 1,..., X k,..., X K of occupation of at least a part of said set of channels, at respective time instants k: 0 < k ≤ K, where K defines a given observation time window, - Performing probabilities calculations so as to determine, for each transmission band, an estimated activation rate τ, on the basis of said measurements X 1,..., X k,..., X K, said estimated activation rate τ corresponding to an occupation rate of a transmission band i by an interferer within said given observation time window, a) For a fixed estimate τ ( t), define the expectation Q ( τ|τ ( t)) such that Q τ | τ t = E Z ∈ Ω | X, τ t log P X, Z | τ b) Find the estimate τ ( t +1) that maximizes Q ( τ|τ ( t)): τ t + 1 = argmax Q τ | τ t, and wherein successive t iterations of steps a) and b) lead to a converging solution giving said estimated activation rate τ ( t +1) for each transmission band.