# Automatic detection of the parasite Trypanosoma cruzi in blood smears using a machine learning approach applied to mobile phone images

Mauro César Cafundó Morais, Diogo Silva, Matheus Marques Milagre, Maykon Tavares De Oliveira, Thaís Pereira, João Santana Silva, Luciano da F. Costa, Paola Minoprio, Roberto Marcondes Cesar Junior, Ricardo Gazzinelli, et al.

► To cite this version:

Mauro César Cafundó Morais, Diogo Silva, Matheus Marques Milagre, Maykon Tavares De Oliveira, Thaís Pereira, et al. Automatic detection of the parasite Trypanosoma cruzi in blood smears using a machine learning approach applied to mobile phone images. PeerJ, 2022, 10, pp.e13470. 10.7717/peerj.13470. hal-04627513

HAL Id: hal-04627513

https://hal.science/hal-04627513

Submitted on 28 Jun 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Automatic detection of the parasite Trypanosoma cruzi in blood smears using a machine learning approach applied to mobile phone images

Mauro César Cafundó Morais, Diogo Silva, Matheus Marques Milagre, Maykon Tavares De Oliveira, Thaís Pereira, João Santana Silva, Luciano da F. Costa, Paola Minoprio, Roberto Marcondes Cesar Junior, Ricardo Gazzinelli, et al.

► To cite this version:

Mauro César Cafundó Morais, Diogo Silva, Matheus Marques Milagre, Maykon Tavares De Oliveira, Thaís Pereira, et al. Automatic detection of the parasite Trypanosoma cruzi in blood smears using a machine learning approach applied to mobile phone images. PeerJ, 2022, 10, pp.e13470. 10.7717/peerj.13470. hal-04627513

HAL Id: hal-04627513

https://hal.science/hal-04627513v1

Submitted on 28 Jun 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Automatic detection of the parasite *Trypanosoma cruzi* in blood smears using a machine learning approach applied to mobile phone images

Mauro César Cafundó Morais<fnref n="1" />,<fnref n="2" />,<fnref n="3" />,*, Diogo Silva<fnref n="3" />,*, Matheus Marques Milagre<fnref n="4" />, Maykon Tavares de Oliveira<fnref n="6" />, Thaís Pereira<fnref n="5" />, João Santana Silva<fnref n="6" />, Luciano da F. Costa<fnref n="7" />, Paola Minoprio<fnref n="2" />, Roberto Marcondes Cesar Junior<fnref n="8" />, Ricardo Gazzinelli<fnref n="5" />, Marta de Lana<fnref n="4" />,<fnref n="9" /> and Helder I. Nakaya<fnref n="1" />,<fnref n="2" />,<fnref n="3" />,<fnref n="10" />

<figure><img src="image_2.png" /></figure>

## ABSTRACT

Chagas disease is a life-threatening illness caused by the parasite *Trypanosoma cruzi*. The diagnosis of the acute form of the disease is performed by trained microscopists who detect parasites in blood smear samples. Since this method requires a dedicated high-resolution camera system attached to the microscope, the diagnostic method is more expensive and often prohibitive for low-income settings. Here, we present a machine learning approach based on a random forest (RF) algorithm for the detection and counting of *T. cruzi* trypomastigotes in mobile phone images. We analyzed micrographs of blood smear samples that were acquired using a mobile device camera capable of capturing images in a resolution of 12 megapixels. We extracted a set of features that describe morphometric parameters (geometry and curvature), as well as color, and texture measurements of 1,314 parasites. The features were divided into train and test sets (4:1) and classified using the RF algorithm. The values of precision, sensitivity, and area under the receiver operating characteristic (ROC) curve of the proposed method were 87.6%, 90.5%, and 0.942, respectively. Automating image analysis acquired with a mobile device is a viable alternative for reducing costs and gaining efficiency in the use of the optical microscope.

<footnote n="1">¹ Hospital Israelita Albert Einstein, São Paulo, Brazil</footnote>

<footnote n="2">² Scientific Platform Pasteur-University of São Paulo (SPPU), Universidade de São Paulo, São Paulo, SP, Brazil</footnote>

<footnote n="3">³ Department of Clinical and Toxicological Analysis, School of Pharmaceutical Sciences, Universidade de São Paulo, Sao Paulo, SP, Brazil</footnote>

<footnote n="4">⁴ Departamento de Análises Clínicas (DEACL), Programa de Pós-graduação em Ciências Farmacêuticas (CiPHARMA), Universidade Federal de Ouro Preto, Ouro Preto, MG, Brazil</footnote>

<footnote n="5">⁵ Laboratório de Imunopatologia, Instituto René Rachou, Fundação Oswaldo Cruz, Universidade Federal de Minas Gerais, Belo Horizonte, MG, Brazil</footnote>

<footnote n="6">⁶ Fiocruz- Bi-Institutional Translational Medicine Project, FIOCRUZ/SP, Ribeirão Preto, SP, Brazil</footnote>

<footnote n="7">⁷ São Carlos Institute of Physics (DFCM- IFSC), Universidade de São Paulo, São Carlos, SP, Brazil</footnote>

<footnote n="8">⁸ Instituto de Matemática e Estatística (IME), Universidade de São Paulo, São Paulo, SP, Brazil</footnote>

<footnote n="9">⁹ Núcleo de Pesquisas em Ciências Biológicas (NUPEB), Universidade Federal de Ouro Preto, Ouro Preto, MG, Brazil</footnote>

<footnote n="10">¹⁰ Center of Research in Inflammatory Diseases (CRID), Universidade de São Paulo, Ribeirão Preto, SP, Brazil</footnote>

<footnote>*These authors contributed equally to this work.</footnote>
# Automatic detection of the parasite *Trypanosoma cruzi* in blood smears using a machine learning approach applied to mobile phone images

Mauro César Cafundó Morais<fnref n="1" />,<fnref n="2" />,<fnref n="3" />,*, Diogo Silva<fnref n="3" />,*, Matheus Marques Milagre<fnref n="4" />, Maykon Tavares de Oliveira<fnref n="6" />, Thaís Pereira<fnref n="5" />, João Santana Silva<fnref n="6" />, Luciano da F. Costa<fnref n="7" />, Paola Minoprio<fnref n="2" />, Roberto Marcondes Cesar Junior<fnref n="8" />, Ricardo Gazzinelli<fnref n="5" />, Marta de Lana<fnref n="4" />,<fnref n="9" /> and Helder I. Nakaya<fnref n="1" />,<fnref n="2" />,<fnref n="3" />,<fnref n="10" />

<figure><img src="image_2.png" /></figure>

## ABSTRACT

Chagas disease is a life-threatening illness caused by the parasite *Trypanosoma cruzi*. The diagnosis of the acute form of the disease is performed by trained microscopists who detect parasites in blood smear samples. Since this method requires a dedicated high-resolution camera system attached to the microscope, the diagnostic method is more expensive and often prohibitive for low-income settings. Here, we present a machine learning approach based on a random forest (RF) algorithm for the detection and counting of *T. cruzi* trypomastigotes in mobile phone images. We analyzed micrographs of blood smear samples that were acquired using a mobile device camera capable of capturing images in a resolution of 12 megapixels. We extracted a set of features that describe morphometric parameters (geometry and curvature), as well as color, and texture measurements of 1,314 parasites. The features were divided into train and test sets (4:1) and classified using the RF algorithm. The values of precision, sensitivity, and area under the receiver operating characteristic (ROC) curve of the proposed method were 87.6%, 90.5%, and 0.942, respectively. Automating image analysis acquired with a mobile device is a viable alternative for reducing costs and gaining efficiency in the use of the optical microscope.

<footnote n="1">¹ Hospital Israelita Albert Einstein, São Paulo, Brazil</footnote>

<footnote n="2">² Scientific Platform Pasteur-University of São Paulo (SPPU), Universidade de São Paulo, São Paulo, SP, Brazil</footnote>

<footnote n="3">³ Department of Clinical and Toxicological Analysis, School of Pharmaceutical Sciences, Universidade de São Paulo, Sao Paulo, SP, Brazil</footnote>

<footnote n="4">⁴ Departamento de Análises Clínicas (DEACL), Programa de Pós-graduação em Ciências Farmacêuticas (CiPHARMA), Universidade Federal de Ouro Preto, Ouro Preto, MG, Brazil</footnote>

<footnote n="5">⁵ Laboratório de Imunopatologia, Instituto René Rachou, Fundação Oswaldo Cruz, Universidade Federal de Minas Gerais, Belo Horizonte, MG, Brazil</footnote>

<footnote n="6">⁶ Fiocruz- Bi-Institutional Translational Medicine Project, FIOCRUZ/SP, Ribeirão Preto, SP, Brazil</footnote>

<footnote n="7">⁷ São Carlos Institute of Physics (DFCM- IFSC), Universidade de São Paulo, São Carlos, SP, Brazil</footnote>

<footnote n="8">⁸ Instituto de Matemática e Estatística (IME), Universidade de São Paulo, São Paulo, SP, Brazil</footnote>

<footnote n="9">⁹ Núcleo de Pesquisas em Ciências Biológicas (NUPEB), Universidade Federal de Ouro Preto, Ouro Preto, MG, Brazil</footnote>

<footnote n="10">¹⁰ Center of Research in Inflammatory Diseases (CRID), Universidade de São Paulo, Ribeirão Preto, SP, Brazil</footnote>

<footnote>*These authors contributed equally to this work.</footnote>
**Subjects** Bioinformatics, Computational Biology, Parasitology, Infectious Diseases, Computational Science

**Keywords** *Trypanosoma cruzi*, Blood trypomastigote, Parasitemia, Machine learning, SVM

INTRODUCTION

Chagas disease is a life-threatening illness caused by infection with the protozoan *Trypanosoma cruzi* (Chagas, 1909). Most of the cases occur when metacyclic trypomastigotes eliminated in the feces or urine of the vector enter the human host. Infection may also develop through oral ingestion of contaminated food, congenital transmission, blood transfusion, transplants of organs, and laboratory accidents (Cancino-Faure et al., 2015; Filigheddu, Górgolas & Ramos, 2017; Luquetti et al., 2015). After penetrating the host cells, the metacyclic trypomastigotes differentiate into amastigote forms in the cytoplasm. Subsequently, the amastigotes multiply themselves by binary fission and transform into trypomastigotes that disrupt the host cell, releasing into the bloodstream. These circulating trypomastigotes may invade other host cells or be ingested by vectors (Lana & Machado, 2017).

The acute phase of Chagas disease is characterized by high parasitemia in the blood (Dias et al., 2016; Luquetti & Schmuñis, 2017). This allows the visualization of bloodstream forms in the blood of infected individuals using a parasitological fresh-blood test, as well as smear and thick drop blood tests (Gomes, Lorena & Luquetti, 2009). Laboratory diagnosis, however, has a few key limitations. First, it should be performed by trained microscopists who observe the parasites. Reliance on professionals with various skills makes the diagnosis prone to errors and heterogeneous. Second, manual search and detection of parasites is a laborious task. This often delays the laboratory result, which further delays the initiation of treatment. Third, methods that improve the search for the parasite in microscopic images, such as attaching a dedicated high-resolution camera system to the microscope, are usually expensive and often prohibitive for low-income settings.

Machine learning (ML) algorithms can assist in the laboratory diagnosis of acute Chagas disease. They can automate the search and detection of the parasites, improving the reproducibility of image analysis. These algorithms have been applied to image detection of *T. cruzi* as well as other parasites that circulate in the blood (Rajaraman et al., 2018; Uc-Cetina, Brito-Loeza & Ruiz-Piña, 2015; Górriz et al., 2018). Previous work relied on an image acquisition system with a dedicated camera attached to the microscope for the detection of *T. cruzi*. This system is important for proper ML training since, it produces homogeneous images regarding the lightness, color, acquisition time, aperture, and resolution. However, these dedicated cameras are often expensive, which makes such system prohibitive to low-income settings. The use of mobile device cameras increases image diversity, which can make ML models even more robust. Models trained in mobile device imaging were developed for the detection of *Plasmodium ssp.* which causes malaria (Rosado et al., 2016; Oliveira et al., 2017; Yang et al., 2020). To date, models for the detection of *T. cruzi* in images acquired with mobile devices have not been developed.
**Subjects** Bioinformatics, Computational Biology, Parasitology, Infectious Diseases, Computational Science

**Keywords** *Trypanosoma cruzi*, Blood trypomastigote, Parasitemia, Machine learning, SVM

## INTRODUCTION

