Skip to main content

The application of drones for mosquito larval habitat identification in rural environments: a practical approach for malaria control?

Abstract

Background

Spatio-temporal trends in mosquito-borne diseases are driven by the locations and seasonality of larval habitat. One method of disease control is to decrease the mosquito population by modifying larval habitat, known as larval source management (LSM). In malaria control, LSM is currently considered impractical in rural areas due to perceived difficulties in identifying target areas. High resolution drone mapping is being considered as a practical solution to address this barrier. In this paper, the authors’ experiences of drone-led larval habitat identification in Malawi were used to assess the feasibility of this approach.

Methods

Drone mapping and larval surveys were conducted in Kasungu district, Malawi between 2018 and 2020. Water bodies and aquatic vegetation were identified in the imagery using manual methods and geographical object-based image analysis (GeoOBIA) and the performances of the classifications were compared. Further, observations were documented on the practical aspects of capturing drone imagery for informing malaria control including cost, time, computing, and skills requirements. Larval sampling sites were characterized by biotic factors visible in drone imagery and generalized linear mixed models were used to determine their association with larval presence.

Results

Imagery covering an area of 8.9 km2 across eight sites was captured. Larval habitat characteristics were successfully identified using GeoOBIA on images captured by a standard camera (median accuracy = 98%) with no notable improvement observed after incorporating data from a near-infrared sensor. This approach however required greater processing time and technical skills compared to manual identification. Larval samples captured from 326 sites confirmed that drone-captured characteristics, including aquatic vegetation presence and type, were significantly associated with larval presence.

Conclusions

This study demonstrates the potential for drone-acquired imagery to support mosquito larval habitat identification in rural, malaria-endemic areas, although technical challenges were identified which may hinder the scale up of this approach. Potential solutions have however been identified, including strengthening linkages with the flourishing drone industry in countries such as Malawi. Further consultations are therefore needed between experts in the fields of drones, image analysis and vector control are needed to develop more detailed guidance on how this technology can be most effectively exploited in malaria control.

Background

Malaria cases in Africa have reduced by over half in the last two decades making transmission more heterogeneous. This has led to a growth of studies applying spatial and temporal analyses to determine where and when remaining transmission foci exist [1], and a focus on how new and existing control methods can be best utilized to reduce this residual transmission [2,3,4].

The geographical spread and extent of malaria transmission is limited by the seasonally-driven mosaic of water bodies available for female mosquitoes in which to lay their eggs. The ecology of preferred breeding grounds for mosquito oviposition vary both within and between species. For example, two of the main sibling species of the Anopheles gambiae complex, Anopheles gambiae and Anopheles arabiensis, tend to be found in transient, sunlit, small pools whereas Anopheles funestus is associated with more permanent, larger vegetated water [5]. At the micro-geographic scale, the presence of mosquito larvae may differ over the course of just a few metres [6,7,8]. Biotic and abiotic factors such as the presence of specific types of vegetation, microbiota, predators, algal density, shade, and water depth influence larval development.

Mosquito larval populations are fixed in space for the duration of their development to adulthood. Typically, eggs hatch into larvae within 2–3 days of oviposition and take 5–10 days to metamorphosize into pupae, although the speed of this process is highly dependent on temperature and the availability of nutrients [9, 10]. One method of controlling diseases transmitted by mosquitoes is to reduce the population by reducing the availability of oviposition sites and/or reduce the likelihood that resulting larvae develop into the adult stage [11]. Larval source management (LSM) involves the environmental, biological or chemical manipulation of the environment in which mosquitoes are present for the purpose of targeting the immature, aquatic stages of the mosquito and hence reducing the adult mosquito population. In the early days of mosquito control, an aggressive approach to searching and removing mosquito breeding sites was successful at reducing (and even eliminating) disease, with historical examples including its use during the construction of the Panama Canal in the early twentieth century, and its role in the elimination of An. gambiae in Brazil by 1940 [12]. In sub-Saharan Africa, LSM was responsible for large reductions in malaria incidence in Zambia copper mines between 1929 and 1949 [13]. LSM is, however, a labour-intensive exercise and following the introduction of IRS by DDT in the 1950s and subsequently the development of ITNs in the 1990s, it fell out of favour as a viable control option, particularly in Africa where the long rainy seasons produce countless sites for Anopheles development [13]. As such, LSM is currently only recommended as a complementary vector control intervention to bed nets and IRS to target residual transmission and as a method of combating insecticide resistance [14]. While its value is acknowledged by the World Health Organization (WHO) and national malaria control programmes (NMCPs) there are several barriers to its widespread implementation.

The primary barrier to implementing LSM is the issue of determining where and when the intervention should be implemented. In rural settings, the WHO recommend the application of LSM in areas where there is high coverage of long-lasting insecticidal nets (LLINs), evidence of outdoor biting and/or insecticide resistance and where larval sources are ‘few, fixed and findable’. Despite the lack of a clear definition of what can be considered ‘few’, ‘fixed’ or ‘findable’, this has led to many considering LSM to be impractical in rural areas with diffuse seasonal larval habitats. The perception of these terms may evolve as technology and processes for implementing LSM advance. This paper focuses on challenging the ‘findable’ component of this trio of conditions.

Geospatial technology is rapidly evolving and what now constitutes as ‘findable’ may switch from less reliance on exhaustive ground-based searches to remotely sensed data. Drone mapping is being touted as at least equivalent (if not superior) to and more cost-effective than mapping larval habitat manually [15, 16] or using remotely sensed satellite imagery [17, 18]. While the latter can cover vast areas in a single day, images are often obscured by clouds and although very high-resolution commercial satellite imagery exists, the resolution (at best 30 cm) is still inferior to that obtained by drones (2–10 cm) with the time of image captured out of the data user’s control.

This paper explores the use of drones as a method for collecting very high resolution (< 10 cm), contemporary imagery of an area for the purposes of identifying larval habitat. The issues addressed include the process of capturing drone imagery (by who, how much, how often), processing the images to extract the required information (what software, image classification methods, computer processing requirements), collecting ‘ground-truth’ data (entomological sampling), and subsequently summarising this information into recommendations that can be used by the control program implementers to guide their LSM activities.

