基于无透镜成像的椭圆细胞检测装置及检测方法与流程

本发明属于医学图像分析技术领域，涉及基于无透镜成像的椭圆细胞检测装置，还涉及椭圆细胞的检测方法。

背景技术：

血细胞形态参数对于疾病诊断、健康监测、新药研制是非常有价值的数据。通常获取这一数据的方法是采集细胞图像进行分析。常用的采集设备时光学显微镜。为了保证成像质量，光学显微镜不得不采用复杂的多透镜组。使得光学显微镜的小型化问题一直受限于透镜组。在面对突发疫情等紧急情况或采集现场检测等即时检测的应用场景中，较高的专业要求和制造成本成为了掣肘，不利于快速部署、低成本大规模部署和即时检测。

2006年，美国加州理工大学的yangchanghuei研究小组首次提出无透镜光流体显微芯片的概念。这一设备的优势在于使用图像传感器不仅降低了设备的制造和使用成本，还提高的便携性。在面对快速部署、低成本大规模部署和即时检测的应用场景中更具有优势。

无透镜设备有两个问题需要解决：1.由于商用图像传感器的像素尺寸与血液细胞十分接近，造成了由无透镜成像系统采集到的细胞图像分辨率较低；2.衍射干扰较小的接触式成像要求成像样品到图像传感器的距离必须缩小到几个可见光波长，对于商用传感器很难实现，不得不采用具有干扰的衍射成像。因此，如何在使用商用传感器的情况下，获得准确的细胞形态参数成为了研究难点。

技术实现要素：

本发明的目的是提供基于无透镜成像的椭圆细胞检测装置，能够针对无透镜成像系统采集细胞图像时存在分辨率低和衍射干扰严重的情况，准确检测出椭圆细胞。

本发明的另一目的是提供椭圆细胞的检测方法，可以准确检测出椭圆细胞。

本发明所采用的第一种技术方案是，基于无透镜成像的椭圆细胞检测装置，包括上下设置的灯光盒和电路板保护盒，灯光盒中由上向下依次设置有位于同一轴线上的led光源、凸透镜、cmos图像传感器，cmos图像传感器位于电路板保护盒上，电路板保护盒中设置有与cmos图像传感器连接的cmos图像传感器电路。

本发明第一种技术方案的特点还在于，

灯光盒与电路板保护盒旋转连接，灯光盒可在电路板保护盒上转动，灯光盒侧壁开设在有开口，电路板保护盒上固定有与灯光盒开口匹配的遮光板。

cmos图像传感器外周设置有样本垫。

样本垫上设置有样本卡槽。

本发明椭圆细胞的检测方法，该方法应用本发明的基于无透镜成像的椭圆细胞检测装置进行椭圆细胞检测，具体按照以下步骤实施：

步骤1、制作待测椭圆细胞样本；

步骤2、将待测椭圆细胞样本置于cmos图像传感器上，并固定于样本卡槽中，旋转灯光盒使其开口与遮光板形成封闭采集室；

步骤3、开启led光源，使led光源发出的光线经凸透镜聚拢垂直射向待测椭圆细胞样本；

步骤4、打开cmos图像传感器，将采集到的图像数据经数据接口导出到cmos图像传感器电路的处理模块上；

步骤5、对导出的图像数据使用霍夫变换算法检测经过衍射干扰后呈现圆形的椭圆细胞的图像的圆心，并以它为中心建立矩形坐标系，对图像进行分割检测细胞衍射环；

步骤6、在步骤5的基础上进行亚像素插值运算得到亚像素图像，并以步骤5中得到的圆心坐标为中心，将亚像素图像从矩形坐标系转换为极坐标系；

步骤7、通过超分辨率算法确定衍射亮环最强光强点到细胞衍射环中心的距离的曲线；

步骤8、根据步骤6的距离曲线，计算出衍射亮环的光强最大值点，表征衍射亮环的位置；

步骤9、以极坐标系参数θ从0度增加至360度，按照设置的步长递增，依次重复步骤7～8，直到得到360度的半径数据，将360度的半径数据组成一个半径序列{x′pk1(θ1)，x′pk2(θ2)，…，x′pkn(θn)}，x′pkn(θn)表示极坐标系参数θ对应的半径值，n为不为0的自然数

步骤10、通过两个半径之和表征衍射环直径，即匹配半径得到衍射环在角度α1,α2,…,αi处的直径长度，i是不为0的自然数：

式(10)中，d1(α1),d2(α2),…,di(αi)是衍射环在角度α1,α2,…,αi处的直径长度；

定义集合a来表示这些数据：

a＝{d1(α1),d2(α2)，…，di(αi)}(11)；

步骤11、进行椭圆细胞的匹配识别，椭圆长轴位置的阴影长度表示为：

dmaj(αmaj)＝min(a)(12)

假设π/4＜αmaj＜3π/4，定义集合b表示为

b＝{d1(β1),d2(β2)，…，di(βi)}(13)

椭圆短轴位置处阴影的长度可以表示为：

dmin(αmin)＝min(b)

定义参数m表示椭圆细胞长轴和短轴之间的差异：

m＝|dmaj(αmaj)-dmin(αmin)|

定义阈值mt，如果m＞mt,那么判断是正常的椭圆细胞，否则为其他细胞。

本发明第二种技术方案的特点还在于，

步骤7具体为，根据图像数据本身携带的光强信息，得到光强和距离衍射发生点的长度的映射关系，衍射发生点即细胞的边缘，这个关系表示在二维坐标系中会出现多个波峰，每一个波峰代表了一个衍射亮环，以波峰的最大值点表征衍射亮环的位置；