Chagas disease is a life-threatening illness caused by infection with the protozoan *Trypanosoma cruzi* (Chagas, 1909). Most of the cases occur when metacyclic trypomastigotes eliminated in the feces or urine of the vector enter the human host. Infection may also develop through oral ingestion of contaminated food, congenital transmission, blood transfusion, transplants of organs, and laboratory accidents (Cancino-Faure et al., 2015; Filigheddu, Górgolas & Ramos, 2017; Luquetti et al., 2015). After penetrating the host cells, the metacyclic trypomastigotes differentiate into amastigote forms in the cytoplasm. Subsequently, the amastigotes multiply themselves by binary fission and transform into trypomastigotes that disrupt the host cell, releasing into the bloodstream. These circulating trypomastigotes may invade other host cells or be ingested by vectors (Lana & Machado, 2017).

The acute phase of Chagas disease is characterized by high parasitemia in the blood (Dias et al., 2016; Luquetti & Schmuñis, 2017). This allows the visualization of bloodstream forms in the blood of infected individuals using a parasitological fresh-blood test, as well as smear and thick drop blood tests (Gomes, Lorena & Luquetti, 2009). Laboratory diagnosis, however, has a few key limitations. First, it should be performed by trained microscopists who observe the parasites. Reliance on professionals with various skills makes the diagnosis prone to errors and heterogeneous. Second, manual search and detection of parasites is a laborious task. This often delays the laboratory result, which further delays the initiation of treatment. Third, methods that improve the search for the parasite in microscopic images, such as attaching a dedicated high-resolution camera system to the microscope, are usually expensive and often prohibitive for low-income settings.

Machine learning (ML) algorithms can assist in the laboratory diagnosis of acute Chagas disease. They can automate the search and detection of the parasites, improving the reproducibility of image analysis. These algorithms have been applied to image detection of *T. cruzi* as well as other parasites that circulate in the blood (Rajaraman et al., 2018; Uc-Cetina, Brito-Loeza & Ruiz-Piña, 2015; Górriz et al., 2018). Previous work relied on an image acquisition system with a dedicated camera attached to the microscope for the detection of *T. cruzi*. This system is important for proper ML training since, it produces homogeneous images regarding the lightness, color, acquisition time, aperture, and resolution. However, these dedicated cameras are often expensive, which makes such system prohibitive to low-income settings. The use of mobile device cameras increases image diversity, which can make ML models even more robust. Models trained in mobile device imaging were developed for the detection of *Plasmodium ssp.* which causes malaria (Rosado et al., 2016; Oliveira et al., 2017; Yang et al., 2020). To date, models for the detection of *T. cruzi* in images acquired with mobile devices have not been developed.
Here, we developed and evaluated a ML algorithm that detects the bloodstream forms of *T. cruzi* based on features of segmented body from the parasites. The extracted features consisted in descriptors of morphology (geometric), color and texture, as well as statistical descriptors (Hu's invariant moments). We tested our method using images acquired with a mobile phone from the acute phase of infection in a murine model. Our model reached 89.5% accuracy in a set of images that were not previously presented in the training process. The proposed method presented an acceptable performance to detect trypomastigotes using mobile phone images of blood smear.

# MATERIALS & METHODS

## Samples analyzed

A total of 33 slides with thin blood smears of Swiss mice experimentally infected with *T. cruzi* Y strain at acute phase of infection were prepared for image annotation and analysis. Sample preparations were obtained from animals of the Laboratory of Chagas disease, Federal University of Ouro Preto where the *T. cruzi* strain was maintained through successive blood passages in mice. The Ethics Committee for the Use of Animals (CEUA) at the Federal University of Ouro Preto, Laboratory of Chagas Disease, Minas Gerais, Brazil provided approval for this research (CEUA No. 2015/50).

## Object segmentation

We standardized the resolution of the images to 768 × 1,024 pixels² before segmenting the parasites. We applied a graph-based segmentation method (Felzenszwalb & Huttenlocher, 2004). In this process, a graph-based representation of the image is defined where pixels corresponds to vertices and neighboring pixels are connected edges. The contours of the regions of interest with the parasites are obtained by selecting the edges between the different regions based on the differences in intensity between the regions, and the difference in intensity between the pixels within each region. As a result, we observed the whole parasite cell body segmented within each region.

Next, we cropped a 100 × 100 pixels² region around each parasite based on manual position annotation. Only the regions of interest with the parasite were selected for processing and feature extraction. This procedure resulted in 1,314 parasites. We selected other segments from the images with features very similar to the parasite using nearest neighbors technique. In other words, we selected objects that do not belong to parasites and clustered under the label “Unknown”. In this way, we were able to obtain the same number of objects (*T. cruzi* or unknown) for the object classification task.

## Feature extraction

After the segmentation of the parasites, we performed the conversion from the RGB color space to CIEL*a*b* color space for object's feature extraction (see the *Supplemental Information* file for details). These features are object descriptors and are classified as geometric (perimeter, area, circularity, thickness ratio, centroid to contour distances, major and minor axis, aspect ratio, etc.); curvature (entropy, bending energy, standard deviation, and variance) (Costa & Cesar Jr, 2009); texture (color co-occurrence matrices,
Here, we developed and evaluated a ML algorithm that detects the bloodstream forms of *T. cruzi* based on features of segmented body from the parasites. The extracted features consisted in descriptors of morphology (geometric), color and texture, as well as statistical descriptors (Hu's invariant moments). We tested our method using images acquired with a mobile phone from the acute phase of infection in a murine model. Our model reached 89.5% accuracy in a set of images that were not previously presented in the training process. The proposed method presented an acceptable performance to detect trypomastigotes using mobile phone images of blood smear.

# MATERIALS & METHODS

## Samples analyzed

A total of 33 slides with thin blood smears of Swiss mice experimentally infected with *T. cruzi* Y strain at acute phase of infection were prepared for image annotation and analysis. Sample preparations were obtained from animals of the Laboratory of Chagas disease, Federal University of Ouro Preto where the *T. cruzi* strain was maintained through successive blood passages in mice. The Ethics Committee for the Use of Animals (CEUA) at the Federal University of Ouro Preto, Laboratory of Chagas Disease, Minas Gerais, Brazil provided approval for this research (CEUA No. 2015/50).

## Object segmentation

We standardized the resolution of the images to 768 × 1,024 pixels² before segmenting the parasites. We applied a graph-based segmentation method (Felzenszwalb & Huttenlocher, 2004). In this process, a graph-based representation of the image is defined where pixels corresponds to vertices and neighboring pixels are connected edges. The contours of the regions of interest with the parasites are obtained by selecting the edges between the different regions based on the differences in intensity between the regions, and the difference in intensity between the pixels within each region. As a result, we observed the whole parasite cell body segmented within each region.

Next, we cropped a 100 × 100 pixels² region around each parasite based on manual position annotation. Only the regions of interest with the parasite were selected for processing and feature extraction. This procedure resulted in 1,314 parasites. We selected other segments from the images with features very similar to the parasite using nearest neighbors technique. In other words, we selected objects that do not belong to parasites and clustered under the label “Unknown”. In this way, we were able to obtain the same number of objects (*T. cruzi* or unknown) for the object classification task.

## Feature extraction

After the segmentation of the parasites, we performed the conversion from the RGB color space to CIEL*a*b* color space for object's feature extraction (see the Supplemental Information file for details). These features are object descriptors and are classified as geometric (perimeter, area, circularity, thickness ratio, centroid to contour distances, major and minor axis, aspect ratio, etc.); curvature (entropy, bending energy, standard deviation, and variance) (Costa & Cesar Jr, 2009); texture (color co-occurrence matrices,
entropy, inverse difference moment, angular second moment, contrast, correlation) color (mean, median, mode, amplitude, and variance) (Gui et al., 2013; Palm, 2004) descriptors and Hu's invariant moments (Huang & Leng, 2010).

## Feature selection

To increase performance, reduce the noise, and avoid overfitting, we applied either Principal Component Analysis (PCA). In this manner, we reduced the number of features used by the algorithm during the training. The dimension reduction with PCA was performed by keeping the 16 highest eigenvalues of the covariance matrix derived from the feature space data. The proportion of variance of the 16 principal components correspond to 95% of the original variance. Therefore, we obtained a new feature space data that consisted of eigenvectors of the principal components. We obtained the final dataset by multiplying the transposed matrix of eigenvectors by the input data matrix centered by the mean.

## Object classification

We applied supervised learning classification with support vector machines (SVM), k-Nearest Neighbours (KNN) and Random Forest (RF) as they presents good generalization with a small dataset (Chen et al., 2020; Hsu, Chang & Lin, 2016; Ben-Hur & Weston, 2010; Cunningham & Delany, 2007; Breiman, 2001).

In the SVM method, each sample is represented by a point in an *n*-dimensional space, where *n* is the number of object features and its values are the coordinates of the point, including its class. The classifier searches for the optimal hyperplane that will be used to split the points belonging to distinct classes. In the KNN method, a sample is classified to the data label which has the most representatives within the *k* nearest neighbors features of the sample. In the RF method, the classification of each object happens through a combination of decision trees. The classification performance of each tree in the ensemble is improved by bootstrapping sampling and aggregation from the training set. The final classification is made by averaging each decision tree probabilistic prediction. In addition, we also performed classification task using an ensemble of these methods based on “soft” voting classifier. The voting classifier combines the results of the different predictors and use the average of predicted probabilities to assign the class labels to each sample.

We used Python's scikit-learn library to train and validate the models (Pedregosa et al., 2011). We performed object feature data standardization before classification since data presented different orders of magnitude. We assessed the classification performance of the methods by comparing the samples in the feature space for a given class with all other samples. After finding the highest accuracy value for each model, we assessed the models in the validation set.

## Statistical analysis

The performance of the constructed models was assessed by Receiver Operating Characteristic (ROC) curve analysis, where the average of the area under the curve (AUC) was calculated for the quantification in both training set and validation sets. We also evaluated the model performance based on sensitivity, specificity, precision, and F1-score.
entropy, inverse difference moment, angular second moment, contrast, correlation) color (mean, median, mode, amplitude, and variance) (Gui et al., 2013; Palm, 2004) descriptors and Hu's invariant moments (Huang & Leng, 2010).

## Feature selection

To increase performance, reduce the noise, and avoid overfitting, we applied either Principal Component Analysis (PCA). In this manner, we reduced the number of features used by the algorithm during the training. The dimension reduction with PCA was performed by keeping the 16 highest eigenvalues of the covariance matrix derived from the feature space data. The proportion of variance of the 16 principal components correspond to 95% of the original variance. Therefore, we obtained a new feature space data that consisted of eigenvectors of the principal components. We obtained the final dataset by multiplying the transposed matrix of eigenvectors by the input data matrix centered by the mean.

## Object classification

We applied supervised learning classification with support vector machines (SVM), k-Nearest Neighbours (KNN) and Random Forest (RF) as they presents good generalization with a small dataset (Chen et al., 2020; Hsu, Chang & Lin, 2016; Ben-Hur & Weston, 2010; Cunningham & Delany, 2007; Breiman, 2001).

In the SVM method, each sample is represented by a point in an *n*-dimensional space, where *n* is the number of object features and its values are the coordinates of the point, including its class. The classifier searches for the optimal hyperplane that will be used to split the points belonging to distinct classes. In the KNN method, a sample is classified to the data label which has the most representatives within the *k* nearest neighbors features of the sample. In the RF method, the classification of each object happens through a combination of decision trees. The classification performance of each tree in the ensemble is improved by bootstrapping sampling and aggregation from the training set. The final classification is made by averaging each decision tree probabilistic prediction. In addition, we also performed classification task using an ensemble of these methods based on “soft” voting classifier. The voting classifier combines the results of the different predictors and use the average of predicted probabilities to assign the class labels to each sample.

We used Python's scikit-learn library to train and validate the models (Pedregosa et al., 2011). We performed object feature data standardization before classification since data presented different orders of magnitude. We assessed the classification performance of the methods by comparing the samples in the feature space for a given class with all other samples. After finding the highest accuracy value for each model, we assessed the models in the validation set.

## Statistical analysis

The performance of the constructed models was assessed by Receiver Operating Characteristic (ROC) curve analysis, where the average of the area under the curve (AUC) was calculated for the quantification in both training set and validation sets. We also evaluated the model performance based on sensitivity, specificity, precision, and F1-score.
<figure><img src="image_3.png" /><figcaption>A. Sample preparation</figcaption></figure>

<figure><img src="image_4.png" /><figcaption>B. Image acquisition</figcaption></figure>

<figure><img src="image_5.png" /><figcaption>C. Parasite labeling and
Image processing</figcaption></figure>

<figure><img src="image_6.png" /><figcaption>E. Machine learning model
training and validation</figcaption></figure>

<figure><img src="image_7.png" /><figcaption>F. Parasite detection</figcaption></figure>

