CN105184853B - 一种基于光流分析的单细胞三维图像生成方法 - Google Patents
一种基于光流分析的单细胞三维图像生成方法 Download PDFInfo
- Publication number
- CN105184853B CN105184853B CN201510500965.1A CN201510500965A CN105184853B CN 105184853 B CN105184853 B CN 105184853B CN 201510500965 A CN201510500965 A CN 201510500965A CN 105184853 B CN105184853 B CN 105184853B
- Authority
- CN
- China
- Prior art keywords
- mrow
- cell
- msub
- dimensional image
- optical flow
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Classifications
-
- C—CHEMISTRY; METALLURGY
- C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
- C12Q—MEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
- C12Q1/00—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions
- C12Q1/02—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions involving viable microorganisms
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N27/00—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means
- G01N27/26—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating electrochemical variables; by using electrolysis or electrophoresis
- G01N27/416—Systems
- G01N27/447—Systems using electrophoresis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Organic Chemistry (AREA)
- Molecular Biology (AREA)
- Immunology (AREA)
- General Health & Medical Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- Proteomics, Peptides & Aminoacids (AREA)
- General Physics & Mathematics (AREA)
- Zoology (AREA)
- Wood Science & Technology (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Microbiology (AREA)
- Computer Graphics (AREA)
- Software Systems (AREA)
- Biotechnology (AREA)
- Biophysics (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Chemical Kinetics & Catalysis (AREA)
- Genetics & Genomics (AREA)
- Theoretical Computer Science (AREA)
- Electrochemistry (AREA)
- Pathology (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
本发明公开一种基于光流分析的单细胞三维图像生成方法。本发明的方法具有如下优点：1)通过显微视觉算法实现的增强光诱导介电泳可操控平台的反馈控制功能；2)采用基于光流场的运动检测方法实现对细胞三维运动的跟踪，从而完成细胞动力学模型的参数估计；3)通过控制单个细胞的旋转，获取细胞二维图像序列，利用最大似然估计方法实现基于可控单细胞的二维图像序列的细胞的三维图像生成技术。本发明的方法，其所用设备简单，成本低，不会产生光毒作用，对样本的要求也较低，提高了可操作性。
Description
技术领域
本发明涉及细胞成像领域，尤其涉及一种基于光流分析的单细胞三维图像生成方法。
背景技术
与早期的成像系统相比，活细胞的3D成像呈现了细胞及其组分的更详细、也更准确的空间视图。技术进步让3D成像成为许多应用的重要工具，如细胞生物学、发育生物学、神经科学以及癌症研究。当前的技术比以往更加准确，能实时给出数据，几乎不需要细胞制备。现有技术的单细胞三维方法主要有以下几种：
1、激光扫描共聚焦显微镜：利用激光束经照明针孔形成点光源对细胞内焦平面的每一点扫描，细胞上的被照射点，在探测针孔处成像，由探测针孔后的光电倍增管(PMT)或冷电耦器件(cCCD)逐点或逐线接收，迅速在计算机监视器屏幕上形成荧光图像。照明针孔与探测针孔相对于物镜焦平面是共轭的，焦平面上的点同时聚焦于照明针孔和发射针孔，焦平面以外的点不会在探测针孔处成像，这样得到的共聚焦图像是细胞的光学横断面。激光扫描共聚焦显微镜通过对同一细胞不同层面的实时扫描成像，进行图像叠加可构成细胞的三维结构图像。
2、白色光衍射断层扫描成像技术：该技术能够为透明样本如活细胞和未标记的细胞成像，基于传统显微镜和白光，在细胞的自然状态下提供高分辨率的3D渲染图像。物镜镜头扫描细胞的整个轴向焦面，产生一叠相位分辨图像，然后通过sparse反卷积算法重建物体的三维结构。能得到350nm的侧向分辨率以及900nm的轴向分辨率。
3、晶格光片显微镜：晶格光片显微镜有两个正交的镜头；一个镜头将光聚焦产生一条非常细的笔状光源，照射在有萤光分子的生物样品，产生荧光；另一个采用宽场成像收集荧光，由面扫描以快速取得3D高清晰度生物影像。依靠空间光调制器同时形成100多条笔状的光束，来增加扫描速度，降低对生物样品的伤害；而且能控制每条光束的距离及形状。
但上述方法均存在不足：
第一种方式需要专用的设备，成本太高；且会产生光毒作用：在激光照射下，许多荧光染料分子会产生单态氧或自由基等细胞毒素，限制扫描时间、激发光强度，以保持样品的活性；标记染料的光漂白：为了获得足够的信噪比必须提高激光的强度；而高强度的激光会使染料在连续扫描过程中迅速褪色。
第二种方式不适用于非半透明样本，对于样本的要求比较高；需要专用设备，成本过高；成像灵活度低。
第三种方式设备成本过高；需要对样本进行荧光处理，降低了可操作性；光毒作用依然存在。
因此，现有技术还有待于改进和发展。
发明内容
鉴于上述现有技术的不足，本发明的目的在于提供一种基于光流分析的单细胞三维图像生成方法，旨在解决现有的细胞三维图像生成方法成本高、有光毒作用、条件苛刻等问题。
本发明的技术方案如下：
一种基于光流分析的单细胞三维图像生成方法，其中，包括步骤：
A、制作光诱导介电泳芯片，所述光诱导介电泳芯片有三层结构组成：有三层结构组成：下层为涂有氢化非晶硅涂层的ITO玻璃，上层是不含涂层的ITO玻璃，在上下两层ITO玻璃之间封装有一个微流体通道，用于注射所需操作的溶液；
B、将细胞和溶液注射到微流体通道，并向上下两层ITO玻璃的电极输入可变频率的交流信号，同时利用入射光照射所述光诱导介电泳芯片，从而在被照射的区域产生非均匀电场；
C、改变交流信号的频率及大小，以控制细胞运动方向，同时采集细胞的图像；
D、对采集的图像进行预处理，然后进行特征提取以及速度计算，最后重构3D细胞图像。
所述的基于光流分析的单细胞三维图像生成方法，其中，所述步骤A中，制作光诱导介电泳芯片的步骤具体包括：
A1、清理ITO玻璃基质；
A2、在ITO玻璃基质上沉积氢化非晶硅涂层；
A3、在氢化非晶硅涂层上涂光刻胶；
A4、在光刻胶上进行板印；
A5、接触腐蚀至ITO玻璃基质；
A6、去除光刻胶；
A7、在ITO玻璃基质上未覆盖氢化非晶硅涂层的区域涂导电粘合剂。
所述的基于光流分析的单细胞三维图像生成方法，其中，所述细胞在非均匀电场中的所受到的平均介电泳力用如下公式描述：
其中FDEP是作用到细胞上的平均介电泳力，R是细胞的半径，εm是细胞所在溶液的介电常数，Erms为所施加交流信号的均方根值，fCM为Clausius-Mossotti因子，在计算平均介电泳力时取该因子的实部Re[fCM]。
所述的基于光流分析的单细胞三维图像生成方法，其中，fCM因子定义如下:
εp*和εm*分别是细胞和溶液的复介电常数。
所述的基于光流分析的单细胞三维图像生成方法，其中，所述复介电常数可表示为:
其中，ε是溶液的介电常数，σ是导电率，ω是所施加交流信号的频率。
所述的基于光流分析的单细胞三维图像生成方法，其中，细胞旋转速度为：
其中E是电场强度，η是溶液的黏稠度，IM[fCM]是Clausius-Mossotti因子的虚部，K为系数。
所述的基于光流分析的单细胞三维图像生成方法，其中，所述预处理包括：高斯滤波处理、亮度调整及模板匹配。
所述的基于光流分析的单细胞三维图像生成方法，其中，所述步骤D中，使用机器学习算法对模型参数进行最大似然估计，按下式最小化的参数取值：
其中Ii+1是细胞旋转图像序列I＝{Ii,i＝1,…,n}中的一帧，M代表细胞上所有点的集合，Ri和Ti分别是旋转矩阵和平移向量,Rimi+Ti是细胞上一点m从i时刻到i+1时刻的三维旋转运动模型，即mi+1，映射是将细胞在i时刻的三维信息投影到细胞旋转图像的某一帧Ii。
有益效果：本发明的方法具有如下优点：1)通过显微视觉算法实现的增强光诱导介电泳可操控平台的反馈控制功能；2)采用基于光流场的运动检测方法实现对细胞三维运动的跟踪，从而完成细胞动力学模型的参数估计；3)通过控制单个细胞的旋转，获取细胞二维图像序列，利用最大似然估计方法实现基于可控单细胞的二维图像序列的细胞的三维图像生成技术。本发明的方法，其所用设备简单，成本低，不会产生光毒作用，对样本的要求也较低，提高了可操作性。
附图说明
图1为本发明一种基于光流分析的单细胞三维图像生成方法较佳实施例的流程图。
图2为本发明中光诱导介电泳平台的结构示意图。
具体实施方式
本发明提供一种基于光流分析的单细胞三维图像生成方法，为使本发明的目的、技术方案及效果更加清楚、明确，以下对本发明进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。
请参阅图1，图1为本发明一种基于光流分析的单细胞三维图像生成方法较佳实施例的流程图，如图所示，其包括步骤：
S100、制作光诱导介电泳芯片(ODEP芯片)，所述光诱导介电泳芯片有三层结构组成：下层为涂有氢化非晶硅涂层的ITO玻璃，上层是不含涂层(即不含氢化非晶硅涂层)的ITO玻璃，在上下两层ITO玻璃之间封装有一个微流体通道，用于注射所需操作的溶液；
S200、向上下两层ITO玻璃的电极输入可变频率的交流信号，同时利用入射光照射所述光诱导介电泳芯片，从而在被照射的区域产生非均匀电场；可先向微流体通道注入细胞和介质(介质即所需操作的溶液，也即细胞所在溶液)。然后输入交流信号。
S300、改变交流信号的频率及大小，以控制细胞运动方向，同时采集细胞的图像；
S400、对采集的图像进行预处理，然后进行特征提取以及速度计算，最后重构3D细胞图像。
进一步，所述的步骤S100中，制作光诱导介电泳芯片的步骤具体包括：
S101、清理ITO玻璃基质；
清理ITO玻璃基质的表面，保证接触面的洁净度。
S102、在ITO玻璃基质上沉积氢化非晶硅涂层(a-Si:H)；
在ITO玻璃基质表面沉积一层氢化非晶硅，厚度为1微米。
S103、在氢化非晶硅涂层上涂光刻胶；
S104、在光刻胶上进行板印；
板印是按照指定图形制作遮盖物，将遮盖物放在光刻胶表面，用紫外线照射遮盖物，没有被遮盖的光刻胶在紫外线作用下溶解，最终得到与遮盖物形状相同的光刻胶层。
S105、接触腐蚀至ITO玻璃基质；具体是用草酸腐蚀制作的芯片表层，以去除没有覆盖光刻胶的氢化非晶硅涂层。
S106、去除光刻胶；即将光刻胶从氢化非晶硅涂层表面去除。
S107、在ITO玻璃基质上未覆盖氢化非晶硅涂层的区域涂导电粘合剂。即在ITO玻璃的表面没有覆盖氢化非晶硅涂层的位置添加一个导电触点。
而上层的ITO玻璃清理干净之后，涂导电粘合剂即可。
在上下两层ITO玻璃之间封装有一个微流体通道(100微米高)，具体是通过PDMS或是双面胶封装出一个微流体通道。
在步骤S200中，如图2所示，首先搭建光诱导介电泳平台。除了步骤S100制作的ODEP芯片20，平台还需要一台光学显微镜10、一台光学投影仪(高分辨率)、一个可编程信号发生电路和主机系统。所述主机系统包括：图像采集模块、显微视觉算法处理模块、生物芯片驱动控制器、虚拟电极生成模块以及显示输出模块。所述图像采集模块用来采集光学显微镜20的图像，并交由显微视觉算法处理模块来进行处理并通过显示输出模块来显示，所述显微视觉算法处理模块还向生物芯片驱动控制器及虚拟电极生成模块发出信号用来控制二者工作。所述生物芯片驱动控制器连接所述可编程信号发生电路来改变信号频率和大小。所述可编程信号发生电路通过电极连接所述ODEP芯片20。所述光学投影仪设置在ODEP芯片20下方，用来对其进行入射光照射。所述虚拟电极生成模块连接所述光学投影仪。
其中光学显微镜参数如下：
尼康CFI60无限远光学系统；
电动对焦，可上下移动(上13mm/下2mm)；
三目镜筒，光分布：目镜/相机100％/0，20％/100％，0/100％；
目镜放大倍率：10x；
聚光器：防水，工作距离：7.2mm；
物镜：20x，高度消色透镜，纳米结晶涂层；
载物台：电动X轴和Y轴，分辨率：0.1微米；
紫外线截止滤光块；
荧光滤波套装：FITC/GFP。
在平台搭建好后，可通过生物芯片驱动控制器向可编程信号发生电路发出信号，然后可编程信号发生电路向上下两层ITO玻璃的电极输入可变频率的交流信号，同时光学投影仪利用入射光照射所述光诱导介电泳芯片，从而在被照射的区域产生非均匀电场。
在所述步骤S300中，通过改变交流信号的频率及大小，来改变细胞所受到的介电泳力的方向与大小，以控制细胞运动方向，同时采集细胞的图像(2D细胞图像)，实现高速操纵微纳米实体。
最后在步骤S400中，对采集的图像进行预处理，然后进行特征提取以及速度计算，最后重构3D细胞图像。
下面先着重介绍下，如何实现由改变交流信号的频率及大小来控制细胞运动方向。
细胞在非均匀电场中的所受到的平均介电泳力可以用如下公式描述：
其中FDEP是作用到细胞上的平均介电泳力，R是细胞的半径，εm是细胞所在溶液的介电常数，Erms为所施加电场(交流信号)的均方根值，fCM为Clausius-Mossotti因子，在计算平均介电泳力时取该因子的实部Re[fCM]，该因子定义如下:
εp*和εm*分别是细胞和溶液的复介电常数，公式2中的复介电常数(包括εp*和εm*)可表示为:
其中，ε是溶液的介电常数，σ是导电率，ω是所加电场(交流信号)的频率。
可以看出fCM是一个和频率相关的可变因子。考虑在施加不同频率的交变电场下，当介电泳力与电场强度变化方向相同时，称为正介电泳现象；当所受到的介电泳力与电场强度变化方向相反，称为负介电泳现象。因而可以通过改变所施加的电场的频率，来改变细胞所受到的介电泳力的方向，达到控制细胞运动方向的目的。
由于生物细胞受到非均匀电场的极化作用而产生偶极矩，根据其介电泳力所产生的转矩与所在介质中受到的摩擦力矩达到平衡，细胞旋转速度为：
其中E是电场强度，η是溶液的黏稠度，IM[fCM]是Clausius-Mossotti因子的虚部，K为系数。根据细胞的旋转速度与细胞的介电常数的关系可以对细胞的介电特性进行估算。
细胞受到的介电泳力强度与方向主要取决于介质与细胞的介电特性，如形状、尺寸与电场频率。本发明利用光诱导介电泳力(ODEP)(当施加某频段，电液动力学的一种主导力)以识别与操纵生物细胞。ODEP芯片由可变频率的交流信号驱动，交流信号通过上下两层ITO玻璃的导电触点输入，此时在溶液层只有一小部分分压，并在溶液层中产生均匀电场。当入射光照射ODEP芯片，a-Si:H的光导率由于电子空穴对数的增多而增加几个数量级。由于入射光区域电阻减小，在溶液层中的分压会大大增大，于是入射光区域的a:Si:H将成为一个有效的虚拟电极产生非均匀电场。这种光诱导的非均匀电场会极化区域内的颗粒产生介电泳力，也就是光诱导介电泳力(ODEP)。通过光学显微镜与主机系统可实现程序化的动态运动，且不需要任何手工界面而实现微纳米实体的自动化捕获、操纵、分离与组装。
对于步骤S400中，细胞三维图像生成的原理如下：当一个极化的物体置于非均匀电场时，在偶极矩的作用下物体会向电场最强或者最弱处运动，方向取决于物体相对于介质的极性。根据细胞在介电泳力场下的动力学模型，通过改变驱动光诱导介电泳生物芯片的交流信号的大小与频率，并配合相应的投射到ODEP芯片上的入射光，控制细胞在光诱导介电泳力的作用下的旋转运动。细胞旋转至不同位置的图像序列I＝{Ii,i＝1,…,n}，n为时刻点，细胞的三维图像是从二维图像序列重构细胞的三维结构。重构过程包括：
一、在获得二维图像后，先进行预处理：
1、先对图像进行高斯滤波处理，滤除掉图像中噪声。高斯滤波的核心公式如下：
其中σ为函数的宽度参数,控制了函数的径向作用范围。
2、然后进行亮度调整处理。滤除噪声后，以背景色为基准，对图像整体亮度作线性变换：
g(x,y)＝c+k(f(x,y)-a) (6)
其中，f(x,y)和g(x,y)分别为图像中某点(x,y)的原始亮度和变换后的亮度，a为背景亮度，k为变换系数，c为亮度补偿。
3、再进行模板匹配处理。对图像中的细胞进行模板匹配，比对同一窗口中的块在上下两帧的差异，计算相关系数ρXY。
其中X和Y分别为上下两帧中的块，Cov(X,Y)为X和Y的协方差，D(X)和D(Y)分别为X和Y的方差。
二、然后进行特征提取。
1、找到相关系数的局部最大值以跟踪峰点，从一个峰点到下一个峰点即代表旋转一圈，根据峰点的索引估算出旋转的圈数。
2、光流法分析像素的运动向量。对于图像上的一个像素点X(x,y)在t时刻的亮度值为I(x,y,t)，u(x,y)和v(x,y)表示(x,y)处光流在x和y方向的运动分量。根据图像序列，计算和令 W＝diag(W(X1),…,W(Xn))，其中n为点的个数，diag()为构造对角矩阵，W为高斯函数：
令V＝(u,v)T，那么其计算公式如下：
V＝(ATW2A)-1ATW2b (9)
三、再进行速度计算。
1、根据得到的自转圈数以及相邻两个图像序列的拍摄间隔计算细胞的自转速度ω，结合细胞的旋转模型，即可得到细胞自转的旋转矩阵K。因此旋转后的细胞上所有点的集合M'＝MK，其中M为细胞上所有点的原始坐标。每个点的运动速度
2、利用光流法得到的细胞二维运动速度，修正上一步得到的三维运动速度。对于某一点a，其三维运动速度为v3D＝(x3D,y3D,z3D)，二维运动速度为v2D＝(x2D,y2D)，计算修正系数其中a和b分别为x轴和y轴修正系数。最终得到的三维运动速度V＝Kv3D。
四、最后进行细胞重构。
具体的，根据细胞的刚体模型，细胞上一点m从时刻tk的位置(xk,yk,zk)经过旋转和平移，运动到时刻tk+1的位置(xk+1,yk+1,zk+1)。设旋转矩阵和平移向量分别是Rk和Tk,则细胞三维旋转运动模型为：
mk+1＝Rkmk+Tk (10)
定义映射其中M代表细胞上所有点的集合，将细胞在i时刻的三维信息投影到细胞旋转图像的某一帧Ii。因此，细胞的三维图像生成归结为模型参数的最大似然估计，使用机器学习算法对模型产生进行最大似然估计，按下式最小化的参数取值。
其中Ii+1是细胞旋转图像序列I＝{Ii,i＝1,…,n}中的一帧，M代表细胞上所有点的集合，Ri和Ti分别是旋转矩阵和平移向量,Rimi+Ti是细胞上一点m从i时刻到i+1时刻的三维旋转运动模型，即mi+1。映射是将细胞在i时刻的三维信息投影到细胞旋转图像的某一帧Ii。
利用训练好的参数，即可根据细胞旋转的二维图像重构映射为当前帧的3D细胞模型。本发明通过光流的方法获得的细胞旋转运动信息将更准确的估计旋转矩阵Rk，进而优化细胞旋转的三维模型。
本发明的方法具有如下优点：1)通过显微视觉算法实现的增强光诱导介电泳可操控平台的反馈控制功能；2)采用基于光流场的运动检测方法实现对细胞三维运动的跟踪，从而完成细胞动力学模型的参数估计；3)通过控制单个细胞的旋转，获取细胞二维图像序列，利用最大似然估计方法实现基于可控单细胞的二维图像序列的细胞的三维图像生成技术。本发明的方法，其所用设备简单，成本低，不会产生光毒作用，对样本的要求也较低，提高了可操作性。
应当理解的是，本发明的应用不限于上述的举例，对本领域普通技术人员来说，可以根据上述说明加以改进或变换，所有这些改进和变换都应属于本发明所附权利要求的保护范围。
Claims (6)
1.一种基于光流分析的单细胞三维图像生成方法，其特征在于，包括步骤：
A、制作光诱导介电泳芯片，所述光诱导介电泳芯片有三层结构组成：有三层结构组成：下层为涂有氢化非晶硅涂层的ITO玻璃，上层是不含涂层的ITO玻璃，在上下两层ITO玻璃之间封装有一个微流体通道，用于注射所需操作的溶液；
B、将细胞和溶液注射到微流体通道，并向上下两层ITO玻璃的电极输入可变频率的交流信号，同时利用入射光照射所述光诱导介电泳芯片，从而在被照射的区域产生非均匀电场；
C、改变交流信号的频率及大小，以控制细胞运动方向，同时采集细胞的图像；
D、对采集的图像进行预处理，然后进行特征提取以及速度计算，最后重构3D细胞图像；
所述步骤A中，制作光诱导介电泳芯片的步骤具体包括：
A1、清理ITO玻璃基质；
A2、在ITO玻璃基质上沉积氢化非晶硅涂层；
A3、在氢化非晶硅涂层上涂光刻胶；
A4、在光刻胶上进行板印，按照指定图形制作遮盖物，将遮盖物放在光刻胶表面，用紫外线照射遮盖物，没有被遮盖的光刻胶在紫外线作用下溶解，最终得到与遮盖物形状相同的光刻胶层；
A5、接触腐蚀至ITO玻璃基质；
A6、去除光刻胶；
A7、在ITO玻璃基质上未覆盖氢化非晶硅涂层的区域涂导电粘合剂；
所述步骤D中，根据细胞的刚体模型，细胞的三维图像生成归结为模型参数的最大似然估计，使用机器学习算法对模型参数进行最大似然估计，按下式最小化的参数取值：
<mrow>
<mi>P</mi>
<mo>=</mo>
<munder>
<mrow>
<mi>arg</mi>
<mi> </mi>
<mi>min</mi>
</mrow>
<mi>p</mi>
</munder>
<munder>
<mo>&Sigma;</mo>
<mi>i</mi>
</munder>
<mo>|</mo>
<mo>|</mo>
<msub>
<mi>I</mi>
<mrow>
<mi>i</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>-</mo>
<mi>f</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>R</mi>
<mi>i</mi>
</msub>
<msub>
<mi>M</mi>
<mi>i</mi>
</msub>
<mo>+</mo>
<msub>
<mi>T</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>|</mo>
<msup>
<mo>|</mo>
<mn>2</mn>
</msup>
<mo>;</mo>
</mrow>
其中Ii+1是细胞旋转图像序列I＝{Ii,i＝1,…,n}中的一帧，M代表细胞上所有点的集合，Ri和Ti分别是旋转矩阵和平移向量,Rimi+Ti是细胞上一点m从i时刻到i+1时刻的三维旋转运动模型，即mi+1，映射是将细胞在i时刻的三维信息投影到细胞旋转图像的某一帧Ii。
2.根据权利要求1所述的基于光流分析的单细胞三维图像生成方法，其特征在于，所述细胞在非均匀电场中的所受到的平均介电泳力用如下公式描述：
<mrow>
<msub>
<mover>
<mi>F</mi>
<mo>&RightArrow;</mo>
</mover>
<mrow>
<mi>D</mi>
<mi>E</mi>
<mi>P</mi>
</mrow>
</msub>
<mo>=</mo>
<mn>2</mn>
<msup>
<mi>&pi;R</mi>
<mn>3</mn>
</msup>
<msub>
<mi>&epsiv;</mi>
<mi>m</mi>
</msub>
<mi>Re</mi>
<mo>&lsqb;</mo>
<msub>
<mi>f</mi>
<mrow>
<mi>C</mi>
<mi>M</mi>
</mrow>
</msub>
<mo>&rsqb;</mo>
<mover>
<mo>&dtri;</mo>
<mo>&RightArrow;</mo>
</mover>
<mrow>
<mo>(</mo>
<msubsup>
<mi>E</mi>
<mrow>
<mi>r</mi>
<mi>m</mi>
<mi>s</mi>
</mrow>
<mn>2</mn>
</msubsup>
<mo>)</mo>
</mrow>
</mrow>
其中FDEP是作用到细胞上的平均介电泳力，R是细胞的半径，εm是细胞所在溶液的介电常数，Erms为所施加交流信号的均方根值，fCM为Clausius-Mossotti因子，在计算平均介电泳力时取该因子的实部Re[fCM]。
3.根据权利要求1所述的基于光流分析的单细胞三维图像生成方法，其特征在于，fCM因子定义如下:
<mrow>
<msub>
<mi>f</mi>
<mrow>
<mi>C</mi>
<mi>M</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msubsup>
<mi>&epsiv;</mi>
<mi>p</mi>
<mo>*</mo>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>&epsiv;</mi>
<mi>m</mi>
<mo>*</mo>
</msubsup>
</mrow>
<mrow>
<msubsup>
<mi>&epsiv;</mi>
<mi>p</mi>
<mo>*</mo>
</msubsup>
<mo>+</mo>
<mn>2</mn>
<msubsup>
<mi>&epsiv;</mi>
<mi>m</mi>
<mo>*</mo>
</msubsup>
</mrow>
</mfrac>
</mrow>
εp*和εm*分别是细胞和溶液的复介电常数。
4.根据权利要求1所述的基于光流分析的单细胞三维图像生成方法，其特征在于，所述复介电常数可表示为:
<mrow>
<msup>
<mi>&epsiv;</mi>
<mo>*</mo>
</msup>
<mo>=</mo>
<mi>&epsiv;</mi>
<mo>-</mo>
<mi>j</mi>
<mfrac>
<mi>&sigma;</mi>
<mi>&omega;</mi>
</mfrac>
</mrow>
其中，ε是溶液的介电常数，σ是导电率，ω是所施加交流信号的频率。
5.根据权利要求1所述的基于光流分析的单细胞三维图像生成方法，其特征在于，细胞旋转速度为：
<mrow>
<mi>R</mi>
<mrow>
<mo>(</mo>
<mi>&Omega;</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mo>-</mo>
<mfrac>
<mrow>
<msub>
<mi>&epsiv;</mi>
<mi>m</mi>
</msub>
<mi>Im</mi>
<mo>&lsqb;</mo>
<msub>
<mi>f</mi>
<mrow>
<mi>C</mi>
<mi>M</mi>
</mrow>
</msub>
<mo>&rsqb;</mo>
<mo>|</mo>
<mi>E</mi>
<msup>
<mo>|</mo>
<mn>2</mn>
</msup>
</mrow>
<mrow>
<mn>2</mn>
<mi>&eta;</mi>
</mrow>
</mfrac>
<mi>K</mi>
</mrow>
其中E是电场强度，η是溶液的黏稠度，IM[fCM]是Clausius-Mossotti因子的虚部，K为系数。
6.根据权利要求1所述的基于光流分析的单细胞三维图像生成方法，其特征在于，所述预处理包括：高斯滤波处理、亮度调整及模板匹配。
Priority Applications (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510500965.1A CN105184853B (zh) | 2015-08-14 | 2015-08-14 | 一种基于光流分析的单细胞三维图像生成方法 |
| PCT/CN2015/088945 WO2017028341A1 (zh) | 2015-08-14 | 2015-09-06 | 一种基于光流分析的单细胞三维图像生成方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510500965.1A CN105184853B (zh) | 2015-08-14 | 2015-08-14 | 一种基于光流分析的单细胞三维图像生成方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105184853A CN105184853A (zh) | 2015-12-23 |
| CN105184853B true CN105184853B (zh) | 2018-04-10 |
Family
ID=54906903
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510500965.1A Active CN105184853B (zh) | 2015-08-14 | 2015-08-14 | 一种基于光流分析的单细胞三维图像生成方法 |
Country Status (2)
| Country | Link |
|---|---|
| CN (1) | CN105184853B (zh) |
| WO (1) | WO2017028341A1 (zh) |
Families Citing this family (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105092679B (zh) * | 2015-08-14 | 2018-07-03 | 深圳大学 | 一种基于光诱导介电泳技术的单细胞控制方法 |
| DE102017201252A1 (de) * | 2017-01-26 | 2018-07-26 | Universität Ulm | Verfahren und Vorrichtung zur Untersuchung von Zellen |
| CN107858289B (zh) * | 2017-12-25 | 2019-06-21 | 上海大学 | 一种细胞划痕芯片、装置及方法 |
| EP3756165A4 (en) * | 2018-04-27 | 2021-10-13 | Hewlett-Packard Development Company, L.P. | OBJECT ROTATION BY NON-UNIFORM NON-ROTATING ELECTRIC FIELD |
| CN109925597B (zh) * | 2019-02-01 | 2023-06-09 | 广州唯思冠电子科技有限公司 | 一种基于衡通仪的细胞呈现方法 |
| CN110314714B (zh) * | 2019-07-09 | 2021-06-25 | 大连海事大学 | 一种基于三维图像特征的细胞活性状态表征监测装置及方法 |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1668913A (zh) * | 2002-05-14 | 2005-09-14 | 阿默森生物科学英国有限公司 | 评价生物膜的方法 |
| CN101135680A (zh) * | 2007-07-13 | 2008-03-05 | 东南大学 | 光诱导介电泳辅助单细胞介电谱自动测试装置及测试方法 |
| CN102037343A (zh) * | 2008-06-12 | 2011-04-27 | 东卡莱罗纳大学 | 用于三维衍射成像的流式细胞仪系统及方法 |
| CN102692416A (zh) * | 2012-06-26 | 2012-09-26 | 天津师范大学 | 基于微操作机器人的胚胎细胞迁徙自动跟踪系统及方法 |
| CN104212706A (zh) * | 2014-09-19 | 2014-12-17 | 成都劲宏科技有限公司 | 一种基于锁相放大技术的细胞微流图像采集系统 |
| CN104567682A (zh) * | 2015-01-14 | 2015-04-29 | 天津大学 | 液态环境下微粒三维位置纳米级分辨力测量方法 |
Family Cites Families (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7686934B2 (en) * | 2005-12-13 | 2010-03-30 | Gamida For Life B.V. | Three dimensional dielectrophoretic separator and methods of use |
| WO2008051248A2 (en) * | 2005-12-21 | 2008-05-02 | Trustees Of The University Of Pennsylvania | System and method for controlling nanoparticles using dielectrophoretic forces |
| WO2010115167A2 (en) * | 2009-04-03 | 2010-10-07 | The Regents Of The University Of California | Methods and devices for sorting cells and other biological particulates |
| CN101697236B (zh) * | 2009-10-21 | 2012-09-19 | 南昌航空大学 | 基于智能优化算法的直线光流场三维重建方法 |
| JP5617532B2 (ja) * | 2010-10-29 | 2014-11-05 | ソニー株式会社 | 誘電サイトメトリ装置及び誘電サイトメトリによる細胞分取方法 |
-
2015
- 2015-08-14 CN CN201510500965.1A patent/CN105184853B/zh active Active
- 2015-09-06 WO PCT/CN2015/088945 patent/WO2017028341A1/zh active Application Filing
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1668913A (zh) * | 2002-05-14 | 2005-09-14 | 阿默森生物科学英国有限公司 | 评价生物膜的方法 |
| CN101135680A (zh) * | 2007-07-13 | 2008-03-05 | 东南大学 | 光诱导介电泳辅助单细胞介电谱自动测试装置及测试方法 |
| CN102037343A (zh) * | 2008-06-12 | 2011-04-27 | 东卡莱罗纳大学 | 用于三维衍射成像的流式细胞仪系统及方法 |
| CN102692416A (zh) * | 2012-06-26 | 2012-09-26 | 天津师范大学 | 基于微操作机器人的胚胎细胞迁徙自动跟踪系统及方法 |
| CN104212706A (zh) * | 2014-09-19 | 2014-12-17 | 成都劲宏科技有限公司 | 一种基于锁相放大技术的细胞微流图像采集系统 |
| CN104567682A (zh) * | 2015-01-14 | 2015-04-29 | 天津大学 | 液态环境下微粒三维位置纳米级分辨力测量方法 |
Non-Patent Citations (3)
| Title |
|---|
| Dielectrophoresis and electrorotation of neurospora slime and murine myeloma cells;Jan Gimsa等;《Biophys.J.Biophysical Society》;19911031;第60卷;第749-760页 * |
| Image-driven cell manipulation;Yen-Heng Lin等;《IEEE Nanotechnology Magazine》;20091231;第3卷(第3期);第6-11页 * |
| Micro Vision Based Cell Motility Analyzing Algorithm by Optically-Induced Dielectrophoresis;Guanglie Zhang等;《Proceedings of the 2012 IEEE International Conference on Robotics and Biomimetics》;20121214;第1779-1781页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| WO2017028341A1 (zh) | 2017-02-23 |
| CN105184853A (zh) | 2015-12-23 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105184853B (zh) | 一种基于光流分析的单细胞三维图像生成方法 | |
| CN105092679B (zh) | 一种基于光诱导介电泳技术的单细胞控制方法 | |
| EP1413911B1 (en) | Method and device for 3 dimensional imaging of suspended micro-objects providing high-resolution microscopy | |
| Backer et al. | A bisected pupil for studying single-molecule orientational dynamics and its application to three-dimensional super-resolution microscopy | |
| Wang et al. | Tracking E. coli runs and tumbles with scattering solutions and digital holographic microscopy | |
| CN102792151B (zh) | 用于2d和3d成像的超分辨率光流体显微镜 | |
| CN107966801A (zh) | 一种基于环形照明的高速傅立叶叠层成像装置及重构方法 | |
| CN105044192B (zh) | 一种基于光诱导介电泳技术的细胞分类方法 | |
| Peng et al. | Algorithms for on-line monitoring of micro spheres in an optical tweezers-based assembly cell | |
| CN103674977B (zh) | 大口径光学元件表面损伤检测装置和相应的检测方法 | |
| CN104441650B (zh) | 立体打印系统以及立体打印方法 | |
| JP2007537729A (ja) | 微粒子及び細胞の操作用の光電子ピンセット | |
| CN103025888B (zh) | 用于检测生物颗粒群的方法 | |
| JP2013517510A (ja) | トモグラフィー光照射野顕微鏡 | |
| CN103809267B (zh) | 投影镜头、使用其的投影装置及光驱动微粒子装置 | |
| Zhang et al. | Data-driven microscopic pose and depth estimation for optical microrobot manipulation | |
| Humbert et al. | Versatile, rapid and robust nano-positioning of single-photon emitters by AFM-nanoxerography | |
| CN104614078A (zh) | 压电摆镜微扫描显微热成像系统 | |
| Kelbauskas et al. | Optical computed tomography for spatially isotropic four-dimensional imaging of live single cells | |
| Hoshino et al. | Electron beam induced fine virtual electrode for mechanical strain microscopy of living cell | |
| Gan et al. | Non-contact massively parallel manipulation of micro-objects by optoelectronic tweezers | |
| Gong et al. | A fully water coupled oblique light-sheet microscope | |
| CN215050149U (zh) | 一种基于微流控芯片的细胞分选与识别系统 | |
| CN106908831A (zh) | 基于微结构阵列的粒子径迹探测装置 | |
| Kim et al. | On-machine measurement for surface flatness of transparent and thin film in laser ablation process |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |