High-altitude platform stations (HAPS) are pivotal in next-generation wireless networks for reducing core network burdens and enabling cost-effective communication. In this article, we propose a spherical stochastic geometry-based analytical framework for the coverage performance evaluation of HAPS networks. Considering the significant influence of directional antenna gain on interference evaluation, we analyze coverage performance under a general channel model that accommodates various beam patterns. Analytical expressions of the uplink and downlink coverage probabilities in cellular and cell-free networks are provided respectively and their accuracy are verified by Monte Carlo simulation. Furthermore, the influences of network-level and physical-level parameters on coverage probability are studied. Finally, several factors that align with the analytical framework of this article are discussed.