<figure><img src="image_8.png" /></figure>

**Figure 1** *T. cruzi* detection image analysis pipeline. (A) Blood smear samples were prepared from mice experimentally infected with *T. cruzi* parasites at the acute infection stage. (B) Images of thin blood smear slides were acquired with a mobile phone camera attached to a microscope ocular lens. (C) Parasite (trypomatigote forms of *T. cruzi*) was segmented by a graph-based algorithm. (D) Images were converted to CIEL*a*b* color space and parasite features were extracted and selected (PCA). (E) Objects feature data were split into training and test sets. Four machine learning models were trained and assessed. (F) Parasites were detected in mobile phone camera images.

Full-size

<figure><img src="image_9.png" /></figure>

DOI: 10.7717/peerj.13470/fig-1

RESULTS

Image analysis pipeline for detecting *T. cruzi*

We developed and tested an image analysis pipeline which was based on a ML model to detect *T. cruzi* in the blood during acute infection (Fig. 1). Initially, the images were prepared from blood samples collected from female Swiss mice acutely infected with the Y strain of *Trypanosoma cruzi*. The collected blood sample smears were stained by the Giemsa method (Vallada, 1999). This technique allows the observation of parasites with oil immersion objectives. Regions of nucleic acid-rich in adenine-thymine bond tend to get darker. On the other hand, regions rich in cytosine-guanine bonds are less prone to embody the Giemsa stain and tend to present clearer stains.

A total of 33 slides with thin blood smears from different animals were stained with Giemsa. About 20 fields of view that correspond to each of the analyzed images were captured from each slide. We manually acquired 674 images from the slides under 100x objective (CFI E Plan Chromat 100x Oil, 1.25 NA/0.23 W.D.) in an optical microscope (Nikon Eclipse E200) with a cell phone camera (Morola Moto G4) attached to the eyepiece (CFI E 10x, F.N. 20 mm) (Fig. 2). The camera was configured with the macro focus and other configurations were set to automatic for acquisition. With these settings, the images were acquired with a resolution of 3,456 × 4,608 pixels², a field of view with a diameter
<figure><img src="image_3.png" /><figcaption>A. Sample preparation</figcaption></figure>

<figure><img src="image_4.png" /><figcaption>B. Image acquisition</figcaption></figure>

<figure><img src="image_5.png" /><figcaption>C. Parasite labeling and
Image processing</figcaption></figure>

<figure><img src="image_6.png" /><figcaption>E. Machine learning model
training and validation</figcaption></figure>

<figure><img src="image_7.png" /><figcaption>F. Parasite detection</figcaption></figure>

<figure><img src="image_8.png" /></figure>

**Figure 1** *T. cruzi* detection image analysis pipeline. (A) Blood smear samples were prepared from mice experimentally infected with *T. cruzi* parasites at the acute infection stage. (B) Images of thin blood smear slides were acquired with a mobile phone camera attached to a microscope ocular lens. (C) Parasite (trypomatigote forms of *T. cruzi*) was segmented by a graph-based algorithm. (D) Images were converted to CIEL*a*b* color space and parasite features were extracted and selected (PCA). (E) Objects feature data were split into training and test sets. Four machine learning models were trained and assessed. (F) Parasites were detected in mobile phone camera images.

Full-size

<figure><img src="image_9.png" /></figure>

DOI: 10.7717/peerj.13470/fig-1

RESULTS

Image analysis pipeline for detecting *T. cruzi*

We developed and tested an image analysis pipeline which was based on a ML model to detect *T. cruzi* in the blood during acute infection (Fig. 1). Initially, the images were prepared from blood samples collected from female Swiss mice acutely infected with the Y strain of *Trypanosoma cruzi*. The collected blood sample smears were stained by the Giemsa method (Vallada, 1999). This technique allows the observation of parasites with oil immersion objectives. Regions of nucleic acid-rich in adenine-thymine bond tend to get darker. On the other hand, regions rich in cytosine-guanine bonds are less prone to embody the Giemsa stain and tend to present clearer stains.

A total of 33 slides with thin blood smears from different animals were stained with Giemsa. About 20 fields of view that correspond to each of the analyzed images were captured from each slide. We manually acquired 674 images from the slides under 100x objective (CFI E Plan Chromat 100x Oil, 1.25 NA/0.23 W.D.) in an optical microscope (Nikon Eclipse E200) with a cell phone camera (Morola Moto G4) attached to the eyepiece (CFI E 10x, F.N. 20 mm) (Fig. 2). The camera was configured with the macro focus and other configurations were set to automatic for acquisition. With these settings, the images were acquired with a resolution of 3,456 × 4,608 pixels², a field of view with a diameter
A.

<figure><img src="image_10.png" /></figure>

<figure><img src="image_11.png" /></figure>

B.

<figure><img src="image_12.png" /></figure>

C.

<figure><img src="image_13.png" /><figcaption>**Figure 2** Mobile phone attached to optical microscope objective lens and image acquisitio. (A) The mobile phone was attached to the microscope ocular lens (eyepiece) with a plastic support device (left). The camera was configured with the macro focus for acquisition. Other configurations were set to automatic. (B) Image of field-of-view from blood smear slide of mice infected with *Trypanosoma cruzi*. The red squares indicate regions with the presence of a parasite. (C) Crop of a field of view with the *T. cruzi* parasite at the center. K, kinetoplast; N, nucleus; Scale bar, 10 µm.</figcaption></figure>

Full-size

<figure><img src="image_14.png" /></figure>

DOI: 10.7717/peerj.13470/fig-2

of 0.2 mm resulted in a pixel lateral size of approximately 0.06 µm. A few images were acquired at 2,448 × 3,264 pixels² due to change in camera configuration. This was useful for classifier test at different input image resolutions. Images were stored in JPEG format, at 100% quality, and file names standardized according to unique identifiers.

Parasites were microscopically identified by two specialist researchers in *T. cruzi*. A total of 1,314 parasites were observed. We then marked the position of each nucleus or cell body of the *T. cruzi* found in the images. The position of the objects of interest (parasite) in the image was obtained by a point-and-click event using an “in-house software”. The image identifier, the pointed object, and the coordinates in X and Y axis information were extracted and stored in a database (Data S1).

To detect the objects of interest, we applied a graph-based approach (Felzenszwalb & Huttenlocher, 2004). Segmented regions with more than 3,000 pixels were considered not to contain the objects of interest. After image segmentation, objects were cropped in a 100 ×
A.

<figure><img src="image_10.png" /></figure>

<figure><img src="image_11.png" /></figure>

B.

<figure><img src="image_12.png" /></figure>

C.

<figure><img src="image_13.png" /><figcaption>**Figure 2** Mobile phone attached to optical microscope objective lens and image acquisitio. (A) The mobile phone was attached to the microscope ocular lens (eyepiece) with a plastic support device (left). The camera was configured with the macro focus for acquisition. Other configurations were set to automatic. (B) Image of field-of-view from blood smear slide of mice infected with *Trypanosoma cruzi*. The red squares indicate regions with the presence of a parasite. (C) Crop of a field of view with the *T. cruzi* parasite at the center. K, kinetoplast; N, nucleus; Scale bar, 10 µm.</figcaption></figure>

Full-size DOI: 10.7717/peerj.13470/fig-2

of 0.2 mm resulted in a pixel lateral size of approximately 0.06 µm. A few images were acquired at 2,448 × 3,264 pixels² due to change in camera configuration. This was useful for classifier test at different input image resolutions. Images were stored in JPEG format, at 100% quality, and file names standardized according to unique identifiers.

Parasites were microscopically identified by two specialist researchers in *T. cruzi*. A total of 1,314 parasites were observed. We then marked the position of each nucleus or cell body of the *T. cruzi* found in the images. The position of the objects of interest (parasite) in the image was obtained by a point-and-click event using an “in-house software”. The image identifier, the pointed object, and the coordinates in X and Y axis information were extracted and stored in a database (Data S1).

To detect the objects of interest, we applied a graph-based approach (Felzenszwalb & Huttenlocher, 2004). Segmented regions with more than 3,000 pixels were considered not to contain the objects of interest. After image segmentation, objects were cropped in a 100 ×
<figure><img src="image_15.png" /><figcaption>**Figure 3** **Object segmentation.** (A) Original image acquired with a mobile phone attached to the microscope. (B) Segmented image with regions highlighted with different colors. Yellow squares indicate the location of the parasite. (C) Segmented parasites in a 100 × 100 pixels². Top row: *T. cruzi* trypomastigotes from original image. Middle row: segmented regions with parasites. Bottom row: segmented parasites within the segmented region of interest highlighted. Only the regions segmented with the parasites were selected for feature extraction.</figcaption></figure>

Full-size

<figure><img src="image_16.png" /></figure>

DOI: 10.7717/peerj.13470/fig-3

100 pixel² sub-region around its X and Y coordinates and labeled into two classes: parasite and unknown (Fig. 3). The segmented regions containing the *T. cruzi* trypomastigote form were labeled as “parasite” (Fig. 4A). The segmented regions that do not contain a parasite or that are over-segmented were labeled as “unknown” (Fig. 4B).

We selected a set of regions labeled as “unknown” to train and validate the classifier method. The regions were selected based on the features values closest to the regions labeled as “parasite” using the nearest neighbors method. In this way we achieved a total of 1,314 segments marked as parasite and the same number of segments marked as unknown (Table 1).

To identify the parasite in the image sub-region, we first extracted features from these regions. These features represent a description of the object's morphology (geometry and curvature), as well as color and texture. We also calculated Hu's invariant moments to capture information regarding shape and intensity regardless of the object's position and size (Table 2). In total, we extracted 49 features from the segmented objects of each class (Data S2).

We then split the region's feature data into two sets based on the number of acquired images. This resulted in a proportion of about 80% of regions for the training set (n = 2181) and 20% for test (n = 447). We also applied principal component analysis (PCA) to reduce the number of features and test whether it may improve the classification performance. The proportion of variance of the 16 principal components corresponds to 95% of the original variance of the data (Fig. S1). Therefore, we used 16 features of the transformed values matrix to train the model.
<figure><img src="image_14.png" /><figcaption>**Figure 3** **Object segmentation.** (A) Original image acquired with a mobile phone attached to the microscope. (B) Segmented image with regions highlighted with different colors. Yellow squares indicate the location of the parasite. (C) Segmented parasites in a 100 × 100 pixels². Top row: *T. cruzi* trypomastigotes from original image. Middle row: segmented regions with parasites. Bottom row: segmented parasites within the segmented region of interest highlighted. Only the regions segmented with the parasites were selected for feature extraction.</figcaption></figure>

Full-size

<figure><img src="image_15.png" /></figure>

DOI: 10.7717/peerj.13470/fig-3

100 pixel² sub-region around its X and Y coordinates and labeled into two classes: parasite and unknown (Fig. 3). The segmented regions containing the *T. cruzi* trypomastigote form were labeled as “parasite” (Fig. 4A). The segmented regions that do not contain a parasite or that are over-segmented were labeled as “unknown” (Fig. 4B).

We selected a set of regions labeled as “unknown” to train and validate the classifier method. The regions were selected based on the features values closest to the regions labeled as “parasite” using the nearest neighbors method. In this way we achieved a total of 1,314 segments marked as parasite and the same number of segments marked as unknown (Table 1).

To identify the parasite in the image sub-region, we first extracted features from these regions. These features represent a description of the object's morphology (geometry and curvature), as well as color and texture. We also calculated Hu's invariant moments to capture information regarding shape and intensity regardless of the object's position and size (Table 2). In total, we extracted 49 features from the segmented objects of each class (Data S2).

We then split the region's feature data into two sets based on the number of acquired images. This resulted in a proportion of about 80% of regions for the training set (n = 2181) and 20% for test (n = 447). We also applied principal component analysis (PCA) to reduce the number of features and test whether it may improve the classification performance. The proportion of variance of the 16 principal components corresponds to 95% of the original variance of the data (Fig. S1). Therefore, we used 16 features of the transformed values matrix to train the model.
<figure><img src="image_17.png" /><figcaption>**Figure 4** **Parasite's segmentation and region labeling.** (A) Example of segmented regions that contain the parasite. The segmented regions containing the *T. cruzi* trypomastigote form were labeled as "parasite". (B) The segmented regions that do not contain a parasite or that are over-segmented were labeled as "unknown".</figcaption></figure>

Full-size

DOI: 10.7717/peerj.13470/fig-4

<table><caption>**Table 1** The number of objects by classes used in the training and test sets.</caption><thead><tr><td>Class</td><td>Training</td><td>Test</td></tr></thead><tbody><tr><td>Parasites</td><td>1103</td><td>211</td></tr><tr><td>Unknown</td><td>1078</td><td>236</td></tr><tr><td>Total</td><td>2181</td><td>447</td></tr></tbody></table>

