CN105427327A - Image registration method of DNA sequencing and apparatus thereof - Google Patents

Abstract

The invention relates to an image registration method of DNA sequencing and an apparatus thereof. The method comprises the following steps of acquiring map information and acquiring a DNA map in a sampling time interval of each map respectively; acquiring a first pixel and a second pixel of the DNA map; calculating a global threshold T of an average gray value of the first pixel and the second pixel; calculating a variance sigma 2 of the first pixel and the second pixel; if the variance is in a preset scope, taking the T as the global threshold to segment the map; acquiring four kinds of fluorescence images; with the map as a reference image, selecting an image transformation model and carrying out registration on the four kinds of fluorescence images so as to obtain registered fluorescence images; finding magnetic bead center points; carrying out fusion processing on the center points so as to acquire a fused magnetic bead center point. Operation time is short and an image registration effect is good. After image identification of a reaction chip is increased, image identification accuracy is increased and determination of a basic group type is accurate. And conditions that an image in a normal map is blurred and the magnetic bead is missed to be checked are avoided.

Description

A kind of method for registering images of DNA sequencing and device

Technical field

The present invention relates to DNA sequencing analysis field, particularly relate to a kind of method for registering images and device of DNA sequencing.

Background technology

In DNA sequencing technical field, integrated operation flow process is described below: DNA sample is by after fragmentation, application is built storehouse reagent and is carried out adding joint, strand is caught, be bonded to microballoon, microemulsion pcr amplification, breakdown of emulsion liquid, obtain the DNA library be based upon on microballoon, application adds model by laies such as the enzymes of library and sequencing reaction needs to the sequence testing chip with micro reaction pool, sequence testing chip and sequencing reagent are mounted on main frame, order-checking program is started according to module number and position by computer for controlling, sequencing reaction is carried out in robotization, the data produced transfer to Data Analysis Computer, after completing order-checking, computation analysis software carries out image procossing, sequence reads, quality analysis, the work such as sequence assembly, finally obtain the sequence information of DNA sample.Micro reaction pool sequence testing chip is the carrier of sequencing reaction, and the DNABeads and the various sequencing reaction enzyme that are loaded with sequencing template are all arranged in the sequence testing chip being carved with micro reaction pool.

After the image acquisition to reaction chip, identification, the judgement to base type is had a strong impact on to the accuracy of image registration.

In view of above-mentioned defect, creator of the present invention obtains this creation finally through long research and practice.

Summary of the invention

The object of the present invention is to provide a kind of method for registering images and device of DNA sequencing, in order to overcome above-mentioned technological deficiency.

For achieving the above object, the invention provides a kind of method for registering images of DNA sequencing, comprising:

Obtain profile information, in the sampling time interval to each collection of illustrative plates, obtain DNA collection of illustrative plates respectively;

Obtain the first pixel and second pixel of described DNA collection of illustrative plates, wherein, the first pixel A is object pixel, and the gray-scale value of the first pixel is more than or equal to initial segmentation threshold value T ₀, sum of all pixels is N; Second pixel B is background pixel, and the gray-scale value of the second pixel is less than initial segmentation threshold value T ₀, sum of all pixels is M; The maximal value of collection of illustrative plates f (i, j) is V _max, minimum value is V _min

Wherein, T ₀=1/2 (V _min+ V _max) (1);

Calculate the global threshold T of the gray average of the first pixel and the second pixel;

$T=1/2(\frac{\underset{f(i,j)\&GreaterEqual;T_0}{\&Sigma;}f(i,j)}{N}+\frac{\underset{f(i,j)<T_0}{\&Sigma;}(i,j)}{M})---\left(2\right);$

Calculate the variances sigma of the first pixel and the second pixel ²

σ ²＝(P _A+P _B)(T-T ₀) ²(3)；

Wherein, the probability of the first pixel is:

$P_A=\underset{k=1}{\overset{M+N}{\&Sigma;}}\frac{N}{M+N}---\left(4\right)$

The probability of the second pixel is:

$P_B=\underset{k=1}{\overset{M+N}{\&Sigma;}}\frac{M}{M+N}---\left(5\right)$

If variance is in preset range, be then that global threshold is split described collection of illustrative plates with T;

