<!DOCTYPE article PUBLIC "-//NLM//DTD JATS (Z39.96) Journal Archiving and Interchange DTD v1.0 20120330//EN" "JATS-archivearticle1.dtd">
<article xmlns:xlink="http://www.w3.org/1999/xlink">
  <front>
    <journal-meta />
    <article-meta>
      <title-group>
        <article-title>Supervised learning and image processing for e cient malaria detection</article-title>
      </title-group>
      <contrib-group>
        <aff id="aff0">
          <label>0</label>
          <institution>University of Cape Town</institution>
          ,
          <addr-line>Cape Town</addr-line>
          ,
          <country country="ZA">South Africa</country>
        </aff>
      </contrib-group>
      <abstract>
        <p>Malaria is a devastating disease that leads to many deaths each year. Currently, most malaria diagnoses are performed manually, which is time consuming. This may result in it taking longer to diagnose patients, especially those in poor and rural areas, motivating the development of automated detection tools. Deep learning approaches, particularly convolutional neural networks (CNNs), have seen success in the existing literature. However, CNNs are computationally expensive and require signi cant amounts of training data, which may limit their real-world viability, especially in poorer and rural communities. Nondeep supervised techniques are largely free from these limitations but have received less attention in the existing literature. This paper differs from existing work using non-deep systems by investigating the use of RFs, adopting a more rigourous testing methodology and conducting a broader exploration of pre-processing techniques. Two non-deep supervised systems are proposed, based on random forests (RFs) and support vector machines (SVMs). The RF system performs better, having achieved an accuracy of 96.29% when tested on 20 000 images, with runtimes of less than two seconds. Testing on a small dataset of images gathered from a di erent source achieves similar performance, suggesting the model may generalise to di erent imaging conditions. The system achieves higher recall than existing non-deep approaches, and its accuracy, recall and precision are within 4% of the highest performing CNN approach.</p>
      </abstract>
      <kwd-group>
        <kwd>Supervised Learning Image Processing Computer-Aided Diagnosis Machine Learning Image Classi cation Feature Extraction Malaria</kwd>
      </kwd-group>
    </article-meta>
  </front>
  <body>
    <sec id="sec-1">
      <title>-</title>
      <p>
        Malaria is a parasitic disease that can have devastating e ects, not only for
individuals who contract it, but also for their families and communities, which
may su er economic harm as a result [
        <xref ref-type="bibr" rid="ref14">14</xref>
        ]. African countries are
disproportionately a ected, with 92% of global infections and 93% of global deaths falling in
the World Health Organisation African region. Malaria also disproportionately
a ects poorer and rural areas where access to diagnosis and treatment is limited.
      </p>
      <p>
        This presents a clear problem as those who are most in need of medical
assistance are less likely to receive it in time. Currently, the gold standard for
malaria diagnosis is manual microscopy performed by an expert. While this is a
reliable method of diagnosis, it is also costly and time consuming [
        <xref ref-type="bibr" rid="ref15">15</xref>
        ]. Part of the
reason poorer and rural areas are worst a ected by malaria may be attributed
to a lack of access to these experts. This motivates the need for low-cost and
reliable automated detection systems that minimise the time burden on medical
experts.
      </p>
      <p>
        Deep learning approaches, speci cally convolutional neural networks (CNNs),
have proven e ective when applied to the problem of malaria diagnosis, with
one example achieving over 99% accuracy [
        <xref ref-type="bibr" rid="ref12">12</xref>
        ]. However, CNNs are reliant on
large sets of training data and signi cant computational power for training.
Since imaging conditions vary across di erent clinics and laboratories, additional
training on local blood sample data may be necessary. As such, limited amounts
of labelled data may hamper widespread adoption. Moreover, the computational
resources required to run CNN models may limit their application in poorer and
rural areas.
      </p>
      <p>Non-deep machine learning approaches have not been thoroughly
investigated in the existing literature. However, these approaches do not su er from
the limitations noted above, and may therefore prove to be more suitable for
use in environments with limited computing resources. This paper presents two
non-deep supervised learning systems, using random forests (RFs) and support
vector machines (SVMs). This paper di ers from existing work using non-deep
models by adopting more rigourous testing on a larger set of data,
investigating RFs as a potential solution and exploring more image ltering and feature
extraction approaches.</p>
      <p>
        Through an initial evaluation process, feature extraction and image ltering
approaches were selected to combine with the RF and SVM models. It was
found that the SVM model operates signi cantly better on Haralick texture
attributes [
        <xref ref-type="bibr" rid="ref4">4</xref>
        ], while the RF model achieves its best performance when combined
with histogram extraction. Both systems showed better performance when paired
with appropriate image ltering. In particular, isolating the saturation channel
of an input image during pre-processing works well with both systems. The SVM
system achieves even better performance when applying additional contrast and
binary thresholding to this isolated channel.
      </p>
      <p>
        The RF system outperforms the SVM system, though both improve upon
