Team:TU Eindhoven/Modeling/Bacterial Cell Counter

From 2014.igem.org

Revision as of 01:05, 18 October 2014 by Mijaafa (Talk | contribs)

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 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 for this image are 9 false positives and 11 false negatives on a total of 93 cell in 126 droplets. Because some droplets had multiple false positives or false negatives the next step to improve this number is selecting to droplets better.
Click here to download the code

Figure 4. Histogram of the example image
iGEM Team TU Eindhoven 2014