As a recent approach for time series analysis, singular spectrum analysis (SSA) has been successfully applied for feature extraction in hyperspectral imaging (HSI), leading to increased accuracy in pixel-based classification tasks. However, one of the main drawbacks of conventional SSA in HSI is the extremely high computational complexity, where each pixel requires individual and complete singular value decomposition (SVD) analyses. To address this issue, a fast implementation of SSA (F-SSA) is proposed for efficient feature extraction in HSI. Rather than applying pixel-based SVD as conventional SSA does, the fast implementation only needs one SVD applied to a representative pixel, i.e., either the median or the mean spectral vector of the HSI hypercube. The result of SVD is employed as a unique transform matrix for all the pixels within the hypercube. As demonstrated in experiments using two well-known publicly available data sets, almost identical results are produced by the fast implementation in terms of accuracy of data classification, using the support vector machine (SVM) classifier. However, the overall computational complexity has been significantly reduced.