the recall achieved by previous non-deep attempts. Moreover, the RF system
provides comparable performance to that o ered by the best existing CNN-based
approach, achieving accuracy, recall and precision within 4% of that system [
        <xref ref-type="bibr" rid="ref12">12</xref>
        ].
These results are achieved when testing on 20 000 images, which is a signi cantly
larger testing set than has been seen in the existing literature. The RF system is
also computationally e cient, with an average runtime of less than two seconds
on a system without a dedicated GPU. The combination of high accuracy and
good computational e ciency suggests that this system may be suitable for
lowresource environments such as rural clinics, where the technology is needed the
most.
      </p>
      <p>Section 2 examines the existing work around automated malaria diagnosis,
while Section 3 discusses the design of the systems presented in this paper.
Section 4 outlines the experimental methodology that was followed. Section 5
discusses the results of the conducted experiments. Finally, Section 6 draws
relevant conclusions and suggests future work.
2</p>
    </sec>
    <sec id="sec-2">
      <title>Related work</title>
      <p>There has been signi cant research interest around the topic of automated
malaria diagnosis, and the existing literature has demonstrated varying levels
of success using a wide variety of supervised learning approaches, including deep
convolutional networks, transfer learning models and non-deep models.</p>
      <p>
        Rajaraman et al. used CNNs to achieve accuracy of 98.6% and later improved
this to 99.5% with an ensemble approach [
        <xref ref-type="bibr" rid="ref11">11</xref>
        ][
        <xref ref-type="bibr" rid="ref12">12</xref>
        ]. The improved model
demonstrated precision of 99.8%. However, these metrics are given at the patient level
rather than the cell level, which is more commonly used in the existing literature.
Unfortunately, the study used a dataset, provided by the U.S. National Library
of Medicine, that consists of images taken under the same conditions, with the
same staining method and from the same archive. As such, it is not clear how
well the model would generalise to images collected under di erent conditions.
The deep convolutional architecture limits the viability of the model being
deployed in areas with scarce computational resources. The ensemble approach of
the later attempt exacerbates this problem, as multiple CNNs must be run to
get the nal result.
      </p>
      <p>
        Transfer learning is another technique that involves using a CNN for feature
extraction before performing classi cation with an external algorithm.
Mehanian et al. [
        <xref ref-type="bibr" rid="ref8">8</xref>
        ] propose a solution for malaria diagnosis using transfer learning
with a feature extraction CNN and non-deep logistic regression classi er. They
achieved sensitivity of 91.6% and speci city of 94.1%. Though the results of their
study are not quite as positive as Rajaraman et al., they also tested the model
on substantially larger dataset that originates from 12 countries. This suggests
that the model developed by Mehanian et al. may be more robust and have
more real-world applicability. While the logistic regression classi er is a
nondeep supervised learning technique, the reliance on the underlying CNN feature
extractor suggests that the real-world application of this approach would still be
limited by its computational complexity.
      </p>
      <p>
        While deep learning approaches have been more prevalent in the existing
literature, some authors have used non-deep supervised learning techniques. Diaz
et al. [
        <xref ref-type="bibr" rid="ref3">3</xref>
        ] are the most successful example, having achieved 94% sensitivity and
99.7% speci city using an SVM based approach. However, the validity of these
results must be questioned, as they were obtained when running the model on
the full dataset used in the study, including those images used for training.
Unfortunately, their study makes use of a private datset, which hinders direct
comparisons. Other non-deep approaches have not seen comparable success, and
many have adopted awed evaluation methodologies [
        <xref ref-type="bibr" rid="ref6">6</xref>
        ] [
        <xref ref-type="bibr" rid="ref17">17</xref>
        ]. As such, they are
not discussed in detail here.
      </p>
    </sec>
    <sec id="sec-3">
      <title>System design</title>
      <p>This section begins by discussing the software used to develop the proposed
systems. The feature extraction and image ltering algorithms that were considered
are then outlined. Subsequently, the nal designs for the systems are presented.
3.1</p>
      <sec id="sec-3-1">
        <title>Software frameworks</title>
        <p>
          The systems proposed in this papers were constructed using Python 3.6 along
with several machine learning and computer vision libraries. Speci cally, the
systems use scikit-learn [
          <xref ref-type="bibr" rid="ref10">10</xref>
          ] for the underlying machine learning models. For
image processing and feature extraction tasks, two libraries were used, namely
Mahotas CV [
          <xref ref-type="bibr" rid="ref2">2</xref>
          ] and OpenCV [
          <xref ref-type="bibr" rid="ref1">1</xref>
          ].
3.2
        </p>
      </sec>
      <sec id="sec-3-2">
        <title>Feature extraction algorithms</title>
        <p>
          Feature extraction algorithms can be applied before input is supplied to machine
learning models in order to reduce the input dimensionality. By doing so, models
can run much faster and, in some cases, achieve higher performance. In this
paper, three approaches were considered:
Hu moments Hu [
          <xref ref-type="bibr" rid="ref5">5</xref>
          ] proposed a set of seven moment invariants, all of which