Obtain CY3, CY5, FAM and TXR tetra-kinds of fluoroscopic images;

With described collection of illustrative plates for benchmark image, selected digital image transformation model, uses described image transform model to carry out registration to described four kinds of fluoroscopic images, obtains the fluoroscopic image after registration;

The magnetic bead central point searched in four kinds of fluoroscopic images after described registration;

Fusion treatment is carried out to described central point, to obtain the magnetic bead central point after fusion.

Further, described image transform model is any one in benchmark image and the four kinds of fluoroscopic images mapping relations on space and gray scale.

Further, the magnetic bead central point searched in four kinds of fluoroscopic images after described registration, specifically comprises:

To kind of the fluoroscopic image identification magnetic bead pixel of four after described registration, if the absolute value of f (i, j) minimum value min (i, j) difference is more than or equal to T ₀, be then identified as magnetic bead, otherwise, be background pixel.

Further, the magnetic bead central point searched in four kinds of fluoroscopic images after described registration, also comprises:

Travel through described magnetic bead pixel, obtain magnetic bead center pixel, wherein, magnetic bead center pixel is, the gray-scale value of current point equals the maximal value of gray scale in the 4*4 pixel region centered by current point, and four neighborhoods of current point are all object pixels.

Further, when obtaining profile information, to signal waveform, within every continuous print K cycle, each cycle selects the sampled point of n predetermined time, at interval of time T ₀sampling should be carried out, serial sampling M time; Can referential and accuracy in order to what ensure sampled data, the time interval Δ t of the n selected within each a cycle sampled point calculates according to the following equation,

$\mathrm{\&Delta;}t=\frac{a\&CenterDot;sin(\mathrm{\&omega;}+\mathrm{\&beta;})\&CenterDot;T}{\mathrm{\&lambda;}}---\left(6\right)$

In formula, Δ t represents the time interval of sampled point, and a is correction factor, its size is determined by sampled point quantity, and ω represents the angular frequency of photographing signals, is determined by the performance of CCD camera, β is initial phase angle, and T represents the time of signal period, and λ represents the peak value of signal waveform.

The present invention also provides a kind of image registration device of DNA sequencing, comprising:

Acquisition module, for obtaining profile information, obtains DNA collection of illustrative plates respectively in the sampling time interval to each collection of illustrative plates;

Computing module, for obtaining the first pixel and second pixel of described DNA collection of illustrative plates, wherein, the first pixel A is object pixel, and the gray-scale value of the first pixel is more than or equal to initial segmentation threshold value T ₀, sum of all pixels is N; Second pixel B is background pixel, and the gray-scale value of the second pixel is less than initial segmentation threshold value T ₀, sum of all pixels is M; The maximal value of collection of illustrative plates f (i, j) is V _max, minimum value is V _min

Wherein, T ₀=1/2 (V _min+ V _max) (1);

Calculate the global threshold T of the gray average of the first pixel and the second pixel;

$T=1/2(\frac{\underset{f(i,j)\&GreaterEqual;T_0}{\&Sigma;}f(i,j)}{N}+\frac{\underset{f(i,j)<T_0}{\&Sigma;}f(i,j)}{M})---\left(2\right);$

Variance computing module, for calculating the variances sigma of the first pixel and the second pixel ²

σ ²＝(P _A+P _B)(T-T ₀) ²(3)；

Wherein, the probability of the first pixel is:

$P_A=\underset{k=1}{\overset{M+N}{\&Sigma;}}\frac{N}{M+N}---\left(4\right);$

The probability of the second pixel is:

$P_B=\underset{k=1}{\overset{M+N}{\&Sigma;}}\frac{M}{M+N}---\left(5\right);$

Segmentation module, for when variance is in preset range, is that global threshold is split described collection of illustrative plates with T;

Also comprise: image registration module, described image registration module comprises:

Fluorescence acquisition module, for obtaining CY3, CY5, FAM and TXR tetra-kinds of fluoroscopic images;

Alignment modules, for described collection of illustrative plates for benchmark image, selected digital image transformation model, uses described image transform model to carry out registration to described four kinds of fluoroscopic images, obtains the fluoroscopic image after registration;

