[发明专利]基于希尔伯特变换的结构照明快速三维彩色显微成像方法有效
申请号: | 202010061033.2 | 申请日: | 2020-01-19 |
公开(公告)号: | CN111123496B | 公开(公告)日: | 2020-11-06 |
发明(设计)人: | 雷铭 | 申请(专利权)人: | 西安交通大学 |
主分类号: | G02B21/06 | 分类号: | G02B21/06;G02B21/36;G06T7/90 |
代理公司: | 北京市诚辉律师事务所 11430 | 代理人: | 范盈 |
地址: | 710049 陕西省西*** | 国省代码: | 陕西;61 |
权利要求书: | 查看更多 | 说明书: | 查看更多 |
摘要: | |||
搜索关键词: | 基于 希尔伯特 变换 结构 照明 快速 三维 彩色 显微 成像 方法 | ||
1.基于希尔伯特变换的结构照明快速三维彩色显微成像方法,其特征在于包括以下步骤:
1)产生结构照明光场:
由单色LED或白光LED照明,通过数字微镜器件DMD产生2个空间方向相同、相位相差π的任意方向的正弦条纹结构照明光场,结构照明光场经过显微物镜照射载物台上的样品;
2)彩色CCD相机采集图像:
对应于两个相位的结构照明光场,彩色CCD相机分别采集样品当前层的两幅彩色二维图像I0(RGB)、Iπ(RGB);
3)图像处理:
3.1)将步骤2)中采集到的两幅不同相位的彩色二维图像I0(RGB)和Iπ(RGB),按照公式(1)计算得到去除条纹的彩色宽场图像Iwide(RGB);
3.2)将步骤2)中采集的两幅不同相位的彩色二维图像I0(RGB)和Iπ(RGB)转换为灰度图像I0gray(x,y)和Iπgray(x,y),再利用公式(2)将I0gray(x,y)和Iπgray(x,y)写为;
其中,Iin表示二维灰度焦面图像,Iout表示离焦背景,m为正弦条纹的调制度,p为条纹的空间频率;
3.3)计算样品当前层的二维灰度焦面图像Iin(x,y);
3.4)将公式(1)得到的彩色宽场图Iwide(RGB)和当前位置的二维灰度焦面图像Iin(x,y)相乘,即可得到彩色光切片图像Iin(RGB):
Iin(RGB)=Iwide(RGB)·Iin(x,y) (8);
3.5)将样品在轴向以一定层间距移动,然后利用彩色CCD相机在每一层采集两幅彩色二维图像I0(RGB)、Iπ(RGB);
3.6)按照步骤3.1)-3.4)对每一层进行重构,得到每一层对应的彩色光切片图像Iin(RGB);
3.7)将所有层的彩色光切片图像Iin(RGB)进行叠加,得到样品的三维彩色图像Iresult(RGB)。
2.根据权利要求1所述的基于希尔伯特变换的结构照明快速三维彩色显微成像方法,其特征在于:
步骤3.3)具体为:
a)将步骤3.2)得到的两幅灰度图像利用公式(3)相减,得到差值图像ΔI(x,y):
其中,C=2msin(π/2)
b)对差值图像ΔI(x,y)二维矩阵中任意一行进行一维操作:
b1)令y=y0,代入公式(3)计算得到ΔI(x,y0),将ΔI(x,y0)代入公式(4)做一维希尔伯特变换:
其中y0表示差值图像二维矩阵中的任意一行,HT{}表示希尔伯特变换算符,表示卷积算符;
b2)对比公式(3)和公式(4)发现,正弦条纹经过希尔伯特变换后产生了π/2的相移,鉴于此,构造一个复数信号IA(x,y0):
IA(x,y0)=ΔI(x,y0)+i·HTx{ΔI(x,y0)} (5)
取复数信号的模值:
|IA(x,y0)|=C·Iin(x,y0) (6)
c)按照步骤b)的操作方法对差值图像二维矩阵所有行进行一维操作,即获得一幅完整的灰度焦面图像Iin(x,y):
3.根据权利要求1或2所述的基于希尔伯特变换的结构照明快速三维彩色显微成像方法,其特征在于:所述步骤1)提到的任意方向为:X方向或Y方向。
该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于西安交通大学,未经西安交通大学许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服】
本文链接:http://www.vipzhuanli.com/pat/books/202010061033.2/1.html,转载请声明来源钻瓜专利网。