Skip to main content
  • Research Article
  • Open access
  • Published:

GPS-based fine-scale mapping surveys for schistosomiasis assessment: a practical introduction and documentation of field implementation



Fine-scale mapping of schistosomiasis to guide micro-targeting of interventions will gain importance in elimination settings, where the heterogeneity of transmission is often pronounced. Novel mobile applications offer new opportunities for disease mapping. We provide a practical introduction and documentation of the strengths and shortcomings of GPS-based household identification and participant recruitment using tablet-based applications for fine-scale schistosomiasis mapping at sub-district level in a remote area in Pemba, Tanzania.


A community-based household survey for urogenital schistosomiasis assessment was conducted from November 2020 until February 2021 in 20 small administrative areas in Pemba. For the survey, 1400 housing structures were prospectively and randomly selected from shapefile data. To identify pre-selected structures and collect survey-related data, field enumerators searched for the houses’ geolocation using the mobile applications Open Data Kit (ODK) and MAPS.ME. The number of inhabited and uninhabited structures, the median distance between the pre-selected and recorded locations, and the dropout rates due to non-participation or non-submission of urine samples of sufficient volume for schistosomiasis testing was assessed.


Among the 1400 randomly selected housing structures, 1396 (99.7%) were identified by the enumerators. The median distance between the pre-selected and recorded structures was 5.4 m. A total of 1098 (78.7%) were residential houses. Among them, 99 (9.0%) were dropped due to continuous absence of residents and 40 (3.6%) households refused to participate. In 797 (83.1%) among the 959 participating households, all eligible household members or all but one provided a urine sample of sufficient volume.


The fine-scale mapping approach using a combination of ODK and an offline navigation application installed on tablet computers allows a very precise identification of housing structures. Dropouts due to non-residential housing structures, absence, non-participation and lack of urine need to be considered in survey designs. Our findings can guide the planning and implementation of future household-based mapping or longitudinal surveys and thus support micro-targeting and follow-up of interventions for schistosomiasis control and elimination in remote areas.

Trial registration ISRCTN, ISCRCTN91431493. Registered 11 February 2020,


Schistosomiasis is among the 20 neglected tropical diseases (NTDs) defined by the World Health Organization (WHO) and is endemic in 78 countries worldwide [1,2,3]. The causative agent of schistosomiasis is a blood fluke of the genus Schistosoma that infects more than 200 million people worldwide [4]. The parasite is transmitted through skin contact with freshwater containing the infectious larval stages, which are released by an aquatic intermediate host snail vector [5]. Over the past decades, great progress was made in schistosomiasis control [6, 7]. In their new roadmap for neglected tropical diseases 2021–2030, the WHO set the global elimination of schistosomiasis as a public health problem and the validated absence of infection in humans in 25 among the 78 endemic countries as targets for 2030 [1].

The sustained implementation of control interventions, including mass drug administration (MDA) as the cornerstone, often supplemented with snail control, educational measures, or improvements in the socio-economic standard and access to clean water and sanitation, has resulted in decreasing schistosomiasis prevalences in many countries [8,9,10,11]. In areas where prevalences decline, the focality and heterogeneity of schistosomiasis becomes more pronounced [8, 12,13,14]. Typically, within a few years of annual or biannual MDA across districts and countries, many communities reach low schistosomiasis prevalence and infection intensity [10, 15, 16]. Some areas, however, remain as pockets of high transmission with persistent or reoccurring high prevalence levels despite intense interventions [8, 12,13,14, 17,18,19,20].

To achieve the WHO elimination goals, the pronounced spatial and temporal heterogeneity in schistosomiasis elimination settings will need to be considered in future intervention planning [8, 14, 20,21,22,23].

Fine-scale mapping of schistosomiasis is considered an essential requirement to move from morbidity control towards interruption of transmission in endemic areas [23]. Detailed fine-scale mapping of infection and disease patterns can help to identify low-risk and high-risk areas of transmission at sub-district level and to micro-target interventions in line with infection levels for optimal treatment and resource allocation [10, 14, 24, 25]. Android tablet-based applications containing or being able to connect to high-resolution maps to geolocate individuals’ residences in remote resource poor areas have been suggested as attractive approaches for fine-scale disease mapping [26, 27].

Here, we provide a practical introduction and documentation of the strengths and shortcomings of Global Position System (GPS)-based household identification and participant recruitment using Android tablet-based applications for fine-scale schistosomiasis mapping in a remote area in Pemba, United Republic of Tanzania.


Study area