Search module, for the magnetic bead central point searched in kind of the fluoroscopic image of four after described registration;

Fusion Module, for carrying out fusion treatment to described central point, to obtain the magnetic bead central point after fusion.

Further, search described in module specifically for:

To spectrum recognition magnetic bead pixel, if the absolute value of f (i, j) minimum value min (i, j) difference is more than or equal to T ₀, be then identified as magnetic bead, otherwise, be background pixel.

Further, search described in module also for:

Travel through described magnetic bead pixel, obtain magnetic bead center pixel, wherein, magnetic bead center pixel is, the gray-scale value of current point equals the maximal value of gray scale in the 4*4 pixel region centered by current point, and four neighborhoods of current point are all object pixels.

Further, described acquisition module is specifically for when obtaining profile information, and to signal waveform, within every continuous print K cycle, each cycle selects the sampled point of n predetermined time, at interval of time T ₀sampling should be carried out, serial sampling M time; Can referential and accuracy in order to what ensure sampled data, the time interval Δ t of the n selected within each a cycle sampled point calculates according to the following equation,

$\mathrm{\&Delta;}t=\frac{a\&CenterDot;sin(\mathrm{\&omega;}+\mathrm{\&beta;})\&CenterDot;T}{\mathrm{\&lambda;}}---\left(6\right)$

In formula, Δ t represents the time interval of sampled point, and a is correction factor, its size is determined by sampled point quantity, and ω represents the angular frequency of photographing signals, is determined by the performance of CCD camera, β is initial phase angle, and T represents the time of signal period, and λ represents the peak value of signal waveform.

The invention provides a kind of method for registering images and device of DNA sequencing, by obtaining the first pixel and the second pixel, calculating the global threshold T of the gray average of the first pixel and the second pixel; And calculate the variances sigma of the first pixel and the second pixel ²if variance is in preset range, be then that global threshold is split described collection of illustrative plates with T.Obtain CY3, CY5, FAM and TXR tetra-kinds of fluoroscopic images; With described collection of illustrative plates for benchmark image, selected digital image transformation model, uses described image transform model to carry out registration to described four kinds of fluoroscopic images, obtains the fluoroscopic image after registration; The magnetic bead central point searched in four kinds of fluoroscopic images after described registration; Fusion treatment is carried out to described central point, to obtain the magnetic bead central point after fusion.Working time is short, effective to image registration, after improving the image recognition to reaction chip, to the accuracy of image recognition, and then the judgement accurately to base type.Avoid in conventional collection of illustrative plates image blurring unclear, the situation of magnetic bead under-enumeration.Further, recognizer is simple, and speed is fast, improves magnetic bead discrimination.

Accompanying drawing explanation

The process flow diagram of the method for registering images of the DNA sequencing that Fig. 1 provides for the embodiment of the present invention one;

The process flow diagram of the method for registering images of the DNA sequencing that Fig. 2 provides for the embodiment of the present invention two;

The structural representation of the image registration device of the DNA sequencing that Fig. 3 provides for the embodiment of the present invention three.

Embodiment

Below in conjunction with accompanying drawing, to above-mentioned being described in more detail with other technical characteristic and advantage of the present invention.

Embodiment one

Refer to shown in Fig. 1, the process flow diagram of the method for registering images of its DNA sequencing provided for the embodiment of the present invention one, comprising:

Step 101, obtains profile information, in the sampling time interval to each collection of illustrative plates, obtains DNA collection of illustrative plates respectively;

Particularly, the acquisition to DNA collection of illustrative plates can be realized in the following manner:

Reactant liquor enters in the reaction chip of DNA sequencer and carries out chemical reaction, produces visible ray; CCD camera, at suitable picture-taking position, is taken pictures to the visible ray information produced in reaction chip and is gathered DNA collection of illustrative plates.

More specifically, when obtaining profile information, to signal waveform, within every continuous print K cycle, each cycle selects the sampled point of n predetermined time, at interval of time T ₀sampling should be carried out, serial sampling M time; Can referential and accuracy in order to what ensure sampled data, the time interval Δ t of the n selected within each a cycle sampled point calculates according to the following equation,

