主权项 |
一种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"><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>N</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="FDA0001128075250000011.GIF" wi="1573" he="164" /></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"><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="FDA0001128075250000012.GIF" wi="1613" he="134" /></maths>第二像素的概率为:<maths num="0003"><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="FDA0001128075250000013.GIF" wi="1509" he="132" /></maths>若方差在预设范围内,则以T为全局阈值对所述图谱进行分割。 |