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. <EOS>