We propose a robust method for estimating the time-varying spectrum of a non-stationary random process. Our approach extends Thomson's powerful multiple window spectrum estimation scheme to the time-frequency and time-scale planes. The method refines previous extensions of Thomson's method through optimally concentrated window and wavelet functions and a statistical test for extracting chirping line components.