## Object classification task presented acceptable performance

We developed a classifier algorithm based on the object features. We observed better performance in the classification task without the feature selection. The models trained on features data selected by PCA were not able to generalize well on the test set, indicating overfitting (Table S1). Without the feature selection, the random forest classifier model presented acceptable performance with an accuracy of 99.7% and area under the ROC curve of 1.0 in the training set, all the while presenting accuracy of 89.5% and AUC of 0.942 in the test set (Table 3). The voting classifier presented accuracy of 93.3% and AUC of 0.978 in the training set, and accuracy of 79.4% and AUC of 0.884 in the test set (Fig. 5).

The voting classifier confusion matrix presented sensitivity and specificity values of 76.8% and 81.8%, respectively (Table 4). The lower performance presented by the
<figure><img src="image_16.png" /><figcaption>**Figure 4** **Parasite's segmentation and region labeling.** (A) Example of segmented regions that contain the parasite. The segmented regions containing the *T. cruzi* trypomastigote form were labeled as "parasite". (B) The segmented regions that do not contain a parasite or that are over-segmented were labeled as "unknown".</figcaption></figure>

Full-size

DOI: 10.7717/peerj.13470/fig-4

<table><caption>**Table 1** The number of objects by classes used in the training and test sets.</caption><thead><tr><td>Class</td><td>Training</td><td>Test</td></tr></thead><tbody><tr><td>Parasites</td><td>1103</td><td>211</td></tr><tr><td>Unknown</td><td>1078</td><td>236</td></tr><tr><td>Total</td><td>2181</td><td>447</td></tr></tbody></table>

## Object classification task presented acceptable performance

We developed a classifier algorithm based on the object features. We observed better performance in the classification task without the feature selection. The models trained on features data selected by PCA were not able to generalize well on the test set, indicating overfitting (Table S1). Without the feature selection, the random forest classifier model presented acceptable performance with an accuracy of 99.7% and area under the ROC curve of 1.0 in the training set, all the while presenting accuracy of 89.5% and AUC of 0.942 in the test set (Table 3). The voting classifier presented accuracy of 93.3% and AUC of 0.978 in the training set, and accuracy of 79.4% and AUC of 0.884 in the test set (Fig. 5).

The voting classifier confusion matrix presented sensitivity and specificity values of 76.8% and 81.8%, respectively (Table 4). The lower performance presented by the
<table><caption>**Table 2** Object feature metrics.</caption><thead><tr><th>Feature</th><th>Description</th><th>References</th></tr></thead><tbody><tr><td><strong>Geometric</strong></td><td></td><td></td></tr><tr><td>Perimeter (P)</td><td>Parametric representation of the contour and its points identified by the coordinates x(t) and y(t)</td><td>Costa & Cesar (2009)</td></tr><tr><td>Area (A)</td><td>Integral of the contour</td><td></td></tr><tr><td>Area and perimeter ratio</td><td>A/P</td><td></td></tr><tr><td>Circularity</td><td>4π A/p²</td><td></td></tr><tr><td>Thickness ratio</td><td>p²/A</td><td></td></tr><tr><td>Centroid</td><td>Given the center of mass M of a contour of complex signal u(n), the centroid coordinates (z₁, z₂) were obtained by the average of all the points in u(n).</td><td></td></tr><tr><td>Centroid to contour maximum distance</td><td>The distance between the centroid and the furthest point on the contour.</td><td></td></tr><tr><td>Centroid to contour minimal distance</td><td>The distance between the centroid and the nearest point on the contour.</td><td></td></tr><tr><td>Centroid to contour average distance</td><td>The average of the distances between the centroid and all points in the contour.</td><td></td></tr><tr><td>Major axis</td><td>Pair of more distant points belonging to the object.</td><td></td></tr><tr><td>Minor axis</td><td>Pair of closest points belonging to the object.</td><td></td></tr><tr><td>Aspect ratio</td><td>Majoraxis/Minoraxis</td><td></td></tr><tr><td>Perimeter and Major axis ratio</td><td>P/Majoraxis</td><td></td></tr><tr><td>Bilateral symmetry</td><td>Bilateral symmetry is given by the proportion of the number of pixels between the intersection of an object and its reflecting shape with respect to the major axis, and the union between those two objects.</td><td>Hu (1962); Huang & Leng (2010)</td></tr><tr><td><strong>Hu's invariant moments<sup>a</sup></strong></td><td></td><td></td></tr><tr><td>φ<sub>1</sub></td><td>η<sub>20</sub> + η<sub>02</sub></td><td></td></tr><tr><td>φ<sub>2</sub></td><td>(η<sub>20</sub> - η<sub>02</sub>)<sup>2</sup> + 4η<sub>11</sub><sup>2</sup></td><td></td></tr><tr><td>φ<sub>3</sub></td><td>(η|30 - 3η<sub>12</sub>)<sup>2</sup> + (3η|21 - μ<sub>03</sub>)<sup>2</sup></td><td></td></tr><tr><td>φ<sub>4</sub></td><td>(η<sub>30</sub> - η<sub>12</sub>)<sup>2</sup> + (η<sub>21</sub> + μ<sub>03</sub>)<sup>2</sup></td><td></td></tr><tr><td>φ<sub>5</sub></td><td>(η<sub>30</sub> - 3η<sub>12</sub>)(η<sub>30</sub> + η<sub>12</sub>)[(η<sub>30</sub> + η<sub>12</sub>) - 3(η<sub>21</sub> + η<sub>03</sub>)<sup>2</sup>] + (3η<sub>21</sub> - η<sub>03</sub>)(η<sub>21</sub> + η<sub>03</sub>)[3(η<sub>30</sub> + η<sub>12</sub>)<sup>2</sup> - (η<sub>21</sub> + η<sub>03</sub>)<sup>2</sup>]</td><td></td></tr><tr><td>φ<sub>6</sub></td><td>(η<sub>20</sub> - η<sub>02</sub>)[(η<sub>30</sub> + η<sub>12</sub>)<sup>2</sup> - (η<sub>21</sub> + η<sub>03</sub>)<sup>2</sup>] + 4η<sub>11</sub>(η<sub>30</sub> + η<sub>12</sub>)(η<sub>21</sub> - η<sub>03</sub>)</td><td></td></tr><tr><td>φ<sub>7</sub></td><td>(3η<sub>21</sub> - η<sub>03</sub>)(η<sub>30</sub> + η<sub>12</sub>)[(η<sub>30</sub> + η<sub>12</sub>)<sup>2</sup> - 3(η<sub>21</sub> - η<sub>03</sub>)<sup>2</sup>] - (η<sub>30</sub> - 3η<sub>12</sub>)(η<sub>21</sub> + η<sub>03</sub>)[3(η<sub>30</sub> + η<sub>12</sub>)<sup>2</sup> - (η<sub>21</sub> + η<sub>03</sub>)<sup>2</sup>]</td><td></td></tr></tbody></table>

(continued on next page)
<table><caption>**Table 2** Object feature metrics.</caption><thead><tr><th>Feature</th><th>Description</th><th>References</th></tr></thead><tbody><tr><td><strong>Geometric</strong></td><td></td><td></td></tr><tr><td>Perimeter (P)</td><td>Parametric representation of the contour and its points identified by the coordinates x(t) and y(t)</td><td rowspan="6">Costa & Cesar (2009)</td></tr><tr><td>Area (A)</td><td>Integral of the contour</td></tr><tr><td>Area and perimeter ratio</td><td>A/P</td></tr><tr><td>Circularity</td><td>4π A/p²</td></tr><tr><td>Thickness ratio</td><td>p²/A</td></tr><tr><td>Centroid</td><td>Given the center of mass M of a contour of complex signal u(n), the centroid coordinates (z₁, z₂) were obtained by the average of all the points in u(n).</td></tr><tr><td>Centroid to contour maximum distance</td><td>The distance between the centroid and the furthest point on the contour.</td><td></td></tr><tr><td>Centroid to contour minimal distance</td><td>The distance between the centroid and the nearest point on the contour.</td><td></td></tr><tr><td>Centroid to contour average distance</td><td>The average of the distances between the centroid and all points in the contour.</td><td></td></tr><tr><td>Major axis</td><td>Pair of more distant points belonging to the object.</td><td></td></tr><tr><td>Minor axis</td><td>Pair of closest points belonging to the object.</td><td></td></tr><tr><td>Aspect ratio</td><td>Majoraxis/Minoraxis</td><td></td></tr><tr><td>Perimeter and Major axis ratio</td><td>P/Majoraxis</td><td></td></tr><tr><td>Bilateral symmetry</td><td>Bilateral symmetry is given by the proportion of the number of pixels between the intersection of an object and its reflecting shape with respect to the major axis, and the union between those two objects.</td><td>Hu (1962); Huang & Leng (2010)</td></tr><tr><td><strong>Hu's invariant moments<sup>a</sup></strong></td><td></td><td></td></tr><tr><td>φ<sub>1</sub></td><td>η<sub>20</sub> + η<sub>02</sub></td><td></td></tr><tr><td>φ<sub>2</sub></td><td>(η<sub>20</sub> - η<sub>02</sub>)<sup>2</sup> + 4η<sub>11</sub><sup>2</sup></td><td></td></tr><tr><td>φ<sub>3</sub></td><td>(η|30 - 3η<sub>12</sub>)<sup>2</sup> + (3η|21 - μ<sub>03</sub>)<sup>2</sup></td><td></td></tr><tr><td>φ<sub>4</sub></td><td>(η<sub>30</sub> - η<sub>12</sub>)<sup>2</sup> + (η<sub>21</sub> + μ<sub>03</sub>)<sup>2</sup></td><td></td></tr><tr><td>φ<sub>5</sub></td><td>(η<sub>30</sub> - 3η<sub>12</sub>)(η<sub>30</sub> + η<sub>12</sub>)[(η<sub>30</sub> + η<sub>12</sub>) - 3(η<sub>21</sub> + η<sub>03</sub>)<sup>2</sup>] + (3η<sub>21</sub> - η<sub>03</sub>)(η<sub>21</sub> + η<sub>03</sub>)[3(η<sub>30</sub> + η<sub>12</sub>)<sup>2</sup> - (η<sub>21</sub> + η<sub>03</sub>)<sup>2</sup>]</td><td></td></tr><tr><td>φ<sub>6</sub></td><td>(η<sub>20</sub> - η<sub>02</sub>)[(η<sub>30</sub> + η<sub>12</sub>)<sup>2</sup> - (η<sub>21</sub> + η<sub>03</sub>)<sup>2</sup>] + 4η<sub>11</sub>(η<sub>30</sub> + η<sub>12</sub>)(η<sub>21</sub> - η<sub>03</sub>)</td><td></td></tr><tr><td>φ<sub>7</sub></td><td>(3η<sub>21</sub> - η<sub>03</sub>)(η<sub>30</sub> + η<sub>12</sub>)[(η<sub>30</sub> + η<sub>12</sub>)<sup>2</sup> - 3(η<sub>21</sub> - η<sub>03</sub>)<sup>2</sup>] - (η<sub>30</sub> - 3η<sub>12</sub>)(η<sub>21</sub> + η<sub>03</sub>)[3(η<sub>30</sub> + η<sub>12</sub>)<sup>2</sup> - (η<sub>21</sub> + η<sub>03</sub>)<sup>2</sup>]</td><td></td></tr></tbody></table>