$\mathrm{\&Delta;}t=\frac{a\&CenterDot;sin(\mathrm{\&omega;}+\mathrm{\&beta;})\&CenterDot;T}{\mathrm{\&lambda;}}---\left(6\right)$

In formula, Δ t represents the time interval of sampled point, and a is correction factor, its size is determined by sampled point quantity, and ω represents the angular frequency of photographing signals, is determined by the performance of CCD camera, β is initial phase angle, and T represents the time of signal period, and λ represents the peak value of signal waveform.

Through above-mentioned formula (6) sampling, when signal amplitude is larger, sample more intensive, sampled data can referential stronger; Sampling according to pre-conditioned sampling to signal data, makes follow-up signal transacting data volume reduce, alleviates numerous and diverse computing of data processing.

Step 102, obtains the first pixel and second pixel of described DNA collection of illustrative plates, and wherein, the first pixel A is object pixel, and the gray-scale value of the first pixel is more than or equal to initial segmentation threshold value T ₀, sum of all pixels is N; Second pixel B is background pixel, and the gray-scale value of the second pixel is less than initial segmentation threshold value T ₀, sum of all pixels is M; The maximal value of collection of illustrative plates f (i, j) is V _max, minimum value is V _min

Wherein, T ₀=1/2 (V _min+ V _max) (1);

Step 103, calculates the global threshold T of the gray average of the first pixel and the second pixel;

$T=1/2(\frac{\underset{f(i,j)\&GreaterEqual;T_0}{\&Sigma;}f(i,j)}{N}+\frac{\underset{f(i,j)<T_0}{\&Sigma;}f(i,j)}{M})---\left(2\right);$

Step 104, calculates the variances sigma of the first pixel and the second pixel ²

σ ²＝(P _A+P _B)(T-T ₀) ²(3)；

Wherein, the probability of the first pixel is:

$P_A=\underset{k=1}{\overset{M+N}{\&Sigma;}}\frac{N}{M+N}---\left(4\right)$

The probability of the second pixel is:

$P_B=\underset{k=1}{\overset{M+N}{\&Sigma;}}\frac{M}{M+N}---\left(5\right)$

Step 105 if variance is in preset range, is then that global threshold is split described collection of illustrative plates with T.

Step 106, obtains CY3, CY5, FAM and TXR tetra-kinds of fluoroscopic images;

Particularly, when CCD camera take pictures the DNA collection of illustrative plates that obtains smudgy time, and when fluoroscopic image is known, by obtaining fluoroscopic image, magnetic bead identification can be carried out.

Step 107, with described collection of illustrative plates for benchmark image, selected digital image transformation model, uses described image transform model to carry out registration to described four kinds of fluoroscopic images, obtains the fluoroscopic image after registration;

Preferably, described image transform model is any one in benchmark image and the four kinds of fluoroscopic images mapping relations on space and gray scale.

Step 108, the magnetic bead central point searched in the fluoroscopic image after described registration;

Step 109, carries out fusion treatment to described central point, to obtain the magnetic bead central point after fusion.

Particularly, under normal circumstances, the piece image of each magnetic bead in CY3, CY5, FAM and TXR tetra-kinds of images is brighter, and other three amplitude ratios are darker.But during different channel acquisition signal, because the time shutter is different with the setting of yield value, occur that in two width images, same magnetic bead is all bright, and the magnetic bead central point found out from two width images is not on same pixel, now needs to kick out of unnecessary magnetic bead central point.

In embodiment, by obtaining the first pixel and the second pixel, calculate the global threshold T of the gray average of the first pixel and the second pixel; And calculate the variances sigma of the first pixel and the second pixel ²if variance is in preset range, be then that global threshold is split described collection of illustrative plates with T.Obtain CY3, CY5, FAM and TXR tetra-kinds of fluoroscopic images; With described collection of illustrative plates for benchmark image, selected digital image transformation model, uses described image transform model to carry out registration to described four kinds of fluoroscopic images, obtains the fluoroscopic image after registration; The magnetic bead central point searched in four kinds of fluoroscopic images after described registration; Fusion treatment is carried out to described central point, to obtain the magnetic bead central point after fusion.Working time is short, effective to image registration, after improving the image recognition to reaction chip, to the accuracy of image recognition, and then the judgement accurately to base type.Avoid in conventional collection of illustrative plates image blurring unclear, the situation of magnetic bead under-enumeration.Further, recognizer is simple, and speed is fast, improves magnetic bead discrimination.

