Schistosomiasis is a parasitic disease in Tanzania affecting over 50% of the population. Current control strategies involve mass drug administration (MDA) campaigns at the district level, which have led to problems of over- and under-treatment in different areas. WHO guidelines have called for more targeted MDA to circumvent these problems, however a scarcity of prevalence data inhibits decision makers from prioritizing sub-district areas for MDA. This study demonstrated how geostatistics can be used to inform planning for targeted MDA. Geostatistical sub-district (ward-level) prevalence estimates were generated through combining a zero-inflated poisson model and kriging approach (regression kriging). To make predictions, the model used prevalence survey data collected in 2021 of 17,400 school children in six regions of Tanzania, along with several open source ecological and socio-demographic variables with known associations with schistosomiasis. The model results show that regression kriging can be used to effectively predict the ward level parasite prevalence of the two species of Schistosoma endemic to the study area. Kriging was found to further improve the regression model fit, with an adjusted R-squared value of 0.51 and 0.32 for intestinal and urogenital schistosomiasis, respectively. Targeted treatment based on model predictions would represent a shift in treatment away from 193 wards estimated to be over-treated to 149 wards that would have been omitted from the district level MDA. Geostatistical models can help to support NTD program efficiency and reduce disease transmission by facilitating WHO recommended targeted MDA treatment through provision of prevalence estimates where data is scarce.