Methods

Drone image capturing

A series of image data capture exercises were conducted within Kasungu district, central Malawi in an area that has been designated by the Government of Malawi, in collaboration with UNICEF, as a ‘humanitarian drone testing corridor’ (Fig. 1). Authorization to conduct these flights was obtained from the Malawi Department of Civil Aviation. Malaria transmission occurs all year round in this area, with parasite prevalence in children between 2 and 10 years old estimated at 19% in 2017 [19]. This transmission is potentially driven by a number of reservoirs (artificial lakes) which provide permanent sources of water within which female Anopheles can lay their eggs [20]. Images were captured over three visits in June 2018 (early dry season), October 2019 (late dry season) and February 2020 (wet season) using two drones, both of which were able to capture images using a standard RGB camera, plus a near-infrared (NIR) camera. Both drones were purchased off-the-shelf from commercial vendors (Table 1). The first was a multirotor (quadcopter) type aimed towards the ‘hobbyist’ market (the DJI Phantom 4 Pro), supplemented by an additional NIR sensor by Sentera [21]. The second was a fixed-wing drone marketed towards the agriculture industry (eBee SQ) which incorporated a Parrot Sequoia multispectral camera.

Fig. 1
figure 1

Locations of sites surveyed within the ‘humanitarian drone testing corridor’, centred on Kasungu town, Central Malawi (inset). Coordinates can be found in Additional file 1

Table 1 Comparison of approximate costs required to capture and process imagery

Drone image processing

Individual images captured by both drones during each mapping mission were stitched together into orthomosaics using a structure-from-motion (SfM) workflow within the commercial image processing software Agisoft Metashape Professional (version 1.4.2). SfM is a photogrammetry approach which uses multiple 2D images overlapping images to construct a 3D landscape in the form of an elevation point cloud, and subsequently geometrically correct and combine the aerial images [22]. A subset of images captured during the wet season were then classified using a geographical object-based image analysis (GeoOBIA), using the LargeScaleMeanShift algorithm within the open source software Orfeo Toolbox (OTB, version 7.1.0), applied within the QGIS environment (version 3.8.1). GeoOBIA involves grouping contiguous pixels into ‘objects’ or ‘segments’ such that each segment is relatively homogenous (within a prespecified threshold) with respect to pixel characteristics. In this instance, pixels were grouped into segments according the values of red, green, blue and elevation, with the latter being estimated using photogrammetric methods within Agisoft Metashape and then rescaled to lie between 0 and 255 to match the scale of the RGB values. Trial and error was used to select the optimal segmentation parameters i.e. the spatial radius, range radius and minimum segment size. The smoothing radius determines how the amount by which the image is smoothed or filtered prior to the segmentation algorithm being implemented, whereas the range radius determines the similarity between pixels for grouping within the same segment. Similarity in this context refers to the Euclidean distance between two pixels.

Supervised classification was then undertaken to assign each segment to one of 12 land cover classes identified in the image: open water, floating aquatic vegetation, emergent aquatic vegetation, submerged aquatic vegetation, trees/bushes, grass, bare soil, iron roof, rusted iron roof, thatched roof, tarmac road, untarmacked roads/paths (Additional file 2). Additional file 3 displays the segment-level characteristics used to train the classification algorithm, which incorporated characteristics related to segment texture (Haralick textural features [23]) in addition to a range of water and vegetation indices. The mean and variance of each of these were used in the classification.

Classification was undertaken using a set of 1800 segments all of which were firstly manually classified by the research team. One-third of the segments (n = 600) were within a 400 m by 400 m area and were used to train the classification algorithm. An additional one-third were in the same 400 m by 400 m area and were used for evaluating the accuracy of the classification within the same geographical area used for training (spatial interpolation), whereas the remaining 600 segments were distributed outside of the area used for training (spatial extrapolation). Classification was undertaken in R (version 3.6.1) using the caret package and the Random Forests (RF) classification algorithm [24]. OTB can be used to undertake the classification process, however R was used to more efficiently interrogate the data to determine the influence of training and testing set selection on classification accuracy. Therefore, whilst a wide range of supervised classification algorithms are available in caret, RF was chosen as this method is available in OTB and other commonly used commercial GeoOBIA software such as eCognition [25, 26]. Manually classified segments within the 400 m by 400 m area were randomly split into training and testing segments. A ten-fold cross validation approach was used to determine the optimal tuning parameters for the RF algorithm, and the resulting model was then applied to the testing segments. Scores were then produced to determine the accuracy of the classification for the interpolated and extrapolated testing segments considering all 12 land cover classes, followed by a reduced classification that only differentiated between surface water (open water and aquatic vegetation) and any other class. These scores included the percentage of classifications that were correct (accuracy), Cohen’s kappa agreement statistic, quantity disagreement, i.e. the difference in the number of segments classified as each category between the manual and automated classification, and allocation disagreement, i.e. difference in locations of segments classified in each category, as overall measures. Quantity and allocation disagreement were calculated using the diffeR package in R [27]. Further, the producer and user accuracy, also known as recall and precision respectively, were calculated for each of the 12 individual classes, and for the reduced classification categories. The classification was then applied to the entire 600 m by 800 m area, and the percentage of the area classed as being covered in surface water was calculated. This process of randomly splitting the segments into training and testing groups was repeated 100 times, and the median and inter-quartile range of the resulting accuracy, kappa agreement and percentage surface water cover were reported.

A manual classification of surface water was also undertaken which involved systematically scanning through the image from left to right and drawing a polygon around each area determined by the assessor to be surface water, regardless of whether it was open water or aquatic vegetation. This process was undertaken independently of the manual generation of training and testing segments, hence polygons were not constrained by the boundaries of the segments. This manual classification was undertaken by two researchers independently to evaluate the consistency of manual classifications. The union of these two areas was then considered as the manually classified surface water layer. As both the derivation of the training and testing segments plus the fully manual classification were undertaken by the same researchers and involve subjectively determining whether a segment (GeoOBIA approach) or area of pixels (manual approach) was classed as surface water, the accuracy of the two approaches cannot be objectively compared against one another. Instead, the level of agreement between the two resulting surface water polygons was evaluated by calculating the area covered by the union of the outputs obtained using the two approaches and calculating the percentage of the union covered by the intersect of the two outputs. A percentage of 100 implies the two outputs are identical whereas a value of zero implies there is complete disagreement between the two. The time taken and the computer resources required for each of these tasks were also recorded. An overview of the image capture, processing and classification procedures is presented in Fig. 2.