Embodiment two

The process flow diagram of the method for registering images of the DNA sequencing that Fig. 2 provides for the embodiment of the present invention two, refers to Fig. 2, and the present embodiment, on the basis of embodiment one, limits further, in the present embodiment,

The magnetic bead central point searched in fluoroscopic image after described registration, specifically comprises:

Step 201, to spectrum recognition magnetic bead pixel, if the absolute value of f (i, j) minimum value min (i, j) difference is more than or equal to T ₀, be then identified as magnetic bead, otherwise, be background pixel.

Preferably, can also comprise:

Step 202, travels through described magnetic bead pixel, and obtain magnetic bead center pixel, wherein, magnetic bead center pixel is, the gray-scale value of current point equals the maximal value of gray scale in the 4*4 pixel region centered by current point, and four neighborhoods of current point are all object pixels.

In embodiment, time smudgy for DNA collection of illustrative plates, and when fluoroscopic image is known, by obtaining fluoroscopic image, magnetic bead identification can be carried out.By obtaining CY3, CY5, FAM and TXR tetra-kinds of fluoroscopic images, with described collection of illustrative plates for benchmark image, registration is carried out to described four kinds of fluoroscopic images, obtain the fluoroscopic image after registration, the magnetic bead central point searched in fluoroscopic image after described registration, carries out fusion treatment to described central point, to obtain the magnetic bead central point after fusion, avoid in conventional collection of illustrative plates image blurring unclear, the situation of magnetic bead under-enumeration.Further, recognizer is simple, and speed is fast, improves magnetic bead discrimination.

Embodiment three

The structural representation of the image registration device of the DNA sequencing that Fig. 3 provides for the embodiment of the present invention three, refers to Fig. 3, and the pattern recognition device of a kind of DNA sequencing that the present embodiment provides, comprising:

Acquisition module 1, for obtaining profile information, obtains DNA collection of illustrative plates respectively in the sampling time interval to each collection of illustrative plates;

Particularly, the acquisition to DNA collection of illustrative plates can be realized in the following manner:

Reactant liquor enters in the reaction chip of DNA sequencer and carries out chemical reaction, produces visible ray; CCD camera, at suitable picture-taking position, is taken pictures to the visible ray information produced in reaction chip and is gathered DNA collection of illustrative plates.

More specifically, described acquisition module 1 is specifically for when obtaining profile information, and to signal waveform, within every continuous print K cycle, each cycle selects the sampled point of n predetermined time, at interval of time T ₀sampling should be carried out, serial sampling M time; Can referential and accuracy in order to what ensure sampled data, the time interval Δ t of the n selected within each a cycle sampled point calculates according to the following equation,

$\mathrm{\&Delta;}t=\frac{a\&CenterDot;sin(\mathrm{\&omega;}+\mathrm{\&beta;})\&CenterDot;T}{\mathrm{\&lambda;}}---\left(6\right)$

In formula, Δ t represents the time interval of sampled point, and a is correction factor, its size is determined by sampled point quantity, and ω represents the angular frequency of photographing signals, is determined by the performance of CCD camera, β is initial phase angle, and T represents the time of signal period, and λ represents the peak value of signal waveform.

Through above-mentioned formula (6) sampling, when signal amplitude is larger, sample more intensive, sampled data can referential stronger; Sampling according to pre-conditioned sampling to signal data, makes follow-up signal transacting data volume reduce, alleviates numerous and diverse computing of data processing.

Computing module 2, for obtaining the first pixel and second pixel of described DNA collection of illustrative plates, wherein, the first pixel A is object pixel, and the gray-scale value of the first pixel is more than or equal to initial segmentation threshold value T ₀, sum of all pixels is N; Second pixel B is background pixel, and the gray-scale value of the second pixel is less than initial segmentation threshold value T ₀, sum of all pixels is M; The maximal value of collection of illustrative plates f (i, j) is V _max, minimum value is V _min

Wherein, T ₀=1/2 (V _min+ V _max) (1);

Calculate the global threshold T of the gray average of the first pixel and the second pixel;

