A time sequence of raw radar data for a region of space is subdivided into a plurality of processing frames. The processing frames are subdivided into a plurality of processing cells and iteratively processed by selecting a single processing cell for processing, transforming the radar data of the processing cell to form transformed radar data in either the time domain or the Fourier domain. The transformed data is converted to a Power Spectrum Density Matrix in the case of the Fourier domain and a Time Space Correlation Matrix in the case of the time domain. This is smoothed and thresholded and then the clutter for the processing cell is estimated. Estimated local non-speckle clutter is estimated and removed from the transformed radar data, with the cleaned transformed radar data converted back to the time domain if required.