A kind of circulating tumor cell detection method
Technical field
The present invention relates to a kind of circulating tumor cell detection methods.
Background technique
Cancer is the second largest assailant of current human beings worldwide's death, is only second to traffic accident.Many experts and scholars are in the world
In the molding reason and method of discrimination early period of research cancer.If cancer can find in time in early stage and take corresponding treatment
Measure is that available effective control reduces dead probability.A kind of current novel cancer cell method of discrimination early period is
It is widely used in --- circulating tumor cell detects (Circulating Tumer Cell, CTC).Usually differentiating process
In, the blood of patient can be extracted and manufacture into sample, put and observe under the microscope.These blood samples are usually by three kinds of dyestuffs
It is dyed, and is observed with corresponding three channels of fluorescence microscope.Electronic fluorescence microscope is commonly to acquire
Tool can guarantee that the image of acquisition covers entire sample using image mosaic technology and Autofocus Technology.But due to lacking
Few effective digital method of discrimination, staff usually will do it artificial detection.But the efficiency of artificial detection is extremely low, and quasi-
True rate is not high.Thus leverage the accuracy and practicability of this kind of method detection.
Summary of the invention
It is a kind of accurately swollen with efficient circulation the technical problem to be solved by the present invention is to be provided for the above-mentioned prior art
Oncocyte detection method.
The present invention solves technical solution used by above-mentioned technical problem:A kind of circulating tumor cell detection method,
It is characterized in that:Include the following steps
Step (1), the fluorescent image that three different channels are acquired by fluorescence microscope, three different channels are corresponding glimmering
Light image is respectively green fluorescence image, red fluorescence images and blue-fluorescence image;
Step (2), by step (1) acquisition green fluorescence image, red fluorescence images and blue-fluorescence image respectively into
Row initialization process:
Step (2-1) obtains green fluorescence image and red fluorescence images and blue first with maximum variance between clusters
Initialization area in fluorescent image where each cell;
Step (2-2), the initialization to all acquisitions in green fluorescence image, red fluorescence images and blue-fluorescence image
Region carries out median filtering, and the convolution kernel size of filtering is N*N, and the value of N is 2~5, respectively obtains filtered green fluorescence
Image and red fluorescence images and blue-fluorescence image;
Step (2-3), using watershed algorithm to filtered green fluorescence image, red fluorescence images and blue-fluorescence
Image is split processing, then records the position of each cell split, area and average brightness, thus
Respectively obtain the position that the cell come out is divided in green fluorescence image, red fluorescence images and blue-fluorescence image, area
And average brightness;
Step (3) will be divided in the resulting green fluorescence image of step (2), red fluorescence images and blue-fluorescence image
The average brightness for all cells for cutting out carries out then finding out position in preceding 10% cell respectively from bright to dark sequence respectively
In the average brightness of the cell of lower critical value, as positive cell brightness discrimination threshold, find out in rear 10% cell positioned at upper
The average brightness of critical value cell, as negative cells brightness discrimination threshold;The judgment criteria of positive cell is the flat of the cell
Equal brightness is more than or equal to positive cell brightness discrimination threshold;The judgment criteria of negative cells be the cell average brightness be less than etc.
In negative cells brightness discrimination threshold;
Step (4), according to the judgment principle of step (3), find the negative cells in green fluorescence image, record green is glimmering
The cell position of all non-negative cells in light image, and rejected;Simultaneously according to institute in the green fluorescence image of record
The cell position for the non-negative cells having rejects the cell of the same position in red fluorescence images and blue-fluorescence image, and
Update the cell distribution in red fluorescence images and blue-fluorescence image;Then, it also according to the judgment principle of step (3), looks for
To the positive cell in blue-fluorescence image, the cell position of all non-positive cells in blue-fluorescence image is recorded, and is given
It rejects, while rejecting the cell of the same position in green fluorescence image and red fluorescence images, and update green fluorescence image
With the cell distribution in red fluorescence images;Finally, being found in red fluorescence images also according to the judgment principle of step (3)
Positive cell, record the position of all non-positive cells in red fluorescence images, and rejected, while rejecting glimmering in green
The cell of same position in light image and blue-fluorescence image, and update the cell in green fluorescence image and blue-fluorescence image
Distribution;Green fluorescence image, red fluorescence images and the blue-fluorescence image of interference cell are eliminated to obtain three width;
Step (5), by step (4) it is obtained eliminate the interference green fluorescence image of cell, red fluorescence images and
Blue-fluorescence image is synthesized, and the standard of synthesis has:
A, the cell existed simultaneously in green fluorescence image, red fluorescence images and blue-fluorescence image is retained;
B, it is retained in the cell that the cell area in blue-fluorescence image is greater than the cell area in red fluorescence images;
The cell for meeting a and b condition simultaneously is retained, these cells are output to as the circulating tumor cell detected
In final examining report.
Compared with the prior art, the advantages of the present invention are as follows:The present invention is recycled by the method for Digital Image Processing
Tumour cell is detected, and to the processing of cell level analysis, is able to solve what traditional micro-image was generated in pixel layer face treatment
A large amount of cavities and noise problem can be removed effectively due to raw brought erroneous detection of making an uproar, while respectively to the figure in three channels
Judgment criteria is defined as progress different disposal can be convenient user, to improve the accuracy of detection, and then improves the effect of detection
Rate.
Detailed description of the invention
Fig. 1 is circulating tumor cell detection method flow chart in the embodiment of the present invention.
Specific embodiment
The present invention will be described in further detail below with reference to the embodiments of the drawings.
The present invention provides a kind of circulating tumor cell detection methods comprising following steps
Step (1), the fluorescent image that three different channels are acquired by fluorescence microscope, three different channels are corresponding glimmering
Light image is respectively green fluorescence image, red fluorescence images and blue-fluorescence image;In the present embodiment, three differences of acquisition
The colour filter module for the fluorescence microscope that the fluorescent image in channel utilizes is different, and corresponding excitation wavelength is 475nm,
555nm and 385nm, and the whole resulting image resolution ratio of lid fragmentation is scanned up to 11500*36100 using 20x scanning objective;
Step (2), by step (1) acquisition green fluorescence image, red fluorescence images and blue-fluorescence image respectively into
Row initialization process:
Step (2-1) obtains green fluorescence image and red fluorescence images and blue first with maximum variance between clusters
Initialization area in fluorescent image where each cell, maximum variance between clusters are also referred to as big law, and this method is routine side
Method;
Step (2-2), the initialization to all acquisitions in green fluorescence image, red fluorescence images and blue-fluorescence image
Region carries out median filtering, and the convolution kernel size of filtering is N*N, and the value of N is 2~5, preferably 3, respectively obtains filtered green
Color fluorescent image and red fluorescence images and blue-fluorescence image;
Step (2-3), using watershed algorithm to filtered green fluorescence image, red fluorescence images and blue-fluorescence
Image is split processing, then records the position of each cell split, area and average brightness, thus
Respectively obtain the position that the cell come out is divided in green fluorescence image, red fluorescence images and blue-fluorescence image, area
And average brightness;
Step (3) will be divided in the resulting green fluorescence image of step (2), red fluorescence images and blue-fluorescence image
The average brightness for all cells for cutting out carries out then finding out position in preceding 10% cell respectively from bright to dark sequence respectively
In the average brightness of the cell of lower critical value, as positive cell brightness discrimination threshold, find out in rear 10% cell positioned at upper
The average brightness of critical value cell, as negative cells brightness discrimination threshold;The judgment criteria of positive cell is the flat of the cell
Equal brightness is more than or equal to positive cell brightness discrimination threshold;The judgment criteria of negative cells be the cell average brightness be less than etc.
In negative cells brightness discrimination threshold;
Step (4), according to the judgment principle of step (3), find the negative cells in green fluorescence image, record green is glimmering
The cell position of all non-negative cells in light image, and rejected;Simultaneously according to institute in the green fluorescence image of record
The cell position for the non-negative cells having rejects the cell of the same position in red fluorescence images and blue-fluorescence image, and
Update the cell distribution in red fluorescence images and blue-fluorescence image;Then, it also according to the judgment principle of step (3), looks for
To the positive cell in blue-fluorescence image, the cell position of all non-positive cells in blue-fluorescence image is recorded, and is given
It rejects, while rejecting the cell of the same position in green fluorescence image and red fluorescence images, and update green fluorescence image
With the cell distribution in red fluorescence images;Finally, being found in red fluorescence images also according to the judgment principle of step (3)
Positive cell, record the position of all non-positive cells in red fluorescence images, and rejected, while rejecting glimmering in green
The cell of same position in light image and blue-fluorescence image, and update the cell in green fluorescence image and blue-fluorescence image
Distribution;Green fluorescence image, red fluorescence images and the blue-fluorescence image of interference cell are eliminated to obtain three width;
Step (5), by step (4) it is obtained eliminate the interference green fluorescence image of cell, red fluorescence images and
Blue-fluorescence image is synthesized, and the standard of synthesis has:
A, the cell existed simultaneously in green fluorescence image, red fluorescence images and blue-fluorescence image is retained;
B, it is retained in the cell that the cell area in blue-fluorescence image is greater than the cell area in red fluorescence images;
The cell for meeting a and b condition simultaneously is retained, these cells are output to as the circulating tumor cell detected
In final examining report.
The present invention can be removed effectively by the detection of cell level due to raw brought erroneous detection of making an uproar, to cell level
Analysis processing, is able to solve a large amount of empty and noise problems that traditional micro-image is generated in pixel layer face treatment, and cellular layer
For surface treatment by each cell as individual, the cell for only meeting specified conditions can just enter number system;Simultaneously respectively
User can be convenient to the picture progress different disposal in three channels and define judgment criteria, so that the present invention has more practicability and standard
True property.