$T=1/2(\frac{\underset{f(i,j)\&GreaterEqual;T_0}{\&Sigma;}f(i,j)}{N}+\frac{\underset{f(i,j)<T_0}{\&Sigma;}f(i,j)}{M})---\left(2\right);$

Variance computing module 3, for calculating the variances sigma of the first pixel and the second pixel ²

σ ²＝(P _A+P _B)(T-T ₀) ²(3)；

Wherein, the probability of the first pixel is:

$P_A=\underset{k=1}{\overset{M+N}{\&Sigma;}}\frac{N}{M+N}---\left(4\right);$

The probability of the second pixel is:

$P_B=\underset{k=1}{\overset{M+N}{\&Sigma;}}\frac{M}{M+N}---\left(5\right);$

Segmentation module 4, for when variance is in preset range, is that global threshold is split described collection of illustrative plates with T.

Image registration module 5, described image registration module 5 comprises:

Fluorescence acquisition module 51, for obtaining CY3, CY5, FAM and TXR tetra-kinds of fluoroscopic images;

Particularly, when CCD camera take pictures the DNA collection of illustrative plates that obtains smudgy time, and when fluoroscopic image is known, by obtaining fluoroscopic image, magnetic bead identification can be carried out.

Alignment modules 52, for described collection of illustrative plates for benchmark image, selected digital image transformation model, uses described image transform model to carry out registration to described four kinds of fluoroscopic images, obtains the fluoroscopic image after registration;

Search module 53, for the magnetic bead central point searched in kind of the fluoroscopic image of four after described registration;

Fusion Module 54, for carrying out fusion treatment to described central point, to obtain the magnetic bead central point after fusion.

Particularly, under normal circumstances, the piece image of each magnetic bead in CY3, CY5, FAM and TXR tetra-kinds of images is brighter, and other three amplitude ratios are darker.But during different channel acquisition signal, because the time shutter is different with the setting of yield value, occur that in two width images, same magnetic bead is all bright, and the magnetic bead central point found out from two width images is not on same pixel, now needs to kick out of unnecessary magnetic bead central point.

In embodiment, by obtaining the first pixel and the second pixel, calculate the global threshold T of the gray average of the first pixel and the second pixel; And calculate the variances sigma of the first pixel and the second pixel ²if variance is in preset range, be then that global threshold is split described collection of illustrative plates with T.Obtain CY3, CY5, FAM and TXR tetra-kinds of fluoroscopic images; With described collection of illustrative plates for benchmark image, selected digital image transformation model, uses described image transform model to carry out registration to described four kinds of fluoroscopic images, obtains the fluoroscopic image after registration; The magnetic bead central point searched in four kinds of fluoroscopic images after described registration; Fusion treatment is carried out to described central point, to obtain the magnetic bead central point after fusion.

Working time is short, effective to image registration, after improving the image recognition to reaction chip, to the accuracy of image recognition, and then the judgement accurately to base type.Avoid in conventional collection of illustrative plates image blurring unclear, the situation of magnetic bead under-enumeration.Further, recognizer is simple, and speed is fast, improves magnetic bead discrimination.

Embodiment four

Please again consult Fig. 3, the present embodiment, on the basis of embodiment three, limits further, in the present embodiment, search module 53 specifically for: to spectrum recognition magnetic bead pixel, if f (i, j) absolute value of minimum value min (i, j) difference is more than or equal to T ₀, be then identified as magnetic bead, otherwise, be background pixel.

Preferably, search module 53 described in can also be used for:

Travel through described magnetic bead pixel, obtain magnetic bead center pixel, wherein, magnetic bead center pixel is, the gray-scale value of current point equals the maximal value of gray scale in the 4*4 pixel region centered by current point, and four neighborhoods of current point are all object pixels.

In embodiment, time smudgy for DNA collection of illustrative plates, and when fluoroscopic image is known, by obtaining fluoroscopic image, magnetic bead identification can be carried out.By obtaining CY3, CY5, FAM and TXR tetra-kinds of fluoroscopic images, with described collection of illustrative plates for benchmark image, registration is carried out to described four kinds of fluoroscopic images, obtain the fluoroscopic image after registration, the magnetic bead central point searched in fluoroscopic image after described registration, carries out fusion treatment to described central point, to obtain the magnetic bead central point after fusion, avoid in conventional collection of illustrative plates image blurring unclear, the situation of magnetic bead under-enumeration.Further, recognizer is simple, and speed is fast, improves magnetic bead discrimination.