(continued on next page)
<table><caption>Table 2 (continued)</caption><thead><tr><th>Feature</th><th>Description</th><th>References</th></tr></thead><tbody><tr><td><strong>Color</strong></td><td></td><td>Burger &amp; Burge (2016)</td></tr><tr><td>Mean</td><td><sup>1</sup>/<sub>MN</sub> &sum;<sub>i=1</sub><sup>M</sup> &sum;<sub>j=1</sub><sup>N</sup> P<sub>ij</sub></td><td></td></tr><tr><td>Median</td><td>P<sub>(n+1)/2</sub></td><td></td></tr><tr><td>Mode</td><td>Pixel value that occurs with greatest frequency</td><td></td></tr><tr><td>Amplitude</td><td>max(p) - min(p)</td><td></td></tr><tr><td>Variance</td><td>&sum;<sub>i=1</sub><sup>M</sup> &sum;<sub>j=1</sub><sup>N</sup> P<sub>ij</sub> &middot; (x<sub>ij</sub> - &mu;)<sup>2</sup></td><td></td></tr><tr><td><strong>Curvature<sup>b</sup></strong></td><td></td><td>Costa &amp; Cesar (2009)</td></tr><tr><td>Bending energy</td><td>B = <sup>1</sup>/<sub>P</sub> &int; k(t)<sup>2</sup> dt</td><td></td></tr><tr><td>Variance</td><td>Var(t) = &sum;<sub>i=0</sub><sup>t</sup> p<sub>i</sub> &middot; (x<sub>i</sub> - &mu;)<sup>2</sup></td><td></td></tr><tr><td>Entropy</td><td>H(t) = - &sum;<sub>i=0</sub><sup>t</sup> p(k(t)) &middot; log(p(k(t)))</td><td></td></tr><tr><td><strong>Color texture<sup>c</sup></strong></td><td></td><td>Gui et al. (2013)</td></tr><tr><td>Entropy (E)</td><td>E = - &sum;<sub>i=1</sub><sup>L</sup> &sum;<sub>j=1</sub><sup>L</sup> p(i,j) log(p(i,j))</td><td></td></tr><tr><td>Angular second moment (ASM)</td><td>ASM = &sum;<sub>i=1</sub><sup>L</sup> &sum;<sub>j=1</sub><sup>L</sup> (p(i,j))<sup>2</sup></td><td></td></tr><tr><td>Contrast (CON)</td><td>CON = &sum;<sub>k=0</sub><sup>L-1</sup> k<sup>2</sup> (&sum;<sub>|i-j|=k</sub> p(i,j))</td><td></td></tr><tr><td>Inverse difference moment (IDM)</td><td>IDM = <sup>&sum;<sub>i=1</sub><sup>L</sup> &sum;<sub>j=1</sub><sup>L</sup> p(i,j)</sup>/<sub>|1+i-j|</sub></td><td></td></tr><tr><td>Correlation (COR)</td><td>COR = <sup>&sum;<sub>i=1</sub><sup>L</sup> &sum;<sub>j=1</sub><sup>L</sup> (i,j) p(i,j) - &mu;<sub>x</sub>&mu;<sub>y</sub></sup>/<sub>&sigma;<sup>2</sup></sub></td><td></td></tr></tbody></table>

Notes.

<footnote><sup>a</sup>Refer to Huang & Leng (2010) for μ and η equations.</footnote>