are invariant to scale, rotation and translation. These features, known as the Hu
moments, have been applied broadly to problems in the eld of computer vision.
For example, Otiniano-Rodriguez et al. [
          <xref ref-type="bibr" rid="ref9">9</xref>
          ] achieved over 90% accuracy using the
Hu moments as input to an SVM classi er for sign language recognition.
Haralick texture attributes Haralick et al. [
          <xref ref-type="bibr" rid="ref4">4</xref>
          ] presented a set of 14
textural features that can be extracted from images to improve image classi cation
accuracy. These features have since been used in a broad range of image
classi cation tasks. For example, Roula et al. [
          <xref ref-type="bibr" rid="ref13">13</xref>
          ] used them for classi cation of
prostatic neoplasia in microscopic images of samples taken by needle biopsy.
Histograms Histograms count the number of pixels that fall into a speci ed
number of intensity bins, for each colour channel. Existing works have used
histograms for various image classi cation tasks. For example, Szummer and
Picard [
          <xref ref-type="bibr" rid="ref16">16</xref>
          ] used histograms as part of a feature set for indoor-outdoor image
classi cation, resulting in 90.3% accuracy.
3.3
        </p>
      </sec>
      <sec id="sec-3-3">
        <title>Image lters</title>
        <p>Applying lters to images before they are passed to machine learning models
may serve to accentuate di erences in visual characteristics between positive
and negative cases. This may, in turn, result in improved model performance.
Five image ltering con gurations were considered for use in this paper:
1. No lters
2. HSV colour space conversion
3. Saturation channel isolation
4. Saturation channel isolation with contrast
5. Saturation channel isolation with contrast and binary thresholding
The SVM system uses the scikit-learn support vector classi er for its underlying
model. This model uses a third degree polynomial kernel function, as this showed
the best performance when compared to radial basis function, linear and fth
degree polynomial kernels. Gamma values between 0.001 and 1 were evaluated,
as well as the automatically generated gamma value provided by scikit-learn.</p>
        <p>A gamma value of 0.1 demonstrated the highest accuracy, but this was only
an improvement of 0.02% over the automatically generated value. The
automatically generated value achieved 0.51% higher recall at the cost of some precision.
Typically higher recall is preferred in the context of initial screening tests, and
so the automatically generated value was selected. Other hyperparameters were
set to the default scikit-learn values.</p>
        <p>The system extracts Haralick texture attributes as features. Initial
evaluations showed that, in combination with the binary thresholding image
ltering approach, this provided much higher accuracy than either Hu moments or
histogram features. Unfortunately, the computational cost of Haralick feature
extraction is much higher than that of these other approaches.
3.5</p>
      </sec>
      <sec id="sec-3-4">
        <title>Random Forest (RF) system</title>
        <p>The RF system uses the scikit-learn random forest classi er for its underlying
model. During hyperparameter tuning, forest sizes ranging from 1 to 250 were
evaluated and peak performance was seen at a value of 100. Maximum depths
ranging from 10 to 1000, as well as unlimited depth, were also evaluated. Minimal
performance improvements were seen above a value of 100, so it was decided
to select this value. Other hyperparameters were set to the default scikit-learn
values.</p>
        <p>Histogram extraction was adopted as the feature extraction approach for
this system. Initial evaluations showed that, in combination with the saturation
channel isolation lter, this provided the best classi cation accuracy and
computational performance. Haralick texture attributes resulted in slightly lower
accuracy and are far more computationally expensive. On the other hand, Hu
moments resulted in poor model performance, though they are not
computationally expensive.</p>
      </sec>
    </sec>
    <sec id="sec-4">
      <title>Experimental methodology</title>
      <p>This section lays out the process by which the research aims were addressed.
Speci cally, two research questions are posed:
1. Can the proposed non-deep systems achieve accuracy within 5% of the top
performing CNN approach?
2. Do the performances of models degrade when applied to di erently sourced
data?</p>
      <p>The details of the datasets and computational equipment used are presented,
followed by explanations of the two experiments that were conducted.
4.1</p>
      <sec id="sec-4-1">
        <title>Evaluation criteria</title>
        <p>Systems are evaluated along three evaluation criteria, namely accuracy, recall
and precision. Accuracy refers to the overall ability of a system to make correct
predictions, while recall refers to the ability of a system to not label positive
examples as negative. In other words, recall is the ability of a system to avoid
false negatives. On the other hand, precision refers to the ability of a system
to not label negative examples as positive or, in other words, to avoid false
positives. Speci city is another metric used by some papers, which measures the
ability of a system to correctly identify negative examples. Qualitatively, this is
quite similar to precision, though it is calculated di erently. Typically studies
will not report both speci city and precision, and so it becomes necessary to
draw comparisons between these metrics.
4.2</p>
      </sec>
      <sec id="sec-4-2">
        <title>Datasets</title>
        <p>
          Two datasets are used to run the experiments detailed in this study. The rst
is a publicly accessible dataset, originating from a study by Rajaraman et al.
[
          <xref ref-type="bibr" rid="ref11">11</xref>
          ] and made available by the U.S. National Library of Medicine (NLM). This
