[发明专利]一种DNA测序的互信息图像配准方法及装置

说明书

技术领域

本发明涉及DNA测序分析领域，尤其涉及一种DNA测序的互信息图像配准方法及装置。

背景技术

在DNA测序技术领域，整体操作流程描述如下：DNA样品通过破碎后，应用建库试剂进行加接头、单链捕获、结合至微球、微乳液PCR扩增、破乳液，获得建立在微球上的DNA文库，应用加样板将文库和测序反应需要的酶等铺放至具有微反应池的测序芯片，测序芯片和测序试剂安装至主机上，通过控制计算机根据模块数量和位置启动测序程序，自动化进行测序反应，产生的数据传输至数据分析计算机，完成测序后应用计算分析软件进行图像处理、序列读出、质量分析、序列拼接等工作，最终得到DNA样本的序列信息。微反应池测序芯片是测序反应的载体，载有测序模板的DNA Beads及各种测序反应用酶均位于刻有微反应池的测序芯片中。

在对反应芯片的图像采集、识别后，对图像配准的精确性严重影响对碱基类型的判定。

鉴于上述缺陷，本发明创作者经过长时间的研究和实践终于获得了本创作。

发明内容

本发明的目的在于提供一种DNA测序的互信息图像配准方法及装置，用以克服上述技术缺陷。

为实现上述目的，本发明提供一种DNA测序的互信息图像配准方法，包括:

获取图谱信息，在对每个图谱的采样时间间隔内分别获取DNA图谱；

获取所述DNA图谱的第一像素和第二像素,其中,第一像素A为目标像素,第一像素的灰度值大于或等于初始分割阈值T₀,像素总数为N；第二像素B为背景像素,第二像素的灰度值小于初始分割阈值T₀,像素总数为M；图谱f(i,j)的最大值为V_max，最小值为V_min

其中，T₀＝1/2(V_min+V_max)(1)；

计算第一像素和第二像素的灰度均值的全局阈值T；

计算第一像素和第二像素的方差σ²

σ²＝(P_A+P_B)(T-T₀)²(3)；

其中，第一像素的概率为：

第二像素的概率为：

若方差在预设范围内，则以T为全局阈值对所述图谱进行分割；

获取CY3、CY5、FAM和TXR四种荧光图像；

以所述图谱为基准图像，对所述四种荧光图像进行互信息配准，得到配准后的荧光图像；

在所述配准后的四种荧光图像中查找的磁珠中心点；

对所述中心点进行融合处理，以获得融合后的磁珠中心点。

进一步地，所述以所述图谱为基准图像，对所述四种荧光图像进行互信息配准，得到配准后的荧光图像；

具体包括：

获取基准图像f₁(x,y)和第一图像f₂(x,y)，设置初始配准参数值(x,y,θ)和初始搜索方向，其中，第一图像为CY3、CY5、FAM和TXR四种荧光图像中的一种；

对所述第一图像f₂(x,y)进行空间几何变换，得到结果图像；

计算所述基准图像f₁(x,y)和所述记过图像的互信息值；

根据最大互信息理论，判断此时的配准参数值是否为最优，若是则输出所述配准参数值，使用所述配准参数值对所述第一图像进行几何变换和重采样，得到配准后的荧光图像；若不是，则修改配准参数值(x,y,θ)，直至所述配准参数值为最优。

进一步地，在所述配准后的四种荧光图像中查找的磁珠中心点，具体包括：

对所述配准后的四种荧光图像识别磁珠像素，如果f(i,j)最小值min(i,j)差值的绝对值大于等于T₀，则识别为磁珠，否则，为背景像素。

进一步地，在所述配准后的四种荧光图像中查找的磁珠中心点，还包括：

遍历所述磁珠像素，获取磁珠中心像素，其中，磁珠中心像素为，当前点的灰度值等于以当前点为中心的4*4像素区域内灰度的最大值，且当前点的四邻域都是目标像素。