发明名称 一种应用于DNA测序的图像处理系统
摘要 本发明涉及一种应用于DNA测序的图像处理系统,所述的图像处理及控制系统获取图谱信息,在对每个图谱Δt的采样时间间隔内分别获取m个图谱,m为偶数;对同一碱基的拍照图谱进行分组,然后通过冗余判定的方式进行判定,并与标准的碱基图谱信息进行比对,A,T,C,G四种碱基信息,确定图谱信息中包含的碱基类型。本发明中冗余判断的方式,能够有效判断图谱采集存在问题的准确位置,方便采取措施,且软件运行的成本较低,可根据准确性任意设置图谱采集频率,适合大规模应用;同时又避免单一的灰度值对碱基类型的误判。
申请公布号 CN105400689A 申请公布日期 2016.03.16
申请号 CN201510872776.7 申请日期 2015.12.02
申请人 北京中科紫鑫科技有限责任公司 发明人 任鲁风;冯玉臣;张睿;陈哲;刘元杰;蔡亦梅;高静;范东雨;王绪敏;王者馥;殷金龙
分类号 C12M1/34(2006.01)I 主分类号 C12M1/34(2006.01)I
代理机构 代理人
主权项 一种应用于DNA测序的图像处理系统,其特征在于,包括设置有反应芯片的反应仓组件、运动平台组件、运动控制单元、设置在运动平台组件上并受运动控制单元控制的CCD相机、平台支撑组件以及控制CCD相机获取图像并对图像进行冗余处理的图像处理及控制系统,其中,所述的反应仓组件,在控制系统的作用下,反应液进入反应芯片中进行化学反应,产生可见光;所述的运动平台组件,其按照图像处理及控制系统的控制指令将CCD相机进给到适当的拍照位置,对反应芯片中产生的可见光信息进行拍照采集;所述的平台支撑组件,其支撑所述的反应仓组件、运动平台组件;所述的图像处理及控制系统获取图谱信息,在对每个图谱Δt的采样时间间隔内分别获取m个图谱,m为偶数;对同一碱基的拍照图谱进行分组,然后通过冗余判定的方式进行判定,并与标准的碱基图谱信息进行比对,A,T,C,G四种碱基信息,确定图谱信息中包含的碱基类型;所述的图像处理及控制系统包括分组单元、数据采集控制单元、比较模块和逻辑控制模块,所述的分组单元按照二位矩阵将上述采集的图谱或图像信息进行分组整理;获取m为12的连续图谱信息,分别每两个相邻近的时间,作为一二维组灰度值的两个维度值,按照标号1‑12,分别标记为6组二维灰度值,其中的第一、二灰度值为第一组,第三、四灰度值为第二组,第五、六灰度值为第三组;第七、八灰度值为第四组,第九、十灰度值为第五组,第十一、十二灰度值为第六组,所述分组单元将采集的灰度值信息按照上述分组将三组二维灰度值矩阵进行整理得到矩阵(i<sub>1</sub>,i<sub>2</sub>)、(i<sub>3</sub>,i<sub>4</sub>)、(i<sub>5</sub>,i<sub>6</sub>)、(i<sub>7</sub>,i<sub>8</sub>)、(i<sub>9</sub>,i<sub>10</sub>)、(i<sub>11</sub>,i<sub>12</sub>);所述比较模块内包括一比较器,其对所有数据采集模块内存储的灰度值采集信号进行运算处理,分别计算得出一重合度P<sub>n</sub>,并将其传输至所述逻辑控制模块内,所述比较模块按照下述公式(1)计算该两组二维灰度值矩阵信号中灰度值的重合度P<sub>21</sub>,<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>P</mi><mn>21</mn></msub><mo>=</mo><mfrac><mrow><msub><mi>&Sigma;</mi><mrow><msub><mi>i</mi><mn>4</mn></msub><mo>,</mo><msub><mi>i</mi><mn>3</mn></msub></mrow></msub><mrow><mo>(</mo><mi>T</mi><mo>(</mo><mrow><msub><mi>i</mi><mn>4</mn></msub><mo>,</mo><msub><mi>i</mi><mn>3</mn></msub></mrow><mo>)</mo><mo>*</mo><mi>I</mi><mo>(</mo><mrow><msub><mi>i</mi><mn>2</mn></msub><mo>+</mo><msub><mi>i</mi><mn>4</mn></msub><mo>,</mo><msub><mi>i</mi><mn>1</mn></msub><mo>+</mo><msub><mi>i</mi><mn>3</mn></msub></mrow><mo>)</mo><mo>)</mo></mrow></mrow><msqrt><mrow><msub><mi>&Sigma;</mi><mrow><msub><mi>i</mi><mn>4</mn></msub><mo>,</mo><msub><mi>i</mi><mn>3</mn></msub></mrow></msub><mi>T</mi><msup><mrow><mo>(</mo><msub><mi>i</mi><mn>4</mn></msub><mo>,</mo><msub><mi>i</mi><mn>3</mn></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>*</mo><msub><mi>&Sigma;</mi><mrow><msub><mi>i</mi><mn>4</mn></msub><mo>,</mo><msub><mi>i</mi><mn>3</mn></msub></mrow></msub><mi>I</mi><msup><mrow><mo>(</mo><msub><mi>i</mi><mn>2</mn></msub><mo>+</mo><msub><mi>i</mi><mn>4</mn></msub><mo>,</mo><msub><mi>i</mi><mn>1</mn></msub><mo>+</mo><msub><mi>i</mi><mn>3</mn></msub><mo>)</mo></mrow><mn>2</mn></msup></mrow></msqrt></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000865497760000011.GIF" wi="1572" he="247" /></maths>式中,P<sub>21</sub>表示每两组灰度值的重合度,i<sub>1</sub>和i<sub>2</sub>分别表示所述第一组二维灰度值矩阵的灰度值,i<sub>1</sub>表示第一灰度值的采样值,i<sub>2</sub>表示第二灰度值的采样值;i<sub>3</sub>和i<sub>4</sub>分别表示所述第二组二维灰度值矩阵的灰度值,i<sub>3</sub>表示第三灰度值的采样值,i<sub>4</sub>表示第四灰度值的采样值;T表示均方差运算,I表示积分运算;所述比较模块将上述重合度两两做差值并取绝对值,计算出结果后,将信号重合度与差值绝对值计算结果传输至所述逻辑控制模块内。
地址 101111 北京市大兴区经济技术开发区科创六街88号亦庄生物医药园3号楼608