dataset is made up of pre-cropped Giemsa-stained blood cell images, with 13 779
classi ed as infected and the same number classi ed as uninfected, for a total
of 27 558 images. This NLM dataset was also used during the model tuning
process.
        </p>
        <p>The second is a private dataset provided to the author by PathCare
Laboratory Services. This dataset was provided as Giemsa-stained blood slide images,
not pre-cropped but with infected cells having been identi ed by a pathologist.
The dataset was manually cropped into a set of 120 cell images, with 60 labelled
as infected and 60 labelled as uninfected.</p>
        <p>Images in each dataset were not of a standardised size, and so it was decided
to resize all images to 50 pixel by 50 pixel squares. These measurements were
decided to avoid artifacts generated from excessive upscaling or stretching in
either axis, with most of the original images being about the same size or larger,
and roughly square in shape.</p>
        <p>Dataset splitting To ensure the validity of results, an iterative random
holdout approach was taken, whereby random subsets of the dataset were selected
at each iteration to form the training and test sets. The test set is not used for
training models so that when testing occurs, that data is unseen by the model.
The metrics achieved for each iteration are aggregated to form the nal result.</p>
      </sec>
      <sec id="sec-4-3">
        <title>Computational equipment</title>
        <p>Development, training and testing were all conducted on the same machine: a
laptop running Ubuntu 18.04.1 LTS. The machine had 16GB of RAM and a four
core Intel i7 CPU, with clock speeds of 2.7GHz. The machine did not have a
dedicated graphics card.
4.4</p>
      </sec>
      <sec id="sec-4-4">
        <title>Experiment design</title>
        <p>Two experiments were conducted on the proposed systems, allowing for
evaluation of various experimental hypotheses.</p>
        <p>Experiment 1: Performance of systems on NLM data Systems were
evaluated on their ability to operate on data that is part of the same dataset
used for training. Results were obtained from ten iterations, each time trained on
5000 images and tested on a disjoint set of 20 000 images. Initial testing indicated
that a training set of 5000 images was su cient for models to converge, so the rest
of the data could be used to form a more extensive testing set. This experiment
was designed to evaluate the ability of the systems to predict infection in images
collected in a similar way to the training data. It was hypothesised that the
systems would outperform existing non-deep attempts, and achieve performance
within 5% of that demonstrated by the top-performing CNN approach.
Experiment 2: Performance of systems on PathCare test data The
systems were put through a nal evaluation to test their ability to operate on
datasets other than the one used for training. Results were obtained by loading
models from the rst experiment, pre-trained on the NLM dataset, and tested
on the full PathCare dataset of 60 infected and 60 uninfected images. This
experiment was designed to evaluate the generalisability of systems trained on a
dataset gathered from a singular source, answering the second research question.
It was hypothesised that the prediction performance on the PathCare dataset
would be worse than that seen in the rst experiment, as it is unlikely that the
systems will generalise successfully without a broad range of collected data.
5</p>
      </sec>
    </sec>
    <sec id="sec-5">
      <title>Results and discussion</title>
      <p>In this section, the results of the two experiments are presented. The implications
of these results are discussed, and possible explanations are laid out. Finally,
limitations of the experimentation process are noted.
5.1</p>
      <sec id="sec-5-1">
        <title>Experiment 1: Performance of tuned models on NLM data</title>
        <p>
          The results observed when testing on the NLM dataset are encouraging. The
SVM system achieved an accuracy of 93.06%, recall of 96.12% and precision
of 90.58%. The best performing existing approach by Diaz et al. [
          <xref ref-type="bibr" rid="ref3">3</xref>
          ] does not
report precision but reports speci city of 99.7% and recall of 94%. While the
SVM system's precision is signi cantly lower than the speci city reported by
Diaz et al., it also demonstrates a 2.12% improvement in recall. The RF system
achieved an accuracy of 96.29%, recall of 96.06% and precision of 96.49%. This
amounts to a 2.06% increase in recall, with a 3.68% lower precision than the
speci city achieved by Diaz et al. See Figure 2 for a graphical representation of
this comparison.
        </p>
        <p>
          Medical experts argued that attaining higher recall while maintaining
comparable overall accuracy is better than attaining higher speci city or precision
in the context of initial screening tests [
          <xref ref-type="bibr" rid="ref7">7</xref>
          ]. This is because it is far more
important to ensure that false negatives do not occur, as these could lead to patients
not being given treatment when needed. Following the high sensitivity
screening test, another test with higher precision may be performed to ensure overall
accuracy. An automated malaria screening system could identify cells with a
high probability of infection and present these to a pathologist for con rmation.
This reduces the number of cells which must be checked by the expert while still
ensuring diagnostic precision. In this sense, the hypothesis that the proposed
systems would improve on existing supervised approaches is con rmed, though
it is acknowledged that the improvement is not observed across all metrics.
        </p>
        <p>The current work also evaluates models on signi cantly more testing data