<footnote><sup>b</sup>The curvature k(t) of a parametric curve c(t) = (x(t), y(t)) was defined as $k(t) = \frac{x'(t)y''(t) - x''(t)y'(t)}{(x'(t))^2 - y'(t)^2}^{1/2}$, where x'(t), y'(t) and x''(t), y''(t) are the first and second derivative of the contour signal x(t) and y(t), respectively.</footnote>

<footnote><sup>c</sup>Texture features were extracted based on the color co-occurrence matrix (CCM).</footnote>
<table><caption>Table 2 (continued)</caption><thead><tr><th>Feature</th><th>Description</th><th>References</th></tr></thead><tbody><tr><td><strong>Color</strong></td><td></td><td>Burger &amp; Burge (2016)</td></tr><tr><td>Mean</td><td><sup>1</sup>/<sub>MN</sub> &sum;<sub>i=1</sub><sup>M</sup> &sum;<sub>j=1</sub><sup>N</sup> P<sub>ij</sub></td><td></td></tr><tr><td>Median</td><td>P<sub>(n+1)/2</sub></td><td></td></tr><tr><td>Mode</td><td>Pixel value that occurs with greatest frequency</td><td></td></tr><tr><td>Amplitude</td><td>max(p) - min(p)</td><td></td></tr><tr><td>Variance</td><td>&sum;<sub>i=1</sub><sup>M</sup> &sum;<sub>j=1</sub><sup>N</sup> P<sub>ij</sub> &middot; (x<sub>ij</sub> - &mu;)<sup>2</sup></td><td></td></tr><tr><td><strong>Curvature<sup>b</sup></strong></td><td></td><td>Costa &amp; Cesar (2009)</td></tr><tr><td>Bending energy</td><td>B = <sup>1</sup>/<sub>P</sub> &int; k(t)<sup>2</sup> dt</td><td></td></tr><tr><td>Variance</td><td>Var(t) = &sum;<sub>i=0</sub><sup>t</sup> p<sub>i</sub> &middot; (x<sub>i</sub> - &mu;)<sup>2</sup></td><td></td></tr><tr><td>Entropy</td><td>H(t) = -&sum;<sub>i=0</sub><sup>t</sup> p(k(t)) &middot; log(p(k(t)))</td><td></td></tr><tr><td><strong>Color texture<sup>c</sup></strong></td><td></td><td>Gui et al. (2013)</td></tr><tr><td>Entropy (E)</td><td>E = -&sum;<sub>i=1</sub><sup>L</sup> &sum;<sub>j=1</sub><sup>L</sup> p(i,j) log(p(i,j))</td><td></td></tr><tr><td>Angular second moment (ASM)</td><td>ASM = &sum;<sub>i=1</sub><sup>L</sup> &sum;<sub>j=1</sub><sup>L</sup> (p(i,j))<sup>2</sup></td><td></td></tr><tr><td>Contrast (CON)</td><td>CON = &sum;<sub>k=0</sub><sup>L-1</sup> k<sup>2</sup> (&sum;<sub>|i-j|=k</sub> p(i,j))</td><td></td></tr><tr><td>Inverse difference moment (IDM)</td><td>IDM = <sup>&sum;<sub>i=1</sub><sup>L</sup> &sum;<sub>j=1</sub><sup>L</sup> p(i,j)</sup>&frasl;<sub>|1+i-j|</sub></td><td></td></tr><tr><td>Correlation (COR)</td><td>COR = <sup>&sum;<sub>i=1</sub><sup>L</sup> &sum;<sub>j=1</sub><sup>L</sup> (i,j)p(i,j) - &mu;<sub>x</sub>&mu;<sub>y</sub></sup>&frasl;<sub>&sigma;<sup>2</sup></sub></td><td></td></tr></tbody></table>

Notes.

<footnote><sup>a</sup>Refer to Huang & Leng (2010) for μ and η equations.</footnote>

<footnote><sup>b</sup>The curvature k(t) of a parametric curve c(t) = (x(t), y(t)) was defined as $k(t) = \frac{x'(t)y''(t) - x''(t)y'(t)}{(x'(t))^2 - y'(t)^2}$. where x'(t), y'(t) and x''(t), y''(t) are the first and second derivative of the contour signal x(t) and y(t), respectively.</footnote>

<footnote><sup>c</sup>Texture features were extracted based on the color co-occurrence matrix (CCM).</footnote>
<table><caption>Table 3 Prediction performance of models on the training and testing sets.</caption><thead><tr><th rowspan="2">Set</th><th rowspan="2">Feature selection</th><th rowspan="2">Model</th><th colspan="6">Metrics (%)</th></tr><tr><th>Sensitivity</th><th>Specificity</th><th>Precision</th><th>Accuracy</th><th>F<sub>1</sub>-score</th><th>AUC</th></tr></thead><tbody><tr><td rowspan="4">Train</td><td rowspan="4">None</td><td>SVM</td><td>67.4</td><td>80.2</td><td>77.7</td><td>73.8</td><td>72.2</td><td>0.797</td></tr><tr><td>KNN</td><td>75.6</td><td>84.2</td><td>83.1</td><td>79.9</td><td>79.2</td><td>0.878</td></tr><tr><td>RF</td><td>99.8</td><td>99.5</td><td>99.5</td><td>99.7</td><td>99.7</td><td>1.0</td></tr><tr><td>Ensemble</td><td>88.6</td><td>92.0</td><td>91.9</td><td>90.3</td><td>90.2</td><td>0.978</td></tr><tr><td rowspan="4">Test</td><td rowspan="4">None</td><td>SVM</td><td>69.7</td><td>75.4</td><td>71.7</td><td>72.7</td><td>70.7</td><td>0.78</td></tr><tr><td>KNN</td><td>69.7</td><td>75.4</td><td>71.7</td><td>72.7</td><td>70.7</td><td>0.759</td></tr><tr><td>RF</td><td>90.5</td><td>88.6</td><td>87.6</td><td>89.5</td><td>89.0</td><td>0.942</td></tr><tr><td>Ensemble</td><td>76.8</td><td>81.8</td><td>79.0</td><td>79.4</td><td>77.9</td><td>0.884</td></tr></tbody></table>

Notes.

AUC, Area under the curve; SVM, Support vector machines; KNN, k-nearest neighbors; RF, random forest.

<figure><img src="image_18.png" /><figcaption>**Figure 5** Model classification performance. (A) Receiver operating characteristic (ROC) curve in the training set. (B) ROC curve in the testing set. AUC, area under the curve.</figcaption></figure>

Full-size DOI: 10.7717/peerj.13470/fig-5

voting classifier is because the ensemble's prediction is based on the average of the prediction probabilities of each classifier. Since the SVM and KNN classifiers presented lower performance individually (both presented sensitivity of 69.7% and specificity of 75.4%, Table 3) the prediction of the ensemble was lower. On the other hand, the RF classifier had a sensitivity of 90.5% and a specificity of 88.6% (Table 5).

The objects found in the images were then marked by the algorithm as parasite. We found two problems in the parasite recognition task. The first is related to the high rate of false positives (Fig. 6A). The regions of the images with leukocytes or high density of red cells showed over-stained areas. These areas were difficult to classify by the algorithm. The second problem was the false negative rate (Fig. 6B). Parasites in regions of the image that presented low contrast or low sharpness (“out of focus”), most commonly found at the
<table><caption>**Table 3** Prediction performance of models on the training and testing sets.</caption><thead><tr><th rowspan="2">Set</th><th rowspan="2">Feature selection</th><th rowspan="2">Model</th><th colspan="6">Metrics (%)</th></tr><tr><th>Sensitivity</th><th>Specificity</th><th>Precision</th><th>Accuracy</th><th>F<sub>1</sub>-score</th><th>AUC</th></tr></thead><tbody><tr><td rowspan="4">Train</td><td rowspan="4">None</td><td>SVM</td><td>67.4</td><td>80.2</td><td>77.7</td><td>73.8</td><td>72.2</td><td>0.797</td></tr><tr><td>KNN</td><td>75.6</td><td>84.2</td><td>83.1</td><td>79.9</td><td>79.2</td><td>0.878</td></tr><tr><td>RF</td><td>99.8</td><td>99.5</td><td>99.5</td><td>99.7</td><td>99.7</td><td>1.0</td></tr><tr><td>Ensemble</td><td>88.6</td><td>92.0</td><td>91.9</td><td>90.3</td><td>90.2</td><td>0.978</td></tr><tr><td rowspan="4">Test</td><td rowspan="4">None</td><td>SVM</td><td>69.7</td><td>75.4</td><td>71.7</td><td>72.7</td><td>70.7</td><td>0.78</td></tr><tr><td>KNN</td><td>69.7</td><td>75.4</td><td>71.7</td><td>72.7</td><td>70.7</td><td>0.759</td></tr><tr><td>RF</td><td>90.5</td><td>88.6</td><td>87.6</td><td>89.5</td><td>89.0</td><td>0.942</td></tr><tr><td>Ensemble</td><td>76.8</td><td>81.8</td><td>79.0</td><td>79.4</td><td>77.9</td><td>0.884</td></tr></tbody></table>

**Notes.**

AUC, Area under the curve; SVM, Support vector machines; KNN, k-nearest neighbors; RF, random forest.

<figure><img src="image_17.png" /><figcaption>**Figure 5** Model classification performance. (A) Receiver operating characteristic (ROC) curve in the training set. (B) ROC curve in the testing set. AUC, area under the curve.</figcaption></figure>

Full-size DOI: 10.7717/peerj.13470/fig-5

voting classifier is because the ensemble's prediction is based on the average of the prediction probabilities of each classifier. Since the SVM and KNN classifiers presented lower performance individually (both presented sensitivity of 69.7% and specificity of 75.4%, Table 3) the prediction of the ensemble was lower. On the other hand, the RF classifier had a sensitivity of 90.5% and a specificity of 88.6% (Table 5).

The objects found in the images were then marked by the algorithm as parasite. We found two problems in the parasite recognition task. The first is related to the high rate of false positives (Fig. 6A). The regions of the images with leukocytes or high density of red cells showed over-stained areas. These areas were difficult to classify by the algorithm. The second problem was the false negative rate (Fig. 6B). Parasites in regions of the image that presented low contrast or low sharpness (“out of focus”), most commonly found at the
<table><caption>**Table 4** Confusion matrix of the voting classifier (ensemble's prediction) in the test set.</caption><thead><tr><th rowspan="2"></th><th rowspan="2"></th><th colspan="2">True label</th></tr><tr><th>Parasite</th><th>Unknown</th></tr></thead><tbody><tr><th rowspan="2">Predicted label</th><th>Parasite</th><td>162</td><td>43</td></tr><tr><th>Unknown</th><td>49</td><td>193</td></tr></tbody></table>

<table><caption>**Table 5** Confusion matrix of the Random Forest classification model in the test set.</caption><thead><tr><th rowspan="2"></th><th rowspan="2"></th><th colspan="2">True label</th></tr><tr><th>Parasite</th><th>Unknown</th></tr></thead><tbody><tr><th rowspan="2">Predicted label</th><th>Parasite</th><td>191</td><td>27</td></tr><tr><th>Unknown</th><td>20</td><td>207</td></tr></tbody></table>

edges of the field of view, were not recognized by the algorithm. This second problem is much more significant, since undiagnosed Chagas disease can put a person's life at risk.

## DISCUSSION AND CONCLUSIONS

In this work, we present an algorithm for automatic detection of the *T. cruzi* parasite in images acquired with a mobile phone device. Our approach involved image segmentation with a graph-based method, extraction of parasite features, selection of the most important features, and classification of these features with a Random Forest model for the detection of the parasite in the image.

The detection of *T. cruzi* in images was previously done using several classification models, such as gaussian discriminant, k-nearest neighbors, AdaBoost + SVM, and convolutional neural networks (CNN) (*Soberanis-Mukul et al., 2013*; *Uc-Cetina, Brito-Loeza & Ruiz-Piña, 2013*; *Uc-Cetina, Brito-Loeza & Ruiz-Piña, 2015*; *Pereira et al., 2020*). Despite these works reported a good performance (sensitivity and specificity >85%), all of them made use of images acquired with a dedicated camera system. Our method obtained a sensitivity of 90.5% and a specificity of 88.6% even though we used images with lower resolution (less than 1 megapixel). Machine learning approaches applied to images obtained from mobile device cameras showed similar performance (sensitivity of 80.5% and specificity of 93.8%) in detecting the malaria agent *Plasmodium* spp. (*Oliveira et al., 2017*; *Rosado et al., 2016*). Therefore, our method was the first to combine machine learning algorithms and low-resolution images to automatically detect *T. cruzi* parasite (Table 6).

We can enhance the classification task by testing other models. Currently, one of the techniques most used in pattern recognition are deep learning approaches (*Acevedo et al., 2019*; *Moen et al., 2019*; *Schmidhuber, 2015*). Deep learning approaches presented better performance in detecting *Plasmodium* spp. (sensitivity of 94.5% and specificity of 96.9%) (*Rajaraman et al., 2018*). However, to build an effective model using these techniques to detect *T. cruzi*, huge data sets are required where the performance of the model increases in logarithmic proportion to the volume of images (*Sun et al., 2017*). Another challenge
<table><caption>**Table 4** Confusion matrix of the voting classifier (ensemble's prediction) in the test set.</caption><thead><tr><th rowspan="2"></th><th rowspan="2"></th><th colspan="2">True label</th></tr><tr><th>Parasite</th><th>Unknown</th></tr></thead><tbody><tr><th rowspan="2">Predicted label</th><th>Parasite</th><td>162</td><td>43</td></tr><tr><th>Unknown</th><td>49</td><td>193</td></tr></tbody></table>

<table><caption>**Table 5** Confusion matrix of the Random Forest classification model in the test set.</caption><thead><tr><th rowspan="2"></th><th rowspan="2"></th><th colspan="2">True label</th></tr><tr><th>Parasite</th><th>Unknown</th></tr></thead><tbody><tr><th rowspan="2">Predicted label</th><th>Parasite</th><td>191</td><td>27</td></tr><tr><th>Unknown</th><td>20</td><td>207</td></tr></tbody></table>

edges of the field of view, were not recognized by the algorithm. This second problem is much more significant, since undiagnosed Chagas disease can put a person's life at risk.

## DISCUSSION AND CONCLUSIONS

In this work, we present an algorithm for automatic detection of the *T. cruzi* parasite in images acquired with a mobile phone device. Our approach involved image segmentation with a graph-based method, extraction of parasite features, selection of the most important features, and classification of these features with a Random Forest model for the detection of the parasite in the image.

The detection of *T. cruzi* in images was previously done using several classification models, such as gaussian discriminant, k-nearest neighbors, AdaBoost + SVM, and convolutional neural networks (CNN) (*Soberanis-Mukul et al., 2013*; *Uc-Cetina, Brito-Loeza & Ruiz-Piña, 2013*; *Uc-Cetina, Brito-Loeza & Ruiz-Piña, 2015*; *Pereira et al., 2020*). Despite these works reported a good performance (sensitivity and specificity >85%), all of them made use of images acquired with a dedicated camera system. Our method obtained a sensitivity of 90.5% and a specificity of 88.6% even though we used images with lower resolution (less than 1 megapixel). Machine learning approaches applied to images obtained from mobile device cameras showed similar performance (sensitivity of 80.5% and specificity of 93.8%) in detecting the malaria agent *Plasmodium* spp. (*Oliveira et al., 2017*; *Rosado et al., 2016*). Therefore, our method was the first to combine machine learning algorithms and low-resolution images to automatically detect *T. cruzi* parasite (Table 6).

We can enhance the classification task by testing other models. Currently, one of the techniques most used in pattern recognition are deep learning approaches (*Acevedo et al., 2019*; *Moen et al., 2019*; *Schmidhuber, 2015*). Deep learning approaches presented better performance in detecting *Plasmodium* spp. (sensitivity of 94.5% and specificity of 96.9%) (*Rajaraman et al., 2018*). However, to build an effective model using these techniques to detect *T. cruzi*, huge data sets are required where the performance of the model increases in logarithmic proportion to the volume of images (*Sun et al., 2017*). Another challenge
A

<figure><img src="image_19.png" /></figure>

B

<figure><img src="image_20.png" /><figcaption>**Figure 6** Sample images of false-positive and false-negative Chagas parasite detection algorithm. (A)</figcaption></figure>

The regions of the images with leukocytes (left) or high density of red cells (right) showed overstained areas that made it difficult to properly classify by the algorithm. (B) Parasites in regions of the image that presented low contrast (left) or low sharpness (right) were not recognized by the algorithm.

Full-size

<figure><img src="image_21.png" /></figure>

DOI: 10.7717/peerj.13470/fig-6

in this kind of study is image acquisition and ground truth annotation. To obtain and annotate this large number of images is particularly difficult, and it is more difficult to apply these techniques in a neglected disease context.

The mobile phone camera mainly affects the outline of objects in the image and the sharpness. The algorithm we present also has its results affected by such characteristics. A higher false-positive rate was observed in regions of the image with leukocytes or high red cell density. A higher false-negative rate was observed in regions of low contrast and sharpness. Therefore, smear quality directly affects classifier performance. It is extremely important to reduce the false-negative rate, since undiagnosed patients can be left without proper treatment and in a life-threatening situation. We recommend evaluating the algorithm on images acquired from samples with different staining time and
A

<figure><img src="image_18.png" /></figure>

B

<figure><img src="image_19.png" /><figcaption>**Figure 6** Sample images of false-positive and false-negative Chagas parasite detection algorithm. (A)</figcaption></figure>

The regions of the images with leukocytes (left) or high density of red cells (right) showed overstained areas that made it difficult to properly classify by the algorithm. (B) Parasites in regions of the image that presented low contrast (left) or low sharpness (right) were not recognized by the algorithm.

Full-size

<figure><img src="image_20.png" /></figure>

DOI: 10.7717/peerj.13470/fig-6

in this kind of study is image acquisition and ground truth annotation. To obtain and annotate this large number of images is particularly difficult, and it is more difficult to apply these techniques in a neglected disease context.

The mobile phone camera mainly affects the outline of objects in the image and the sharpness. The algorithm we present also has its results affected by such characteristics. A higher false-positive rate was observed in regions of the image with leukocytes or high red cell density. A higher false-negative rate was observed in regions of low contrast and sharpness. Therefore, smear quality directly affects classifier performance. It is extremely important to reduce the false-negative rate, since undiagnosed patients can be left without proper treatment and in a life-threatening situation. We recommend evaluating the algorithm on images acquired from samples with different staining time and
<table><caption>**Table 6** Comparison of the results of our algorithm with other published studies.</caption><thead><tr><th>Reference</th><th>Parasite</th><th>Image capture device</th><th>ML Model</th><th>Sensitivity (%)</th><th>Specifity (%)</th></tr></thead><tbody><tr><td>The present work</td><td>T. cruzi</td><td>Mobile phone camera</td><td>RF</td><td>90.5</td><td>88.6</td></tr><tr><td>Uc-Cetina, Brito-Loeza & Ruiz-Piña (2013)</td><td>T. cruzi</td><td>Dedicated camera</td><td>Gaussian discriminant</td><td>98.3</td><td>84.4</td></tr><tr><td>Uc-Cetina, Brito-Loeza & Ruiz-Piña (2015)</td><td>T. cruzi</td><td>Dedicated camera</td><td>AdaBoost + SVM</td><td>100</td><td>93.2</td></tr><tr><td>Pereira et al. (2020)</td><td>T. cruzi</td><td>Dedicated camera</td><td>Convolutional Neural Network</td><td>97.6</td><td>95.2</td></tr><tr><td>Soberanis-Mukul et al. (2013)</td><td>T. cruzi</td><td>Dedicated camera</td><td>KNN</td><td>98</td><td>85</td></tr><tr><td>Savkare & Narote (2012)</td><td>Plasmodium spp.</td><td>Dedicated camera</td><td>SVM</td><td>96.3</td><td>99.1</td></tr><tr><td>Yang et al. (2020)</td><td>Plasmodium spp.</td><td>Mobile phone camera</td><td>Convolutional Neural Network</td><td>92.6</td><td>94.3</td></tr><tr><td>Rajaraman et al. (2018)</td><td>Plasmodium spp.</td><td>Mobile phone camera</td><td>Convolutional Neural Network</td><td>94.5</td><td>96.9</td></tr><tr><td>Rosado et al. (2016)</td><td>Plasmodium spp.</td><td>Mobile phone camera</td><td>SVM</td><td>80.5</td><td>93.8</td></tr><tr><td>Oliveira et al. (2017)</td><td>Plasmodium spp.</td><td>Mobile phone camera</td><td>Adaboost</td><td>59</td><td>95</td></tr></tbody></table>

**Notes.**

RF, random forest; SVM, support vector machine; KNN, k-nearest neighbours.
<table><caption>**Table 6** Comparison of the results of our algorithm with other published studies.</caption><thead><tr><th>Reference</th><th>Parasite</th><th>Image capture device</th><th>ML Model</th><th>Sensitivity (%)</th><th>Specifity (%)</th></tr></thead><tbody><tr><td>The present work</td><td>T. cruzi</td><td>Mobile phone camera</td><td>RF</td><td>90.5</td><td>88.6</td></tr><tr><td>Uc-Cetina, Brito-Loeza & Ruiz-Piña (2013)</td><td>T. cruzi</td><td>Dedicated camera</td><td>Gaussian discriminant</td><td>98.3</td><td>84.4</td></tr><tr><td>Uc-Cetina, Brito-Loeza & Ruiz-Piña (2015)</td><td>T. cruzi</td><td>Dedicated camera</td><td>AdaBoost + SVM</td><td>100</td><td>93.2</td></tr><tr><td>Pereira et al. (2020)</td><td>T. cruzi</td><td>Dedicated camera</td><td>Convolutional Neural Network</td><td>97.6</td><td>95.2</td></tr><tr><td>Soberanis-Mukul et al. (2013)</td><td>T. cruzi</td><td>Dedicated camera</td><td>KNN</td><td>98</td><td>85</td></tr><tr><td>Savkare & Narote (2012)</td><td>Plasmodium spp.</td><td>Dedicated camera</td><td>SVM</td><td>96.3</td><td>99.1</td></tr><tr><td>Yang et al. (2020)</td><td>Plasmodium spp.</td><td>Mobile phone camera</td><td>Convolutional Neural Network</td><td>92.6</td><td>94.3</td></tr><tr><td>Rajaraman et al. (2018)</td><td>Plasmodium spp.</td><td>Mobile phone camera</td><td>Convolutional Neural Network</td><td>94.5</td><td>96.9</td></tr><tr><td>Rosado et al. (2016)</td><td>Plasmodium spp.</td><td>Mobile phone camera</td><td>SVM</td><td>80.5</td><td>93.8</td></tr><tr><td>Oliveira et al. (2017)</td><td>Plasmodium spp.</td><td>Mobile phone camera</td><td>Adaboost</td><td>59</td><td>95</td></tr></tbody></table>

**Notes.**

RF, random forest; SVM, support vector machine; KNN, k-nearest neighbours.
dye concentration. Such an assessment can further validate the robustness of the algorithm and identify optimal sample preparation.

In summary, our results demonstrate that the proposed algorithm can detect trypomastigote forms of *T. cruzi* in images acquired with a mobile device attached to a microscope. Automating image analysis acquired with a mobile device is a viable alternative for reducing costs and gaining efficiency in the use of the optical microscope. We hope that this algorithm can serve as a tool for early diagnosis of Chagas disease.

# ADDITIONAL INFORMATION AND DECLARATIONS

## Funding

This work was supported by the São Paulo Research Foundation (FAPESP, grant numbers 2020/12017-9 to Mauro César Cafundó Morais, 2018/14933-2 to Helder Nakaya, 2015/22308 to Luciano da F. Costa) and National Council for Research (CNPq grant n. 307085/2018-0). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

## Grant Disclosures

The following grant information was disclosed by the authors:
São Paulo Research Foundation (FAPESP): 2020/12017-9, 2018/14933-2, 2015/22308.
National Council for Research (CNPq): 307085/2018-0.

## Competing Interests

Helder I. Nakaya is an Academic Editor for PeerJ.

## Author Contributions

* Mauro César Cafundó Morais conceived and designed the experiments, performed the experiments, analyzed the data, prepared figures and/or tables, authored or reviewed drafts of the paper, and approved the final draft.

* Diogo Silva conceived and designed the experiments, performed the experiments, analyzed the data, prepared figures and/or tables, authored or reviewed drafts of the paper, and approved the final draft.

* Matheus Marques Milagre conceived and designed the experiments, performed the experiments, analyzed the data, authored or reviewed drafts of the paper, and approved the final draft.

* Maykon Tavares de Oliveira conceived and designed the experiments, performed the experiments, authored or reviewed drafts of the paper, and approved the final draft.

* Thaís Pereira performed the experiments, authored or reviewed drafts of the paper, and approved the final draft.

* João Santana Silva conceived and designed the experiments, analyzed the data, authored or reviewed drafts of the paper, and approved the final draft.

* Luciano da F. Costa conceived and designed the experiments, analyzed the data, authored or reviewed drafts of the paper, and approved the final draft.
dye concentration. Such an assessment can further validate the robustness of the algorithm and identify optimal sample preparation.

In summary, our results demonstrate that the proposed algorithm can detect trypomastigote forms of *T. cruzi* in images acquired with a mobile device attached to a microscope. Automating image analysis acquired with a mobile device is a viable alternative for reducing costs and gaining efficiency in the use of the optical microscope. We hope that this algorithm can serve as a tool for early diagnosis of Chagas disease.

# ADDITIONAL INFORMATION AND DECLARATIONS

## Funding

This work was supported by the São Paulo Research Foundation (FAPESP, grant numbers 2020/12017-9 to Mauro César Cafundó Morais, 2018/14933-2 to Helder Nakaya, 2015/22308 to Luciano da F. Costa) and National Council for Research (CNPq grant n. 307085/2018-0). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

## Grant Disclosures

The following grant information was disclosed by the authors:
São Paulo Research Foundation (FAPESP): 2020/12017-9, 2018/14933-2, 2015/22308.
National Council for Research (CNPq): 307085/2018-0.

## Competing Interests

Helder I. Nakaya is an Academic Editor for PeerJ.

## Author Contributions

* Mauro César Cafundó Morais conceived and designed the experiments, performed the experiments, analyzed the data, prepared figures and/or tables, authored or reviewed drafts of the paper, and approved the final draft.

* Diogo Silva conceived and designed the experiments, performed the experiments, analyzed the data, prepared figures and/or tables, authored or reviewed drafts of the paper, and approved the final draft.

* Matheus Marques Milagre conceived and designed the experiments, performed the experiments, analyzed the data, authored or reviewed drafts of the paper, and approved the final draft.

* Maykon Tavares de Oliveira conceived and designed the experiments, performed the experiments, authored or reviewed drafts of the paper, and approved the final draft.

* Thaís Pereira performed the experiments, authored or reviewed drafts of the paper, and approved the final draft.

* João Santana Silva conceived and designed the experiments, analyzed the data, authored or reviewed drafts of the paper, and approved the final draft.

* Luciano da F. Costa conceived and designed the experiments, analyzed the data, authored or reviewed drafts of the paper, and approved the final draft.
* Paola Minoprio conceived and designed the experiments, analyzed the data, authored or reviewed drafts of the paper, and approved the final draft.

* Roberto Marcondes Cesar Junior conceived and designed the experiments, analyzed the data, authored or reviewed drafts of the paper, and approved the final draft.

* Ricardo Gazzinelli conceived and designed the experiments, authored or reviewed drafts of the paper, and approved the final draft.

* Marta de Lana conceived and designed the experiments, analyzed the data, authored or reviewed drafts of the paper, and approved the final draft.

* Helder I. Nakaya conceived and designed the experiments, analyzed the data, authored or reviewed drafts of the paper, and approved the final draft.

## Animal Ethics

The following information was supplied relating to ethical approvals (i.e., approving body and any reference numbers):

The Ethics Committee for the Use of Animals (CEUA) at the Federal University of Ouro Preto, Laboratory of Chagas Disease, Minas Gerais, Brazil provided approval for this research (CEUA No. 2015/50).

## Data Availability

The following information was supplied regarding data availability:

The image data files are available at Zenodo: https://zenodo.org/record/5123062.

The image annotations and features extracted from the parasite's objects are available in the Supplementary Files.

The code used for model development are available at Bitbucket: https://bitbucket.org/dmatos88/jmire2/src/master/.

The code to test the algorithm is available at GitHub: https://github.com/csbl-br/chagas_detection/.

## Supplemental Information

Supplemental information for this article can be found online at http://dx.doi.org/10.7717/peerj.13470#supplemental-information.

## REFERENCES

**Acevedo A, Alférez S, Merino A, Puigví L, Rodellar J. 2019.** Recognition of peripheral blood cell images using convolutional neural networks. *Computer Methods and Programs in Biomedicine* **180**:105020 DOI 10.1016/j.cmpb.2019.105020.

**Ben-Hur A, Weston J. 2010.** A user's guide to support vector machines. In: Carugo O, Eisenhaber F, eds. *Data mining techniques for the life sciences*. Totowa, NJ, USA: Humana Press, 223–239.

**Breiman L. 2001.** Random forests. *Machine Learning* **45**:5–32
DOI 10.1023/A:1010933404324.

**Burger W, Burge MJ. 2016.** *Digital image processing: an algorithmic introduction using Java*. London: Springer London.
* Paola Minoprio conceived and designed the experiments, analyzed the data, authored or reviewed drafts of the paper, and approved the final draft.

* Roberto Marcondes Cesar Junior conceived and designed the experiments, analyzed the data, authored or reviewed drafts of the paper, and approved the final draft.

* Ricardo Gazzinelli conceived and designed the experiments, authored or reviewed drafts of the paper, and approved the final draft.

* Marta de Lana conceived and designed the experiments, analyzed the data, authored or reviewed drafts of the paper, and approved the final draft.

* Helder I. Nakaya conceived and designed the experiments, analyzed the data, authored or reviewed drafts of the paper, and approved the final draft.

## Animal Ethics

The following information was supplied relating to ethical approvals (i.e., approving body and any reference numbers):

The Ethics Committee for the Use of Animals (CEUA) at the Federal University of Ouro Preto, Laboratory of Chagas Disease, Minas Gerais, Brazil provided approval for this research (CEUA No. 2015/50).

## Data Availability

The following information was supplied regarding data availability:

The image data files are available at Zenodo: https://zenodo.org/record/5123062.

The image annotations and features extracted from the parasite's objects are available in the Supplementary Files.

The code used for model development are available at Bitbucket: https://bitbucket.org/dmatos88/jmire2/src/master/.

The code to test the algorithm is available at GitHub: https://github.com/csbl-br/chagas_detection/.

## Supplemental Information

Supplemental information for this article can be found online at http://dx.doi.org/10.7717/peerj.13470#supplemental-information.

## REFERENCES

**Acevedo A, Alférez S, Merino A, Puigví L, Rodellar J. 2019.** Recognition of peripheral blood cell images using convolutional neural networks. *Computer Methods and Programs in Biomedicine* **180**:105020 DOI 10.1016/j.cmpb.2019.105020.

**Ben-Hur A, Weston J. 2010.** A user's guide to support vector machines. In: Carugo O, Eisenhaber F, eds. *Data mining techniques for the life sciences*. Totowa, NJ, USA: Humana Press, 223–239.

**Breiman L. 2001.** Random forests. *Machine Learning* **45**:5–32
DOI 10.1023/A:1010933404324.

**Burger W, Burge MJ. 2016.** *Digital image processing: an algorithmic introduction using Java*. London: Springer London.
Cancino-Faure B, Fisa R, Riera C, Bula I, Girona-Llobera E, Jimenez-Marco T. 2015. Evidence of meaningful levels of Trypanosoma cruzi in platelet concentrates from seropositive blood donors. *Transfusion* 55:1249–1255 DOI 10.1111/trf.12989.

Chagas C. 1909. Nova tripanozomiaze humana: estudos sobre a morfolojia e o ci-
clo evolutivo do Schizotrypanum cruzi n. gen. n. sp. ajente etiolojico de nova
entidade morbida do homem. *Memórias Do Instituto Oswaldo Cruz* 1:159-218
DOI 10.1590/S0074-02761909000200008.

Chen RC, Dewi C, Huang SW, Caraka RE. 2020. Selecting critical features for data classification based on machine learning methods. *Journal of Big Data* 7:52 DOI 10.1186/s40537-020-00327-4.

Cunningham P, Delany SJ. 2007. K-Nearest neighbour classifiers. *Multiple Classifier Systems* 34:1–17.

Dias JCP, Ramos Jr AN, Gontijo ED, Luquetti A, Shikanai-Yasuda MA, Coura JR, Torres RM, Melo JR da C, de Almeida EA, de Oliveira Jr W, Silveira AC, de Rezende JM, Pinto FS, Ferreira AW, Rassi A, Filho AAF, de Sousa AS, Correia D, Jansen AM, Andrade GMQ, Britto CFP de C, Pinto AYN, Rassi Jr A, Campos DE, Abad-Franch F, Santos SE, Chiari E, Hasslocher-Moreno AM, Moreira EF, Marques DSO, Silva EL, Marin-Neto JA, Galvão LM da C, Xavier SS, Valente SA da S, Carvalho NB, Cardoso AV, Silva RA, da Costa VM, Vivaldini SM, Oliveira SM, Valente V da C, Lima MM, Alves RV. 2016. 2nd Brazilian consensus on Chagas disease, 2015. *Revista da Sociedade Brasileira de Medicina Tropical* 49(Suppl 1):3–60 DOI 10.1590/0037-8682-0505-2016.

Costa da FL, Cesar Jr RM. 2009. *Shape classification and analysis: theory and practice.*
Second edition. Boca Raton, FL, USA: CRC Press.

Felzenszwalb PF, Huttenlocher DP. 2004. Efficient graph-based image segmentation.
*International Journal of Computer Vision* 59:167–181
DOI 10.1023/B:VISI.0000022288.19776.77.

Filigheddu MT, Górgolas M, Ramos JM. 2017. Enfermedad de Chagas de transmisión oral. *Medical Clínica* 148:125-131.

Gomes YM, Lorena VM, Luquetti AO. 2009. Diagnosis of Chagas disease: what has been achieved? What remains to be done with regard to diagnosis and follow up studies? *Memórias do Instituto Oswaldo Cruz* 104:115-121 DOI 10.1590/S0074-02762009000900017.

Górriz M, Aparico A, Raventós B, Vilaplana V, Sayrol E, López-Codina D. 2018. Leishmaniasis parasite segmentation and classification using deep learning. *ArXiv preprint*. arXiv:1812.11586 [cs] 10945.

Gui W, Liu J, Yang C, Chen N, Liao X. 2013. Color co-occurrence matrix based froth image texture extraction for mineral flotation. *Minerals Engineering* 46-47:60-67.

Hsu C-W, Chang C-C, Lin C-J. 2016. A practical guide to support vector classification 16.

Hu M-K. 1962. Visual pattern recognition by moment invariants. *Information Theory, IRE Transactions* 8:179–187.
Cancino-Faure B, Fisa R, Riera C, Bula I, Girona-Llobera E, Jimenez-Marco T. 2015. Evidence of meaningful levels of Trypanosoma cruzi in platelet concentrates from seropositive blood donors. *Transfusion* **55**:1249–1255 DOI **10.1111/trf.12989**.

Chagas C. 1909. Nova tripanozomiaze humana: estudos sobre a morfolojia e o ci-
clo evolutivo do Schizotrypanum cruzi n. gen. n. sp. ajente etiolojico de nova
entidade morbida do homem. *Memórias Do Instituto Oswaldo Cruz* **1**:159-218
DOI **10.1590/S0074-02761909000200008**.

Chen RC, Dewi C, Huang SW, Caraka RE. 2020. Selecting critical features for data classification based on machine learning methods. *Journal of Big Data* **7**:52 DOI **10.1186/s40537-020-00327-4**.

Cunningham P, Delany SJ. 2007. K-Nearest neighbour classifiers. *Multiple Classifier Systems* **34**:1–17.

Dias JCP, Ramos Jr AN, Gontijo ED, Luquetti A, Shikanai-Yasuda MA, Coura JR, Torres RM, Melo JR da C, de Almeida EA, de Oliveira Jr W, Silveira AC, de Rezende JM, Pinto FS, Ferreira AW, Rassi A, Filho AAF, de Sousa AS, Correia D, Jansen AM, Andrade GMQ, Britto CFP de C, Pinto AYN, Rassi Jr A, Campos DE, Abad-Franch F, Santos SE, Chiari E, Hasslocher-Moreno AM, Moreira EF, Marques DSO, Silva EL, Marin-Neto JA, Galvão LM da C, Xavier SS, Valente SA da S, Carvalho NB, Cardoso AV, Silva RA, da Costa VM, Vivaldini SM, Oliveira SM, Valente V da C, Lima MM, Alves RV. 2016. 2nd Brazilian consensus on Chagas disease, 2015. *Revista da Sociedade Brasileira de Medicina Tropical* **49**(Suppl 1):3–60 DOI **10.1590/0037-8682-0505-2016**.

Costa da FL, Cesar Jr RM. 2009. *Shape classification and analysis: theory and practice.*
Second edition. Boca Raton, FL, USA: CRC Press.

Felzenszwalb PF, Huttenlocher DP. 2004. Efficient graph-based image segmentation.
*International Journal of Computer Vision* **59**:167–181
DOI **10.1023/B:VISI.0000022288.19776.77**.

Filigheddu MT, Górgolas M, Ramos JM. 2017. Enfermedad de Chagas de transmisión oral. *Medical Clínica* **148**:125-131.

Gomes YM, Lorena VM, Luquetti AO. 2009. Diagnosis of Chagas disease: what has been achieved? What remains to be done with regard to diagnosis and follow up studies? *Memórias do Instituto Oswaldo Cruz* **104**:115-121 DOI **10.1590/S0074-02762009000900017**.

Górriz M, Aparico A, Raventós B, Vilaplana V, Sayrol E, López-Codina D. 2018. Leishmaniasis parasite segmentation and classification using deep learning. ArXiv preprint. arXiv:1812.11586 [cs] 10945.

Gui W, Liu J, Yang C, Chen N, Liao X. 2013. Color co-occurrence matrix based froth image texture extraction for mineral flotation. *Minerals Engineering* **46-47**:60–67.

Hsu C-W, Chang C-C, Lin C-J. 2016. A practical guide to support vector classification 16.

Hu M-K. 1962. Visual pattern recognition by moment invariants. *Information Theory, IRE Transactions* **8**:179–187.
**Huang Z, Leng J. 2010.** Analysis of Hu's moment invariants on image scaling and rotation. In: *2010 2nd international conference on computer engineering and technology*, V7-476-V7-480.

**Lana M, Machado EMM. 2017.** Biology of *Trypanosoma cruzi* and biological diversity. In: *American trypanosomiasis Chagas disease: one hundred years of research*. 2nd edition. London, England: Elsevier.

**Luquetti AO, Schmuñis GA. 2017.** 29 - Diagnosis of *Trypanosoma cruzi* infection. In: *Teleria J, Tibayrenc M, eds. American Trypanosomiasis chagas disease*. Second edition. London: Elsevier, 687-730 DOI 10.1016/B978-0-12-801029-7.00030-7.

**Luquetti AO, Tavares SB do N, Siriano L da R, Oliveira RA de, Campos DE, Morais CA de, Oliveira EC de. 2015.** Congenital transmission of Trypanosoma cruzi in central Brazil. A study of 1, 211 individuals born to infected mothers. *Memórias do Instituto Oswaldo Cruz* 110:369-376 DOI 10.1590/0074-02760140410.

**Moen E, Bannon D, Kudo T, Graf W, Covert M, Valen DV. 2019.** Deep learning for cellular image analysis. *Nature Methods* 16:1233-1246 DOI 10.1038/s41592-019-0403-1.

**Oliveira AD, Prats C, Espasa M, Zarzuela Serrat F, Montañola Sales C, Silgado A, Codina DL, Arruda ME, Prat JGI, Albuquerque J. 2017.** The malaria system microapp: a new, mobile device-based tool for malaria diagnosis. *JMIR Research Protocols* 6:e70 DOI 10.2196/resprot.6758.

**Palm C. 2004.** Color texture classification by integrative co-occurrence matrices. *Pattern Recognition* 37:965-976 DOI 10.1016/j.patcog.2003.09.010.

**Pedregosa F, Varouquaux G, Gramfort A, Michel V, Thirion B. 2011.** Scikit-learn: machine learning in python. *Journal of Machine Learning Research* 12:2825-2830.

**Pereira AS, Pyrrho AS, Vanzan DF, Mazza LO, Gomes JGRC. 2020.** Deep Convolutional Neural Network applied to Chagas Disease Parasitemia Assessment, in Anais do 14. In: *Congresso Brasileiro de Inteligência Computacional 1-8 (ABRICOM, 2020)*. doi:10.21528/CBIC, volume 21. 9-119.

**Rajaraman S, Antani SK, Poostchi M, Silamut K, Hossain MA, Maude RJ, Jaeger S, Thoma GR. 2018.** Pre-trained convolutional neural networks as feature extractors toward improved malaria parasite detection in thin blood smear images. *PeerJ* 6:e4568 DOI 10.7717/peerj.4568.

**Rosado L, da Costa JMC, Elias D, Cardoso JS. 2016.** Automated detection of malaria parasites on thick blood smears via mobile devices. *Procedia Computer Science* 90:138-144 DOI 10.1016/j.procs.2016.07.024.

**Savkare SS, Narote SP. 2012.** Automatic system for classification of erythrocytes infected with malaria and identification of parasite's life stage. *Procedia Technology* 6:405-410 DOI 10.1016/j.protcy.2012.10.048.

**Schmidhuber J. 2015.** Deep learning in neural networks: an overview. *Neural Networks* 61:85-117 DOI 10.1016/j.neunet.2014.09.003.

**Soberanis-Mukul R, Uc-Cetina V, Brito-Loeza C, Ruiz-Piña H. 2013.** An automatic algorithm for the detection of Trypanosoma cruzi parasites in blood
**Huang Z, Leng J. 2010.** Analysis of Hu's moment invariants on image scaling and rotation. In: *2010 2nd international conference on computer engineering and technology*, V7-476-V7-480.

**Lana M, Machado EMM. 2017.** Biology of *Trypanosoma cruzi* and biological diversity. In: *American trypanosomiasis Chagas disease: one hundred years of research*. 2nd edition. London, England: Elsevier.

**Luquetti AO, Schmuñis GA. 2017.** 29 - Diagnosis of *Trypanosoma cruzi* infection. In: *Teleria J, Tibayrenc M, eds. American Trypanosomiasis chagas disease*. Second edition. London: Elsevier, 687-730 DOI 10.1016/B978-0-12-801029-7.00030-7.

**Luquetti AO, Tavares SB do N, Siriano L da R, Oliveira RA de, Campos DE, Morais CA de, Oliveira EC de. 2015.** Congenital transmission of Trypanosoma cruzi in central Brazil. A study of 1, 211 individuals born to infected mothers. *Memórias do Instituto Oswaldo Cruz* 110:369-376 DOI 10.1590/0074-02760140410.

**Moen E, Bannon D, Kudo T, Graf W, Covert M, Valen DV. 2019.** Deep learning for cellular image analysis. *Nature Methods* 16:1233-1246 DOI 10.1038/s41592-019-0403-1.

**Oliveira AD, Prats C, Espasa M, Zarzuela Serrat F, Montañola Sales C, Silgado A, Codina DL, Arruda ME, Prat JGI, Albuquerque J. 2017.** The malaria system microapp: a new, mobile device-based tool for malaria diagnosis. *JMIR Research Protocols* 6:e70 DOI 10.2196/resprot.6758.

**Palm C. 2004.** Color texture classification by integrative co-occurrence matrices. *Pattern Recognition* 37:965-976 DOI 10.1016/j.patcog.2003.09.010.

**Pedregosa F, Varouquaux G, Gramfort A, Michel V, Thirion B. 2011.** Scikit-learn: machine learning in python. *Journal of Machine Learning Research* 12:2825-2830.

**Pereira AS, Pyrrho AS, Vanzan DF, Mazza LO, Gomes JGRC. 2020.** Deep Convolutional Neural Network applied to Chagas Disease Parasitemia Assessment, in Anais do 14. In: *Congresso Brasileiro de Inteligência Computacional 1-8 (ABRICOM, 2020)*. doi:10.21528/CBIC, volume 21. 9-119.

**Rajaraman S, Antani SK, Poostchi M, Silamut K, Hossain MA, Maude RJ, Jaeger S, Thoma GR. 2018.** Pre-trained convolutional neural networks as feature extractors toward improved malaria parasite detection in thin blood smear images. *PeerJ* 6:e4568 DOI 10.7717/peerj.4568.

**Rosado L, da Costa JMC, Elias D, Cardoso JS. 2016.** Automated detection of malaria parasites on thick blood smears via mobile devices. *Procedia Computer Science* 90:138-144 DOI 10.1016/j.procs.2016.07.024.

**Savkare SS, Narote SP. 2012.** Automatic system for classification of erythrocytes infected with malaria and identification of parasite's life stage. *Procedia Technology* 6:405-410 DOI 10.1016/j.protcy.2012.10.048.

**Schmidhuber J. 2015.** Deep learning in neural networks: an overview. *Neural Networks* 61:85-117 DOI 10.1016/j.neunet.2014.09.003.

**Soberanis-Mukul R, Uc-Cetina V, Brito-Loeza C, Ruiz-Piña H. 2013.** An automatic algorithm for the detection of Trypanosoma cruzi parasites in blood
sample images. *Computer Methods and Programs in Biomedicine* **112**:633–639
DOI 10.1016/j.cmpb.2013.07.013.

**Sun C, Shrivastava A, Singh S, Gupta A. 2017.** Revisiting unreasonable effectiveness of data in deep learning era. ArXiv preprint. arXiv:170702968 Cs.

**Uc-Cetina V, Brito-Loeza C, Ruiz-Piña H. 2013.** Chagas parasites detection through Gaussian discriminant analysis. *Abstraction & Application* **8**:6–17.

**Uc-Cetina V, Brito-Loeza C, Ruiz-Piña H. 2015.** Chagas parasite detection in blood images using adaboost. *Computational and Mathematical Methods* **2015**:1–13.

**Vallada EP. 1999.** *Manual de Técnicas Hematológicas*. São Paulo: Atheneu.

**Yang F, Poostchi M, Yu H, Zhou Z, Silamut K, Yu J, Maude RJ, Jaeger S, Antani S. 2020.** Deep learning for smartphone-based malaria parasite detection in thick blood smears. *IEEE Journal of Biomedical and Health Informatics* **24**:1427–1438
DOI 10.1109/JBHI.2019.2939121.
sample images. *Computer Methods and Programs in Biomedicine* **112**:633–639
DOI 10.1016/j.cmpb.2013.07.013.

**Sun C, Shrivastava A, Singh S, Gupta A. 2017.** Revisiting unreasonable effectiveness of data in deep learning era. ArXiv preprint. arXiv:170702968 Cs.

**Uc-Cetina V, Brito-Loeza C, Ruiz-Piña H. 2013.** Chagas parasites detection through Gaussian discriminant analysis. *Abstraction & Application* **8**:6–17.

**Uc-Cetina V, Brito-Loeza C, Ruiz-Piña H. 2015.** Chagas parasite detection in blood images using adaboost. *Computational and Mathematical Methods* **2015**:1–13.

**Vallada EP. 1999.** *Manual de Técnicas Hematológicas*. São Paulo: Atheneu.

**Yang F, Poostchi M, Yu H, Zhou Z, Silamut K, Yu J, Maude RJ, Jaeger S, Antani S. 2020.** Deep learning for smartphone-based malaria parasite detection in thick blood smears. *IEEE Journal of Biomedical and Health Informatics* **24**:1427–1438
DOI 10.1109/JBHI.2019.2939121.