Photoacoustic tomography (PAT) has emerged as a promising biomedical imaging technique. The combination of optical contrast and ultrasound spatial resolution in photoacoustic tomography overcomes the limitations of optical scattering, enabling clear imaging of tissue structures. However, achieving high-resolution photoacoustic images typically requires a large number of sensor detection elements for sufficient angular coverage. This demand for extensive data acquisition and processing raises concerns about efficiency and system complexity. While sparse sampling strategies can improve efficiency, preserving detailed structural information becomes challenging with a minimal number of detectors. To address the challenges of sparse sampling, compressed sensing (CS) techniques have been successfully applied for image reconstructions in 2D and 3D photoacoustic embodiments. In this context, we propose a joint graph regularization group sparse dictionary and total variational regularization (GRGS-TV) algorithm based on our previous work of a group sparse dictionary. It preserves structured information and geometric relationships among dictionary atoms. Moreover, TV regularization effectively preserves edge structures while exhibiting a certain degree of robustness and flexibility. Numerical simulations and in vivo experiments on mice validate the effectiveness of this method in improving photoacoustic image quality and suppressing artifacts. Comparative evaluations against other algorithms show enhanced performance in terms of image reconstruction evaluation indices. This innovative approach holds promise for advancing photoacoustic imaging in biomedical research and clinical diagnostics.