than Diaz et al.: 20 000 blood cells compared to 12 557. Moreover, the
parasitemia of the testing data used by Diaz et al. is reported as 5.6%, amounting to
approximately 703 infected blood cells. On the other hand, the testing data used
in this experiment was balanced, with 10 000 infected cells, possibly resulting in
a more accurate re ection of the real-world performance that can be expected
from the system.</p>
        <p>
          As expected, the proposed systems did not achieve better results than the
CNN approaches detailed by Rajaraman et al. [
          <xref ref-type="bibr" rid="ref11">11</xref>
          ][
          <xref ref-type="bibr" rid="ref12">12</xref>
          ]. However, the RF model's
metrics all fell within 4%, and the SVM achieved recall within 4%. These results,
particularly those of the RF system, con rm the hypothesis that the proposed
systems would show performance within 5% of the best CNN-based approach.
Mehanian et al. [
          <xref ref-type="bibr" rid="ref8">8</xref>
          ] achieved recall of 91.6% and speci city of 94.1% with their
transfer learning model. Both the SVM and RF systems improve signi cantly on
this recall, and the RF also shows higher precision than the speci city reported
by Mehanian et al. The testing set used in this paper was signi cantly larger
than those used by both Mehanian et al. and Rajaraman et al., though that used
by the former was more diverse, with images taken under di erent conditions
and originating from 12 countries. See Figure 3 for a graphical representation of
this comparison.
        </p>
        <p>Papers in the existing literature have not reported training or testing time
metrics, so it is not possible to provide a comparison here. However, comparing
the RF and SVM systems proposed in this paper, it is clear that the RF system
is signi cantly less computationally expensive. The testing time of the SVM
system is a factor of 236 greater than that of the RF system, while the training
time is a factor of 53 greater. These runtime metrics suggest that the RF system
may be more suitable for deployment in situations where computational power
is limited. Table 1 shows the full results of this experiment.
Both the SVM and RF systems saw slight improvements in accuracy when run
on the PathCare dataset. While the precision of both decreased, by 1.03% for
the SVM system and 2.74% for the RF system, each achieved 100% recall. It
must be noted that these results are likely overly optimistic due to the small
size of the PathCare dataset, however this encourages further evaluation of the
systems on a larger PathCare dataset.</p>
        <p>The runtime metrics of this experiment con rm the observation made in the
rst experiment: The RF system is more computationally e cient than the SVM
system. While train times are not evaluated, as the systems were pre-trained,
the test time of the SVM system is a factor of 128 greater than that of the RF
system. Table 2 shows the full results of this experiment.
5.3</p>
      </sec>
      <sec id="sec-5-2">
        <title>Limitations</title>
        <p>The small size of the PathCare dataset limits the acceptability of the results
of the second experiment. While these results were very positive as an initial
evaluation, further experimentation on a larger set of PathCare data would be
necessary to con rm them.</p>
        <p>The data that was made available for this research was only classi ed in two
classes: as parasitised or non-parasitised. However, in reality there are multiple
species of the malaria parasite, each with various life stages. Because the data
was only classi ed in this binary manner, it is unclear whether the models
generalise to provide similar performance for all species and life stage combinations.
Unfortunately, it is time-consuming to manually label data in this manner and
the labelling must be conducted by an expert pathologist. For this reason, it
is di cult to acquire large enough datasets to allow for adequate training and
testing of models.</p>
        <p>
          Both datasets were made up of images of Giemsa-stained blood cells. It is not
clear whether the systems would maintain similar performance when di erent
staining methods are employed. Giemsa staining is the method recommended
by the World Health Organisation [
          <xref ref-type="bibr" rid="ref18">18</xref>
          ], but some clinics in lower-income areas
may vary from this. Large datasets of blood cells not stained with Giemsa are
di cult to acquire, as it is uncommon for major pathology practices to use other
staining methods.
        </p>
        <p>We tested a range of values for each hyperparameter during model tuning,
however, due to limited computational resources, a fairly coarse search was
performed. It is possible that a ner grid search or more sophisticated intelligent
search techniques may nd better con gurations.</p>
      </sec>
    </sec>
    <sec id="sec-6">
      <title>6 Conclusions</title>
      <p>The experiments presented in this paper demonstrate that non-deep supervised
learning techniques may be used as an alternative to popular deep learning
approaches for malaria classi cation. Several conclusions can be drawn from the
results of these experiments.</p>
      <p>Firstly, the proposed RF and SVM systems outperform existing non-deep
supervised approaches in terms of recall, which is arguably the most important
metric by which medical screening tests can be judged. Secondly, the RF system
appears to be a more suitable solution, as it achieves better accuracy and is
far more computationally e cient than the SVM system. Thirdly, the proposed
RF system is able to achieve recall, precision and accuracy within 4% of that
reported by the best existing CNN approach. This indicates that the system may
be a viable alternative in situations where the high computational power required
by CNN systems is not possible, such as in rural clinics. Fourthly, initial results
when run on a small dataset gathered from a di erent source to the training data
seem to indicate that the promising performance demonstrated by the proposed
systems may generalise to various imaging conditions. However, further work is
necessary to con rm this.</p>
      <p>Finally, it is noted that feature extraction and image processing can signi
cantly impact both the computational e ciency and classi cation performance
of systems. Extracting histogram features is computationally cheap and greatly
improves the accuracy of the RF model, but the SVM model saw better
classi cation performance when operating on Haralick texture attributes. However
these attributes are computationally expensive to extract, limiting the viability
of the SVM system. Hu moments performed the worst for both SVM and RF
models, suggesting they may be less suitable for the task of malaria detection.</p>
      <p>
        On the whole, it is shown that non-deep supervised learning techniques have
great promise for reducing the burden on medical professionals in performing
malaria diagnosis. This may, in turn, result in much faster times to diagnosis,
and allow quicker intervention, which is described by the WHO as the most
important factor in preventing severe cases and deaths from occurring [
        <xref ref-type="bibr" rid="ref19">19</xref>
        ]. The
higher computational e ciency of non-deep systems, such as those presented in
this paper, may allow for more widespread adoption, especially in areas where
extensive computational resources are not available. Thus, the value added by
these systems may have a signi cant impact on rural and poor communities,
where access to medical experts is typically limited.
6.1
      </p>
      <sec id="sec-6-1">
        <title>Future work</title>
        <p>The systems proposed in this paper show promise as a computationally cheap
alternative to CNN-based systems, with a lower training data requirement. As
such, future work to develop a fully integrated diagnosis system is warranted.
Such a system should include the full computational pipeline of automated blood
cell image cropping, ltering, feature extraction and classi cation.</p>
        <p>
          Rajaraman et al. [
          <xref ref-type="bibr" rid="ref12">12</xref>
          ] improved the performance of their initial CNN approach
by adopting an ensemble strategy. Similarly, future work may produce better
results by using an ensemble of high performing non-deep models, such as those
proposed in this paper.
7
        </p>
      </sec>
    </sec>
    <sec id="sec-7">
      <title>Ethics</title>
      <p>This paper involved the use of both a publicly-available dataset from the U.S.
National Library of Medicine, as well as a private dataset acquired from PathCare
Laboratory Services. Both datasets are de-identi ed, containing no personal or
demographic information on the patients corresponding to each blood cell
image. Ethics clearance was granted by both the University of Cape Town's Science
Faculty Research Ethics Committee and the PathCare Research Committee.
8</p>
    </sec>
    <sec id="sec-8">
      <title>Acknowledgements</title>
      <p>Thanks are extended to PathCare Laboratory services for their contribution of
blood sample data.</p>
    </sec>
  </body>
  <back>
    <ref-list>
      <ref id="ref1">
        <mixed-citation>
          1.
          <string-name>
            <surname>Bradski</surname>
            ,
            <given-names>G.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Kaehler</surname>
            ,
            <given-names>A.</given-names>
          </string-name>
          :
          <article-title>Learning OpenCV: Computer vision with the OpenCV library. "</article-title>
          <string-name>
            <surname>O'Reilly Media</surname>
          </string-name>
          ,
          <source>Inc."</source>
          (
          <year>2008</year>
          )
        </mixed-citation>
      </ref>
      <ref id="ref2">
        <mixed-citation>
          2.
          <string-name>
            <surname>Coelho</surname>
            ,
            <given-names>L.P.</given-names>
          </string-name>
          : Mahotas:
          <article-title>Open source software for scriptable computer vision</article-title>
          .
          <source>arXiv preprint arXiv:1211.4907</source>
          (
          <year>2012</year>
          )
        </mixed-citation>
      </ref>
      <ref id="ref3">
        <mixed-citation>
          3.
          <string-name>
            <given-names>D</given-names>
            <surname>az</surname>
          </string-name>
          , G.,
          <string-name>
            <surname>Gonzalez</surname>
            ,
            <given-names>F.A.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Romero</surname>
          </string-name>
          , E.:
          <article-title>A semi-automatic method for quanti cation and classi cation of erythrocytes infected with malaria parasites in microscopic images</article-title>
          .
          <source>Journal of Biomedical Informatics</source>
          <volume>42</volume>
          (
          <issue>2</issue>
          ),
          <volume>296</volume>
          {
          <fpage>307</fpage>
          (
          <year>2009</year>
          )
        </mixed-citation>
      </ref>
      <ref id="ref4">
        <mixed-citation>
          4.
          <string-name>
            <surname>Haralick</surname>
            ,
            <given-names>R.M.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Shanmugam</surname>
            ,
            <given-names>K.</given-names>
          </string-name>
          , et al.:
          <article-title>Textural features for image classi cation</article-title>
          .
          <source>IEEE Transactions on systems, man, and cybernetics (6)</source>
          ,
          <volume>610</volume>
          {
          <fpage>621</fpage>
          (
          <year>1973</year>
          )
        </mixed-citation>
      </ref>
      <ref id="ref5">
        <mixed-citation>
          5.
          <string-name>
            <surname>Hu</surname>
            ,
            <given-names>M.K.</given-names>
          </string-name>
          :
          <article-title>Visual pattern recognition by moment invariants</article-title>
          .
          <source>IRE transactions on information theory 8</source>
          (
          <issue>2</issue>
          ),
          <volume>179</volume>
          {
          <fpage>187</fpage>
          (
          <year>1962</year>
          )
        </mixed-citation>
      </ref>
      <ref id="ref6">
        <mixed-citation>
          6.
          <string-name>
            <surname>Khan</surname>
            ,
            <given-names>N.A.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Pervaz</surname>
            ,
            <given-names>H.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Latif</surname>
            ,
            <given-names>A.K.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Musharraf</surname>
            ,
            <given-names>A.</given-names>
          </string-name>
          , et al.:
          <article-title>Unsupervised identi cation of malaria parasites using computer vision</article-title>
          .
          <source>In: 2014 11th International Joint Conference on Computer Science and Software Engineering (JCSSE)</source>
          . pp.
          <volume>263</volume>
          {
          <fpage>267</fpage>
          .
          <string-name>
            <surname>IEEE</surname>
          </string-name>
          (
          <year>2014</year>
          )
        </mixed-citation>
      </ref>
      <ref id="ref7">
        <mixed-citation>
          7.
          <string-name>
            <surname>Lalkhen</surname>
            ,
            <given-names>A.G.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>McCluskey</surname>
            ,
            <given-names>A.:</given-names>
          </string-name>
          <article-title>Clinical tests: sensitivity and speci city</article-title>
          .
          <source>Continuing Education in Anaesthesia Critical Care &amp; Pain</source>
          <volume>8</volume>
          (
          <issue>6</issue>
          ),
          <volume>221</volume>
          {
          <fpage>223</fpage>
          (
          <year>2008</year>
          )
        </mixed-citation>
      </ref>
      <ref id="ref8">
        <mixed-citation>
          8.
          <string-name>
            <surname>Mehanian</surname>
            ,
            <given-names>C.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Jaiswal</surname>
            ,
            <given-names>M.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Delahunt</surname>
            ,
            <given-names>C.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Thompson</surname>
            ,
            <given-names>C.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Horning</surname>
            ,
            <given-names>M.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Hu</surname>
            ,
            <given-names>L.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Ostbye</surname>
            ,
            <given-names>T.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>McGuire</surname>
            ,
            <given-names>S.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Mehanian</surname>
            ,
            <given-names>M.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Champlin</surname>
            ,
            <given-names>C.</given-names>
          </string-name>
          , et al.:
          <article-title>Computer-automated malaria diagnosis and quantitation using convolutional neural networks</article-title>
          .
          <source>In: Proceedings of the IEEE International Conference on Computer Vision</source>
          . pp.
          <volume>116</volume>
          {
          <issue>125</issue>
          (
          <year>2017</year>
          )
        </mixed-citation>
      </ref>
      <ref id="ref9">
        <mixed-citation>
          9.
          <string-name>
            <surname>Otiniano-Rodr guez</surname>
            ,
            <given-names>K.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Camara-Chavez</surname>
            ,
            <given-names>G.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Menotti</surname>
            ,
            <given-names>D.</given-names>
          </string-name>
          :
          <article-title>Hu and zernike moments for sign language recognition</article-title>
          .
          <source>In: Proceedings of international conference on image processing</source>
          ,
          <source>computer vision, and pattern recognition</source>
          . pp.
          <volume>1</volume>
          {
          <issue>5</issue>
          (
          <year>2012</year>
          )
        </mixed-citation>
      </ref>
      <ref id="ref10">
        <mixed-citation>
          10.
          <string-name>
            <surname>Pedregosa</surname>
            ,
            <given-names>F.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Varoquaux</surname>
            ,
            <given-names>G.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Gramfort</surname>
            ,
            <given-names>A.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Michel</surname>
            ,
            <given-names>V.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Thirion</surname>
            ,
            <given-names>B.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Grisel</surname>
            ,
            <given-names>O.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Blondel</surname>
            ,
            <given-names>M.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Prettenhofer</surname>
            ,
            <given-names>P.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Weiss</surname>
            ,
            <given-names>R.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Dubourg</surname>
            ,
            <given-names>V.</given-names>
          </string-name>
          , et al.:
          <article-title>Scikit-learn: Machine learning in python</article-title>
          .
          <source>Journal of machine learning research 12(Oct)</source>
          ,
          <volume>2825</volume>
          {
          <fpage>2830</fpage>
          (
          <year>2011</year>
          )
        </mixed-citation>
      </ref>
      <ref id="ref11">
        <mixed-citation>
          11.
          <string-name>
            <surname>Rajaraman</surname>
            ,
            <given-names>S.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Antani</surname>
            ,
            <given-names>S.K.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Poostchi</surname>
            ,
            <given-names>M.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Silamut</surname>
            ,
            <given-names>K.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Hossain</surname>
            ,
            <given-names>M.A.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Maude</surname>
            ,
            <given-names>R.J.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Jaeger</surname>
            ,
            <given-names>S.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Thoma</surname>
            ,
            <given-names>G.R.</given-names>
          </string-name>
          :
          <article-title>Pre-trained convolutional neural networks as feature extractors toward improved malaria parasite detection in thin blood smear images</article-title>
          .
          <source>PeerJ</source>
          <volume>6</volume>
          ,
          <issue>e4568</issue>
          (
          <year>2018</year>
          )
        </mixed-citation>
      </ref>
      <ref id="ref12">
        <mixed-citation>
          12.
          <string-name>
            <surname>Rajaraman</surname>
            ,
            <given-names>S.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Jaeger</surname>
            ,
            <given-names>S.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Antani</surname>
            ,
            <given-names>S.K.</given-names>
          </string-name>
          :
          <article-title>Performance evaluation of deep neural ensembles toward malaria parasite detection in thin-blood smear images</article-title>
          .
          <source>PeerJ</source>
          <volume>7</volume>
          ,
          <issue>e6977</issue>
          (
          <year>2019</year>
          )
        </mixed-citation>
      </ref>
      <ref id="ref13">
        <mixed-citation>
          13.
          <string-name>
            <surname>Roula</surname>
            ,
            <given-names>M.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Diamond</surname>
            ,
            <given-names>J.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Bouridane</surname>
            ,
            <given-names>A.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Miller</surname>
            ,
            <given-names>P.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Amira</surname>
            ,
            <given-names>A.</given-names>
          </string-name>
          :
          <article-title>A multispectral computer vision system for automatic grading of prostatic neoplasia</article-title>
          .
          <source>In: Proceedings IEEE International Symposium on Biomedical Imaging</source>
          . pp.
          <volume>193</volume>
          {
          <fpage>196</fpage>
          .
          <string-name>
            <surname>IEEE</surname>
          </string-name>
          (
          <year>2002</year>
          )
        </mixed-citation>
      </ref>
      <ref id="ref14">
        <mixed-citation>
          14.
          <string-name>
            <surname>Russell</surname>
            ,
            <given-names>S.:</given-names>
          </string-name>
          <article-title>The economic burden of illness for households in developing countries: a review of studies focusing on malaria, tuberculosis, and human immunode ciency virus/acquired immunode ciency syndrome</article-title>
          .
          <source>The American journal of tropical medicine and hygiene 71</source>
          ,
          <volume>147</volume>
          {
          <fpage>155</fpage>
          (
          <year>2004</year>
          )
        </mixed-citation>
      </ref>
      <ref id="ref15">
        <mixed-citation>
          15.
          <string-name>
            <surname>Shillcutt</surname>
            ,
            <given-names>S.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Morel</surname>
            ,
            <given-names>C.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Goodman</surname>
            ,
            <given-names>C.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Coleman</surname>
            ,
            <given-names>P.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Bell</surname>
            ,
            <given-names>D.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Whitty</surname>
            ,
            <given-names>C.J.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Mills</surname>
            ,
            <given-names>A.</given-names>
          </string-name>
          :
          <article-title>Cost-e ectiveness of malaria diagnostic methods in sub-saharan africa in an era of combination therapy</article-title>
          . Bulletin of the World Health Organization pp.
          <volume>101</volume>
          {
          <issue>110</issue>
          (
          <year>2008</year>
          )
        </mixed-citation>
      </ref>
      <ref id="ref16">
        <mixed-citation>
          16.
          <string-name>
            <surname>Szummer</surname>
            ,
            <given-names>M.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Picard</surname>
          </string-name>
          , R.W.:
          <article-title>Indoor-outdoor image classi cation</article-title>
          .
          <source>In: Proceedings 1998 IEEE International Workshop on Content-Based Access of Image and Video Database</source>
          . pp.
          <volume>42</volume>
          {
          <fpage>51</fpage>
          .
          <string-name>
            <surname>IEEE</surname>
          </string-name>
          (
          <year>1998</year>
          )
        </mixed-citation>
      </ref>
      <ref id="ref17">
        <mixed-citation>
          17.
          <string-name>
            <surname>Tek</surname>
            ,
            <given-names>F.B.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Dempster</surname>
            ,
            <given-names>A.G.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Kale</surname>
            ,
            <given-names>I.</given-names>
          </string-name>
          :
          <article-title>Parasite detection and identi cation for automated thin blood lm malaria diagnosis</article-title>
          .
          <source>Computer vision and image understanding 114(1)</source>
          ,
          <volume>21</volume>
          {
          <fpage>32</fpage>
          (
          <year>2010</year>
          )
        </mixed-citation>
      </ref>
      <ref id="ref18">
        <mixed-citation>
          18. WHO:
          <article-title>Giemsa staining of malaria blood lms (</article-title>
          <year>2018</year>
          )
        </mixed-citation>
      </ref>
      <ref id="ref19">
        <mixed-citation>
          19. WHO:
          <source>World malaria report 2018</source>
          . World Health Organization (
          <year>2018</year>
          )
        </mixed-citation>
      </ref>
    </ref-list>
  </back>
</article>