Fig. 2
figure 2

Processes undertaken to identify larval habitat from drone imagery

Entomological sampling

Larval surveys were conducted concurrently to the drone image capture during each field visit (Fig. 1). Permission to collect these data were obtained from the Kasungu District Council, Kasungu Water Board and private landowners. Sites were selected as part of a larger study exploring the role of permanent water bodies such as local reservoirs on dry season transmission, with drone imagery being captured to monitor the changes in the landscape as the season transitioned from dry to wet. Due to delays in fieldwork as a result of Covid-19, this full study is still ongoing, and results will be published elsewhere. Preliminary work to determine the presence of mosquito larvae in the area in both wet and dry season involved three field visits. The first two field visits were undertaken during the dry season, with larval sampling being undertaken at regular intervals around the periphery of reservoirs only, with surface water being scarce elsewhere. Reservoirs were selected purposively based on their proximity to Kasungu town (for accessibility) and their proximity to human settlements. The third visit was conducted during the wet season, and sampling was focused around one of the reservoirs sampled in the previous dry season. During this period both temporary and permanent water bodies were present and larval surveys were undertaken in a selection of these sites that had been identified using drone imagery captured the previous day. A subset of sites was sampled on four consecutive days to determine their consistency with respect to larval presence.

At each site, the presence and number of larvae were recorded using 10 repeated dips of the surface water, categorised by stage (L1/L2 or L3/L4) and either anopheline or culicine. To characterise malaria vectors in Kasungu as part of the study’s broader efforts to understand transmission in the area, all anopheline larvae were raised to adult stage and identified morphologically to species [28]. The location, description and photographs of each site were recorded using an Android Smartphone and Open Data Kit (ODK). These photographs were later used by the research team to manually evaluate each site according to the amount (none, 0–1/3, 1/3−2/3, > 2/3 of the area) and type of vegetation present plus turbidity (turbid or clear). A generalized linear mixed model was then fitted to the resulting presence/absence data using R (v 3.6.3) to predict the likelihood that larvae were present from biotic site information obtained via drone imagery (vegetation type, coverage, turbidity) which were included as fixed effects, with sample area being included as a random effect. Season (wet/dry) was included as an additional fixed effect in the model. Models were fitted to presence/absence data for any larvae irrespective of stage or genus, and for late stage larvae only as characteristics of habitat containing late stage larvae are considered by WHO to be of greater importance than early stage [11]. The productivity of the sampled sites with respect to the number of early or late stage larvae collected was also considered.

Results

Image capturing

During the three sampling periods, images of 10 distinct areas in Kasungu were captured, covering an area of 8.9 km2. The two drones significantly differ in relation to operational costs, equipment and software requirements and usage. Tables 1 and 2 describe the primary differences in relation to initial costs and operational usage respectively. The fixed wing drone (eBeeSQ) had a greater initial cost than the multirotor Phantom 4 Pro due to it being inclusive of a NIR sensor, costing approximately £7000 (inclusive of an educational discount), compared to £3300 for a standard (RGB sensor) Phantom 4 Pro drone on which a NIR sensor was retrofitted. The eBeeSQ also required a high-spec laptop (£1000+) on which to run the software required to plan and conduct missions, whereas the Phantom 4 Pro was operated using free apps installed on GPS-enabled Android or iOS smartphone or tablet devices. On an operational level, the primary differences are between the flight times per battery, and the ease of use (Table 2). While, overall the Phantom 4 Pro is easier to use due to the small amount of open space required for take-off and landing, its limited battery life means that to cover a relatively modest area of 1 km2, 2–3 individual flights are needed depending on whether a fixed launch site is used or whether this is adapted to minimize flight time. This comes at a cost of both time and money, particularly given each battery comes at a price of £150. The eBeeSQ fixed wing drone requires less energy to fly and therefore batteries last approximately twice as long (up to 1 h in comparison to 30 min theoretical flight time, 45 min in comparison to 22 min practical flight time) than the Phantom 4 Pro. Therefore, while the time required to cover the same area is longer, this area can be comfortably covered using fewer batteries and a single launch site, meaning that in practice the process is more efficient. For example, the eBeeSQ takes 66 min (two batteries) to cover a square with an area of 1 km2 compared to 38 min (three batteries) with the Phantom 4 Pro when flying at 120 m above sea level (asl) with an 80% overlap in captured images. The fixed-wing drone is however more difficult to operate than the Phantom 4 Pro, requires a larger space for take-off and landing, and therefore cannot be used in more densely vegetated areas.

Table 2 Practical and operational differences between the drones used in this study

Image processing and classification

Agisoft Metashape was used to process all images captured. Images captured by the Phantom 4 Pro over Malangano in the wet season were used to demonstrate the processing and classification process (Fig. 1). Flying at 120 m above surface level, with an 70% overlap in images, a total of 782 individual images (6.2 GB) covered an area of 1.77 km2. Processing these images in Agisoft Metashape in order to produce an orthomosaic of the area and an accompanying digital surface model took a total of 250 min using a computer with an Intel Core i7-6700 processor, 32 GB RAM, resulting in an orthomosaic with a spatial resolution of 3 cm (file size = 4.2 GB). A subset of the image covering an area of 0.48 km2 (800 m by 600 m, file size = 1.9 GB) was then selected for classification (Fig. 3).

Fig. 3
figure 3

Image captured by the Phantom 4 Pro in Kasungu in February 2020 covering 800 m by 600 m, with each grid representing 200 m by 200 m. Grids 4, 5, 7 and 8 were used for training the classification algorithm, and an assessment of its accuracy was made using features both within this area (interpolation) and in the surrounding grids (extrapolation)

