[发明专利]基于主动轮廓模型和细胞神经网络的生物芯片分析方法

摘要

本发明公开了一种基于主动轮廓模型和细胞神经网络的生物芯片分析方法，包括以下步骤：对矩形样点采用改进的Hough变换进行倾斜校正，而对于圆形样点则采用改进的Radon变换；基于投影法对样点进行初次定位，并生成优化后的网格；然后基于邻域搜索自适应地调整网格，对样点进行二次精确定位；根据贪婪算法优化主动轮廓模型，利用CNN将样点按信号强弱进行分类；将Multiple snakes与CNN结合，CNN先学习强信号样点snake的收敛行为，再指导弱信号样点snake的收敛，最终实现样点的合理分割；提取并输出微阵列样点信号数据。本发明解决了生物芯片图像倾斜校正、形状不规则样点及弱信号样点分割困难等问题，实现了生物芯片样点的自动识别，适合大规模的生物芯片样点快速分析。

主权项

一种基于主动轮廓模型和细胞神经网络的生物芯片分析方法，其特征在于，该方法为：1)若原始生物芯片图像未发生倾斜，即图像准直，则进入4）；若图像发生倾斜，即倾斜角度在[‑5°,5°]范围内，则进入2）；若图像倾角超出[‑5°,5°]范围，则认为芯片图像的质量不符合要求，不进行后续的分析；2)将生物芯片图像二值化，计算二值的生物芯片图像的倾斜角度，得到其相反数为最佳的倾斜校正角度；3)利用所述最佳的倾斜校正角度和图像几何空间变换的原理，对倾斜的原始生物芯片图像进行旋转校正并灰度化，得到准直的芯片灰度图像；4)对所述准直的芯片灰度图像进行投影计算，得到准直图像在水平和竖直两个方向上的投影波形PH(x),PV(y)，将投影波形进行开闭级联运算操作和自适应二值化处理，得到二值投影方波；5)二值投影方波PH(x)的波峰中心对应着每一行样点的中心，在所述中心划分网格线，确定样点行在水平方向的坐标；二值投影方波PH(x)的波谷中心对应着样点行间空隙的中心位置，在所述中心位置划分网格线，将每个样点分割到各自的子区域中；运用上述同样的方法处理二值投影方波PV(y)；6)读取上述二值投影方波的各个波峰和波谷的宽度，以及波峰中心与波谷中心的位置；运用数学统计学原理，计算出贴近真实的平均波峰宽度和平均波谷宽度，以及平均波峰中心距；7)结合数学统计学原理，通过判断每条网格线对应的方波的波峰宽度和中心距，发现冗余网格线并予以剔除；对剩余的网格线重新进行计算，添补缺失的网格线；对网格线的位置与数量进行调整，生成优化后的网格线，完成对样点的初次定位；8)将步骤6）获得的平均波峰宽度用作二次精确定位网格的初始直径，对每个样点生成一个二次精确定位初始网格；对于圆形样点的芯片图像，网格设计成圆形；对于矩形样点的芯片图像，网格设计成矩形；9)基于邻域搜索和自适应调整的原理，对样点进行二次精确定位：运用邻域搜索算法，对上述产生的二次精确定位网格的中心位置逐个进行调整；通过比较网格内部样点信号数据，确定尺寸调整方向，再按照直径调整准则进行调整；使得网格最终能将样点套住，信号区域完全落在网格内部，同时套住的背景区域像素数目达到最小值；10)根据贪婪算法优化主动轮廓模型，利用细胞神经网络CNN将芯片样点按信号强弱进行分类；将多主动轮廓模型Multiple snakes与CNN结合，CNN先学习强信号芯片样点snake的收敛行为，并得到强信号样点的snake轮廓；再指导弱信号芯片样点snake的收敛，得到弱信号芯片样点的snake轮廓，最终实现芯片样点信号区域的分割，提取芯片样点信号区域内灰度强度作为后续处理的数据；11)计算每个样点信号区域内部的灰度总值、平均值和方差，统计得到芯片样点的信号值；计算样点snake轮廓外部环形区域的灰度强度信息，统计出每个样点的背景值，最后输出生物芯片样点分析数据。