Condition-based monitoring of rotating machines requires robust features for accurate fault diagnosis, which is indeed directly linked to the quality of the features extracted from the signals. This is especially true for vibration data, whose quasi-stationary nature implies that the quality of frequency domain extracted features depends on the signal-to-noise ratio (SNR) condition, operating condition variations, and data segmentation. This paper presents a novel statistical spectral analysis, which leads to highly robust fault diagnosis with poor SNR conditions, different time-window segmentation, and different operating conditions. The amplitudes of spectral contents of the quasi-stationary time vibration signals are sorted and transformed into statistical spectral images. The sort operation leads to the knowledge of the empirical cumulative distribution function (ECDF) of the amplitudes of each frequency band. The ECDF provides a robust statistical information of the distribution of the amplitude under different SNR and operating conditions. Statistical metrics have been adopted for fault classification, by using the ECDFs obtained from the spectral images as fault features. By applying simple statistical metrics, it is possible to achieve fault diagnosis without classifier training, saving both time and computational costs. The proposed algorithm has been tested using a vibration data benchmark: comparison with state-of-the-art fault diagnosis algorithms shows promising results.