Deep learning (DL) exhibits commendable performance in HSI classification because of its powerful feature expression ability. Siamese neural network further improves the performance of DL models by learning similarities withinclass and differences between-class from sample pairs. However, there are still some limitations in siamese neural network. On the one hand, siamese neural network usually needs large number of negative pair samples in training process, leading to computing overhead. On the other hand, current models may lack interpretability because of complex network structure. To overcome the above limitations, we propose a spectral-spatial siamese neural network with Bag-of-Features (S3BoF) for HSI classification. Firstly, we use a siamese neural network with 3D and 2D convolutions to extract the spectral-spatial features. Secondly, we introduce stop-gradient operation and prediction head structure to make the siamese neural network work without negative pair samples, thus reducing computational burden. Thirdly, a Bag-of-Features (BoF) learning module is introduced to enhance the model interpretability and feature representation. Finally, a symmetric loss and a cross entropy loss are respectively used for contrastive learning and classification. Experiments results on four common hyperspectral data sets indicated that S3BoF performs better than the other traditional and stateof-the-art deep learning HSI classification methods in terms of classification accuracy and generalization performance, with improvements in terms of OA around 1.40%-30.01%, 0.27%-