Analyzing the spectral and spatial characteristics of Hyperspectral Imaging (HSI) in a three-dimensional space is a challenging task. Recently, there have been developments in 3D feature extraction methods based on tensor decomposition, which allow for the effective utilization of both global and local information in HSI. These methods also explore the inherent low-rank properties of HSI through tensor decomposition. In this paper, we propose a new approach called variable randomized T-product decomposition (Vrt-SVD), which is a variation of Tensor Singular Spectral Analysis. The goal of this approach is to improve the efficiency of tensor methods for feature extraction and reduce artifacts of image processing. By using a randomized algorithm based on the variable t-SVD, we are able to capture both global and local spatial and spectral information in HSI efficiently, which enables us to explore its low-rank characteristics. To evaluate the effectiveness of the extracted features, we use a Support Vector Machine (SVM) classifier to assess the accuracy of image classification. By conducting numerous numerical experiments, we provide strong evidence to show that the proposed method outperforms several advanced feature extraction techniques.