Accuracy and reliability are determining factors when selecting where to place a global navigation satellite system (GNSS) receiver for landslide monitoring and early warning in real time. Previous studies have revealed that surface occlusion models (SOMs) are beneficial to reduce the multipath error of satellite observation when the satellite line of sight is near the ground and improve the accuracy and reliability of GNSS monitoring. This paper investigates fitting functions to determine appropriate SOMs for GNSS stations located in different spatial environments. Based on their similarity to terrain patterns, Gaussian, Fourier, and sum of sine functions were chosen to build SOMs for the rover sites with different occlusions; the performance of fitting accuracy and application on the positioning were analyzed. The results indicate that the three functions present a similar performance considering the goodness of fit, but the Fourier and sum of sine functions can significantly improve the efficacy of SOMs and enhance performance in positioning for landslide monitoring. Therefore, these two functions can effectively support accurate and reliable GNSS landslide monitoring in real time.