Above-mentioned detailed description is illustrating for one of them possible embodiments of the present invention, and this embodiment is also not used to limit the scope of the claims of the present invention, and the equivalence that all the present invention of disengaging do is implemented or changed, and all should be contained in the scope of technical solution of the present invention.

Claims (9)

1. a method for registering images for DNA sequencing, is characterized in that, comprising:

Obtain profile information, in the sampling time interval to each collection of illustrative plates, obtain DNA collection of illustrative plates respectively;

Obtain the first pixel and second pixel of described DNA collection of illustrative plates, wherein, the first pixel A is object pixel, and the gray-scale value of the first pixel is more than or equal to initial segmentation threshold value T ₀, sum of all pixels is N; Second pixel B is background pixel, and the gray-scale value of the second pixel is less than initial segmentation threshold value T ₀, sum of all pixels is M; The maximal value of collection of illustrative plates f (i, j) is V _max, minimum value is V _min

Wherein, T ₀=1/2 (V _min+ V _max) (1);

Calculate the global threshold T of the gray average of the first pixel and the second pixel;

$T=1/2(\frac{\underset{f(i,j)\&GreaterEqual;T_0}{\&Sigma;}f(i,j)}{N}+\frac{\underset{f(i,j)<T_0}{\&Sigma;}f(i,j)}{M})---\left(2\right);$

Calculate the variances sigma of the first pixel and the second pixel ²

σ ²＝(P _A+P _B)(T-T ₀) ²(3)；

Wherein, the probability of the first pixel is:

$P_A=\underset{k=1}{\overset{M+N}{\&Sigma;}}\frac{N}{M+N}---\left(4\right)$

The probability of the second pixel is:

$P_B=\underset{k=1}{\overset{M+N}{\&Sigma;}}\frac{M}{M+N}---\left(5\right)$

If variance is in preset range, be then that global threshold is split described collection of illustrative plates with T;

Obtain CY3, CY5, FAM and TXR tetra-kinds of fluoroscopic images;

With described collection of illustrative plates for benchmark image, selected digital image transformation model, uses described image transform model to carry out registration to described four kinds of fluoroscopic images, obtains the fluoroscopic image after registration;

The magnetic bead central point searched in four kinds of fluoroscopic images after described registration;

Fusion treatment is carried out to described central point, to obtain the magnetic bead central point after fusion.

2. the method for registering images of DNA sequencing according to claim 1, is characterized in that,

Described image transform model is any one in benchmark image and the four kinds of fluoroscopic images mapping relations on space and gray scale.

3. the method for registering images of DNA sequencing according to claim 2, is characterized in that, the magnetic bead central point searched in four kinds of fluoroscopic images after described registration, specifically comprises:

To kind of the fluoroscopic image identification magnetic bead pixel of four after described registration, if the absolute value of f (i, j) minimum value min (i, j) difference is more than or equal to T ₀, be then identified as magnetic bead, otherwise, be background pixel.

4. the method for registering images of DNA sequencing according to claim 3, is characterized in that, the magnetic bead central point searched in four kinds of fluoroscopic images after described registration, also comprises:

Travel through described magnetic bead pixel, obtain magnetic bead center pixel, wherein, magnetic bead center pixel is, the gray-scale value of current point equals the maximal value of gray scale in the 4*4 pixel region centered by current point, and four neighborhoods of current point are all object pixels.

5. the method for registering images of DNA sequencing according to claim 1, is characterized in that, when obtaining profile information, to signal waveform, within every continuous print K cycle, each cycle selects the sampled point of n predetermined time, at interval of time T ₀sampling should be carried out, serial sampling M time; Can referential and accuracy in order to what ensure sampled data, the time interval Δ t of the n selected within each a cycle sampled point calculates according to the following equation,

$\mathrm{\&Delta;}t=\frac{a\&CenterDot;sin(\mathrm{\&omega;}+\mathrm{\&beta;})\&CenterDot;T}{\mathrm{\&lambda;}}---\left(6\right)$