A set of 1800 training and testing segments were then generated for the 12 identified land classes (Additional file 2), plus an additional category representing areas that were in shadow. This process took approximately six hours. The study area was partitioned into cells of 200 m by 200 m, labelled as cells 1–8 (Fig. 3), and the training and testing segments were proportionally distributed throughout the cells as follows: two thirds (1200) of segments were within cells 4, 5, 7 and 8 covering an area of 400 m by 400 m, referred to in this paper as the internal segments. One third of segments (600) were within the remaining cells (1–3, 6, 9, 10–12), referred to in this paper as external segments. Segmentation was performed using Orfeo Toolbox (OTB) functions within the QGIS environment. Figure 4 demonstrates the impact of varying values of the spatial and range radius on the resulting segmentation and the time taken to perform this segmentation over a 100 m by 100 m area using the computer specifications previously specified (see "Methods"). While increasing the spatial radius provided a more adequate balance between over-segmentation (single discrete features of interest being split into many segments) and under-segmentation (multiple discrete features of interest being grouped into a single segment), this came at the price of substantially increasing the processing time. Additional processing time is required to calculate the segment-level summaries (mean, variance) of each of the variables being used to classify the imagery, with processing time increasing as the number of segments increases.

Fig. 4
figure 4

Examples of the segmentation process under different values for spatial radius s (0,10, 30), range radius r (25, 50) with a minimum segment size of 100. Time t corresponds to the time taken in seconds to segment a 100 m by 100 m image with a spatial resolution of 3 cm using the LargeScaleMeanShift algorithm in Orfeo Toolbox. This process includes calculating the mean and variance of the RGB and elevation values for each segment

The segmentation process was then applied to the entire 800 m by 600 m area using the parameters 10 (spatial radius), 25 (range radius) and 200 (minimum segment size), creating close to 800,000 segments. The total processing time, which includes calculating the segment-level mean and variance of the RGB and elevation values, was 24.5 h with an additional 15 h taken to calculate the mean and variance of each of the additional variables under consideration (Additional file 2). Two classifications were then undertaken, one of which included the NIR-derived variables and one of which did not.

The resulting accuracy assessments of these classifications obtained without using NIR-derived variables are presented in Table 3, with class-level producer and user accuracy across all 100 samples presented in Fig. 5. For interpolated areas, user and producer accuracy was generally high (> 80%) for all classes. In extrapolated areas there was a reduction in accuracy across all classes, most notably in the class representing bare soil. A representation of the classified output from one area of Kasungu excluding NIR-derived variables is shown in Fig. 6, with Additional file 4 presenting the corresponding contingency table for interpolated areas only. A representation of the classified output obtained using NIR-derived variables is presented in Additional file 5. The corresponding variable importance plots are available in Additional file 6. The results are very similar for both models fitted with and without the NIR-derived variables. The variable making the greatest contribution to the classification model in both cases is the mean elevation, with mean red, blue, green and brightness also important. While the NIR-derived variables NDVI and SAVI make the greatest contribution to the classification algorithm in the second model (Additional file 6), Table 3 indicates that the inclusion of NIR-derived variables does not make any significant impact on classification accuracy, with all interpolated accuracy measures being similar when computed using a classification which includes and excludes the NIR-derived variables. For example, the overall median interpolated accuracy obtained using NIR-derived variables is marginally lower (0.907) than that obtained without using NIR-derived variables (0.912), whereas quantity disagreement increases to 22 from 21 when accounting for NIR-derived variables. There is a clear drop in all accuracy measures when considering data from the extrapolation area. Median accuracy reduces to 0.761 and 0.798 when considering overall accuracy without and with NIR-derived variables respectively whereas quantity disagreement increases to 70 and 63 respectively. This reduction in classification quality is less pronounced when considering surface water (open water or aquatic vegetation) accuracy alone compared with trying to distinguish between all 12 land cover classes.

Table 3 Summaries of classification accuracy
Fig. 5
figure 5

Boxplots of class-level user and producer accuracy for the classifications undertaken without using NIR-derived variables

Fig. 6
figure 6

Example of a classification obtained for the entire study area using the random forests algorithm without including NIR-derived variables (left), with a more detailed view of a smaller area comparing the original image (top right) with the classified image (bottom right)

Summaries are for all 12 classes (overall) and for surface water (including open water and aquatic vegetation) versus all other classes for GeoOBIA obtained with and without NIR-derived variables. Measures including accuracy (proportion of segments correctly classified), Cohen’s kappa, (0 = disagreement, 1 = complete agreement), quantity disagreement (different in the number of segments in each class as determined by the manual assessment and the GeoOBIA) and allocation disagreement (difference in the locations of segments in each class between the manual and GeoOBIA).

There is a close agreement in the percentage of the 400 m by 600 m area that is covered in surface water obtained by fitting the model with and without NIR-derived variables (without NIR: median = 18.3%; with NIR: median = 16.6%), with more variability observed in the NIR-inclusive models (Table 3). A fully manual classification of the image, independently undertaken by two researchers, resulted in a larger percentage of the area being identified as surface water (21.2 and 20.1% by researchers 1 and 2 respectively), with the union of these two outputs (‘manual surface water classification’) covering 22.0% of the area. This process of manual classification took approximately 2 h to complete. Using the classifications presented in Fig. 6 and Additional file 5 as an example, the intersect of the manual and GeoOBIA-derived classifications covered 67.2% (without NIR) and 61.2% (with NIR) of the union of the classifications, covering 17.1 and 15.3% of the study area respectively. Without using NIR-derived variables, 19.1% of the union was covered by the manual classification alone, with the remaining 13.7% being covered by the GeoOBIA classification along. Similarly, using NIR-derived variables, 26.8% of the union was covered by the manual classification alone, with the remaining 12.0% being covered by the GeoOBIA classification. Without extensive contemporary ground-based surveys it is not possible to determine which classification approach is more accurate, however it is clear that there are substantial differences between aerial image-based manual and automated classifications.

Entomological surveys

