Method for determining the positions of a plurality of objects in a digital image

Patent metadata

Field	Value
Publication number	US-10733707-B2
Application number	US-201616060887-A
Country	US
Kind code	B2
Filing date	Dec 9, 2016
Priority date	Dec 10, 2015
Publication date	Aug 4, 2020
Grant date	Aug 4, 2020

Abstract

The present invention relates to a method and a system for determining the positions of a plurality of objects in a digital image by discriminating true positive positions of the plurality of objects from false positive candidate positions of the plurality of objects. In particular, the invention relates to a method for determining the positions of a plurality of objects in a digital image by discriminating true positive positions of the plurality of objects from false positive candidate positions of the plurality of objects, the plurality of objects being configured to receive molecules comprising genetic information.

First claim

The invention claimed is:

1. A computer implemented method for determining positions of a plurality of objects in a digital image, the plurality of objects being configured to receive molecules comprising genetic information, the method comprising:
determining a plurality of pixels in the digital image, each pixel having a local maximum intensity value with respect to intensity values of a plurality of neighboring pixels of the respective pixel;
determining a variance of an intensity value of each pixel having a respective local maxim urn intensity value with respect to the intensity value of the respective pixel and the intensity value of each of the plurality of the respective neighboring pixels;
weighting the intensity value of each of the plurality of pixels having a respective local maximum intensity value by the respective determined variance of the intensity value;
creating a histogram comprising the respective weighted intensity of each pixel of the plurality of pixels each having a respective local maximum intensity value;
determining a plurality of local minimums and a global maximum of the histogram;
selecting a local minimum of the plurality of local minimums as a threshold value for discriminating true positive positions of the plurality of objects from false positive positions of the plurality of objects, wherein the local minimum is selected under consideration of a position of the selected local minimum of the histogram relative to a position of the global maximum; and
using the threshold value for discriminating the true positive positions of the plurality of objects from false-positive candidate positions of the plurality of objects.

2. The method according to claim 1 , further comprising:
determining a sub-pixel position of each object by parabola fitting for the determined pixels.

3. The method according to claim 1 , further comprising:
removing weighted pixels having a weighted intensity value below a first threshold or above a second threshold.

4. The method according to claim 1 , further comprising:
scaling weighted pixels with a logarithmic scale.

5. The method according to claim 1 , wherein at least one object of the plurality of objects is a bead configured to receive DNA or RNA.

6. The method according to claim 1 , further comprising:
removing at least one pixel of the digital image having an intensity below a first threshold or above a second threshold.

7. The method according to claim 1 , further comprising:
removing uneven illumination in the digital image.

8. The method according to claim 7 , wherein the removing uneven illumination comprises:
applying a low pass filter to a copy of the digital image;
dividing the digital image by the low pass filtered copy of the digital image;
dividing the digital image into a plurality of N tiles;
determining a first threshold intensity value and a second threshold intensity value for each tile of the plurality of N tiles;
creating a first matrix comprising the first threshold intensity values for the N tiles, and a second matrix comprising the second threshold intensity values for the N tiles;
applying a median filter to the first matrix; and a median filter to the second matrix;
upscaling the first matrix and the second matrix, the upscaled first matrix and the upscaled second matrix having an original size of the digital image; and
normalizing each intensity value of each pixel of the digital image under consideration of the corresponding value of the upscaled first matrix and of the corresponding value of the upscaled second matrix.

9. The method according to claim 8 , further comprising:
applying a low pass filter to the first matrix, and a low pass filter to the second matrix.

10. The method according to claim 8 , wherein the upscaling is performed using bicubic interpolation.

11. The method according to claim 8 , wherein the normalizing each intensity value of each pixel of the digital image is performed by a linear transformation.

12. A non-transitory computer readable media having instructions stored thereon that, when executed by at least one computing device, cause the at least one computing device to perform operations comprising:
determining a plurality of pixels in a digital image; each pixel having a local maximum intensity value with respect to intensity values of a plurality of neighboring pixels of the respective pixel;
determining a variance of an intensity value of each pixel having a respective local maximum intensity value with respect to the intensity value of the respective pixel and the intensity value of each of the plurality of the respective neighboring pixels;
weighting the intensity value of each of the plurality of pixels having a respective local maximum intensity value by the respective determined variance of the intensity value;
creating a histogram comprising the respective weighted intensity of each pixel of the plurality of pixels each having a respective local maximum intensity value;
determining a plurality of local minimums and a global maximum of the histogram;
selecting a local minimum of the plurality of local minimums as a threshold value for discriminating true positive positions of the plurality of objects from false positive positions of the plurality of objects, wherein the local minimum is selected under consideration of a position of the selected local minimum of the histogram relative to a position of the global maximum; and
using the threshold value for discriminating the true positive positions of the plurality of objects from false-positive candidate positions of the plurality of objects.

13. The non-transitory computer readable media according to claim 12 , wherein the operations further comprise:
removing uneven illumination in the digital image.

14. The non-transitory computer readable media according to claim 13 , wherein the removing uneven illumination comprises:
applying a low pass filter to a copy of the digital image;
dividing the digital image by the low pass filtered copy of the digital image;
dividing the digital image into a plurality of N tiles;
determining a first threshold intensity value and a second threshold intensity value for each tile of the plurality of N tiles;
creating a first matrix comprising the first threshold intensity values for the N tiles, and a second matrix comprising the second threshold intensity values for the N tiles;
applying a median filter to the first matrix, and a median filter to the second matrix;
upscaling the first matrix the second matrix, the upscaled first matrix and the upscaled second matrix having an original size of the digital image; and
normalizing each intensity value of each pixel of the digital image under consideration of the corresponding value of the upscaled first matrix and of the corresponding value of the upscaled second matrix.

15. The non-transitory computer readable media according to claim 14 , wherein the operations further comprise:
applying a low pass filter to the first matrix, and a low pass filter to the second matrix.

16. A system for determining positions of a plurality of objects in a digital image, the plurality of objects being configured to receive molecules comprising genetic information, wherein the system comprises:
a memory; and
at least one processor coupled to the memory and configured to:
determine a plurality of pixels in the digital image, each pixel having a local maximum intensity value with respect to intensity values of a plurality of neighboring pixels of the respective pixel;
determine a variance of an intensity value o

Assignees

• Qiagen Gmbh

Inventors

• Unterer Thomas
• Schuster Guido
• Lohel Maiko
• Labitzke Björn
• Zerfass Thorsten

Classifications

• G06T7/0012 Primary

  Biomedical image inspection · CPC title

• G06T2207/10064

  Fluorescence image · CPC title

• G06T2207/30072

  Microarray; Biochip, DNA array; Well plate · CPC title

• G06T2207/20021

  Dividing image into blocks, subimages or windows · CPC title

• G06T7/11

  Region-based segmentation · CPC title