The fine-scale mapping approach described here is part of the cross-sectional baseline survey of the “SchistoBreak” study, a multi-year project that aims to investigate new tools and strategies for breaking schistosomiasis transmission [28]. Pemba and Unguja are the main islands of the Zanzibar Archipelago, a semi-autonomous part of the United Republic of Tanzania. Pemba, where our study is implemented, is divided into four districts: Micheweni, Wete, Chake Chake and Mkoani [29]. In total, the districts in Pemba contain 129 small administrative areas, which are called shehias [30]. The study area of the SchistoBreak project, where the here reported fine-scale mapping survey is part of, covers a total of 20 shehias in the north of Pemba. Among them, 12 shehias are located in Micheweni and eight in Wete district, respectively. The average population size in our study shehias was ~ 3900 individuals per shehia in 2012, with a minimum of ~ 2200 and a maximum of ~ 8600 inhabitants [31]. As indicated in Fig. 1, shehias in Pemba are mostly connected by tarmacked main roads. However, some rough road connections exist. Housing structures within the shehias are grouped in clusters, which are usually aligned along the main roads or located near unpaved streets or footpaths. Housing structures are mostly not characterized or identifiable by street names or house numbers. Urogenital schistosomiasis is endemic in the north of Pemba [32, 33]. Due to intense control and elimination efforts implemented across the Zanzibar islands over the past decades, the overall S. haematobium prevalence had reached a very low level of 3.4% in schoolchildren in 2020, constituted by many low-prevalence areas and a few remaining hotspot areas of transmission [8].

Fig. 1
figure 1

