Bioluminescence tomography (BLT) is one kind of noninvasive optical molecular imaging technology, widely used to study molecular activities and disease progression inside live animals. By combining the optical propagation model and inversion algorithm, BLT enables three-dimensional imaging and quantitative analysis of light sources within organisms. However, challenges like light scattering and absorption in tissues, and the complexity of biological structures, significantly impact the accuracy of BLT reconstructions. Here, we propose a dictionary learning method based on K-sparse approximation and Orthogonal Procrustes analysis (KSAOPA). KSAOPA uses an iterative alternating optimization strategy, enhancing solution sparsity with k-coefficients Lipschitzian mappings for sparsity(K-LIMAPS) in the sparse coding stage, and reducing errors with Orthogonal Procrustes analysis in the dictionary update stage, leading to stable and precise reconstructions. We assessed the method performance through simulations and in vivo experiments, which showed that KSAOPA excels in localization accuracy, morphological recovery, and in vivo applicability compared to other methods.