Based on recent results of applying graph signal processing (GSP) to narrowband angle of arrival estimation for uniform linear arrays, we generalise the analysis to the case of arrays with elements placed arbitrarily in three dimensional space. We comment on the selection of the adjacency matrix, analyse how this new approach compares to the multiple signal classification (MUSIC) algorithm, and provide an efficient implementation. We demonstrate that the GSP approach can perform as well as the MUSIC algorithm in terms of accuracy and computational cost. Simulations indicate that the proposed GSP approach avoids the severe performance degradation with which MUSIC is associated at low signal to noise ratios.