During three separate field visits (June 2018, October 2019, and February 2020) a total of 326 larval sites were sampled (available through the Figshare repository [29]). During the dry periods these samples were focused along the shorelines of larger permanent water bodies (296 sites), with a mixture of 30 temporary and permanent sites surveyed during the wet season (Table 4). Both anopheline and culicine larvae were found throughout the area during each sampling period (56% of sites sampled), with the lowest proportion of positive sites found in the late dry season (76% in June 2018, 16% in October 2019, 70% in February 2020). No clear sympatry was observed between anopheline and culicine larvae in this study. For example, of the 321 sites where late stage larvae data were recorded (excluding five sites with missing data), larvae were observed in 31% (101) of samples with only 7% (24) sites containing both anophelines and culicines. In the area surrounding the Malangano site (Fig. 3) in February 2020, 177 out of 297 anopheline specimens were morphologically identified to species level, finding a predominance of An. gambiae sensu lato (s.l.) (87.6%, 155/177) followed by Anopheles coustani (8.5%, 15/177) and very few Anopheles pretoriensis (2.3%, 4/177) and An. funestus (1.7%, 3/177).

Table 4 Summaries of the larval sampling sites by presence/absence of mosquito larvae found

At each site, GPS coordinates were recorded using the ODK app, photographs were taken using a smartphone, and aerial imagery was captured (Fig. 7). Samples were taken within approximately one metre of where the researcher stood to record the coordinate, however, as GPS coordinates have an accuracy of approximately three metres it was not possible to pinpoint precisely where the samples were taken within the aerial imagery.

Fig. 7
figure 7

Examples of sampling sites where anopheline larvae were found. The top row indicates the precise GPS location captured using ODK (yellow circle), the expected sampling area based on these coordinates (1 m radius), and the expected accuracy of the coordinates (3 m radius), overlaid on top of the drone imagery. The middle row presents the classified imagery for these sites and the bottom row contains photographs of each site taken at the time of sampling

Using the photographs, each site was characterized according to presence/absence of larvae and sample site characteristics including dominant vegetation type (none, floating, submerged, emerging), vegetation cover (none, < 1/3, 1/3−2/3, > 2/3) and turbidity (Additional file 8). Vegetation was present in most sites sampled (284/326, 87%). Of these, 64% (183) contained emergent vegetation, 20% (57) contained submerged and 15% (44) contained floating vegetation (Table 4). Vegetation cover varied evenly across sites, with 32% having low (0–1/3) coverage, 35% having moderate (1/3−2/3 coverage) and 33% having high (> 2/3) coverage. There was an interaction between vegetation type and coverage, such that sites with floating vegetation rarely had high vegetation coverage (Additional file 9). With regards to turbidity, 40% (129) of sites were classed as turbid whereas the remaining 60% (197) were clear.

Due to the strong interaction was observed between vegetation type and coverage, when fitting the GLMMs to the presence/absence data these variables were not included in the model simultaneously, but rather explored which of the two resulted in the best fitting model with regards to Akaike Information Criterion (AIC). After counting for the effect of sampling period and site, there was a strong association between the presence of vegetation and the likelihood of any larvae (log-OR = 1.44, 95% CI [0.34, 2.67], p = 0.01), however accounting for vegetation coverage or vegetation type did not improve the model further. When considering Anopheles L3 and L4 larvae only (Additional file 7), the model was improved when vegetation type was considered such that larvae were more likely to be present when emerging (log-OR = 1.14, 95% CI [− 0.05, 2.51], p = 0.07) or submerged (log-OR = 1.90, 95% CI [0.51, 3.44], p = 0.01) vegetation were available, compared to sites with no vegetation. With regards to productivity, while there was variability in the abundance of larvae sampled per site (145 sites, min = 1, median = 4, max = 56), there were insufficient high productivity sites to formally explore any trends in their characteristics.

During the wet season, 10 sites were repeatedly sampled over four consecutive days, with larvae consistently observed in four sites and no larvae being found on at least one day in the remaining six sites. Due to changes in the environment it was difficult to resample the same locations across larger time scales. Temporary surface water observed in the wet season dried up even after just a few days without rain and shorelines of permanent water bodies varied substantially both between seasons and between the same season over consecutive years (Fig. 8). For example, images captured later in the dry season (October) in 2019 were wetter than those captured in the early dry season (June) in 2018 at the Chitete reservoir.

Fig. 8
figure 8

Comparisons of aerial images captured at different seasonal time points. Left images display a comparison between consecutive dry (Oct 2019) and wet (Feb 2020) seasons around the Malangano dam. Right images display comparisons between dry seasons over 2 consecutive years (June 2018, Oct 2019) around the Chitete dam

Discussion

This study has established that established that drone-captured imagery can be used to accurately identify environmental characteristics associated with larval habitat in an area of land surrounding a manmade reservoir in Kasungu district, Malawi. Whilst this is a useful demonstration which provides valuable information on the theoretical aspects of mapping larval habitats using drone imagery, the primary objective for conducting this study was to explore the feasibility of this technology being used as part of a malaria control programmes’ toolkit. This discussion section therefore focuses on the potential benefits and practical bottlenecks of the methods applied in this paper with respect to drone image capture and drone image processing, plus the continued role of ground-based larval surveys.

Image capturing

