ABSTRACT Unsupervised hyperspectral image (HSI) band selection methods have been attracting ever-increasing attention. However, the local structural features captured by most of the existing methods suffer from a certain degree of risk of feature ambiguity. Moreover, these methods typically have difficulty preserving the latent manifold structure of the HSI in selected subsets of the bands. To address these problems, firstly, the graph approximate representation learning (GARL) model is proposed in this article by incorporating graph-regularized sparse coding into self-representation learning to preserve the local manifold structure information of HSI in the selected subset of bands. Meanwhile, to alleviate the ambiguity of the local graph structures, we design a local-global feature relationship expression model to construct the latent long-distance contextual connectivity (LDCC) graphs between the obtained local graph structures. Then, to preserve HSI’s local and non-local manifold structure information as much as possible in the selected subset of bands, a multigraph approximate representation learning (MGARL) model is proposed in this article by incorporating the obtained LDCC graph into the proposed GARL model. Next, we design a solution method to solve the proposed formulation. Finally, extensive experimental results on four datasets reveal the promising performance of MGARL over some state-of-the-art methods.