Research | Open | Published:
Prediction of high-risk areas for visceral leishmaniasis using socioeconomic indicators and remote sensing data
International Journal of Health Geographicsvolume 13, Article number: 13 (2014)
Spatial heterogeneity in the incidence of visceral leishmaniasis (VL) is an important aspect to be considered in planning control actions for the disease. The objective of this study was to predict areas at high risk for visceral leishmaniasis (VL) based on socioeconomic indicators and remote sensing data. We applied classification and regression trees to develop and validate prediction models. Performance of the models was assessed by means of sensitivity, specificity and area under the ROC curve. The model developed was able to discriminate 15 subsets of census tracts (CT) with different probabilities of containing CT with high risk of VL occurrence. The model presented, respectively, in the validation and learning samples, sensitivity of 79% and 52%, specificity of 75% and 66%, and area under the ROC curve of 83% and 66%. Considering the complex network of factors involved in the occurrence of VL in urban areas, the results of this study showed that the development of a predictive model for VL might be feasible and useful for guiding interventions against the disease, but it is still a challenge as demonstrated by the unsatisfactory predictive performance of the model developed.
Since the 1980s, visceral leishmaniasis (VL) had its epidemiological profile modified in Brazil; no longer being characterized as a predominantly rural disease, but established in the urban environment as well .
The introduction, propagation and dissemination of VL in urban settings is associated to multiple and complex conditions, such as environmental changes due to migration movements, disorderly occupation of city’s outskirts, high population density, and inadequate living conditions [2, 3].
The spatial distribution of urban VL is markedly heterogeneous, which may lead to a substantial increase in transmission levels . In this situation, focusing intervention on high risk areas might be an efficient strategy to reduce transmission rates .
The objective of the present study is to characterize and predict high risk areas for the occurrence of VL in Teresina, Piauí, based on socioeconomic indicators and environmental data, obtained through remote sensing.
Area of study
Teresina, capital of Piauí state, in the northeast region of Brazil, is located between the confluence of the rivers Parnaíba e Poti, at a 5°5′ south latitude and 42°48′ west longitude. Teresina has a population of 793,915 inhabitants and population density of 444.2 inhabitants/km . It has a tropical sub-humid climate, with high registered temperatures throughout the year, varying from 22°C to 40°C. Rainy season is mainly from January until April.
This is an ecologic study, in which the units of analysis are the 430 and 653 urban census tracts (CT) of Teresina for the years of 1991 and 2000, respectively.
Data and variables
The average annual incidence rates of VL in the CT were calculated for the periods of 1993–1996 and 2000–2006. Census tracts with rates above the 3rd quartile (>0.62 cases per 1000 person-years for the period of 1993–1996 and >0.24 cases per 1000 person-years for the period of 2001–2006) were classified as “high risk” and the ones that did not meet the criteria as “low risk” (Figure 1).
Socioeconomic and demographic variables were obtained from the Demographic Census of 1991 and 2000. The 22 indicators analyzed were categorized according to their quantile distribution, 19 of them were categorized by quartiles, two by the median and one by tertiles (Table 1).
The environmental classification was obtained through the extraction of land coverage characteristics of Teresina using a Landsat 5 Thematic Mapper (TM) scene of August 1990 and a Landsat 5 Thematic Mapper (TM) scene of June 2003.
Satellite image classification
Bands 3 (red) and 4 (near-infrared) were used for the image classification, due to larger spectral differences .
The image classification, performed with the software Definiens Developer 7.0, encompassed two basic steps: multiresolution segmentation and algorithm classification using fuzzy and Boolean logics. This approach uses a combination of spectral, textural, and contextual/topologic information .
The rules used to segment and classify the image from 1990 were applied to the image from 2003.
The environmental indicators for both analyzed periods were built by means of calculating the proportion of each thematic class in each census tract (Table 2). This calculation was performed using the program LEGAL in the software SPRING (National Institute of Spatial Research - INPE).
The data from the first period (1993–1996) was used as a learning or sample for the development of the predictive model and the second period (2001–2006) as a validation sample.
The CART algorithm (Classification and Regression Trees) was used to attain the predictive model for CTs at high risk for VL occurrence . Generally CART generates a very large tree which shows a minimum number of classification errors, but it is a model excessively adjusted to the data with a limited capacity of generalization. Therefore, the tree has to be reduced (“trimmed”). To get a smaller tree (smaller number of terminal nodes, or “leaves”) it was initially established a restriction of a minimum of 20 observations to be included before a split was attempted and a minimum of 10 observations in a terminal node. Subsequently, by means of the graphic inspection of the relation between the number of terminal nodes of the tree and gains on the classification homogeneity, a tree with 15 nodes was concluded to be most adequate.
To calculate accuracy measures, the cutoff point of 25% to the probability of finding high risk census tracts on the terminal nodes of the tree was used, which simultaneously maximized the sensitivity and specificity, with no important decrease of the global accuracy.
CART models were implemented on the software Splus 4.5. The evaluatiton of the performance of the models was done based on the calculations of the sensitivity, specificity, accuracy and area under ROC curve indicators; estimated in the software Stata 11.0.
For the period of 1993–1996, more than 70% of the CT reported VL cases, having maximum and median incidence of 4.54 and 0.29 cases per 1000 inhabitants respectively. On the second period analyzed, around 50% of the CT reported VL cases, with maximum and median incidence of 5.91 and zero cases per 1000 inhabitants respectively. The satellite image classification results indicate the expansion of the urban area towards the outskirts of the city, where there was more vegetation coverage (Figure 2).The model developed was capable of discriminating 15 CT sets (Figure 3), with different probabilities of containing CTs at high risk for VL occurrence. This figure shows the 15 sets of CT corresponding to the terminal nodes (in red) and inside each of them the probabilities (P) of the existence of high risk CT and the number of CT (N).
The subset with the lower probability of containing CTs at high risk for VL occurrence (3.9%) included those CTs with the percentage of literate heads of the household higher than the median (>64.2%) and with the income ratio smaller than the median (≤2.49%). The subset with the highest probability of containing CTs at high risk for VL occurrence (92%) encompassed CTs with percentage of literate heads of the household below than the median (≤64.2%), with a larger area covered with dense vegetation, with a percentage of household with up to 3 residents above the third quartile (>31.6%). The other 13 subsets presented the probability of containing CTs at high risk for VL occurrence varying from 5.6% to 82%.
Table 3 shows the model’s performance on the learning (1993–1996) and validation (2001–2006) samples in terms of area under the ROC curve, sensitivity, specificity, accuracy, positive and negative predictive value. The model with 15 terminal nodes (Figure 3) had 79% sensitivity, 74% specificity, 75% global accuracy and 83% area under ROC curve. When applied on the validation sample with the same prediction cutoff (25%), it presented a 52% sensitivity, 66% specificity, 62% global accuracy, and 66% area under the ROC curve (Table 3). ROC curves figures for the validation and learning samples are provided in Additional file 1.
The results of this study reinforce the notion that the spatial heterogeneity found in the occurrence of the disease is directly related to the living conditions of the population and environmental characteristics of household neighborhood . However, though these relations were possible to be characterized in this study, the validation sample results were not fully satisfactory, indicating that the prediction of areas at high risk for VL incidence is a more complex challenge than the simple identification of associations between environmental and socioeconomic factors and disease incidence, which has been already shown in Brazil and in the world [8–14].
A series of potential explanations for the deficient predictive ability of the model developed could be identified. First, the time interval between both periods analyzed implied in substantial modifications on the dimension and structure of the geographical area under study, including the incorporation on the validation sample of areas that were considered rural in the period of 1993–1996. Second, the data used to derive the predictive model refers to an epidemic period, while the validation sample entails endemic years [12, 15]. It is reasonable to suppose that in epidemic situations, due to the typical large magnitude of transmission during these periods, the disease spreads more largely in the geographic space, affecting population subsets that could have eventually been spared in endemic periods. However, having in mind that the VL was only introduced in Teresina in the beginning of the 1980s triggering two epidemics (1981–1985 and 1993–1996), there is still no set of complete data for a typical endemic period that not of 2001–2006. An alternative would be to restrict the analysis strictly to the years between 2001 and 2006, but that would bring difficulties related to the small number of annual cases observed in this period.
Third, the unsatisfactory performance of the model could be due to positive quantitative variations of the socioeconomic indicators, even if the distribution quartiles were used as cutoff points. Since the 1980s in Brazil there has been an improvement on the social and economic indicators. For example, a historical series analysis since the mid-1970s shows a substantial and unequivocal fall of inequality from 2001 to 2004, this last year having the smallest income inequality of the period analyzed . In Piauí the situation is not different, having occurred an increase of 12% of the average income (1985–2006), contributing to the increase in the human development index from 0.57 to 0.70 (1991–2005) and to the reduction of the illiteracy rate from 41.7% to 23.3% (1991–2007) [17, 18]. A possible alternative to minimize this problem would be employing a categorization of the indicators by strata of homogeneous areas according to living conditions [19, 20].
The VL transmission niches in urban environments not only present a heterogeneous distribution, but also constitute areas with varied landscape and epidemiologic characteristics, where distinct forms of occupation and soil coverage implicate in ecological and social processes which result in huge differences of magnitude on the incidence of the disease.
Therefore, it is possible to infer that the process of establishing and dissemination of VL in the urban environment of Teresina, having a markedly heterogeneous spatial distribution, results from the socio-territorial organization. Nevertheless, we believe that part of the predictive deficiencies of the model is owing to the lack of a better definition of the territory as a spatial-social-environmental unit favorable to VL occurrence.
The improvement on the predictive performance of area at high risk for VL could be attained with the use of more environmental indicators extracted from the satellite images of medium resolution [21, 22] or with the use of remote sensing images with high spatial resolution . The satellite images used in this study come from the Landsat Thematic Mapper 5 sensors, which have spatial resolution of 30 meters, in other words, it does not allow the discrimination of elements on the earth’s surface for areas smaller than 900 m2. However, VL studies that used medium spatial resolution were able to identify some elements related to vegetation coverage, soil use and patterns of urban occupation associated to the risk of the disease [12, 22, 24, 25], but the difficulties and limitations of these images to properly characterize the local features that determine the pattern of transmission of VL on the urban context are evident . In future studies, high spatial resolution images should be used to better define the local environmental features related to VL occurrence .
At last, considering the complexity of the factors involved in the VL dissemination in the urban environment, the results of this study demonstrate that the occurrence of VL on the outskirts of Teresina is intensely related to socioeconomic and environmental problems, arising from the urban expansion process and from the changes in the vector habitat, due to the environmental imbalance caused by deforestation and land occupation with lack of adequate urban infrastructure . In this perspective, focusing interventions in these considered high risk areas could be a useful strategy to improve the effectiveness of control measures while decreasing operational costs.
Ministério da Saúde: Manual de Vigilância e Controle da Leishmaniose Visceral. 2006, Brasília, DF: Série A. Normas e Manuais Técnicos
Costa CHN, Tapety CMM, Werneck GL: Controle da leishmaniose visceral em meio urbano: estudo de intervenção randomizado fatorial. Rev Soc Bras Med Trop. 2007, 40 (4): 415-419. 10.1590/S0037-86822007000400009.
Werneck GL: Forum: geographic spread and urbanization of visceral leishmaniasis in Brazil. Introduction Cad Saúde Pública. 2008, 24 (12): 2937-2940.
Koopman JS, Simon CP, Riolo CP: When to control endemic infections by focusing on high-risk groups. Epidemiology. 2005, 16: 621-627. 10.1097/01.ede.0000172133.46385.18.
Zha Y, Gao J, Ni S: Use of normalized difference build-up index in automatically mapping urban areas from TM imagery. Int J Remote Sens. 2003, 24 (17): 583-594.
Jacquin A, Misakova L, Gay M: A hybrid object-based classification approach for mapping urban sprawl in periurban environment. Landsc Urban Plan. 2008, 84: 152-165. 10.1016/j.landurbplan.2007.07.006.
Breiman L, Friedman J, Olshen R, Stone C: Classification and Regression Trees. 1984, Wadsworth & Brooks/Cole: Pacific Grove, CA
Thompson RA, Wellington de Oliveira Lima J, Maguire JH, Braud DH, Scholl DT: Climatic and demographic determinants of American visceral leishmaniasis in northeastern Brazil using remote sensing technology for environmental categorization of rain and region influences on leishmaniasis. Am J Trop Med Hyg. 2002, 67 (6): 648-655.
Costa CHN, Werneck GL, Rodrigues L, Santos MV, Araújo IB, Moura LS, Moreira S, Gomes RB, Lima SS: Household structure and urban services: neglected targets in the control of visceral leishmaniasis. Ann Trop Med Parasitol. 2005, 99 (3): 229-236. 10.1179/136485905X28018.
Sudhakar S, Srinivas T, Palit A, Kar SK, Battacharya SK: Mapping of risk prone areas of kala-azar (Visceral leishmaniasis) in parts of Bihar State, India: an RS and GIS approach. J Vector Borne Dis. 2006, 43 (3): 115-122.
Oliveira CD, Diez-Roux A, César CC, Proietti FA: A case–control study of microenvironmental risk factors for urban visceral leishmaniasis in a large city in Brazil, 1999–2000. Rev Panam Salud Publica. 2006, 20 (6): 369-376. 10.1590/S1020-49892006001100002.
Cerbino Neto J, Werneck GL, Costa CH: Factors associated with the incidence of urban visceral leishmaniasis: an ecological study in Teresina, Piauí State, Brazil. Cad Saude Publica. 2009, 25 (7): 1543-1551. 10.1590/S0102-311X2009000700012.
Boelaert M, Meheus F, Sanchez A, Singh SP, Vanlerberghe V, Picado A, Meessen B, Sundar S: The poorest of the poor: a poverty appraisal of households affected by visceral leishmaniasis in Bihar, India. Trop Med Int Health. 2009, 14 (6): 639-644. 10.1111/j.1365-3156.2009.02279.x.
Bhunia GS, Kumar V, Kumar AJ, Das P, Kesari S: The use of remote sensing in the identification of the eco-environmental factors associated with the risk of human visceral leishmaniasis (kala-azar) on the Gangetic plain, in north-eastern India. Ann Trop Med Parasitol. 2010, 104 (1): 35-53. 10.1179/136485910X12607012373678.
Werneck GL, Pereira TJCF, Farias GC, Silva FO, Chaves FC, Gouvêa MV, Costa CHNC, Carvalho FAA: Avaliação da efetividade das estratégias de controle da leishmaniose visceral na cidade de Teresina, estado do Piauí, Brasil: resultados do inquérito inicial – 2004. Epidemiologia e Serviços de Saúde. 2008, 17 (2): 87-96.
Soares SSD: Distribuição de Renda no Brasil de 1976 a 2004 com Ênfase no Período Entre 2001 e 2004. 2006, Brasília: Ipea. Texto para Discussão no 1166, fev
Fundação CEPRO (Centro de Pesquisas Econômicas e Sociais do Piauí): Piauí em Números. 2011, Teresina: Fundação CEPRO, 8
PNUD (Programa das Nações Unidas para o Desenvolvimento): Atlas do desenvolvimento humano no brasil. 2003, Disponível em: http://www.pnud.org.br/atlas/ Disponível em:
Lacerda JT, Calvo MCM, Freitas SFT: Diferenciais intra-urbanos no município de florianópolis, santa catarina, brasil: potencial de uso para o planejamento em saúde. Cad Saude Publica. 2002, 18 (5): 1331-1338. 10.1590/S0102-311X2002000500026.
Gamarra CJ, Valente JG, Azevedo e Silva G: Magnitude da mortalidade por câncer do colo do útero na Região Nordeste do Brasil e fatores socioeconômicos. Rev Panam Salud Publica. 2010, 28 (2): 100-106.
Van Benthem BHB, Vanwambeke SO, Khantikul N, Burghoorn-Maas C, Panart K, Oskam L, Lambin EF, Somboon P: Spatial patterns of and risk factors for seropositivity for dengue infection. Am J Trop Med Hyg. 2005, 72 (2): 201-208.
Correia VRM, Monteiro AMV, Carvalho MS, Werneck GL: Uma aplicação do sensoriamento remoto para a investigação de endemias urbanas. Cad Saude Publica. 2007, 23 (5): 1015-1028. 10.1590/S0102-311X2007000500004.
Leonardi F, Almeida CN, Fonseca LMG, Camargo FF: Avaliação Comparativa entre Classificação Supervisionada por Regiões e Orientada a Objeto para Imagens de Alta Resolução Espacial: Cbers 2B-HRC e QuickBird. Anais XIV Simpósio Brasileiro de Sensoriamento Remoto. 2009, Natal, Brasil: INPE, 981-988. 25–30 abril
Werneck GL, Maguire JH: Spatial modeling using mixed models: an ecologic study of visceral leishmaniasis in Teresina, Piaui State, Brazil. Cad Saude Publica. 2002, 18: 633-637. 10.1590/S0102-311X2002000300012.
Werneck GL, Costa CHN, Walker AM, David JR, Wand M, Maguire JH: Multilevel modelling of the incidence of visceral leishmaniasis in Teresina, Brazil. Epidemiol Infect. 2007, 135: 195-201. 10.1017/S0950268806006881.
Souza IM, Alves CD, Almeida CM, Pinho CMD: Caracterização socioeconômica do espaço residencial construído utilizando imagens de alta resolução espacial e análise orientada a objeto. Geografia. 2007, 16 (1): 119-142.
The research in this paper was funded by CNPq (160571/2011-1, 202088/2012-0 and 306267/2010-1).
The authors declare that they have no competing interests.
ASA participated in the article’s conceptualization and conducted the literature review, structured the database, analyzed and interpreted the compiled data, and wrote the article. GLW participated in the article’s conceptualization and contributed to the analysis and interpretation of the results and helped write the article. Both authors read and approved the final manuscript.