Image capture using drones inevitably leads to technical and skills-based challenges and a few of these have been highlighted here in the context of searching for water bodies in a rural setting. Aside from hardware and software issues, it was noted that flight experience was a key requirement to determine optimal flight times as neither the rotor or fixed wing drone could be flown in wet or windy conditions, and the study team experienced the impact of extreme weather on the hardware with multiple occasions of over-heating on warm days. This emphasized the need for extensive drone piloting training by the operator. The country’s drone regulations also need to be taken into careful consideration. In Malawi, data capture was facilitated by the relationship between UNICEF and the Department of Civil Aviation and a toolkit is currently being developed to outline the procedures that need to be followed by those wishing to fly drones for non-commercial purposes (https://www.updwg.org/wp-content/uploads/2019/12/Malawi-RPA-Toolkit-2019_Dec.-Final.pdf). While regulations vary by country, national civil aviation authorities are also requiring drone pilots to obtain accredited qualifications and seek appropriate permissions before using drones for research or humanitarian purposes. Training courses which cover both the operational and regulatory aspects of drone flying are currently quite sparse in sub-Saharan Africa and this may require future pilots to travel outside of their own country to gain the necessary experience. Should a malaria control implementer wish to use drone imagery within their programmes, they may therefore incur significant expense both in purchasing the equipment and training their staff. A solution to this would be to outsource the image capturing to qualified drone pilots operating in the area.

An additional bottleneck is the availability of hardware within the country of operation. While it may be possible to purchase off-the-shelf drones in-country, should any technical issues arise, obtaining part replacements or repairs becomes problematic and expensive. Investments are therefore being made in ‘home-grown’ drones, to support local economies, decrease the cost of equipment, and make repairs much more easily accessible. In Malawi for example, MicroMek (https://www.micromek.net/) manufacture the low-cost fixed wing drone known as EcoSoar [30], for both transporting goods and capturing imagery.

Image processing

Processing drone imagery to create the orthomosaics is time-consuming, requires a high-spec computer and a large capacity for data storage. Therefore, to use this imagery in the field, an NMCP would require people skilled in both image capture and processing, plus access to the relevant software. These skills are not usually taught as part of standard drone pilot training, however this may change as the potential for using drone technology for humanitarian purposes is increasingly realized. For example, the African Drone and Data Academy was launched in January 2020 in Malawi to build capacity in both drone piloting and drone image processing and analysis [31].

Image classification is appealing because once the algorithm has been trained, it can simply be applied to any additional imagery captured without any or only a little additional data being required. In this analysis there was a decline in classification accuracy in areas within very close proximity to that used to train the algorithm and noted that even in this small area there were important land cover classes in the extrapolation area that did not appear in the training area e.g. red algae in the water. This challenge is likely to be exacerbated when considering areas further apart, or data collected at different time points. As more data are collected these limitations may be overcome, however in the short term, the effort required by the end-user e.g. an individual NMCP to train a classification algorithm may outweigh its benefits. In this demonstration, a geographical object-based classification approach was implemented which generated the segments and computing segment-level characteristics prior to training and applying the classification algorithm. The segment-generating process can be very time-consuming depending on the values of the segmentation parameters, the size of the area being classified and the segmentation algorithm being used [32]. While other classification techniques such as pixel-based classification may be quicker to perform, an object-based approach is the most appropriate for very high-resolution images such as that generated by drones [33]. The cost and benefits of accuracy against processing time therefore need to be considered should an NMCP wish to perform image classification in-house. The role of additional sensors in the image classification process is also unclear. In this analysis, comparisons were made between classification accuracy using imagery captured from a standard camera only (RGB) and additional imagery captured by a much more expensive NIR sensor. While the performance metrics indicated very little difference in the accuracy obtained using the two approaches, a more extensive investigation would need to be undertaken over a more environmentally diverse area before any conclusions can be drawn on the benefits to incorporating this additional technology. Further sensors such as shortwave infrared (SWIR) and thermal infrared may further improve the surface water classification accuracy [34].

A more practical solution to ‘automated’ image classification may be to persevere with the less efficient, but lower skilled task of manual classification. This task is, however, not without its drawbacks, as human error can easily miss small areas, or misclassify water containing a lot of aquatic vegetation as land and vice versa. These latter ‘missed’ areas are of significance, as Anopheles mosquitoes are generally found in water containing vegetation. The fact that there was a 10% discrepancy in the manual classification undertaken by two independent researchers, both of whom were familiar with the study area, demonstrates the fallibility in this method.

As with the drone image capture, an alternative is to outsource these activities to an organization which specializes in image processing and classification. Additionally, cloud-based computing services such as DroneDeploy’s Map Engine [35] and Google Earth Engine [36] could be used as these allow individuals/groups to harness the power of remote servers to manage and manipulate the data. This approach could facilitate the development of more automated habitat classification approaches i.e. using data from other organizations, previous field or professional expertise in remote sensing to develop classification algorithms that do not require the use of bespoke training data. The TropWet tool developed by Hardy, Oakes, and Ettritch [37] is a demonstration of this in which satellite imagery (Landsat, 30 m resolution) is automatically classified for a user-specified area and time period using a Google Earth Engine interface.

There are still practical challenges with these approaches, particularly relating to the upload of large image files to enable these processes to be undertaken remotely, however these may be preferable to the more technical challenge of managing the data in-house.

Entomological survey

Larval surveys are an important part of the process of LSM both to confirm the species of mosquitoes found in the area, to characterize the types of surface water where larvae are likely to be found, and to monitor the progress of any subsequent intervention. Larval surveys are however a time-consuming process, particularly when undertaken during the wet season during which areas become inaccessible following heavy rains. The role of drones in LSM is not to completely remove the need for larval surveys, but to help differentiate between water bodies with respect to their potential as larval habitat and/or to differentiate sites according to their potential larval productivity.

It was noted that vegetation coverage was important when considering presence/absence of late stage Anopheles larvae, with coverage correlated with the type of vegetation found i.e. coverage of floating vegetation was likely to be less than that of emergent or submerged vegetation. A full understanding of the larval ecology of the local individual malaria vectors would greatly assist a targeted LSM approach aided by drone-imagery support. In south-eastern Tanzania, a basic characterization of An. funestus larval habitats provides support that this species occupies small spring-fed pools, permanent natural ponds and slow-moving waters each of which fall under the ‘few, fixed and findable’ paradigm [5]. In a recent study in Southern Malawi [6], An. arabiensis was the dominant species, with high densities being found in aquatic habitats surrounded by bare soil. A species-specific approach to identifying larval habitat using drone imagery may therefore be required, with imagery captured throughout the year to better understand the temporal dynamics of larval habitat and thereby optimize the impact of any potential intervention. These images could further be used to monitor the progress of LSM campaigns with, for example, a more accurate estimates of LSM coverage and demonstrable changes in the landscape because of habitat removal/modification. Further entomological surveillance remains pivotal to establish where and when LSM should be deployed and measure the impact of the intervention on malaria transmission potential.

Conclusions

This study demonstrates the potential for drone imagery to be used as a tool to support the identification of mosquito larval habitat in rural areas where malaria is endemic. While this technology has the capacity to complement the more labour-intensive approach of identifying larval habitat from the ground, there are technical challenges to overcome before it can be smoothly integrated into malaria control activities. The authors believe that outsourcing the capturing and processing of drone imagery to private companies with the equipment and skills necessary to extract the required information is a more practical approach to developing equivalent skills in house. These services are becoming increasingly available in other sectors such as agriculture, forestry and environmental monitoring and there are promising developments in the African drone sector to support this local capacity. It is however important to emphasize that drone imagery should not be used to completely replace larval surveys. Instead, this technology could provide Additional information which may help to reduce the time spent finding locations to be sampled, monitor environmental changes over time and help to guide the frequency and scale of any LSM intervention, ultimately increasing its potential for success. Further consultations between experts and stakeholders in the fields of drones, image analysis and vector control are needed to develop more detailed guidance on how this technology can be most effectively exploited.

Availability of data and materials

The entomological dataset generated and analysed during the current study is available in the Figshare repository: https://figshare.com/s/25102b084f56f41a1ca8.

All drone imagery referenced in this paper is available from OpenAerialMap: (https://map.openaerialmap.org/#/33.489675521850586,-13.050633215031182,13/user/5f1fe6f357ddda00054a0647?_k=mmgy4y).

References

  1. Stresman G, Bousema T, Cook J. Malaria hotspots: is there epidemiological evidence for fine-scale spatial targeting of interventions? Trends Parasitol. 2019;35:822–34.

    Article  Google Scholar 

  2. Bousema T, Stresman G, Baidjoe AY, et al. The impact of hotspot-targeted interventions on malaria transmission in Rachuonyo South District in the Western Kenyan Highlands: a cluster-randomized controlled trial. PLoS Med. 2016;13:e1001993.

    Article  Google Scholar 

  3. Hsiang MS, Ntuku H, Roberts KW, Dufour MSK, Whittemore B, Tambo M, et al. Effectiveness of reactive focal mass drug administration and reactive focal vector control to reduce malaria transmission in the low malaria-endemic setting of Namibia: a cluster-randomised controlled, open-label, two-by-two factorial design trial. Lancet. 2020;395:1361–73.

    Article  CAS  Google Scholar 

  4. Sy O, Niang EHA, Diallo A, Ndiaye A, Konaté L, Ba EHCC, et al. Evaluation of the effectiveness of a targeted community-based IRS approach for malaria elimination in an area of low malaria transmission of the central-western Senegal. Parasite Epidemiol Control. 2019;6:e00109.

    Article  Google Scholar 

  5. Nambunga IH, Ngowo HS, Mapua SA, Hape EE, Msugupakulya BJ, Msaky DS, et al. Aquatic habitats of the malaria vector Anopheles funestus in rural south-eastern Tanzania. Malar J. 2020;19:219.

    Article  CAS  Google Scholar 

  6. Gowelo SA, Chirombo J, Koenraadt CJM, Mzilahowa T, Berg H, Takken W, et al. Characterisation of anopheline larval habitats in southern Malawi. Acta Trop. 2020;210:105558.

    Article  Google Scholar 

  7. Eneh LK, Fillinger U, Borg Karlson AK, Kuttuva Rajarao G, Lindh J. Anopheles arabiensis oviposition site selection in response to habitat persistence and associated physicochemical parameters, bacteria and volatile profiles. Med Vet Entomol. 2019;33:56–67.

    Article  CAS  Google Scholar 

  8. Musiime AK, Smith DL, Kilama M, Geoffrey O, Kyagamba P, Rek J, et al. Identification and characterization of immature Anopheles and culicines (Diptera: Culicidae) at three sites of varying malaria transmission intensities in Uganda. Malar J. 2020;19:221.

    Article  CAS  Google Scholar 

  9. Beck-Johnson LM, Nelson WA, Paaijmans KP, Read AF, Thomas MB, Bjørnstad ON. The effect of temperature on Anopheles mosquito population dynamics and the potential for malaria transmission. PLoS ONE. 2013;8:79276.

    Article  Google Scholar 

  10. Vantaux A, Ouattarra I, Lefèvre T, Dabiré KR. Effects of larvicidal and larval nutritional stresses on Anopheles gambiae development, survival and competence for Plasmodium falciparum. Parasit Vectors. 2016;9:226.

    Article  Google Scholar 

  11. WHO. Larval source management—a supplementary measure for malaria vector control. An operational manual. Geneva, World Health Organization. 2013. http://www.who.int/malaria/publications/atoz/9789241505604/en/.

  12. Tusting LS, Thwing J, Sinclair D, Fillinger U, Gimnig J, Bonner KE, et al. Mosquito larval source management for controlling malaria. Cochrane Database Syst Rev. 2013;CD008923.

  13. Fillinger U, Lindsay SW. Larval source management for malaria control in Africa: myths and reality. Malar J. 2011;10:353.

    Article  Google Scholar 

  14. Killeen GF. Characterizing, controlling and eliminating residual malaria transmission. Malar J. 2014;13:330.

    Article  Google Scholar 

  15. Hardy A, Makame M, Cross D, Majambere S, Msellem M. Using low-cost drones to map malaria vector habitats. Parasit Vectors. 2017;10:29.

    Article  Google Scholar 

  16. Carrasco-Escobar G, Manrique E, Ruiz-Cabrejos J, Saavedra M, Alava F, Bickersmith S, et al. High-accuracy detection of malaria vector larval habitats using drone-based multispectral imagery. PLoS Negl Trop Dis. 2019;13:e0007105.

    Article  Google Scholar 

  17. Tokarz R, Novak RJ. Spatial–temporal distribution of Anopheles larval habitats in Uganda using GIS/remote sensing technologies. Malar J. 2018;17:420.

    Article  Google Scholar 

  18. Hardy A, Ettritch G, Cross D, Bunting P, Liywalii F, Sakala J, et al. Automatic detection of open and vegetated water bodies using sentinel 1 to map African malaria vector mosquito breeding habitats. Remote Sens. 2019;11:593.

    Article  Google Scholar 

  19. Chipeta MG, Giorgi E, Mategula D, Macharia PM, Ligomba C, Munyenyembe A, et al. Geostatistical analysis of Malawi’s changing malaria transmission from 2010 to 2017. Wellcome Open Res. 2019;4:57.

    Article  Google Scholar 

  20. Kibret S, Lautze J, McCartney M, Nhamo L, Wilson GG. Malaria and large dams in sub-Saharan Africa: future impacts in a changing climate. Malar J. 2016;15:448.

    Article  Google Scholar 

  21. Sentera. 2020. https://sentera.com/sensors/.

  22. Anderson K, Westoby MJ, James MR. Low-budget topographic surveying comes of age: structure from motion photogrammetry in geography and the geosciences. Prog Phys Geogr Earth Environ. 2019;43:163–73.

    Article  Google Scholar 

  23. Haralick RM, Dinstein I, Shanmugam K. Textural features for image classification. IEEE Trans Syst Man Cybern. 1973;SMC-3:610–21.

    Article  Google Scholar 

  24. Kuhn M. Building predictive models in R using the caret package. J Stat Softw. 2008;28:1–26.

    Article  Google Scholar 

  25. De Luca G, Silva N, Cerasoli JM, Araújo S, Campos J, Di Fazio J. S, et al. Object-based land cover classification of cork oak woodlands using UAV imagery and Orfeo ToolBox. Remote Sens. 2019;11:1238.

    Article  Google Scholar 

  26. Ma L, Fu T, Blaschke T, Li M, Tiede D, Zhou Z, et al. Evaluation of feature selection methods for object-based land cover mapping of unmanned aerial vehicle imagery using random forest and support vector machine classifiers. ISPRS Int J Geo-Inf. 2017;6:51.

    Article  CAS  Google Scholar 

  27. Pontius RG, Millones M. Death to Kappa: Birth of quantity disagreement and allocation disagreement for accuracy assessment. Int J Remote Sens. 2011;32:4407–29.

    Article  Google Scholar 

  28. Coetzee M. Key to the females of Afrotropical Anopheles mosquitoes (Diptera: Culicidae). Malar J. 2020;19:70.

    Article  Google Scholar 

  29. Stanton MC. Figshare dataset: larval sampling data, Kasungu. Figshare. 2020. https://figshare.com/s/25102b084f56f41a1ca8.

  30. Standridge Z. Design and development of low-cost multi-function UAV suitable for production and operation in low resource environments. Thesis, MSc Aerospace Engineering, Virginia Polytechnic Institute and State University; 2018.

  31. UNICEF. African Drone and Data Academy. 2020.  https://www.unicef.org/malawi/african-drone-and-data-academy-malawi.

  32. Hossain MD, Chen D. Segmentation for Object-Based Image Analysis (OBIA): a review of algorithms and challenges from remote sensing perspective. J Photogramm Remote Sens. 2019;150:115–34.

    Article  Google Scholar 

  33. Pande-Chhetri R, Abd-Elrahman A, Liu T, Morton J, Wilhelm VL. Object-based classification of wetland vegetation using very high-resolution unmanned air system imagery. Eur J Remote Sens. 2017;50:564–76.

    Article  Google Scholar 

  34. Khanh Ni TN, Tin HC, Thach VT, Jamet C, Saizen I. Mapping submerged aquatic vegetation along the central Vietnamese coast using multi-source remote sensing. Int J Geo-Inf. 2020;9:395.

    Article  Google Scholar 

  35. DroneDeploy. Introducing map engine. 2018. https://medium.com/aerial-acuity/introducing-map-engine-cd3ef93bc730.

  36. Mutanga O, Kumar L. Google earth engine applications. Remote Sens. 2019;11:591.

    Article  Google Scholar 

  37. Hardy A, Oakes G, Ettritch G. Tropical wetland (TropWet) mapping tool: the automatic detection of open and vegetated waterbodies in Google Earth engine for tropical wetlands. Remote Sens. 2020;12:1182.

    Article  Google Scholar 

Download references

Acknowledgements

Not applicable.

Funding

This work was funded by awards from the Medical Research Council (MR/M014975/1) and the Wellcome Trust (215184/A/19/Z). These funders had no role in the design of this study, data collection, analyses, interpretation of the data, or the writing of the manuscript.

Author information

Authors and Affiliations

Authors

Contributions

MCS, CMJ and RHS conceived the ideas and designed methodology; all authors collected the data; MCS, PK and KZ analysed the data; MCS and CMJ led the writing of the manuscript. All authors contributed critically to the drafts and gave final approval for publication. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Michelle C. Stanton.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1.

Locations of reservoirs around which the drone image capture and entomological surveys were focused. The area covered represents the size of the area captured in the drone imagery.

Additional file 2.

Number of training and testing segments that were manually created for each land class.

Additional file 3.

A description of the variables that were derived from the images captured by the drones.

Additional file 4.

An example of a contingency table and individual class accuracy summaries obtained using a classification obtained without NIR-derived variables, applied to the interpolated area only.

Additional file 5.

Example of a classification obtained using the random forests algorithm including NIR-derived variables

Additional file 6.

Variable importance plot corresponding to the classification presented in Fig. 5 (without NIR) and Additional file 5: Fig. S1

Additional file 7.

Summaries of the larval sampling sites by presence/absence of late stage mosquito larvae found.

Additional file 8.

Examples of the vegetation types observed in the sampling sites – (from L to R) none, floating, submerged, emerging.

Additional file 9.

Summary of larval sampling sites by presence/absence of larvae, vegetation type and vegetation cover.

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 http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) 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

Stanton, M.C., Kalonde, P., Zembere, K. et al. The application of drones for mosquito larval habitat identification in rural environments: a practical approach for malaria control?. Malar J 20, 244 (2021). https://0-doi-org.brum.beds.ac.uk/10.1186/s12936-021-03759-2

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://0-doi-org.brum.beds.ac.uk/10.1186/s12936-021-03759-2

Keywords