Gathering very accurate spatially explicit data related to the distribution of mean annual precipitation is required when laying the groundwork for the prevention and mitigation of water-related disasters. In this study, four Bayesian maximum entropy (BME) models were compared to estimate the spatial distribution of mean annual precipitation of the selected areas. Meteorological data from 48 meteorological stations were used, and spatial correlations between three meteorological factors and two topological factors were analyzed to improve the mapping results including annual precipitation, average temperature, average water vapor pressure, elevation, and distance to coastline. Some missing annual precipitation data were estimated based on their historical probability distribution and were assimilated as soft data in the BME method. Based on this, the univariate BME, multivariate BME, univariate BME with soft data, and multivariate BME with soft data analysis methods were compared. The estimation accuracy was assessed by cross-validation with the mean error (ME), mean absolute error (MAE), and root mean square error (RMSE). The results showed that multivariate BME with soft data outperformed the other methods, indicating that adding the spatial correlations between multivariate factors and soft data can help improve the estimation performance.