A method of enhancing a noisy frequency-domain channel estimate H in, including the steps of transforming the noisy frequency-domain channel estimate H in into a time-domain channel estimate h that is virtually over-sampled by an integer factor m; dividing the time-domain channel estimate h into m de-interleaved sub-vectors h i; suppressing noisy taps from each de-interleaved sub-vectors h i to form m noise suppressed sub-vectors h ' i; selecting a noise suppressed sub-vector h ' i 0 associated with a highest signal-to-noise ratio γ i; converting the selected noise suppressed sub-vector h ' i 0 into a noise suppressed preliminary frequency-domain channel estimate H '; and phase-adjusting the taps of the noise suppressed preliminary frequency-domain channel estimate H ' to form an enhanced frequency-domain channel estimate H out, wherein the suppressing step includes determining a respective position of a peak power tap for each interleaved sub-vector h i; setting (S15+S16) a respective noise threshold Tε i for each sub-vector h i; suppressing taps with a power that falls below the respective noise threshold Tε i and are positioned within a predetermined channel length L from the determined respective position as well as taps that are not positioned within the channel length L from the determined respective position, while leaving the peak power tap unchanged, wherein the dividing step includes forming m de-interleaved sub-vectors h i in accordance with: h i = h 0, i, h 1, i,..., h n, i,..., h N − 1, i, i = − max ⌈ m 2 ⌉ − 1,0,..., ⌊ m 2 ⌋ where h n, i = h mod i + nm, mN, n = 0,1,..., N − 1, N is the total number of taps in the frequency-domain channel estimate H in, ┌ ┐is the ceiling function, └ ┘is the floor function, wherein the transforming step includes: zero-padding the noisy frequency-domain channel estimate H in with ( m -1) ·N zeros, where N is the total number of taps in the noisy frequency-domain channel estimate H in; and performing an m·N point inverse discrete Fourier transformation on the zero-padded noisy frequency-domain channel estimate H in or includes: performing an N point inverse discrete Fourier transformation on the noisy frequency-domain channel estimate H in, where N is the total number of taps in the noisy frequency-domain channel estimate H in, to form a noisy time-domain channel estimate; padding m -1 zeros after each sample of the noisy time-domain channel estimate; and lowpass filtering the zero-padded noisy time-domain channel estimate. zero-padding the noisy frequency-domain channel estimate H in with ( m -1) ·N zeros, where N is the total number of taps in the noisy frequency-domain channel estimate H in; and performing an m·N point inverse discrete Fourier transformation on the zero-padded noisy frequency-domain channel estimate H in or includes: performing an N point inverse discrete Fourier transformation on the noisy frequency-domain channel estimate H in, where N is the total number of taps in the noisy frequency-domain channel estimate H in, to form a noisy time-domain channel estimate; padding m -1 zeros after each sample of the noisy time-domain channel estimate; and lowpass filtering the zero-padded noisy time-domain channel estimate.