Team:TU Eindhoven/Modeling/Bacterial Cell Counter

From 2014.igem.org

iGEM Team TU Eindhoven 2014

iGEM Team TU Eindhoven 2014

Figure 1. The used sample image

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.

Figure 2. The results after EdgeDetect

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.

Figure 3. The results of clustering shown as blue circles, layered onto the original image

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 corresponding output of the program from the sample image shown in Figure 1 can be found in Figure 4. The analysis of the sample image included 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. In the future, droplet selection has to be optimized in order to increase accuracy of the program.
Click here to download the code

Figure 4. The output of the program is a histogram showing how many droplets contained a certain number of cells
iGEM Team TU Eindhoven 2014