主权项 |
一种DNA测序的图像配准方法,其特征在于,包括:获取图谱信息,在对每个图谱的采样时间间隔内分别获取DNA图谱;获取所述DNA图谱的第一像素和第二像素,其中,第一像素A为目标像素,第一像素的灰度值大于或等于初始分割阈值T<sub>0</sub>,像素总数为N;第二像素B为背景像素,第二像素的灰度值小于初始分割阈值T<sub>0</sub>,像素总数为M;图谱f(i,j)的最大值为V<sub>max</sub>,最小值为V<sub>min</sub>其中,T<sub>0</sub>=1/2(V<sub>min</sub>+V<sub>max</sub>) (1);计算第一像素和第二像素的灰度均值的全局阈值T;<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>T</mi><mo>=</mo><mn>1</mn><mo>/</mo><mn>2</mn><mrow><mo>(</mo><mfrac><mrow><munder><mo>Σ</mo><mrow><mi>f</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>≥</mo><msub><mi>T</mi><mn>0</mn></msub></mrow></munder><mi>f</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow></mrow><mi>N</mi></mfrac><mo>+</mo><mfrac><mrow><munder><mo>Σ</mo><mrow><mi>f</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo><</mo><msub><mi>T</mi><mn>0</mn></msub></mrow></munder><mi>f</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow></mrow><mi>M</mi></mfrac><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow><mo>;</mo></mrow>]]></math><img file="FDA0000873446250000011.GIF" wi="1580" he="175" /></maths>计算第一像素和第二像素的方差σ<sup>2</sup>σ<sup>2</sup>=(P<sub>A</sub>+P<sub>B</sub>)(T‑T<sub>0</sub>)<sup>2</sup> (3);其中,第一像素的概率为:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>P</mi><mi>A</mi></msub><mo>=</mo><munderover><mo>Σ</mo><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mrow><mi>M</mi><mo>+</mo><mi>N</mi></mrow></munderover><mfrac><mi>N</mi><mrow><mi>M</mi><mo>+</mo><mi>N</mi></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000873446250000012.GIF" wi="1630" he="134" /></maths>第二像素的概率为:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>P</mi><mi>B</mi></msub><mo>=</mo><munderover><mo>Σ</mo><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mrow><mi>M</mi><mo>+</mo><mi>N</mi></mrow></munderover><mfrac><mi>M</mi><mrow><mi>M</mi><mo>+</mo><mi>N</mi></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000873446250000013.GIF" wi="1692" he="134" /></maths>若方差在预设范围内,则以T为全局阈值对所述图谱进行分割;获取CY3、CY5、FAM和TXR四种荧光图像;以所述图谱为基准图像,选定图像变换模型,使用所述图像变换模型对所述四种荧光图像进行配准,得到配准后的荧光图像;在所述配准后的四种荧光图像中查找的磁珠中心点;对所述中心点进行融合处理,以获得融合后的磁珠中心点。 |