Study area of the SchistoBreak project in Pemba, Tanzania. The exemplary section shows the distribution of housing structures and main roads. The image base map (United Republic of Tanzania—Subnational administrative boundaries) was downloaded from OCHA services ( The data source is: Tanzania National Bureau of Statistics/UN OCHA ROSA. The data are published under the following license: Creative Commons Attribution for Intergovernmental Organisations [CC BY-IGO; (]. Additionally, we received written permission to use and adapt the data from the UN Office for the Coordination of Humanitarian Affairs (OCHA). Additional shapefiles for the map (shehia boundaries and housing structures) were provided by the Zanzibar Commission for Lands to the Zanzibar Neglected Diseases Program

Design of cross-sectional household survey

In the baseline survey of the SchistoBreak study, we aimed to assess the apparent prevalence of S. haematobium infections in the 20 study shehias in the north of Pemba. To gain a clear picture of the schistosomiasis distribution we conducted a cross-sectional community-based household survey from November 2020 until February 2021, based on randomly selected housing structures. For fine-scale mapping, we aimed to include 50 households with an average of five household members per shehia, and hence a total of 250 individuals per shehia. Details of the sample size calculations are provided in our published study protocol [28]. Eligible to participate in the household survey were all individuals aged 4 years or older living in the study shehias. Accounting for an estimated 30% dropout due to housing structures not being residential houses or refusal of household members to participate in the study, we randomized 70 housing structures per shehia. Our survey team consisted of seven field-enumerators. Each of them visited 10 housing structures per shehia. Shehias were visited for 3 subsequent days to allow enough time for revisiting houses in case the inhabitants were absent on Day 1 and/or for delayed urine submission on Day 2 or Day 3. Each enumerator was equipped with a tablet computer for navigation to their pre-selected housing structures using the offline navigation map application © 2020 MAPS.ME ( and for electronic capturing of participant data using the data collection software © 2020 Open Data Kit (ODK;

Randomization of housing structures

For the randomization of houses, we used shapefile data of housing structures that were provided by the Zanzibar Commission for Lands to the Zanzibar Neglected Diseases Program. Shape file data could have also been obtained from © OpenStreetMap contributors (OSM). The OSM house shapefile data are published under the Open Database License ( Housing structures are indicated as polygons in shapefiles. For randomization of housing structures, for navigating to the randomized structures, and for conducting the questionnaire, we used the statistical software R version 3.5.1 (, the geospatial analysis software Esri ArcMap version 10.6.1, the offline map application MAPS.ME, and the data collection software ODK, respectively.

For randomization of the 70 housing structures, we first imported the shapefile polygon data into ArcMap and used its feature-to-point tool to calculate the centroid points of all housing polygons. As the information about the location of the centroids was provided in meters, we converted the meters into latitude and longitude coordinates as this is required by MAPS.ME. Subsequently, we exported the centroids given as shapefile data from ArcMap into R and randomly selected 70 houses per shehia.

Stratification of housing structures per shehia and enumerator

After randomly selecting 70 houses per shehia in R, we imported the centroids of randomized housing structures back into ArcMap. With the selection-by-lasso tool of ArcMap, we stratified the centroids in equal-sized groups of 10 per each of the seven enumerators. Hereby, we took care that the 10 housing structures for each enumerator were as close in distance as possible to avoid the need for far-distance walking. We saved each group of housing structures per enumerator as a separate layer and shapefile and subsequently exported all single shapefiles to R. In R, we added a new variable “interviewer” and assigned a unique value to each group of housing structures per enumerator. With the new interviewer variable, we ensured to keep the group of housing structures in the next step, in which we combined all single shapefiles of housing groups into one large dataset. Subsequently, we saved the merged data containing the shehia name, the geolocation of all 70 randomized housing structures per shehia, and the group variable for seven groups with 10 housing structures each per shehia as csv-files and as shapefiles for display.

Import of grouped housing structures per enumerator into ODK

Once the randomization process was completed, each housing structure was given a unique identifier code (with the variable name: "houseID"), which together with the shehia name, and latitude and longitude coordinates, was exported from the csv-files into an ODK excel sheet. In order to limit the choices displayed in ODK Collect to only those that belonged to the selected enumerator's group of housing structures, we used the “choice_filter” option of ODK. An example for using the choice_filter option with regards to the shehia and the enumerator is presented in Table 1, line 1–3 and Table 2, line 1–12. The example in Tables 1 and 2 shows a model for two shehias and two enumerators, each taking care of one group of housing structures per shehia.

Table 1 Customizing the ODK excel sheet for data collection and combined use with MAPS.ME: the survey sheet
Table 2 Customizing the ODK excel sheet for data collection and combined use with MAPS.ME: the choices sheet

Combination of MAPS.ME and ODK for wayfinding to selected housing structures

The enumerators used the navigator app MAPS.ME to identify and find the way to each randomized housing structure in each shehia. We selected this app as it works offline and downloadable maps contain the housing structures and street data published on OSM. ODK and MAPS.ME apps were installed on the mobile devices (Samsung Galaxy Tab A tablets; Samsung Electronics, Seoul, South Korea) that were used by the enumerators. To connect ODK with MAPS.ME, we inserted into the ODK excel sheet the information presented in Table 1, line 3–8, and Table 2, line 5–15. As indicated in Table 1, line 3–5, the calculation of the geolocation is based on the houseID that is selected by the enumerator and presented in Table 2, line 5–12. Table 1, line 6 and the corresponding lines 13–15 in Table 2 serve to select different transport modes, e.g. walking, going by vehicle or by bicycle so that transport ways and transport time can be calculated by MAPS.ME. By opening the ODK questionnaire, choosing a houseID and selecting the transport mode on the terminal mobile device, a button saying “Click to open” appears as indicated in Table 1, line 7. By clicking the appearing button, the prior installed navigator app MAPS.ME opens and the geolocation of the selected houseID is automatically set as destination point and displays the way to go from the point of being to the housing structure of interest.

Community-based data collection

To allow enough time for household identification, participant recruitment, questionnaire interview and urine collection, a shehia was visited by the field-enumerators for three subsequent days. On Day 1, most or all of the 10 pre-selected housing structures per enumerator were identified in the shehia community, and the geolocation and type of housing structure was recorded. In case the house was inhabited, household members were invited to participate in the study. Once a present adult household member agreed to participate, a questionnaire interview was conducted and instructions for urine collection from all household members aged ≥ 4 years were provided. One urine collection container per eligible household member was distributed and labelled with a unique picture (e.g. with a boat, star, cat or smiley) linked to the name of the participant recorded on a paper sheet to avoid confusion of urine samples. On Day 2, written informed consent forms for each participant, either signed by the participant or, in case of children, by their parent or legal guardian, were collected by the enumerator together with the urine samples from the participants. For all children 12–17 years old, an additional assent form signed by the adolescents themselves was also collected. Day 3 served as a mop-up day, when urine samples that were not submitted on Day 1 or Day 2 were collected.

Data management and statistical analysis

All data were collected with ODK and transferred to the ODK Central Server hosted at Swiss Tropical and Public Health Institute (Swiss TPH). All data were cleaned and analyzed with R version 3.5.1. Anonymized data are available as Additional files 1, 2. To document the strengths and shortcomings of GPS-based household identification and participant recruitment for fine-scale schistosomiasis mapping in a remote area in Africa (and in our case Pemba island), we evaluated the results of the community-based household survey once the data collection was completed.

First, we calculated the number of identified housing structures that were no residential houses but other structures such as buildings under construction, schools or shops.

Second, to assess the concordance of the location of the pre-selected housing structure with the location of the housing structures identified in the shehia communities by the enumerator, we determined the median, minimum and maximum distance between the initially randomized centroids and the recorded geolocations across all shehias and per shehia.

Third, we determined the number of residential houses where inhabitants were not at home at any day of the visits or where all inhabitants refused to participate in the study. The dropout rate of housing structures was calculated by dividing the number of not included housing structures by the total number of initially selected housing structures.

Fourth, for households that agreed to participate in the study, we assessed the total number of household members per inhabited house and the number of individuals eligible to participate in the study. The dropout rate of individuals who were eligible but did not provide a urine sample of sufficient volume for schistosomiasis testing was calculated by dividing the number of participants who submitted a urine sample by the total number of eligible individuals, stratified by sex and age category (i.e. adults versus children, defined as individuals below the age of 18 years).

Finally, we determined the number of urine samples received from each household and calculated the percentage of households where all participants provided a urine sample of sufficient volume for testing.


Housing structures and participating households

Among the 1400 housing structures that were initially randomly selected to be surveyed in the 20 shehias, the enumerators were able to locate and identify 1396 housing structures (99.7%) (Fig. 2). Four housing structures were not identified due to inaccessibility. Among the 1396 housing structures correctly identified, 298 (21.3%) were no residential houses: 99 were buildings under constructions, 57 were broken housing structures, 22 were Islamic schools, 21 were mosques, 18 were shops, 15 were animal stalls and 66 were other buildings, such as public schools and offices. Among the 1098 residential houses, in 99 (9.0%) houses, inhabitants were absent during any of the visits of the enumerator and in 40 (3.6%) houses, inhabitants did not agree to participate in the household survey. Hence, the overall dropout rate of housing structures was 31.5% (441/1400).

Fig. 2
figure 2

Housing structures and households participating in the fine-scale mapping survey in Pemba, Tanzania

Distance between centroids and recorded geolocation of housing structures

As indicated in Fig. 3, the median distance between the randomized centroids of the housing structures and the geolocation recorded by the enumerator when arriving at the housing structure was 5.4 m (range: 0.2–940.3 m). While in the six shehias surveyed in the early study period (weeks 1–4) the median distance was 7.4 m (range: 0.3–940.3 m), in the 14 shehias surveyed in the later study period (week 4–12) the median distance was 4.8 m (range: 0.2–54.4 m).

Fig. 3
figure 3

Distance between pre-selected and recorded geolocations of housing structures in the mapping survey in Pemba, Tanzania, per shehia. The horizontal middle lines in each boxplot indicate the median distance per shehia, the boxes indicate the lower and upper quartile, the vertical lines indicate the minimum and maximum excluding outliers, and the dots indicate outliers (more than 1.5 times the interquartile range below the first or above the 3rd quartile). For shehia 1, 2 and 6 the maximum outliers are indicated separately. The red dashed line indicates the median distance between pre-selected and recorded geolocations across all 20 shehias

Among the 1396 identified housing structures, four recorded geolocations were more than 100 m away from the pre-selected geopoints, taken by four different enumerators. In general, for all seven enumerators there was no major difference in the distribution of distances between the preselected and recorded geopoints (Fig. 4).

Fig. 4
figure 4

Distance between pre-selected and recorded geolocations of housing structures in the mapping survey in Pemba, Tanzania, per enumerator

Study participants

In total, 5340 individuals aged 0–99 years lived in the 959 participating households. In average, one household consisted of six individuals (range: 1–16 inhabitants).

Among the 5340 individuals, 4599 (86.1%) were eligible to participate in the study since their age was ≥ 4 years. Among them, 2324 (50.5%) were female and 2275 (49.5%) were male, and 2364 (51.4%) individuals were adults aged ≥ 18 years and 2233 (48.6%) were children aged 4–17 years. Two participants did not report their age. A total of 713 participants did not provide a urine sample, resulting in a dropout rate of 15.5% (713/4599). The dropout rate was higher in males (19.8%, 450/2275) than females (11.3%, 263/2323) and in adults (18.6%, 440/2364) compared with children (12.2%, 272/2233).

The number of participating households per shehia ranged from 40 to 56. In six among the 20 study shehias, the number of 50 households we aimed to include was achieved or exceeded (Fig. 5A). The number of participating individuals per shehia ranged from 171 to 319. In five shehias, the number of 250 participants we aimed for was achieved or exceeded (Fig. 5B).

Fig. 5
figure 5

Number of households (A) and individuals (B) participating in the mapping survey in Pemba, Tanzania, per shehia

Among the 959 participating households, in 598 (62.4%) households all eligible members provided a urine sample of sufficient volume (≥ 10 ml) for laboratory examinations, in 199 (20.8%) households all inhabitants but one individual provided a sufficiently large urine sample, in 81 (8.4%) households all inhabitants but two provided a sufficiently large urine sample, and in 81 (8.4%) households all inhabitants but three or more individuals provided a sufficiently large urine sample (Fig. 6).

Fig. 6
figure 6

Provision of urine samples. Percentage of households with study participants providing a urine sample of sufficient volume (≥ 10 ml) for schistosomiasis testing in the fine-scale mapping survey of the SchistoBreak study in Pemba, Tanzania


In areas where schistosomiasis control progresses towards elimination, heterogeneity of transmission becomes more pronounced and fine-scale mapping of infection patterns is needed for effective micro-targeting of interventions at sub-district level [14, 23, 24, 28]. Novel tablet-based tools containing offline high-resolution maps offer new opportunities for disease mapping and household identification, particularly in remote areas where the distribution of housing structures do not follow an order such as street names or house numbers and where internet connectivity is poor [26].

Here, we aimed to provide a practical introduction and documentation of the strengths and shortcomings of tablet application-based household identification and participant recruitment for fine-scale schistosomiasis mapping in a rural area in Pemba, United Republic of Tanzania.

In our baseline mapping survey, almost all (99.7%) pre-selected housing structures were identified at the spot, in the communities, and located with a median distance of 5.4 m between the randomized centroids of the housing structures and the geolocation recorded by the enumerator. The computer-based randomization and offline navigation app-based geolocation of housing structures worked hence very well, even in remote villages in Pemba where “houses” may be small in size and located closely to each other and where internet connectivity is often limited or completely absent. The correct and exact identification of prospectively randomized housing structures is important for the accurate implementation of a community-based cross-sectional survey and also enables an easy follow-up of individuals or households for treatment or in longitudinal studies. Moreover, highly accurate small-scale spatial data can help to link household participant infection data with possible transmission sites of schistosomiasis, and hence to directly identify or model high-risk zones for infection and transmission that require special attention and targeted interventions. The location data can also be used to assess the distance to the next primary health care unit. A study conducted in rural Zambia showed that the preciseness of geographical cluster locations is crucial to link individuals to the correct health provider, such as primary health care units [34].

The dropout rate due to non-residential housing structures in our survey was 21.3%, absence of inhabitants was 9.0%, refused participation of households was 3.6%, and inability to provide a sufficiently large urine sample among all study participants was 15.5%. A very similar dropout rate of 24.1% due to uninhabited houses was observed in a household study based on spatial randomization from rural and urban regions of Cameroon [35]. The absence of inhabitants in our study was, however, considerably lower than in a comparable GPS-based household survey conducted in Haiti where the absence of inhabitants was 16.4% [36]. This lower rate in our study might be explained by the fact that our enumerators visited the residential houses in a shehia for a second day if no individuals were at home the first day. The refusal dropout rate of 3.6% was, however, the same in our study as in the GPS-based household study conducted in Haiti [36]. A considerably higher dropout rate of 29.0% due to absence and/or inability to provide a urine and stool sample was observed in a household survey in Yemen, which was conducted to assess the local schistosomiasis prevalence [37]. A reason for the higher dropout in this study might be that not only urine but also stool samples were collected, which is more embarrassing to produce, collect and submit. In summary, our results show that our originally estimated overall dropout rate of 30% due to housing structures not being residential houses or refusal of household members to participate in the study (versus the actually observed dropout rate of 31.5%) was realistic and can serve as a guiding value for future survey design planning and sample size calculations, provided that 2 or 3 days are planned for participant recruitment and sample collection. In case households are only visited on a single day, higher dropout rates need to be considered in schistosomiasis mapping surveys, due to the potential absence of household members, or their inability or unwillingness to provide a urine or stool sample straight away. This being said, for S. haematobium diagnosis, urine samples are ideally collected in the late morning or early afternoon, when egg excretion is highest [38]. If urine samples can be produced and kept for the enumerator at any time of the day or night, this might impact on the diagnostic results and hamper accurate disease mapping.

With our practical introduction to and documentation of geolocation-based household identification and participant recruitment, we provide a rigorous basis for future survey planning. Particular strengths of our suggested approach are that it is low-cost, precise and easy to implement in the field as the final ODK questionnaire on the mobile device includes an automated link to the offline navigation system MAPS.ME. Due to the fully offline procedure it can also be used in remote areas without access to stable internet connections. Second, if data of housing structures are available on OSM, the fine-scale mapping approach is easy to replicate in other countries and for other purposes, particularly with the hands-on explanation provided in this article.

However, there are also important limitations. First, the approach relies on the availability of up-to-date (housing) data uploaded on OSM. In Pemba, new housing structures appear and disappear within a short time, or might be “under construction” for several years as has also been observed in other environments [36]. Without current data that include all new housing structures, individuals living in these housing structures will not be included in the randomization and missed for disease mapping. Without current data that mark or deleted demolished housing structures, the dropout rate will increase. To complement a higher dropout of housing structures if such invalid housing structures still appear in the data, the sample size could be increased or replacement-housing structures could be included in the survey, as done in other studies [35, 39]. Second, our approach relies on the accessibility of the housing structures. For example, four housing structures could not be included in our survey due to inaccessibility, mainly as they were located in a military area. The data that provided the basis for the randomization of housing structures did not show the military area and therefore the area was not excluded prior to randomization. Hence, the success of our mapping approach is very much depending on the availability of up-to-date high-resolution data, clearly indicating zones with restricted access.


We show that the fine-scale mapping approach with the combination of ODK and an offline navigation application installed on tablet computers allows a very precise identification of pre-randomized housing structures. Dropouts due to non-residential houses, absence, non-participation and lack of urine need to be considered in sample size calculations. Our results can guide the planning and implementation of future household-based mapping or longitudinal surveys and thus support micro-targeting and follow-up of interventions for schistosomiasis control and elimination in remote areas.

Availability of data and materials

Most relevant data are within the manuscript and its additional information files. Geopoints of the housing structures’ centroids and geolocations recorded by the field enumerators cannot be shared to preserve spatial confidentiality.



Ethikkommission Nordwest- und Zentralschweiz


Mass drug administration


Neglected tropical disease


Open Data Kit


Open Street Map

Swiss TPH:

Swiss Tropical and Public Health Institute


World Health Organization


Zanzibar Health Research Institute


  1. WHO. Ending the neglect to attain the Sustainable Development Goals—a road map for neglected tropical diseases 2021–2030. Geneva: World Health Organization; 2020.

    Google Scholar 

  2. WHO. Schistosomiasis key facts 2020. Accessed 15 Dec 2021.

  3. GBD 2017 DALYs and HALE Collaborators. Global, regional, and national disability-adjusted life-years (DALYs) for 359 diseases and injuries and healthy life expectancy (HALE) for 195 countries and territories, 1990–2017: a systematic analysis for the Global Burden of Disease Study 2017. Lancet. 2018;392(10159):1859–922.

    Article  Google Scholar 

  4. Steinmann P, Keiser J, Bos R, Tanner M, Utzinger J. Schistosomiasis and water resources development: systematic review, meta-analysis, and estimates of people at risk. Lancet Infect Dis. 2006;6(7):411–25.

    Article  Google Scholar 

  5. Colley DG, Bustinduy AL, Secor WE, King CH. Human schistosomiasis. Lancet. 2014;383(9936):2253–64.

    Article  Google Scholar 

  6. Deol AK, Fleming FM, Calvo-Urbano B, Walker M, Bucumi V, Gnandou I, et al. Schistosomiasis—assessing progress towards the 2020 and 2025 global goals. N Engl J Med. 2019;381(26):2519–28.

    Article  Google Scholar 

  7. Webster JP, Molyneux DH, Hotez PJ, Fenwick A. The contribution of mass drug administration to global health: past, present and future. Philos Trans R Soc Lond B Biol Sci. 2014;369(1645):20130434.

    Article  Google Scholar 

  8. Trippler L, Ame SM, Hattendorf J, Juma S, Abubakar S, Ali SM, et al. Impact of seven years of mass drug administration and recrudescence of Schistosoma haematobium infections after one year of treatment gap in Zanzibar: repeated cross-sectional studies. PLoS Negl Trop Dis. 2021;15(2):e0009127.

    Article  CAS  Google Scholar 

  9. Rollinson D, Knopp S, Levitz S, Stothard JR, Tchuem Tchuenté LA, Garba A, et al. Time to set the agenda for schistosomiasis elimination. Acta Trop. 2013;128(2):423–40.

    Article  Google Scholar 

  10. King CH, Kittur N, Binder S, Campbell CH, N’Goran EK, Meite A, et al. Impact of different mass drug administration strategies for gaining and sustaining control of Schistosoma mansoni and Schistosoma haematobium infection in Africa. Am J Trop Med Hyg. 2020;103(1_Suppl):14–23.

    Article  Google Scholar 

  11. Zhou Y, Chen Y, Jiang Q. History of human schistosomiasis (bilharziasis) in China: from discovery to elimination. Acta Parasitol. 2021;66(3):760–9.

    Article  Google Scholar 

  12. Santos IGA, Bezerra LP, Cirilo TM, Silva LO, Machado JPV, Lima PD, et al. New epidemiological profile of schistosomiasis from an area of low prevalence in Brazil. Rev Soc Bras Med Trop. 2020;53:e20200335.

    Article  Google Scholar 

  13. Souza MDR, de Jesus DMS, Santos AHC, Lima S, Dos Santos A, Tavares DDS, et al. Risk clusters of Schistosoma mansoni infection in an endemic state of Brazil: space-time modelling and association with socio-economic and environmental factors. Trans R Soc Trop Med Hyg. 2021.

    Article  PubMed  PubMed Central  Google Scholar 

  14. Fornace KM, Fronterrè C, Fleming FM, Simpson H, Zoure H, Rebollo M, et al. Evaluating survey designs for targeting preventive chemotherapy against Schistosoma haematobium and Schistosoma mansoni across sub-Saharan Africa: a geostatistical analysis and modelling study. Parasit Vectors. 2020;13(1):555.

    Article  Google Scholar 

  15. Shen Y, Wiegand RE, Olsen A, King CH, Kittur N, Binder S, et al. Five-year impact of different multi-year mass drug administration strategies on childhood Schistosoma mansoni-associated morbidity: a combined analysis from the Schistosomiasis Consortium for Operational Research and Evaluation cohort studies in the Lake Victoria regions of Kenya and Tanzania. Am J Trop Med Hyg. 2019;101(6):1336–44.

    Article  Google Scholar 

  16. Phillips AE, Tohon Z, Dhanani NA, Sofo B, Gnandou I, Sidikou B, et al. Evaluating the impact of biannual school-based and community-wide treatment on urogenital schistosomiasis in Niger. Parasit Vectors. 2020;13(1):557.

    Article  Google Scholar 

  17. Camara Y, Sanneh B, Joof E, Sanyang AM, Sambou SM, Sey AP, et al. Mapping survey of schistosomiasis and soil-transmitted helminthiases towards mass drug administration in The Gambia. PLoS Negl Trop Dis. 2021;15(7):e0009462-e.

    Article  Google Scholar 

  18. Mawa PA, Kincaid-Smith J, Tukahebwa EM, Webster JP, Wilson S. Schistosomiasis morbidity hotspots: roles of the human host, the parasite and their interface in the development of severe morbidity. Front Immunol. 2021;12:635869.

    Article  CAS  Google Scholar 

  19. Wiegand RE, Mwinzi PNM, Montgomery SP, Chan YL, Andiego K, Omedo M, et al. A persistent hotspot of Schistosoma mansoni infection in a five-year randomized trial of praziquantel preventative chemotherapy strategies. J Infect Dis. 2017;216(11):1425–33.

    Article  Google Scholar 

  20. Knopp S, Ame SM, Person B, Hattendorf J, Rabone M, Juma S, et al. A 5-year intervention study on elimination of urogenital schistosomiasis in Zanzibar: parasitological results of annual cross-sectional surveys. PLoS Negl Trop Dis. 2019;13(5):e0007268.

    Article  Google Scholar 

  21. Menghistu HT, Mersha TT, Shumuye NA, Woldie BM, Redda YT, Hadush B, et al. Neglected tropical zoonotic diseases in Tigray region, northern Ethiopia: spatial distribution and trend analysis of rabies, tuberculosis, schistosomiasis, and visceral leishmaniasis in humans. Zoonoses Public Health. 2021;68(7):823–33.

    Article  Google Scholar 

  22. Hatherell H-A, Simpson H, Baggaley RF, Hollingsworth TD, Pullan RL. Sustainable surveillance of neglected tropical diseases for the post-elimination era. Clin Infect Dis. 2021;72(Suppl 3):S210–6.

    Article  Google Scholar 

  23. Tchuem Tchuenté LA, Stothard JR, Rollinson D, Reinhard-Rupp J. Precision mapping: an innovative tool and way forward to shrink the map, better target interventions, and accelerate toward the elimination of schistosomiasis. PLoS Negl Trop Dis. 2018;12(8):e0006563.

    Article  Google Scholar 

  24. Soares Magalhães RJ, Salamat MS, Leonardo L, Gray DJ, Carabin H, Halton K, et al. Geographical distribution of human Schistosoma japonicum infection in The Philippines: tools to support disease control and further elimination. Int J Parasitol. 2014;44(13):977–84.

    Article  Google Scholar 

  25. Monnier N, Barth-Jaeggi T, Knopp S, Steinmann P. Core components, concepts and strategies for parasitic and vector-borne disease elimination with a focus on schistosomiasis: a landscape analysis. PLoS Negl Trop Dis. 2020;14(10):e0008837.

    Article  CAS  Google Scholar 

  26. Fornace KM, Surendra H, Abidin TR, Reyes R, Macalinao MLM, Stresman G, et al. Use of mobile technology-based participatory mapping approaches to geolocate health facility attendees for disease surveillance in low resource settings. Int J Health Geogr. 2018;17(1):21.

    Article  Google Scholar 

  27. Wagenaar BH, Augusto O, Ásbjörnsdóttir K, Akullian A, Manaca N, Chale F, et al. Developing a representative community health survey sampling frame using open-source remote satellite imagery in Mozambique. Int J Health Geogr. 2018;17(1):37.

    Article  Google Scholar 

  28. Trippler L, Hattendorf J, Ali SM, Ame SM, Juma S, Kabole F, et al. Novel tools and strategies for breaking schistosomiasis transmission: study protocol for an intervention study. BMC Infect Dis. 2021;21(1):1024.

    Article  Google Scholar 

  29. Office of Chief Government Statistician Zanzibar. Zanzibar statistical abstract 2019. 2020.

  30. Office of Chief Government Statistician Zanzibar. Zanzibar in figures 2019. 2020.

  31. Office of Chief Government Statistician Zanzibar. 2012 Population and housing census: population distribution by administrative areas. 2013.

  32. Hatz C, Savioli L, Mayombana C, Dhunputh J, Kisumku UM, Tanner M. Measurement of schistosomiasis-related morbidity at community level in areas of different endemicity. Bull World Health Organ. 1990;68(6):777–87.

    CAS  PubMed  PubMed Central  Google Scholar 

  33. Savioli L, Mott KE. Urinary schistosomiasis on Pemba Island: low-cost diagnosis for control in a primary health care setting. Parasitol Today. 1989;5(10):333–7.

    Article  CAS  Google Scholar 

  34. Carter ED, Munos MK. Impact of imprecise household location on effective coverage estimates generated through linking household and health provider data by geographic proximity: a simulation study. Int J Health Geogr. 2021;20(1):38.

    Article  Google Scholar 

  35. Alcoba G, Ochoa C, Babo Martins S, Ruiz de Castañeda R, Bolon I, Wanda F, et al. Novel transdisciplinary methodology for cross-sectional analysis of snakebite epidemiology at national scale. PLoS Negl Trop Dis. 2021;15(2):e0009023.

    Article  Google Scholar 

  36. Wampler PJ, Rediske RR, Molla AR. Using ArcMap, Google Earth, and Global Positioning Systems to select and locate random households in rural Haiti. Int J Health Geogr. 2013;12:3.

    Article  Google Scholar 

  37. Sady H, Al-Mekhlafi HM, Mahdy MAK, Lim YAL, Mahmud R, Surin J. Prevalence and associated factors of schistosomiasis among children in Yemen: implications for an effective control programme. PLoS Negl Trop Dis. 2013;7(8):e2377.

    Article  Google Scholar 

  38. Scott JA, Stimmel CM. The regularity of egg output of Schistosoma haematobium. Tex Rep Biol Med. 1956;14(4):440–58.

    CAS  PubMed  Google Scholar 

  39. Kondo MC, Bream KDW, Barg FK, Branas CC. A random spatial sampling method in a rural developing nation. BMC Public Health. 2014;14:338.

    Article  Google Scholar 

Download references


We thank all household members who participated in our household survey from November 2020 until February 2021 and the field-enumerators in Pemba who conducted the household survey. Aurelio Di Pasquale is acknowledged for his repeated support with the ODK software and server at Swiss TPH. We are grateful for the information kindly provided by Mario Andrés Gronert Álvarez, which helped us to link our ODK questionnaires with MAPS.ME. Furthermore, we thank the Mapping Section, Department of Survey and Urban Planning, Zanzibar Commission for Lands for the shapefiles containing information about shehia borders and housing structures in the North of Pemba.


Funding for the study has been obtained from the Swiss National Science Foundation (SNSF) via the PRIMA grant (PR00P3_179753) of Stefanie Knopp. This funding source had no role in the design of this study and will not have any role during its execution, analyses, interpretation of the data, or decision to submit results.

Author information

Authors and Affiliations



LT, SaMA, JH and SK conceptualized the SchistoBreak study and SaMA and SK initiated the SchistoBreak study. LT performed the data curation, formal analysis and data visualization for the research article. SK acquired the funding for the SchistoBreak study. MNA contributed significantly to the data collection and implementation of the household survey. LT, JH and SK conceptualized the statistical methods of the study. LT, ShMA, SaMA, FK and SK were responsible for the management of the research activity planning. LT conducted the randomization of housing structures and the connection of ODK with MAPS.ME for the household survey. LT, ShMA, SaMA and SK mentored the core team and MNA, ShMA, SaMA and FK advised for local conditions. LT and SK have drafted the first version of the research article. All authors have made substantial contributions to the study, and have reviewed and approved the submitted version of the manuscript.

Corresponding author

Correspondence to Stefanie Knopp.

Ethics declarations

Ethics approval and consent to participate

The study protocol has been waived by the ethics committee of Switzerland (Ethikkommission Nordwest- und Zentralschweiz; EKNZ) on October 23, 2019 (Req-2019-00951) and has been approved by the ethics committee of Zanzibar (Zanzibar Health Research Institute; ZAHRI) on December 13, 2019 (ZAHREC/02/November/2019/16). All study participants submitted an informed consent form signed by a parent or legal guardian in case of participating children and adolescents, or signed by the participant in case of participating adults. The study was prospectively registered at ISRCTN (ISCRCTN91431493).

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Supplementary Information

Additional file 1.

Supplementary data that support the findings of this study from Nov 2020 until Feb 2021.

Additional file 2.

Supplementary data dictionary that explains the data collected from Nov 2020 until Feb 2021.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Trippler, L., Ali, M.N., Ame, S.M. et al. GPS-based fine-scale mapping surveys for schistosomiasis assessment: a practical introduction and documentation of field implementation. Infect Dis Poverty 11, 8 (2022).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: