Content-based histopathology image retrieval using CometCloud

  • Xin Qi1, 2Email author,

    Affiliated with

    • Daihou Wang2, 3,

      Affiliated with

      • Ivan Rodero3,

        Affiliated with

        • Javier Diaz-Montes3,

          Affiliated with

          • Rebekah H Gensure1,

            Affiliated with

            • Fuyong Xing5,

              Affiliated with

              • Hua Zhong1,

                Affiliated with

                • Lauri Goodell1,

                  Affiliated with

                  • Manish Parashar3,

                    Affiliated with

                    • David J Foran1, 2, 4 and

                      Affiliated with

                      • Lin Yang5

                        Affiliated with

                        BMC Bioinformatics201415:287

                        DOI: 10.1186/1471-2105-15-287

                        Received: 26 November 2013

                        Accepted: 12 August 2014

                        Published: 26 August 2014



                        The development of digital imaging technology is creating extraordinary levels of accuracy that provide support for improved reliability in different aspects of the image analysis, such as content-based image retrieval, image segmentation, and classification. This has dramatically increased the volume and rate at which data are generated. Together these facts make querying and sharing non-trivial and render centralized solutions unfeasible. Moreover, in many cases this data is often distributed and must be shared across multiple institutions requiring decentralized solutions. In this context, a new generation of data/information driven applications must be developed to take advantage of the national advanced cyber-infrastructure (ACI) which enable investigators to seamlessly and securely interact with information/data which is distributed across geographically disparate resources. This paper presents the development and evaluation of a novel content-based image retrieval (CBIR) framework. The methods were tested extensively using both peripheral blood smears and renal glomeruli specimens. The datasets and performance were evaluated by two pathologists to determine the concordance.


                        The CBIR algorithms that were developed can reliably retrieve the candidate image patches exhibiting intensity and morphological characteristics that are most similar to a given query image. The methods described in this paper are able to reliably discriminate among subtle staining differences and spatial pattern distributions. By integrating a newly developed dual-similarity relevance feedback module into the CBIR framework, the CBIR results were improved substantially. By aggregating the computational power of high performance computing (HPC) and cloud resources, we demonstrated that the method can be successfully executed in minutes on the Cloud compared to weeks using standard computers.


                        In this paper, we present a set of newly developed CBIR algorithms and validate them using two different pathology applications, which are regularly evaluated in the practice of pathology. Comparative experimental results demonstrate excellent performance throughout the course of a set of systematic studies. Additionally, we present and evaluate a framework to enable the execution of these algorithms across distributed resources. We show how parallel searching of content-wise similar images in the dataset significantly reduces the overall computational time to ensure the practical utility of the proposed CBIR algorithms.


                        Histopathology Digital pathology Content-based image retrieval High performance computing


                        A growing number of leading institutions now routinely utilize digital imaging technologies to support investigative research and routine diagnostic procedures. The exponential rate at which images and videos are being generated has resulted in a significant need for efficient content-based image retrieval (CBIR) methods, which allow one to quickly characterize and locate images in large collections based upon the features of a given query image. CBIR has been one of the most active research areas in a wide spectrum of imaging informatics fields over the past few decades [113]. Several domains stand to benefit from the use of CBIR including cinematography, education, investigative basic and clinical research, and the practice of medicine. CBIR has been successfully utilized in applications spanning radiology [4, 11, 14, 15], pathology [9, 1618], dermatology [19, 20], and cytology [2123].

                        There have been several successful CBIR systems that have been developed for medical applications since the 1980’s. Several approaches utilize simple features such as color histograms [24], shape [4, 22], texture [6, 25], or fuzzy features [7] to characterize the content of images while allowing higher level diagnostic abstractions based on systematic queries [4, 2527]. The recent adoption and popularity of case-based reasoning [28] and evidence-based medicine [29] has created a compelling need for more reliable image retrieval strategies to support diagnostic decisions. In fact, a number of state-of-the-art CBIR systems [4, 9, 1113, 15, 16, 25, 3032] have been designed to support the processing of queries across imaging modalities.

                        With the advent of whole-slide imaging technology, the size and scale of image-based data has grown tremendously, making it impractical to perform matching operations across an entire image dataset using traditional methods. To meet this challenge, a new family of strategies are being developed, which enable investigators to perform sub-region searching to automatically identify image patches that exhibit patterns that are consistent with a given query patch. In practice, this approach makes it possible to select a region or object of interest within a digitized specimen as a query while the algorithm systematically identifies regions exhibiting similar characteristics in either the same specimen or across disparate specimens. The results can then be used to draw comparisons among patient samples in order to make informed decisions regarding likely prognoses and most appropriate treatment regimens.

                        To perform a region-of-interest (ROI) query, Vu et al. [33] presented a Sam Match framework-based similarity model. The use of a part-based approach was later reported in [34] to solve the CBIR problem by synthesizing a DoG detector, and a local hashing table search algorithm. The primary limitation of this approach, however, was the time cost of the large number of features that need to be computed. Intra-expansion and inter-expansion strategies were later developed to boost the hash-based search quality based on a bag-of-features model which could more accurately represent the images. Recently, a structured visual search method was developed to perform CBIR in medical image datasets [35]. The primary advantage of this framework is that it is flexible and can be quickly extended to other modalities.

                        Most CBIR algorithms rely on content localization, feature extraction, and user feedback steps [57, 25, 27, 3640]. The retrieved results are then ranked by some criteria, such as appearance similarity or diagnostic relevance, which can also serve as a measure of the practical usability of the algorithm. Typically the retrieved images only include those cases with the most similar appearance to a given query image whereas introducing relevance feedback [4147] to CBIR provides a practical means for addressing the semantic gap between visual and semantic similarity.

                        Large-scale image retrieval applications are generally computationally expensive. In this paper, we present the use of the CometCloud [48, 49] to execute CBIR in a parallel fashion on multiple high performance computing (HPC) and cloud resources as a means for reducing computational time significantly. CometCloud is an autonomic cloud framework that allows dynamic, on-demand federation of distributed infrastructures. It also provides an effective programming platform that supports MapReduce, Workflow, and Master-Worker/BOT models making it possible for investigators to quickly develop applications that can run across the federated resources [4953]. The algorithm that our team developed exploits the parallelism of CBIR by combining the HPC assets at Rutgers University with external cloud resources. Moreover, our solution uses cloud abstractions to federate resources elastically to achieve acceleration, while hiding infrastructure and deployment details. In this way, the CBIR algorithm can be made available as accessible services to end users.

                        The contributions of this paper are: 1) a novel CBIR algorithm based on a newly developed coarse-to-fine searching criteria which is coupled with a novel feature called hierarchical annular histogram (HAH); 2) a CBIR refinement schema based on dual-similarity relevance feedback; and 3) a reliable parallel implementation of the CBIR algorithm based on Cloud computing.


                        Research design

                        After discussing the needs and requirements of pathologists from their perspective, the CBIR study is designed to quickly and accurately find images exhibiting similar morphologic and staining characteristics throughout a single or collection of imaged specimens. Our team specifically choose to use Giemsa stained peripheral blood smear and hematoxylin and eosin (H&E) stained renal glomeruli datasets to systematically test the algorithms since these are two routine use case scenarios that our clinical colleagues indicated might benefit from the proposed technology. Leukocytes are often differentiated based on traditional morphological characteristics, however the subtle visible differences exhibited by some lymphomas and leukemias result in a significant number of false negative during routine screenings. In many cases, the diagnosis is only rendered after conducting immunophenotyping and a range of other molecular or cytogenetic studies. The additional studies are expensive, time consuming, and usually require fresh tissues that may not be readily available [54]. Pre-transplantation biopsies of kidney grafts have become a routine means for selecting organs which are suitable for transplantation from marginal donors. The main histopathology characteristics that are routinely evaluated by pathologists are percentage of glomerulosclerosis, interstitial fibrosis, and degree of vascular pathology [55]. The central incentive for developing the CBIR algorithms is to determine a reliable means for assisting pathologists when they are called upon to render diagnostic decisions based on whole-slide scanned specimens.

                        In this paper, we present a novel content-based image retrieval (CBIR) algorithm that is systematically tested on both imaged Giemsa stained peripheral blood smears and digitized H&E stained renal glomeruli specimens. Because of the intense computational requirements of the algorithms, our team systematically investigate the use of high performance computing solutions based on CometCloud to distribute the tasks of performing CBIR to significantly reduce the overall running time. The details of datasets, the relevant CBIR algorithms, and the CometCloud implementation of the methods are explained in detail in the following sections.

                        In the case of Giemsa stained peripheral blood smear datasets, the algorithms operate on a given query patch to quickly and reliably detect other leukocytes of the same class throughout the imaged specimen in support of diagnostic decisions. These hematopathology datasets were acquired using a 20× objective to provide a gross overview of the specimen while also supplying sufficient resolution to distinguish among different classes of leukocytes. The dataset consisted of 925 imaged blood smears (1000×1000 pixels). In the case of the H&E stained renal glomeruli datasets, the algorithms are used to process any given query patch to discriminate necrotic glomeruli and normal glomeruli throughout imaged kidney tissue specimens. In these experiments, our team cropped 32 images (5024×3504 pixels) from within eight whole-slide renal specimens using a 20× objective.

                        Quality control of all datasets was conducted by an experienced pathologist (Dr. Zhong) whereas query image patches and ground-truth classification were determined by two pathologists (Dr. Zhong and Dr. Goodell). The retrieved results were evaluated by both pathologists through a completely independent and blinded process. During the peripheral blood smear experiments, pathologists were asked to assign each leukocyte retrieved using the CBIR algorithm to either the relevant or non-relevant class as a means for judging the appropriateness of each returned patch. In all, there were five different classes of leukocytes used in the studies. During the renal glomeruli studies, either a relevant or non-relevant assignment was made to judge the performance of the algorithms in distinguishing between necrotic glomeruli and normal glomeruli.

                        The CBIR algorithms consist of four major steps: 1) regions of interest (ROIs) localization, 2) hierarchical three-stage searching, 3) retrieval refinement based on dual-similarity relevance feedback, and 4) high performance computing using CometCloud [48]. Figure 1 illustrates the actual workflow of the process.
                        Figure 1

                        Workflow of the proposed CBIR algorithm.

                        Step 1: regions of interest localization

                        The first step is to locate the regions of interest (ROIs) throughout the imaged specimens by excluding the background regions from the candidate objects. Using color-decomposition and morphology [56] based preprocessing, the algorithm identifies application-specific ROIs. These regions serve as candidate searching regions in the subsequent stages of hierarchical searching. Candidate image patches are generated using a sliding window approach with an overlapping ratio within the range of [ 50%,90%].

                        Step 2: hierarchical three-stage searching

                        The hierarchical three-stage searching method includes: coarse searching, fine searching, and mean-shift clustering.

                        Coarse searching

                        Let Q represents a query image patch and P serves the candidate image patches. Each patch is divided into consecutive concentric rectangular bin regions (termed as rings) as shown in Figure 2(a-b). As the number of rings, r, increases, more detailed image characteristics are captured and while the computational time increases accordingly. r is determined based on cross-validation. Figure 2(b) illustrates the process of coarse searching. Given a query image patch, the algorithm computes local features from the innermost ring. Based on a similarity measure between candidate image patches, P, and the query image, Q, retrieved image patches, P, are ranked from high to low, and only the top 50% ranked candidates are reserved at each step. This procedure continues until the outermost ring is reached. This cascade structure significantly reduces the computational time, as 50% of the image patches are eliminated in the very first stage of processing by simply evaluating features in the innermost ring.
                        Figure 2

                        An illustration of the hierarchical searching framework: (a) region of interest, (b) coarse searching step, and (c) fine searching step.

                        Fine searching

                        After the coarse searching stage has been completed, each rectangular annular ring from both the query and candidate patches are equally subdivided into eight segments, and local features are calculated in each segment. The final candidates are chosen based on a similarity measure of a concatenated feature vector corresponding to the eight segments. Figure 2(c) illustrates the process of the fine searching. This stage is designed to capture the spatial configuration of the local features. Due to the limited number of candidates passing through the coarse searching stage, the computational time for completing this stage is dramatically reduced.

                        Mean-shift clustering

                        In order to assemble the final retrieval results, mean-shift (MS) clustering [57] is applied to the top ranked candidate patches, which have survived both the coarse and fine searching stages. The bandwidth b for the mean-shift clustering is calculated as , where w is the width of the query image and h is the height of the query image. In this way, the final CBIR results are obtained.

                        HAH Feature and feature comparison

                        HAH feature

                        To implement the hierarchical searching framework, we develop a hierarchical annular histogram (HAH). The intensity color histograms of consecutive concentric rectangular rings are calculated and concatenated together to form a coarse searching feature vector, H c =(h 1,h 2,…,h r ), where h i is the intensity color histogram of the ith ring, i∈ [ 1,r] and r is the number of rings selected for the HAH feature. For fine searching, each rectangular annular ring is equally divided into eight segments, and the color histogram is calculated from each segment sequentially and then concatenated together to form the fine searching feature vector, H f =(h 1,1,…,h 1,8,h 2,1,…,h 2,8,…,h r,1,…,h r,8), where h i,j is the intensity color histogram of the ith ring within the jth segment, j∈ [ 1,8]. Here superscript c represents coarse searching and f represents fine searching. Throughout the CBIR study, we use Euclidean distance as the similarity measure. The distance D i , between the ith candidate patch v i and the query patch q in coarse searching and fine searching are defined as and , respectively:

                        where . Here H c (q c ),H f (q f ) are the feature vector of query image during coarse searching and fine searching stages, respectively, and are the feature vector of the ith candidate patch in the coarse searching and fine searching stages, respectively.

                        Figure 3(a) and (b) illustrate the calculation of the HAH from the innermost rectangle and the fourth ring from the center. Figure 3(c) and (d) show an example of two image patches with similar traditional color histogram (d), but completely different HAH (c). This demonstrates the capacity of the HAH to differentiate among image patches exhibiting similar total color distributions, but different spatial configurations.
                        Figure 3

                        An illustration of HAH calculation. (a) Color histogram of the central ring. (b) Color histogram of the fourth ring from the center. An example of two patches with (c) different HAH, but (d) similar color histogram of the entire image.

                        In order to compare the performance of the HAH feature in CBIR, the Gabor wavelet feature [58] and co-occurrence texture feature [59, 60] were compared with the HAH feature with respect to both speed and accuracy using both imaged peripheral blood smear and renal glomeruli datasets. For the purpose of the studies, precision and recall were used to measure the performance of the CBIR algorithm. Precision is defined as the ratio between the number of retrieved relevant images and the total number of retrieved images. Recall is defined as the ratio between the number of retrieved relevant images and the total number of relevant images in the datasets.

                        The Gabor wavelet feature

                        The Gabor wavelet feature is used to describe the image patterns at a range of different directions and scales. Throughout the experiments, we utilize a Gabor filter with 8 directions and 5 scales, (M=5,N=8), and the mean value and standard deviation of each filtered image are concatenated to form a feature vector: f=(μ 1,1,σ 1,1,μ 1,2,σ 1,2,…,μ 5,8,σ 5,8), in which μ m,n and σ m,n represent the mean value and standard deviation of the filtered image using Gabor filter at the mth scale and nth direction, m∈ [ 1,M],n∈ [ 1,N]. The distance D i between the ith candidate patch v i , and the query patch q, is defined as

                        where .

                        COOC texture feature

                        Co-occurrence (COOC) matrices, also called spatial gray-level dependence matrices, were first proposed by Haralick et al. [59, 60]. COOC matrices are calculated from an estimation of the second-order joint conditional probability of the image intensity with various distances and four specific orientations (00, 450, 900, 1350). COOC texture feature using the COOC matrices quantifies the distribution of gray-level values within an image. For the feature comparison experiment, COOC texture feature including contrast, correlation, energy, and homogeneity [60], is calculated from the COOC matrices within the candidate ROIs and the query image. The distance, D i , between the ith candidate patch v i , and the query patch q, is defined as

                        where , and F = {contrast, correlation, energy, homogeneity}.

                        Stage 3: CBIR retrieval refinement using a dual-similarity relevance feedback

                        Relevance feedback is an interactive procedure which is used to refine the initial retrieval results. Upon completion of the initial retrieval, top ranked retrieval images were reviewed by two pathologists with consensus to label them as relevant or non-relevant as previously described. These responses are used as users’ feedback to re-rank the retrieval results accordingly.

                        Two types of similarities are used in the above retrieval and feedback procedure: similarity in visual appearance as measured by image feature distance and similarity in semantic category as measured as relevant or non-relevant. Current relevance feedback algorithms typically only consider the second similarity. In our algorithm, we develop a dual-similarity schema that combines both types of similarity measures. This is achieved by rebuilding the initial distributions of training samples in an on-line manner.

                        For each top ranked retrieved image, a 256×3×r dimension feature vector is constructed, where r is the number of rings defined in the hierarchical searching process. Dimension reduction using principal component analysis (PCA) is applied to the original HAH feature space, and the top principal components accounting for 90% of the total variance are used as inputs for the following relevance feedback procedure.

                        Adaboost [61] is utilized to train an ensemble classifier composed of a set of weak learners. Given a training dataset, a strong classifier is built as a weighted sum of weak learners by minimizing the misclassification errors. Define weight W i , to be measured by a normalized Euclidean distance D i , representing the image appearance similarity between a pair of retrieved image and the original query. The initial distribution of the training samples is recalculated to update the classifier to place more weights on the visually similar cases following the relevance feedback step. The algorithm is summarized as follows.

                        Step 4: accelerating CBIR using CometCloud

                        Due to the data-independence property of the CBIR algorithm, we can formulate our problem as a set of heterogeneous and independent or loosely couple tasks. In this way, we can parallelize and solve our problem using the aggregated computational power of distributed resources. Our team has designed and developed a framework that enables the execution of CBIR across distributed, federated resources. Our framework uses cloud abstractions to present the underlying infrastructure as a single elastic pool of resources regardless of their physical location or specific particularities. In this way, computational resources are dynamically provisioned on-demand to meet the application’s requirements. These resources can be high performance computing grids, clouds, or supercomputers. In the current application, the framework is built on top of CometCloud [48]. CometCloud is purposely chosen for this application since it enables dynamic and on-demand federation of advanced cyber-infrastructures (ACIs). It also provides a flexible application programming interface (API), for developing applications that can take advantage of federated ACIs. Furthermore, it provides fault-tolerance in the resulting infrastructure.

                        The framework used to run the CBIR algorithm across federated resources is implemented using the master/worker paradigm. In this scenario, the CBIR software serves as a computational engine, while CometCloud orchestrates the entire execution. The master/worker model is suitable for problems with a large pool of independent tasks, where both the tasks and the resources are heterogeneous. Using this approach, the master component generates tasks, collects results, and verifies that tasks are properly executed. Each task contains the description of the images to be processed. All tasks are automatically placed in the CometCloud-managed distributed task space for execution. Workers are dedicated to carry out tasks pulled from the CometCloud task space and send results back to the master.

                        The implementation that we have presented has several important and highly desirable properties. From the user’s perspective, the framework creates a cloud abstraction on top of the resources that hides infrastructure details and offers the CBIR software as a readily accessible service. In this way, one can query the database using different algorithms via a simple interface without consideration of how and where queries are executed. On the other hand, from the developer’s perspective, the integration of the existing CBIR software with the CometCloud framework does not require any adjustments on the application side. Additionally, the resulting framework completely operates within the limits of the end-user space. This means that it is possible to aggregate computational resources without special privileges, which is very important when using external resources.

                        Results and discussion

                        CBIR results and feature comparison

                        A dual-processor system based on Intel Xeon E5530@2.4 GHz with 24 GB RAM and 64-bit operating system was used for the CBIR study. Initial CBIR results using two pathology image datasets and different feature comparison are presented below. Figure 4 shows an example of the CBIR three-stage hierarchical searching results using one neutrophil as a query image in a peripheral blood smear dataset acquired using 20× magnification objective. Green box labeled regions represent the candidate patches that are similar to the query image patch. Figure 5 shows CBIR results using different classes of leukocytes as query images, including basophil, eosinophil, lymphocyte, monocyte, and neutrophil, respectively. Green box labeled regions represent the candidate patches that are similar to the query image patch. Each box has a number to indicate the ranking order of every candidate patch in the dataset. Figure 6 shows an example of CBIR results for a necrotic glomeruli query image using a testing dataset containing multi-scale regions at 1/2, 1, 2, 3, and 4 times of the original size of the query image. Red box labeled regions indicate the query image. Blue box labeled regions represent the healthy glomeruli for comparison. Green box labeled regions represent the top-ranked 10% of retrieval patches of the 32 randomly selected regions (5024×3504 pixels) cropped from whole-slide scanned images.
                        Figure 4

                        An illustration of results of the three-stage CBIR searching using one neutrophil as a query image from peripheral blood smears acquired at 20× objective, in which green box labeled regions represent the candidate patches.

                        Figure 5

                        CBIR results using different classes of leukocytes as query images, including basophil, eosinophil, lymphocyte, monocyte, and neutrophil, respectively. Here green box labeled regions represent the candidate patches that are similar to the query image patch. Each box has a number to indicate the ranking order of every candidate patch in the dataset. The original sizes of the images were adjusted to fit in the figure.

                        Figure 6

                        An example of top 10% CBIR results for a necrotic glomerulus query image. Red box labeled regions indicate the query image. Blue box labeled regions represent the healthy glomeruli for comparison. Green box labeled regions denote the top 10% ranked retrieved patches, which include multiple scaled regions at 1/2, 1, 2, 3, and 4 times of the original size of the query image. The original sizes of the images were adjusted to fit in the figure.

                        By varying the number of rings ∈ [ 2,3,5,10,15] in the hierarchical searching, the performance of CBIR is summarized as follows. For imaged peripheral blood smears, all five classes of leukocytes were correctly retrieved using three inner rings of the HAH. For imaged renal glomeruli, as the number of rings increased to 10, all necrotic glomeruli were correctly retrieved. With an increase of the number of the rings, the computational time also increased. The number of rings was shown to be dependent upon the complexity of the dataset.

                        For local feature comparison, image retrieval was performed on the same datasets with the same query images using HAH, Gabor wavelet, and COOC texture features. Figure 7(a) and (b) show precision-recall curves and average of feature calculation times using peripheral blood smear images, respectively. Figure 7(c) and (d) show precision-recall curves and average of feature calculation times using renal glomeruli images, respectively. The area under a curve (AUC) value of each feature for peripheral blood smear images and renal glomeruli images are shown in Figure 7(a) and (c), respectively. The average of feature computation times are shown in Figure 7(b) and (d). Based on these experiments, it is clear that HAH feature outperforms Gabor wavelet and COOC texture features with respect to both speed and accuracy.
                        Figure 7

                        Local feature comparison using HAH, Gabor wavelet and COOC texture features. (a) Precision-recall curves of CBIR results using HAH, Gabor wavelet, and COOC texture features on peripheral blood smears. (b) Average of feature calculation times per image patch using HAH, Gabor wavelet, and COOC texture features on peripheral blood smears. (c) Precision-recall curves of CBIR results using HAH, Gabor wavelet, and COOC texture features on renal glomeruli images. (d) Average of feature calculation times per image patch using HAH, Gabor wavelet, and COOC texture features on renal glomeruli images.

                        Validation of relevance feedback

                        To evaluate the performance of the dual-similarity relevance feedback algorithm, both peripheral blood smear and multi-scale renal datasets were used. Table 1 summarizes the numbers of relevant/non-relevant images within initial top retrieved 100 images for peripheral blood smear and renal glomeruli datasets, which were labeled by two pathologists with consensus. In general, the percentages of basophils and eosinophils in a given specimen are quite small (e. g., less than 1% and 4% in our dataset as shown in Table 2). In addition, they can be accurately retrieved as we show in Table 1. Due to this reason, only neutrophils, monocytes, and lymphocytes were utilized for relevance feedback analysis. In those experiments, we applied relevance feedback on the first 100 initial retrieved image patches because this number was sufficient to retrieve all similar cases in the datasets.

                        The original query images, initial top retrieval results, and re-ranked results after relevance feedback are showed in Figures 8 and 9 for blood smear and renal datasets. In both figures, image patches with red rectangles represent the incorrect results (negative examples), and the blue ones represent the correct results (positive examples), which were re-assigned to higher ranking after relevance feedback. For the retrieval results of leukocyte image datasets, the ranking of many correct patches were increased from their initial ranking after relevance feedback. Relevance feedback corrected for 5/6 of the incorrect retrieval patches and increased the ranking for 7 patches from the lower ranking (with initial ranking between 41 and 100) in the neutrophil dataset. This procedure also amended all 10 incorrect patches, and increased ranking for 23 patches in the monocyte dataset. This procedure eliminated all 4 incorrect patches, and increased ranking for 35 patches in the lymphocyte dataset. For the renal dataset, the relevance feedback procedure successfully increased the ranking for all of the 9 correct patches of multi-scale renal dataset shown in Figure 9.

                        Ten-fold cross-validation was applied to evaluate the performance of the proposed dual-similarity relevance feedback with receiver operating characteristic (ROC) curves for both peripheral blood smear and renal datasets. The ROC curves after applying relevance feedback on the peripheral blood smear and multi-scale renal datasets are shown in Figure 10.

                        Another measures of performance for the proposed relevance feedback are the recall rate and processing speed. The relevance feedback (RF) calculation time includes feature vector dimension reduction and Adaboost classifier training. The numbers of training samples were 20, 50, and 90, and the training samples were randomly selected from the datasets. Based on Figure 11, the values of area under recall curves increased as the number of training samples increased for three leukocytes ((a) neutrophil, (b) monocyte, and (c) lymphocyte), and (d) renal glomeruli. The recall rate after RF for neutrophils (a) using 20 training samples was no better than the result before RF. This was because the original retrieval process already provided a good performance. As the value of area under recall curve before RF was already 76.902, which was much higher than the rest of cases ((b) monocyte, (c) lymphocyte, and (d) renal glomeruli). In this specific case, there was no significant improvement using RF in a small training set (e.g., 20 training samples). However, RF significantly improved the recall rate in larger training sets (e. g., 50 and 90 training samples). In general, the values of area under recall curves were significantly increased after RF with the number of training samples increased.
                        Table 1

                        Numbers of relevant/non-relevant images within top 100 initially retrieved images for peripheral blood smear and renal glomeruli datasets, which were labeled by two pathologists with an agreement


                        # of relevant images

                        # of non-relevant images
















                        Renal tissue



                        Table 2

                        Percentage of various leukocytes in adults approximately

                        Various leukocytes

                        From %

                        To %
















                        Figure 8

                        Top ranked patches before and after relevance feedback of three classes of leukocytes ((a) neutrophil, (b) monocyte, and (c) lymphocyte). Patches with red rectangles represent the incorrect results (negative examples), and blue rectangles denote the correct results (positive examples), which were re-assigned to higher rankings through the relevance feedback process. The original sizes of the images were adjusted to fit in the figure.

                        Figure 9

                        Top ranked patches before and after relevance feedback of the renal glomeruli dataset. Patches with red rectangles represent the incorrect results (negative examples), and blue rectangles represent the correct results (positive examples), which were re-assigned to higher rankings through the relevance feedback process.

                        Figure 10

                        The ROC curves of the dual-similarity relevance feedback using the peripheral blood smear image dataset (neutrophil, monocyte, and lymphocyte), and the renal glomeruli dataset.

                        Figure 11

                        The recall curves after relevance feedback (RF) and their calculation times using the peripheral blood smear image dataset ((a) neutrophil, (b) monocyte, and (c) lymphocyte)), and (d) the renal glomeruli dataset. The numbers of training samples were 20, 50 and 90.

                        Acceleration of CBIR using CometCloud

                        We conducted experiments to test the performance of CBIR using CometCloud. For HAH, we evaluated two leukocytes query images against a dataset of 925 peripheral blood smear images. In the case of CBIR using multi-scale image candidate patches, we evaluated two different renal glomeruli query images against a dataset of 32 renal images. All the experiments were repeated three times to obtain average results.

                        During the experiments, the input data were initially located on a single site, the required files were transferred as needed. However, once a file was transferred to a remote site, it was locally staged to minimize the amount of data transferred across sites, especially when multiple tasks require the same input data. To address this issue, a pull model was used where workers request tasks when they become idle. In this way, the workload was uniformly distributed across all workers to address the load imbalance.

                        To accommodate the CBIR algorithms, we federated various resources including HPC clusters and clouds. In particular, we federated a HPC cluster at Rutgers (a Dell Power Edge system with 256 cores in 8-core nodes - “Dell” hereafter), a SMP machine at Rutgers (64 cores - “Snake” hereafter), and 40 large instances from OpenStack [62] (“FutureGrid”, hereafter), which is a cloud similar to Amazon EC2. Currently we are exploiting the inherent task parallelism of the problem, which means that we can divide the algorithm into smaller sub-modules and execute each module independently. This provides a linear scalability as long as we have more tasks than computational cores.

                        Figure 12 presents a summary of the execution time of the proposed hierarchical searching algorithm using two representative peripheral blood smears and a multi-scale renal glomeruli dataset while varying the parameters, respectively. The results illustrate average values, including error bars showing their associated variabilities. Please note that the Y-axes in the sub-figures represent different scales. The figure also demonstrates the execution time of each stage and the time required to transfer the images for processing. Since the image transfer time represents a small fraction of the total execution time (i.e., from a few seconds to a 2–3 minutes depending on the configuration), in our current implementation we copy the images sequentially from a central repository. The execution time varies depending on the algorithm we used, the query and dataset images, and the configuration (e.g., 90% overlapping takes longer than 50% overlapping). The fraction of time spent on each stage of the hierarchical searching is shown in Figure 12.

                        Figure 13 compares the execution time of different configurations using a single system and federated cyber-infrastructure. We observe an average acceleration of 70-fold with a maximum of 96-fold. This is achieved by elastically using multiple resources as discussed below. Figure 14 shows the contribution of the FutureGrid cloud to the execution of the multi-scale algorithm. Cloud resources significantly accelerate the execution of the algorithm. During stages with lower parallelism (e.g., last minutes of the execution), computation can be performed using local HPC resources and cloud resources can be released to reduce operational costs.

                        The variability of the execution time of different tasks is shown in Figures 15 and 16. Figure 15 shows the average task execution time and variability using different configurations. The variability of task execution time is heterogeneous and depends on the configurations and the machine. In general, the longer the execution takes, the larger the variability. Figure 16 shows that the execution time of individual task is relatively heterogeneous. It also demonstrates that the distribution of tasks among different federated resources depends on the number of cores available in each platform (e.g., one of the cores, snake, runs only a few tasks). The results show that the parallelization of CBIR at the image level can dramatically reduce the overall computational time.
                        Figure 12

                        The execution time of hierarchical searching process using (a) peripheral blood smear dataset, and (b) renal glomeruli dataset, with different combinations of number of the HAH rings and the percentage of overlapping.

                        Figure 13

                        The execution time of sequential and federated infrastructure using peripheral blood smear dataset and renal glomeruli dataset with different combinations of number of the rings and the percentage of overlapping. Here the Y-axis is in a logarithmic scale.

                        Figure 14

                        The number of completed tasks over time when testing the CBIR algorithm using the renal glomeruli dataset. The area under “FutureGrid” represents the contribution from the cloud resources.

                        Figure 15

                        The average task execution time per platform using the peripheral blood smear and renal glomeruli datasets with different combinations of number of the rings and the percentage of overlapping. Here “FutureGrid” is abbreviated as “fg”.

                        Figure 16

                        The execution time per task using (a) the peripheral blood smear dataset and (b) the renal glomeruli dataset.


                        In this paper, we present a set of newly developed CBIR algorithms and demonstrate its application on two different pathology applications, which are regularly evaluated in the practice of pathology. The experimental results suggest that the proposed CBIR algorithm using sequential HAH searching follows a progression which parallels to the same logical steps as ever invoked when physicians review digital pathology images. During the review process, the pathologist typically begins by first identifying gross locations of potential regions of interest (coarse searching in the proposed algorithm) before executing the more refined stages (fine searching in the proposed algorithm) to examine the detailed morphometric characteristics.

                        For the peripheral blood smear study, we tested performance using a range of different leukocytes and experimentally showed the reliable performance of the CBIR algorithm. The success of the proposed CBIR algorithm in identifying neutrophils suggests further exploration of the HAH feature in detecting abnormal or hypersegmented neutrophils, which are indicators of megaloblastic anemia and potential risk of gastric cancer. Similarly, a pathologist’s assessment of normal vs. diseased glomeruli in renal biopsies is often used as an indicator of overall kidney health, such as, the determination of graft function from pre-transplantation biopsies [55]. Assisted by the proposed CBIR algorithm, physicians and researchers can quickly review a digital biopsy to evaluate the proportion of ischemic or necrotic glomeruli within a given field to quickly assess whether an incoming specimen is suitable for transplantation or not. This type of review can have multiple applications, such as, determining whether a rejection of the organ might occur by identifying areas of focal and segmental glomerulosclerosis [63]. Currently, our algorithm requires some external feedback to optimize the search. We are exploring different ways of automatizing this process by applying machine learning techniques. On the other hand, although the proposed hierarchical searching has significantly improved the retrieval speed, it is still a computational demanding procedure. Therefore, we are exploring new ways of exploiting parallelism to speed-up this process.

                        We present a generalizable cloud-enabled CBIR algorithm that can be extended to a wide variety of applications. Because of the computational requirements needed for retrieving whole-slide scanned images, we explore the use of federated high performance computing (HPC) cyber-infrastructures and clouds using CometCloud. Comparative results of HPC versus standard computation time demonstrate that the CBIR process can be dramatically accelerated, from weeks to minutes, making real-time clinical practice feasible. Moreover, the proposed framework hides infrastructure and deployment details and offers end-users the CBIR functionality in a readily accessible manner. We are currently working on improving the utilization of resources by exploit the particular capabilities and capacities of each heterogeneous resource, e.g., switching between the usage of the original CBIR implementation in MATLAB (The MathWorks, Natick, MA) when licenses are available or a parallel implementation using graphic processing unit (GPU) and many-core architectures in cases where resources with accelerators are available.



                        This research was funded, in part, by grants from the National Institutes of Health through contract 5R01CA156386-09 and 5R01CA161375-03 from the National Cancer Institute; and contracts 2R01LM009239-05 and 5R01LM011119-03 from the National Library of Medicine. Additional support was provided by a gift from the IBM International Foundation. This work is sponsored by the National Science Foundation (NSF) Office of Industrial Innovation and Partnerships Industry/University Cooperative Research Center (I/UCRC) Program under award 0758596. This work used resources of FutureGrid, which is supported in part by the NSF under Grant No. 0910812. This work is supported in part by the NSF under grants ACI-1339036, IIP-0758566, DMS-0835436, CNS-1305375, ACI-1310283, and by IBM via OCR and Faculty awards. CometCloud was developed as part of the NSF Cloud and Autonomic Computing Center at Rutgers University. The project is also partially supported by the National Center for Research Resources and the National Center for Advancing Translational Sciences, National Institutes of Health, through Grant UL1TR000117 (or TL1 TR000115 or KL2 TR000116). The content is solely the responsibility of the authors and does not necessarily represent the official views of the NIH. Finally we would like to thank Dr. Wenjin Chen for preparing the data for these studies.

                        Authors’ Affiliations

                        Department of Pathology and Laboratory Medicine, Rutger Robert Wood Johnson Medical School
                        Center for Biomedical Imaging and Informatics, Rutgers Cancer Institute of New Jersey
                        Rutgers Discovery Informatics Institute and NSF Cloud and Autonomic Computing Center, Department of Electrical and Computer Engineering, Rutgers University
                        Department of Radiology, Rutgers - Robert Wood Johnson Medical School
                        Division of Biomedical Informatics, Department of Biostatistics, Department of Computer Science, University of Kentucky


                        1. Gudivada VN, Raghavan VV: Content-based image retrieval systems . Computer 1995,28(9):18–22. 10.1109/2.410145View Article
                        2. Flickener M, Sawhney H, Niblack W, Ashley J, Huang Q, Dom B, Gorkani M, Hafner J, Lee D, Petkovic D, Steele D, Yanker P: Query by image and video content: the QBIC system . Computer 1995,28(9):23–32. 10.1109/2.410146View Article
                        3. Smith JR, Chang SF: Visualseek: a fully automated content-based image query system . Proceedings of the Fourth ACM International conference on Multimedia. 1996, 87–98.View Article
                        4. Tagare HD, Jaffe CC, Duncan J: Medical image databases: a content-based retrieval approach . J Am Med Inform Assoc 1997, 4:184–198. 10.1136/jamia.1997.0040184PubMedPubMed CentralView Article
                        5. Smeulders AWM, Worring M, Santini S, Gupta A, Jainh R: Content-based image retrieval at the end of the early years . IEEE Trans Pattern Anal Mach Intell 2000, 22:1349–1380. 10.1109/34.895972View Article
                        6. Wang J, Li J, Wiederhold G: Simplicity: semantics-sensitive integrated matching for picture libraries . IEEE Trans Pattern Anal Mach Intell 2001, 23:947–963. 10.1109/34.955109View Article
                        7. Chen Y, Wang J: A region-based fuzzy feature matching approach to content-based image retrieval . IEEE Trans Pattern Anal Mach Intell 2002, 24:1252–1267. 10.1109/TPAMI.2002.1033216View Article
                        8. Chang E, Goh K, Sychay G, Wu G: Cbsa: content-based soft annotation for multimodal image retrieval using bayes point machines . IEEE Trans Circ Syst Video Tech 2003, 13:26–38. 10.1109/TCSVT.2002.808079View Article
                        9. Zheng L, Wetzel AW, Gilbertson J, Becich MJ: Design and analysis of a content-based pathology image retrieval system . IEEE Trans Inf Technol Biomed 2003,7(4):245–255.
                        10. Muller H, Michoux N, Bandon D, Geissbuhler A: A review of content-based image retrieval systems in medical applications - clinical benefits and future directions . Int J Med Inform 2004, 73:1–23. 10.1016/j.ijmedinf.2003.11.024PubMedView Article
                        11. Lehmann TM, Guld MO, Deselaers T, Keysers D, Schubert H, Spitzer K, Ney H, Wein BB: Automatic categorization of medical images for content-based retrieval and data mining . Comput Med Imaging Graph 2005, 29:143–155. 10.1016/j.compmedimag.2004.09.010PubMedView Article
                        12. Lam M, Disney T, Pham M, Raicu D, Furst J, Susomboon R: Content-based image retrieval for pulmonary computed tomography nodule images . Proc. SPIE Medical Imaging 2007, 6516:65160N-1–65160N-12. 10.1117/12.710297View Article
                        13. Rahman MM, Antani SK, Thoma GR: A learning-based similarity fusion and filtering approach for biomedical image retrieval using SVM classification and relevance feedback . IEEE Trans Inf Technol Biomed 2011,15(4):640–646.PubMedView Article
                        14. Thies C, Malik A, Keysers D, Kohnen M, Fischer B, Lehmann TM: Hierarchical feature clustering for content-based retrieval in medical image databases . Proc SPIE 2003, 5032:598–608. 10.1117/12.481393View Article
                        15. El-Naqa I, Yang Y, Galatsanos NP, Nishikawa RM, Wernick MN: A similarity learning approach to content-based image retrieval: application to digital mammography . IEEE Trans Med Imaging 2004, 23:1233–1244. 10.1109/TMI.2004.834601PubMedView Article
                        16. Akakin HC, Gurcan MN: Content-based microscopic image retrieval system for multi-image queries . IEEE Trans Inf Technol Biomed 2012, 16:758–769.PubMedPubMed CentralView Article
                        17. Zhang Q, Izquierdo E: Histology image retrieval in optimized multifeature spaces . IEEE J Biomed Health Informatics 2013, 17:240–249.View Article
                        18. Tang HL, Hanka R, Ip HH: Histological image retrieval based on semantic content analysis . IEEE Trans Inf Technol Biomed 2003, 7:26–36. 10.1109/TITB.2003.808500PubMedView Article
                        19. Schmidt-Saugenon P, Guillod J, Thiran JP: Towards a computer-aided diagnosis system for pigmented skin lesions . Comput Med Imag Graph 2003, 27:65–78. 10.1016/S0895-6111(02)00048-4View Article
                        20. Sbober A, Eccher C, Blanzieri E, Bauer P, Cristifolini M, Zumiani G, Forti S: A multiple classifier system for early melanoma diagnosis . Artifical Intel Med 2003, 27:29–44. 10.1016/S0933-3657(02)00087-8View Article
                        21. Meyer F: Automatic screening of cytological specimens . Comput Vis Graph Image Process 1986, 35:356–369. 10.1016/0734-189X(86)90005-8View Article
                        22. Mattie ME, L. Staib ES, Tagare HD, Duncan J, Miller PL: Content-based cell image retrieval using automated feature extraction . J Am Med Informatics Assoc 2000, 7:404–415. 10.1136/jamia.2000.0070404View Article
                        23. Beretti S, Bimbo AD, Pala P: Content-based retrieval of 3D cellular structures . 2001 IEEE International Conference on Multimedia and Expo, ICME’01. 2001, 226–226.
                        24. Pentland A, Picard RW, Sclaroff S: Photobook: tools for content-based manipulation of image databases . Int J Comput Vis 1996, 18:233–245. 10.1007/BF00123143View Article
                        25. Lehman TM, Guld MO, Thies C, Fischer B, Spitzer K, Keysers D, Ney H, Kohnen M, Schubert H, Wein BB: Content-based image retrieval in medical applications . Methods Inf Med 2004, 4:354–360.
                        26. Cox IJ, Miller ML, Omohundro SM, Yianilos PN: Target testing and the PicHunter Bayesian multimedia retrieval system . Digital Libraries, Proceedings of the Third Forum on Research and Technology Advances in IEEE. 1996, 66–75.View Article
                        27. Carson C, Belongies S, Greenspan H, Malik J: Region-based image querying . Proceedings of IEEE Workshop on Content-Based Access of Image and Video Libraries. 1997, 42–49.View Article
                        28. LeBozec C, Jaulent MC, Zapletal E, Degoulet P: Unified modeling language and design of a case-based retrieval system in medical imaging . Proceedings of the Annual Symposium of the American Society for Medical Informatics. 1998, 887–891.
                        29. Bui AAT, Taira RK, Dionision JDN, Aberle DR, El-Saden S, Kangarloo H: Evidence-based radiology . Acad Radiol 2002, 9:662–669. 10.1016/S1076-6332(03)80311-2PubMedView Article
                        30. Kong J, Cooper LAD, Wang F, Gutman DA, Gao J, Chisolm C, Sharma A, Pan T, Meir EGV, Kurc TM, Moreno CS, Saltz JH, Brat DJ: Integrative, multimodal analysis of glioblastoma using TCGA molecular data, pathology images, and clinical outcomes . IEEE Trans Biomed Eng 2011, 58:3469–3474.PubMedPubMed CentralView Article
                        31. Cavallaro A, Graf F, Kriegel H, Schubert M, Thoma M: Region of interest queries in CT scans . In Proceedings of the 12th International Conference On Advances in Spatial and Temporal Databases.. Springer Berlin Heidelberg; 2011:56–73.View Article
                        32. Naik J, Doyle S, Basavanhally A, Ganesan S, Feldman MD, Tomaszewski JE, Madabhushi A: A boosted distance metric: application to content based image retrieval and classification of digitized histopathology . Proc SPIE Med Imaging 2009, 7260:1–4.
                        33. Vu K, Hua KA, Tavanapong W: Image retrieval based on regions of interest . IEEE Trans Knowl Data Eng 2003,15(4):1045–1049. 10.1109/TKDE.2003.1209021View Article
                        34. Ke Y, Sukthankar R, Huston L: Efficient near-duplicate detection and subimage retrieval . ACM Multimedia 2004, 4:869–876.
                        35. Simonyan K, Zisserman A, Criminisi A: Immediate structured visual search for medical images . In Medical Image Computing and Computer-Assisted Intervention-MICCAI2011.. Springer Berlin Heidelberg; 2011:288–296.
                        36. Pfund T, Marchang-Maillet S: Dynamic multimedia annotation tool . Internet Imaging III 2002, 4672:206–224.
                        37. Squire DM, Muller W, Muller H, Pun T: Content-based query of image databases: inspirations from text retrieval . Pattern Recogn Lett 2000, 21:1193–1198. 10.1016/S0167-8655(00)00081-7View Article
                        38. Sclaroff S, Taycher L, Cascia ML: Imagerover: a content-based browser for the world wide web . Proceedings of IEEE Workshop on Content-Based Access of Image and Video Libraries. 1997, 2–9.View Article
                        39. Ortega M, Rui Y, Chakrabarti K, Porkaew K, Mehrotra S, Huang TS: Supporting ranked boolean similarity queries in MARS . IEEE Trans. Knowledge Data Eng 1998, 10:905–925. 10.1109/69.738357View Article
                        40. Kuo WJ, Chang RF, Lee CC, Moon WK, Chen DR: Retrieval technique for the diagnosis of solid breast tumors on sonogram . Ultrasound Med Biol 2002, 28:903–909. 10.1016/S0301-5629(02)00541-0PubMedView Article
                        41. Ishikawa Y, RSubramanya R, Faloutsos C: MindReader: querying database through multiple examples . Comput Sci Dep 1998, 551:218–227.
                        42. Porkaew K, Chakrabarti L, Mehrota S: Query refinement for multimedia similarity retrieval in mars . Proceedings of Seventh ACM International Conference on Multimedia (Part 1) 1999, 235–238.View Article
                        43. Liu D, Hua K, Vu K, Yu N: Fast query point movement techniques for large cbir systems . IEEE Trans Knowledge Data Eng 2009, 21:729–743.View Article
                        44. Tieu K, Viola P: Boosting image retrieval . IEEE Comput Soc Conf Comput Vis Pattern Recogn 2000, 1:228–235.
                        45. Tao D, Tang X: Random sampling based SVM for relevance feedback image retrieval . IEEE Comput Soc Conf Comput Vis Pattern Recogn 2004, 2:647–652.
                        46. Yang L, Jin R, Mummert L, Sukthankar R, Goode A, Zheng B, Hoi S, Satyanarayanan M: A boosting framework for visuality-preserving distance metric learning and its application to medical image retrieval . IEEE Trans Pattern Anal Mach Intel 2010, 32:30–44.View Article
                        47. Su J, Huang W, Yu P, Tseng V: Efficient relevance feedback for content-based image retrieval by mining user navigation patterns . IEEE Trans Knowledge Data Eng 2011, 23:360–372.View Article
                        48. The Cloud and Autonomic Computing Center at Rutgers University [http://​nsfcac.​rutgers.​edu/​CometCloud/​]
                        49. Kim H, el-Khamra Y, Jha S, Rodero I, Parashar M: Autonomic management of application workflow on hybrid computing infrastructure . Sci Program 2011, 19:75–89.
                        50. Kim H, Chaudhari S, Parashar M, Marty C: Online risk analysis on the cloud . 9th IEEE/ACM International Symposium on Cluster Computing and the Grid, Cluster Computing and the Grid, CCGRID’09. 2009, 484–489.
                        51. Kim H, el-Khamra Y, Jha S, Parashar M: An autonomic approach to integrated hpc grid and cloud usage . 5th IEEE International Conference on e-Science, e-Science’09. 2009, 366–373.
                        52. Kim H, Parashar M, Foran DJ, Yang L: Investigating the use of cloudbursts for high-throughput medical image registration . 10th IEEE/ACM International Conference on Grid Computing, Grid Computing. 2009, 34–41.
                        53. Parashar M, AbdelBaky M, Rodero I, Devarakonda A: Cloud paradigms and practices for computational and data-enabled science and engineering . Comput Sci Eng 2013,15(4):10–18.View Article
                        54. Yang L, Tuzel O, Chen W, Meer P, Salaru G, Goodell LA, Foran DJ: Pathminer: a web-based tool for computer-assisted diagnostics in pathology . IEEE Trans Inf Technol Biomed 2009, 13:291–299.PubMedPubMed CentralView Article
                        55. Martul VE, Barreiro VA: Importance of kidney biopsy in graft selection . Transplant Proc 2003, 35:1658–1660. 10.1016/S0041-1345(03)00573-6View Article
                        56. Qi X, Xing F, Foran DJ, Yang L: Robust segmentation of overlapping cells in histopathology specimens using parallel seed detection and repulsive level set . IEEE Trans Biomed Eng 2012, 59:754–765.PubMedView Article
                        57. Comaniciu D, Meer P: Mean shift: a robust approach toward feature space analysis . IEEE Trans Pattern Anal Mach Intell 2002, 24:603–619. 10.1109/34.1000236View Article
                        58. Daugman J: Complete discrete 2-D Gabor transforms by neural networks for image analysis and compression . IEEE Trans Acous Speech Signal Process 1988, 36:1169–1179. 10.1109/29.1644View Article
                        59. Haralick RM, Shanmugam K, Dinstein I: Textural features for image classification . IEEE Trans Syst Man Cybern 1973, 3:610–621.View Article
                        60. Haralick RM: Statistical and structural approaches to texture . Proc SPIE 1979, 67:786–804.
                        61. Ratsch G, Onoda T, Muller KR: Soft margins for Adaboost . Mach Learn 2001, 42:287–320. 10.1023/A:1007618119488View Article
                        62. OpenStack: Open source software for building private and public clouds . [http://​www.​openstack.​org]
                        63. Shimizu A, Higo S, Fujita E, Mii A, Kaneko T: focal segmental glomerulosclerosis after renal transplantation . Clin Transplant 2011, 25:6–14.PubMedView Article


                        © Qi et al.; licensee BioMed Central Ltd. 2014

                        This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://​creativecommons.​org/​licenses/​by/​2.​0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://​creativecommons.​org/​publicdomain/​zero/​1.​0/​) applies to the data made available in this article, unless otherwise stated.