Team:TU Eindhoven/Modeling/Bacterial Cell Counter
From 2014.igem.org
Bacterial Cell Counter
Once photos of the microfluidics results have been made, the droplets and the cells inside the droplets have to be counted. In order to have reliable results, this has to be done on a large scale. A computer program has been designed to perform the counting on large numbers of photos. To give a better visualization of the process, a step by step analysis of the image shown in Figure 1 is given below.
The program has the following steps: firstly, to find and count the droplets, then to find and count the cells inside of the droplets and finally to create a histogram of the results. By adding up the histograms of multiple images one get a impression of how the cells are divided over the droplets.
Droplet Detection
To detect the cells EdgeDetect, a function of Mathematica using gradient methods, is used. It is then followed by a dilation in order to make the edges clearer. The results of these actions can be seen in Figure 2.
This results in the clustering of black pixels. The program then looks for the clusters bigger than 2000 pixels and smaller than 6000 pixels. In Figure 3, these clusters are visualized with blue circles.
Parameters for Droplets Detection
The parameters of this function are chosen in a way that prioritizes decrease in false positives over the decrease in false negatives. This is because a false positive means a non-existing droplet and thus false data, while false negative only lowers the sample size, which can be increased by analyzing more images.
Bacterial Cell Detection
The program is currently at the point where the next step is to count the number of cells inside the droplets. Onto the clusters from Figure 2, an Erosion function with value 1 is performed. This makes sure that cells connected to the borders of the droplet are disconnected from this border. The program deletes border components, so that it only finds cells that are not connected to the outside of the droplet. This cancels out outside noise. The program uses EnclosingComponentCount instead of count to get the results, because a lot of cells became empty circles.
The results of the sample image shown Figure 1 are 9 false positives and 11 false negatives, on a total of 93 cell in 126 droplets. Some droplets had multiple false positives or false negatives. The next step optimize droplet selection in order to improve these results.
Click here to download the code