根据径向超分辨率和周向超分辨率得到衍射亮环最强光强点到细胞衍射环中心的距离的曲线：径向超分辨率是找到衍射亮环对应波峰的最快上升点和最快下降点，以这两个点的斜率的交点来修正波峰最大值点的位置；周向超分辨率，通过高斯平滑进行降噪，进而得到衍射亮环修正后的表征位置；通过衍射亮环的表征位置计算出衍射亮环到细胞衍射环的距离，得到衍射亮环最强光强点到细胞衍射环中心的距离的曲线。

本发明的有益效果是：

本发明基于无透镜成像的椭圆细胞检测装置，结构稳定，装置简易；置面向快速部署、低成本大规模部署和即时检测等应用场景；具有成本低，便携性高的优点；能够针对无透镜成像系统采集细胞图像时存在分辨率低和衍射干扰严重的情况，准确检测出椭圆细胞，实现了椭圆细胞检测装置的小型化。

本发明椭圆细胞的检测方法，简单清晰，可操作性强，提高了椭圆形细胞测量的自动化程度，降低了操作的人工成本。

附图说明

图1是本发明基于无透镜成像的椭圆细胞检测装置的结构示意图。

图中，1.灯光盒，2.led光源，3.凸透镜，4.样本卡槽，5.遮光板，6.样本垫，7.cmos图像传感器，8.电路板保护盒，9.样本垫。

具体实施方式

下面结合附图和具体实施方式对本发明进行详细说明。

本发明基于无透镜成像的椭圆细胞检测装置，如图1所示，包括上下设置的灯光盒1和电路板保护盒8，灯光盒1与电路板保护盒8旋转连接，灯光盒1可在电路板保护盒8上转动，灯光盒1侧壁开设在有开口，电路板保护盒8上固定有与灯光盒1开口匹配的遮光板5，灯光盒1中由上向下依次设置有位于同一轴线上的led光源2、凸透镜3、cmos图像传感器7，cmos图像传感器7外周设置有样本垫9，样本垫9上设置有样本卡槽4，cmos图像传感器7位于电路板保护盒8上，电路板保护盒8中设置有与cmos图像传感器7连接的cmos图像传感器电路。

本发明椭圆细胞的检测方法，该方法应用并本发明的基于无透镜成像的椭圆细胞检测装置进行椭圆细胞检测，具体按照以下步骤实施：

步骤1、制作待测椭圆细胞样本6。

步骤2、将待测椭圆细胞样本6置于cmos图像传感器7上，并固定于样本卡槽4中，旋转灯光盒1使其开口与遮光板形成封闭采集室。

步骤3、开启led光源2，使led光源2发出的光线经凸透镜3聚拢垂直射向待测椭圆细胞样本6。

步骤4、打开cmos图像传感器7，将采集到的图像数据经数据接口导出到cmos图像传感器电路的处理模块上。

步骤5、对导出的图像数据使用霍夫变换算法检测血细胞圆心，并依据血细胞圆心对图像进行分割，检测血细胞衍射环。

步骤6、在步骤5的基础上进行亚像素插值运算得到亚像素图像，并以步骤5中得到的圆心坐标为中心，将亚像素图像从矩形坐标系转换为极坐标系。

步骤7、通过超分辨率算法确定衍射亮环最强光强点到细胞衍射环中心的距离的曲线；

步骤7具体为，根据图像数据本身携带的光强信息，得到光强和距离衍射发生点的长度的映射关系，衍射发生点即细胞的边缘，这个关系表示在二维坐标系中会出现多个波峰，每一个波峰代表了一个衍射亮环，以波峰的最大值点表征衍射亮环的位置；

根据径向超分辨率和周向超分辨率得到衍射亮环最强光强点到细胞衍射环中心的距离的曲线：径向超分辨率是找到衍射亮环对应波峰的最快上升点和最快下降点，以这两个点的斜率的交点来修正波峰最大值点的位置；周向超分辨率，通过高斯平滑进行降噪，进而得到衍射亮环修正后的表征位置；通过衍射亮环的表征位置计算出衍射亮环到细胞衍射环的距离，得到衍射亮环最强光强点到细胞衍射环中心的距离的曲线。

步骤8、根据步骤6的距离曲线，计算出衍射亮环的光强最大值点，表征衍射亮环的位置。

步骤9、以极坐标系参数θ从0度增加至360度，按照设置的步长递增，如每次增加1度，依次重复步骤7～8，直到得到360度的半径数据，将360度的半径数据组成一个半径序列{x′pk1(θ1)，x′pk2(θ2)，…，x′pkn(θn)}，x′pkn(θn)表示极坐标系参数θ对应的半径值，n为不为0的自然数。

步骤10、通过两个半径之和表征衍射环直径，匹配半径得到衍射环在角度α1,α2,…,αi处的直径长度，i是不为0的自然数：

式(10)中，d1(α1),d2(α2),…,di(αi)是衍射环在角度α1,α2,…,αi处的直径长度；

定义集合a来表示这些数据：

a＝{d1(α1),d2(α2)，…，di(αi)}(11)。

步骤11、进行椭圆细胞的匹配识别，椭圆长轴位置的阴影长度表示为：

dmaj(αmaj)＝min(a)(12)

假设π/4＜αmaj＜3π/4，定义集合b表示为

b＝{d1(β1),d2(β2)，…，di(βi)}(13)

椭圆短轴位置处阴影的长度可以表示为：

dmin(αmin)＝min(b)

定义参数m表示椭圆细胞长轴和短轴之间的差异：

m＝|dmaj(αmaj)-dmin(αmin)|

定义阈值mt，如果m＞mt,那么判断是正常的椭圆细胞，否则为其他细胞。