In formula, Δ t represents the time interval of sampled point, and a is correction factor, its size is determined by sampled point quantity, and ω represents the angular frequency of photographing signals, is determined by the performance of CCD camera, β is initial phase angle, and T represents the time of signal period, and λ represents the peak value of signal waveform.

6. an image registration device for DNA sequencing, is characterized in that, comprising:

Acquisition module, for obtaining profile information, obtains DNA collection of illustrative plates respectively in the sampling time interval to each collection of illustrative plates;

Computing module, for obtaining the first pixel and second pixel of described DNA collection of illustrative plates, wherein, the first pixel A is object pixel, and the gray-scale value of the first pixel is more than or equal to initial segmentation threshold value T ₀, sum of all pixels is N; Second pixel B is background pixel, and the gray-scale value of the second pixel is less than initial segmentation threshold value T ₀, sum of all pixels is M; The maximal value of collection of illustrative plates f (i, j) is V _max, minimum value is V _min

Wherein, T ₀=1/2 (V _min+ V _max) (1);

Calculate the global threshold T of the gray average of the first pixel and the second pixel;

$T=1/2(\frac{\underset{f(i,j)\&GreaterEqual;T_0}{\&Sigma;}f(i,j)}{N}+\frac{\underset{f(i,j)<T_0}{\&Sigma;}f(i,j)}{M})---\left(2\right);$

Variance computing module, for calculating the variances sigma of the first pixel and the second pixel ²

σ ²＝(P _A+P _B)(T-T ₀) ²(3)；

Wherein, the probability of the first pixel is:

$P_A=\underset{k=1}{\overset{M+N}{\&Sigma;}}\frac{N}{M+N}---\left(4\right);$

The probability of the second pixel is:

$P_B=\underset{k=1}{\overset{M+N}{\&Sigma;}}\frac{M}{M+N}---\left(5\right);$

Segmentation module, for when variance is in preset range, is that global threshold is split described collection of illustrative plates with T;

Also comprise: image registration module, described image registration module comprises:

Fluorescence acquisition module, for obtaining CY3, CY5, FAM and TXR tetra-kinds of fluoroscopic images;

Alignment modules, for described collection of illustrative plates for benchmark image, selected digital image transformation model, uses described image transform model to carry out registration to described four kinds of fluoroscopic images, obtains the fluoroscopic image after registration;

Search module, for the magnetic bead central point searched in kind of the fluoroscopic image of four after described registration;

Fusion Module, for carrying out fusion treatment to described central point, to obtain the magnetic bead central point after fusion.

7. the image registration device of DNA sequencing according to claim 6, is characterized in that, described in search module specifically for:

To spectrum recognition magnetic bead pixel, if the absolute value of f (i, j) minimum value min (i, j) difference is more than or equal to T ₀, be then identified as magnetic bead, otherwise, be background pixel.

8. the image registration device of DNA sequencing according to claim 7, is characterized in that, described in search module also for:

Travel through described magnetic bead pixel, obtain magnetic bead center pixel, wherein, magnetic bead center pixel is, the gray-scale value of current point equals the maximal value of gray scale in the 4*4 pixel region centered by current point, and four neighborhoods of current point are all object pixels.

9. the image registration device of DNA sequencing according to claim 6, is characterized in that, described acquisition module is specifically for when obtaining profile information, to signal waveform, within every continuous print K cycle, each cycle selects the sampled point of n predetermined time, at interval of time T ₀sampling should be carried out, serial sampling M time; Can referential and accuracy in order to what ensure sampled data, the time interval Δ t of the n selected within each a cycle sampled point calculates according to the following equation,

$\mathrm{\&Delta;}t=\frac{a\&CenterDot;sin(\mathrm{\&omega;}+\mathrm{\&beta;})\&CenterDot;T}{\mathrm{\&lambda;}}---\left(6\right)$

In formula, Δ t represents the time interval of sampled point, and a is correction factor, its size is determined by sampled point quantity, and ω represents the angular frequency of photographing signals, is determined by the performance of CCD camera, β is initial phase angle, and T represents the time of signal period, and λ represents the peak value of signal waveform.