CN106097269B - 基于计算鬼成像的高分辨率微视觉系统及获取图像的方法 - Google Patents

基于计算鬼成像的高分辨率微视觉系统及获取图像的方法 Download PDF

Info

Publication number
CN106097269B
CN106097269B CN201610409497.1A CN201610409497A CN106097269B CN 106097269 B CN106097269 B CN 106097269B CN 201610409497 A CN201610409497 A CN 201610409497A CN 106097269 B CN106097269 B CN 106097269B
Authority
CN
China
Prior art keywords
spatial light
light modulator
ccd camera
laser beam
computer
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
Application number
CN201610409497.1A
Other languages
English (en)
Other versions
CN106097269A (zh
Inventor
张宪民
吴衡
李海
甘金强
詹镇辉
何振亚
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
South China University of Technology SCUT
Original Assignee
South China University of Technology SCUT
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by South China University of Technology SCUT filed Critical South China University of Technology SCUT
Priority to CN201610409497.1A priority Critical patent/CN106097269B/zh
Publication of CN106097269A publication Critical patent/CN106097269A/zh
Application granted granted Critical
Publication of CN106097269B publication Critical patent/CN106097269B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G06T5/92
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T3/00Geometric image transformation in the plane of the image
    • G06T3/40Scaling the whole image or part thereof
    • G06T3/4053Super resolution, i.e. output image resolution higher than sensor resolution
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10056Microscopic image

Abstract

本发明公开基于计算鬼成像的高分辨率微视觉系统及获取图像的方法,系统在光路上依次包括激光光源、第一光阑、激光扩束镜、准直透镜、第二光阑、起偏器、空间光调制器、检偏器、第三光阑、反光镜、分束镜、会聚透镜、CCD相机;还包括位于分束镜另一光路上的精密定位载物台;精密定位载物台与计算机连接,计算机还分别与空间光调制器与CCD相机连接,计算机通过计算鬼成像技术获取高分辨率图像。本发明结构简单紧凑,由于采用光场强度关联测量恢复物体信息,消除了经典光学系统成像畸变的问题,可以获得高准确度和对比度的图像。本发明非常有利于微视觉系统的设计及鬼成像技术的研究。

Description

基于计算鬼成像的高分辨率微视觉系统及获取图像的方法
技术领域
本发明涉及计算机微视觉领域,具体涉及一种基于计算鬼成像的高分辨率微视觉系统。
背景技术
计算机微视觉系统是集光学显微镜、视觉成像和计算机视觉技术为一体的可实现实时、可视化检测的测量平台。微视觉系统的构成主要包括光学显微镜、光源、摄像机、图像采集卡、精密定位载物台等硬件以及图像处理软件。其原理是通过显微镜和成像设备(CCD摄像机、图像采集卡等)把被测对象的图像采集到计算机,接着运用图像处理技术、计算机视觉或人工智能等技术对采集到的图像进行处理、识别等操作,从而完成微视觉系统所要求的任务。这种微视觉系统在微观测量、成像等领域有着广泛的应用。
随着科技的不断发展,人们对微观世界的研究已经由微米级、亚微米级进入到纳米级阶段,光学分辨极限的限制逐渐凸现出来,极大地限制了计算机微视觉技术的进一步应用。虽然目前已有多种超分辨成像技术被实现,但是这些方法的成功实施往往需要特殊的条件及硬件支持,使用范围受到一定的限制。面对这些问题,各种新的解决方案被提出来。一方面,直接发展新的显微成像技术,例如借助短波长电子的扫描电子显微镜,原子力显微镜,荧光显微镜等等;另一方面,研究可以突破衍射极限的光学显微成像方法,如量子成像,热光鬼成像,结构光成像技术等等。
鬼成像(ghost imaging),是一种利用双光子符合探测恢复待测物体空间信息的新型成像技术。传统的光学是基于光场的一阶关联(强度与位相)获得信息,而鬼成像利用的光场的二阶或者高阶关联,并结合符合测量技术获取成像信息。鬼成像可以实现非定域成像,无透镜成像,抗大气湍流成像等成像方案,引起了人们的广泛关注。经典成像系统的分辨率受到光学衍射极限的制约,而鬼成像技术具有超越经典分辨率极限的能力,尤其是近年来出现的计算鬼成像,对鬼成像技术的实用化有巨大的推进作用。因此将计算鬼成像技术应用到微视觉系统具有十分重要的意义。
发明内容
针对经典微视觉系统的分辨率受到光学衍射极限制约的问题,本发明提供了一种基于计算鬼成像的高分辨率微视觉系统。该系统结构紧凑,安装方便,抗干扰能力强,结合计算鬼成像技术,能够突破经典光学系统的衍射极限,使系统分辨率不受透镜孔径尺寸限制,成像分辨率和对比度高。
本发明的目的通过如下技术方案实现。
一种基于计算鬼成像的高分辨率微视觉系统,在光路上依次包括激光光源、第一光阑、激光扩束镜、准直透镜、第二光阑、起偏器、空间光调制器、检偏器、第三光阑、反光镜、分束镜、会聚透镜、CCD相机;还包括位于分束镜另一光路上的精密定位载物台;精密定位载物台与计算机连接,计算机还分别与空间光调制器与CCD相机连接,计算机通过计算鬼成像技术获取高分辨率图像。
所述的利用计算鬼成像技术获取图像,其实现步骤如下:
1.利用空间光调制器对激光光强进行调制。
通过计算机生成K张M×M的随机散斑图,散斑图的中心是一个N×N有效散斑区域,且N≤M,有效散斑区域的周围区域为白色,然后将散斑图转换为全息图并存储到硬盘。取一张全息图加载到空间光调制器上,并调节激光源、光阑、激光扩束镜以及准直透镜,使激光扩束镜产生的光斑可以完全覆盖加载到空间光调制器上有效全息区域(对应散斑图上的有效散斑区域)。通过不断的加载新的全息图,可以实现对激光光强的调制。
2.利用CCD相机采集物体表面的光强强度变化。
调节精密控制载物台,确保被测物体的目标区域完全被激光束覆盖。接着调整会聚透镜以及CCD相机的位置,使CCD相机能够接收到被测物体目标区域的反射光。然后通过计算机控制空间光调制器与CCD相机同步工作,即空间光调制器每加载一副全息图,CCD相机就立即拍下被测物体目标区域的光强变化,并将对应的图片保存。将得到的图片的所有像素的灰度值累加,记为Bi,i表示测量的次数,由此可得到测试光路的光强涨落信息。
3.通过计算机模拟获取参考光路光强涨落信息。
激光束未经空间光调制器调制时,在空间光调制器平面处的场强为Es(xs,ys);激光束经过空间光调制器调制后,其场强为
Eo(x,y)=Es(xs,ys)Em(x,y)
上式中Em(x,y)表示用于调制的场强。
激光束经过空间光调制器调制后,在参考光路CCD相机处的场强为
上式中(x,y),(xs,ys)分别表示CCD相机平面、空间光调制器平面的直角坐标;Dr表示CCD相机到空间光调制器的距离;λ为激光的波长;Eo(x,y)表示激光束经空间光调制器调制后的场强。
由上式可以得到CCD相机处的光强为
Ir(x,y)=Er(x,y)Er *(x,y)
4.进行强度关联运算,获取被测物体图像。
将2和3中得到的测试光路和参考光路的光强涨落信息进行关联,即
上式中N表示测量次数。
对G(x,y)进行归一化,即
Gfinal(x,y)=G(x,y)/max(G(x,y))
max(G(x,y))表示取G(x,y)中的最大值。
由上式可以得到被测物体目标区域的图像信息。
采用上述技术方案后,可以设计出一个结构紧凑、安装方便的微视觉系统。结合计算鬼成像技术,能够突破经典光学系统的衍射极限,获得高分辨率和对比度的图像。
与现有技术相比,本发明的有益效果及优点:本发明设计出一个结构简单、安装方便、易于操作的微视觉系统。通过采用基于光场强度关联测量恢复物体信息的鬼成像技术,可以克服经典光学系统普遍存在的衍射极限问题,实现高分辨率成像。同时,由于采用计算鬼成像技术,相比于传统的鬼成像技术,系统的结构得到了简化,实用性更强。另外,由于采用光场强度关联测量恢复物体信息,消除了经典光学系统成像畸变的问题,可以获得高准确度和对比度的图像。本发明非常有利于微视觉系统的设计及鬼成像技术的研究。
附图说明
图1为实施方式中微视觉系统的组成示意图。
图2为实施方式中计算机模拟散斑的示意图。
图3为实施方式中激光束与散斑位置关系的示意图。
具体实施方式
下面结合附图和实施例对本发明的内容进行详细的描述,但本发明的实际应用形式并不仅限于下述的实施例。
如图1所示,本发明提供了一种基于计算鬼成像的高分辨率微视觉系统,该系统由激光光源101、光阑(102,105,109)、激光扩束镜103、准直透镜104、起偏器106、空间光调制器107、检偏器108、反光镜110、会聚透镜111、CCD相机112、分束镜113、计算机114、精密定位载物台115组成。所述的系统利用计算鬼成像技术获取高分辨率图像。
所述的利用计算鬼成像技术获取图像,其实现步骤如下:
1.利用空间光调制器对激光光强进行调制。
通过计算机生成8000张900×900的随机散斑图,散斑图的中心是一个360×360有效散斑区域201,有效散斑区域的周围区域202为白色,如图2所示。然后将散斑图转换为全息图并存储到硬盘。取一张全息图加载到空间光调制器107上,并调节激光源101、光阑(102, 105)、激光扩束镜103以及准直透镜104,使激光扩束镜产生的光斑301可以完全覆盖加载到空间光调制器上有效全息区域(对应散斑图上的有效散斑区域),如图3所示。通过不断的加载新的全息图,可以实现对激光光强的调制。
2.利用CCD相机采集物体表面的光强强度变化。
调节精密定位载物台115,确保被测物体的目标区域完全被激光束覆盖。接着调整会聚透镜111以及CCD相机112的位置,使CCD相机能够接收到被测物体目标区域的反射光。然后通过计算机控制空间光调制器与CCD相机同步工作,即空间光调制器每加载一副全息图,CCD相机就立即拍下被测物体目标区域的光强变化,并将对应的图片保存。将所得到的图片的所有像素的灰度值累加,记为Bi,i表示测量的次数,由此可得到测试光路的光强涨落信息。
3.通过计算机模拟获取参考光路光强涨落信息。
激光束未经空间光调制器调制时,在空间光调制器平面处的场强为Es(xs,ys);激光束经过空间光调制器调制后,其场强为
Eo(x,y)=Es(xs,ys)Em(x,y)
上式中Em(x,y)表示用于调制的场强;
激光束经过空间光调制器调制后,在参考光路CCD相机处的场强为
上式中(x,y),(xs,ys)分别表示CCD相机平面、空间光调制器平面的直角坐标; Dr=800cm表示CCD相机到空间光调制器的距离;λ=635nm为激光的波长;Es(xs,ys) 表示激光束在空间光调制器平面处的场强。
由上式可以得到CCD相机平面处的光强为
Ir(x,y)=Er(x,y)Er *(x,y)
上式中Er(x,y)表示激光束经过空间光调制器调制后在参考光路CCD相机处的场强, Er *(x,y)表示Er(x,y)的共轭函数。
4.进行强度关联运算,获取被测物体图像。
将2和3中得到的测试光路和参考光路的光强涨落信息进行关联,即
上式中 需结合上标进行定义说明表示第i次运算得到的参考光路CCD相机处的光强。
对G(x,y)进行归一化,即
Gfinal(x,y)=G(x,y)/max(G(x,y))
max(G(x,y))表示取G(x,y)中的最大值。
由上式可以得到被测物体目标区域的图像信息,即Gfinal(x,y)。
图1为实施方式中微视觉系统的组成示意图。其中包括激光光源101、光阑(102,105, 109)、激光扩束镜103、准直透镜104、起偏器106、空间光调制器107、检偏器108、反光镜110、会聚透镜111、CCD相机112、分束镜113、计算机114、精密定位载物台115。
图2为实施方式中计算机模拟散斑的示意图,其中201为所述的有效散斑区域,202为有效散斑区域的周围区域。
图3为实施方式中激光束与散斑位置关系的示意图,其中301为激光扩束镜产生的光斑。
本发明结合计算鬼成像技术,能够突破经典光学系统的衍射极限,获得高分辨率和对比度的图像。

Claims (2)

1.一种基于计算鬼成像的高分辨率微视觉系统获取图像的方法,所述高分辨率微视觉系统在光路上依次包括激光光源、第一光阑、激光扩束镜、准直透镜、第二光阑、起偏器、空间光调制器、检偏器、第三光阑、反光镜、分束镜、会聚透镜、CCD相机;还包括位于分束镜另一光路上的精密定位载物台;其特征是包括如下步骤:
(1)利用空间光调制器对激光光强进行调制;
通过计算机生成K张M×M的随机散斑图,K、M、N为正整数,散斑图的中心是一个N×N有效散斑区域,且N≤M,有效散斑区域的周围区域为白色,然后将散斑图转换为全息图并存储到硬盘;取一张全息图加载到空间光调制器上,并调节激光源、第一~第三光阑、激光扩束镜以及准直透镜,使激光扩束镜产生的光斑能完全覆盖加载到空间光调制器上有效全息区域即对应散斑图上的有效散斑区域;通过不断的加载新的全息图,能实现对激光光强的调制;
(2)利用CCD相机采集物体表面的光强强度变化;
调节精密定位载物台,使被测物体的目标区域完全被激光束覆盖;接着调整会聚透镜以及CCD相机的位置,使CCD相机能够接收到被测物体目标区域的反射光;然后通过计算机控制空间光调制器与CCD相机同步工作,即空间光调制器每加载一副全息图,CCD相机就立即拍下被测物体目标区域的光强变化,并将对应的图片保存;将得到的图片的所有像素的灰度值累加,记为Bi,i表示测量的次数,由此可得到测试光路的光强涨落信息;
(3)通过计算机模拟获取参考光路光强涨落信息;激光束未经空间光调制器调制时,在空间光调制器平面处的场强为Es(xs,ys);激光束经过空间光调制器调制后,其场强为
Eo(x,y)=Es(xs,ys)Em(x,y)
上式中Em(x,y)表示用于调制的场强;
激光束经过空间光调制器调制后,在参考光路CCD相机处的场强为
上式中(x,y),(xs,ys)分别表示CCD相机平面、空间光调制器平面的直角坐标;Dr表示CCD相机到空间光调制器的距离;λ为激光的波长;Eo(x,y)表示激光束经空间光调制器调制后的场强;
由上式可以得到CCD相机处的光强为
Ir(x,y)=Er(x,y)Er *(x,y)
上式中Er(x,y)表示激光束经过空间光调制器调制后在参考光路CCD相机平面处的场强,Er *(x,y)表示Er(x,y)的共轭函数;
(4)进行强度关联运算,获取被测物体图像;
将(2)和(3)中得到的测试光路和参考光路的光强涨落信息进行关联,即
上式中N表示测量次数;表示第i次测量运算得到的参考光路CCD相机处的光强;
对G(x,y)进行归一化,即
Gfinal(x,y)=G(x,y)/max(G(x,y))
max(G(x,y))表示取G(x,y)中的最大值;
由上式可以得到被测物体目标区域的图像信息,即Gfinal(x,y)。
2.根据权利要求1所述的基于计算鬼成像的高分辨率微视觉系统获取图像的方法,其特征在于所述高分辨率微视觉系统还包括计算机,精密定位载物台与计算机连接,计算机还分别与空间光调制器与CCD相机连接,计算机通过计算鬼成像技术获取高分辨率图像。
CN201610409497.1A 2016-06-11 2016-06-11 基于计算鬼成像的高分辨率微视觉系统及获取图像的方法 Active CN106097269B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610409497.1A CN106097269B (zh) 2016-06-11 2016-06-11 基于计算鬼成像的高分辨率微视觉系统及获取图像的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610409497.1A CN106097269B (zh) 2016-06-11 2016-06-11 基于计算鬼成像的高分辨率微视觉系统及获取图像的方法

Publications (2)

Publication Number Publication Date
CN106097269A CN106097269A (zh) 2016-11-09
CN106097269B true CN106097269B (zh) 2019-04-09

Family

ID=57227979

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610409497.1A Active CN106097269B (zh) 2016-06-11 2016-06-11 基于计算鬼成像的高分辨率微视觉系统及获取图像的方法

Country Status (1)

Country Link
CN (1) CN106097269B (zh)

Families Citing this family (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106483105B (zh) * 2016-11-18 2023-06-27 华南理工大学 基于强度关联成像的透射式微视觉系统及图像获取方法
CN106772310B (zh) * 2016-12-30 2019-05-07 南京理工大学 一种基于计算鬼成像的目标鉴定装置及其方法
CN107328365B (zh) * 2017-06-01 2019-05-14 华南理工大学 一种柔顺定位平台位移的鲁棒测量方法
CN108901118A (zh) * 2018-07-04 2018-11-27 北京大学 一种激光离子加速器中透明靶体定位系统及其定位方法
CN109151251A (zh) * 2018-08-16 2019-01-04 上海理工大学 基于深度学习和关联成像的双密钥图像传输方法
CN109520619B (zh) * 2018-11-26 2021-03-02 中国科学院上海光学精密机械研究所 基于非瑞利散斑场的关联成像光谱相机及其成像方法
CN109859135B (zh) * 2019-01-31 2021-05-07 北京邮电大学 一种应用于关联成像的图像增强处理方法
CN110230995B (zh) * 2019-05-10 2021-01-05 首都师范大学 一种基于鬼成像的感兴趣区域成像装置
CN110244471B (zh) * 2019-06-18 2021-07-16 华东交通大学 基于柯勒照明的计算鬼成像系统及方法
CN110208294A (zh) * 2019-06-18 2019-09-06 华东交通大学 基于柯勒照明的单像素显微成像方法及系统
CN112802145A (zh) * 2021-01-27 2021-05-14 四川大学 一种基于深度学习的彩色计算鬼成像方法
CN115330756B (zh) * 2022-10-11 2023-02-28 天津恒宇医疗科技有限公司 基于光影特征的oct图像中导丝识别方法及识别系统

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101251655A (zh) * 2008-03-21 2008-08-27 北京理工大学 一种实现光束轨道角动量态叠加和调制的装置
CN101556386A (zh) * 2009-05-15 2009-10-14 哈尔滨理工大学 一种液晶空间光调制器多参数的干涉式双成像测量装置
CN101726868A (zh) * 2009-12-15 2010-06-09 北京理工大学 一种实现光束轨道角动量态复用编码的方法和装置
CN103954367A (zh) * 2014-04-18 2014-07-30 河南科技大学 一种测量分数阶光学涡旋拓扑荷值的装置及其测量方法
US8811768B2 (en) * 2007-12-06 2014-08-19 The United States Of America As Represented By The Secretary Of The Army Image enhancement system and method
CN105444896A (zh) * 2015-11-30 2016-03-30 河南科技大学 基于六角星孔衍射的涡旋光束拓扑荷值测量方法
CN206224605U (zh) * 2016-06-11 2017-06-06 华南理工大学 基于计算鬼成像的高分辨率微视觉系统

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8811768B2 (en) * 2007-12-06 2014-08-19 The United States Of America As Represented By The Secretary Of The Army Image enhancement system and method
CN101251655A (zh) * 2008-03-21 2008-08-27 北京理工大学 一种实现光束轨道角动量态叠加和调制的装置
CN101556386A (zh) * 2009-05-15 2009-10-14 哈尔滨理工大学 一种液晶空间光调制器多参数的干涉式双成像测量装置
CN101726868A (zh) * 2009-12-15 2010-06-09 北京理工大学 一种实现光束轨道角动量态复用编码的方法和装置
CN103954367A (zh) * 2014-04-18 2014-07-30 河南科技大学 一种测量分数阶光学涡旋拓扑荷值的装置及其测量方法
CN105444896A (zh) * 2015-11-30 2016-03-30 河南科技大学 基于六角星孔衍射的涡旋光束拓扑荷值测量方法
CN206224605U (zh) * 2016-06-11 2017-06-06 华南理工大学 基于计算鬼成像的高分辨率微视觉系统

Also Published As

Publication number Publication date
CN106097269A (zh) 2016-11-09

Similar Documents

Publication Publication Date Title
CN106097269B (zh) 基于计算鬼成像的高分辨率微视觉系统及获取图像的方法
CN106483105B (zh) 基于强度关联成像的透射式微视觉系统及图像获取方法
US8570650B2 (en) Method and system for fast three-dimensional structured-illumination-microscopy imaging
CN105659143B (zh) 孔径扫描傅立叶重叠关联成像
US7977625B2 (en) Method and assembly for optical reproduction with depth discrimination
CN105784653B (zh) 一种宽场超分辨荧光显微成像装置
US8675062B2 (en) Shape measuring device, observation device, and image processing method
CN107490562A (zh) 利用波面整形器的超高速三维折射率影像拍摄和荧光结构光照明显微镜系统及其使用方法
Pushkina et al. Superresolution linear optical imaging in the far field
CN104111590B (zh) 基于复合涡旋双瓣聚焦光斑的激光直写装置
CN104089573B (zh) 基于正交偏振光的多通道白光共路干涉显微层析系统
CN109712072A (zh) 基于全内反射的条纹照明傅里叶域迭代更新超分辨显微成像方法
CN116183568B (zh) 一种三维结构光照明超分辨显微成像的高保真重构的方法和装置
WO2024051079A1 (zh) 一种主动结构光照明的超分辨显微成像方法及系统
CN108895986A (zh) 基于条纹成像投影的显微三维形貌测量装置
CN110494794B (zh) 伪散斑图案生成装置、伪散斑图案生成方法、观察装置和观察方法
US11430144B2 (en) Device and process for the contemporary capture of standard images and plenoptic images via correlation plenoptic imaging
CN108254340A (zh) 基于线偏振调制的扫描显微镜
CN100565142C (zh) 光束测量装置
Chen et al. Enhancing speckle statistics for imaging inside scattering media
Wu et al. Autofocusing algorithm for pixel-super-resolved lensfree on-chip microscopy
KR101888924B1 (ko) 디지털 마이크로미러 소자와 시간 복합 구조화 조명을 이용한 구조화 조명 현미경 시스템 및 그 동작 방법
CN206224605U (zh) 基于计算鬼成像的高分辨率微视觉系统
CN109974578A (zh) 一种基于双液晶空间光调制器的涡旋数字全息显微系统
CN105467571B (zh) 一种白光同轴干涉相移显微镜系统和显微成像方法

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