CN108088653B - 共焦显微镜模式像差矫正方法 - Google Patents
共焦显微镜模式像差矫正方法 Download PDFInfo
- Publication number
- CN108088653B CN108088653B CN201711242335.4A CN201711242335A CN108088653B CN 108088653 B CN108088653 B CN 108088653B CN 201711242335 A CN201711242335 A CN 201711242335A CN 108088653 B CN108088653 B CN 108088653B
- Authority
- CN
- China
- Prior art keywords
- aberration
- coefficient
- image
- confocal microscope
- confocal
- 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
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M11/00—Testing of optical apparatus; Testing structures by optical methods not otherwise provided for
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/01—Arrangements or apparatus for facilitating the optical investigation
Landscapes
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Microscoopes, Condenser (AREA)
Abstract
共焦显微镜模式像差矫正方法,属于自适应光学和共焦显微成像技术领域,本发明为了解决现有光学系统的装配误差和光学元件的面形偏差会造成像差,导致焦斑歪曲和分辨率降低的问题。该方法依次加载不同幅值第4‑11阶Zernike项偏置像差,计算相应图像灰度方差函数值,以此利用质心法计算相应像差系数,比较得到像差系数小于0.7rad清零,其余项预先矫正,再重新线性计算像差系数,最后利用两次系数之和进行像差矫正。本发明共焦显微镜模式像差矫正方法,通过空间光调制器反向补偿由光学系统的装配误差与光学元件的面形偏差引起的波前畸变,有效克服共焦显微系统的像差影响,提升其成像质量与分辨率。
Description
技术领域
本发明涉及一种模式像差矫正方法,具体涉及一种利用共焦系统提升共焦成像质量及分辨率的矫正方法,属于自适应光学和共焦显微成像技术领域。
背景技术
共焦显微技术具有高分辨率和可三维深度成像的优点,其高分辨的特性需要焦斑尺寸接近衍射极限附近以获得最佳理论分辨率。然而光学系统的装配误差和光学元件的面形偏差会造成像差,导致焦斑歪曲和分辨率降低。需要动态的纠正这些光学系统所产生的像差,提升成像质量。目前还没有简洁实用的应用于共焦显微系统的像差矫正方法。
发明内容
本发明的目的是提供共焦显微镜模式像差矫正方法,以解决现有光学系统的装配误差和光学元件的面形偏差会造成像差,导致焦斑歪曲和分辨率降低的问题。
所述共焦显微镜模式像差矫正方法所依托的共焦系统,其组成包括激光器、双面45°反光镜、分束镜、半波片、空间光调制器、反射镜、光阑、偏振分束镜、XY扫描振镜、扫描透镜、管镜、物镜、1/4波片、收集透镜、光纤小孔和光电倍增管,激光器发出的激光经双面45°反光镜的一侧反光镜面反射后途经分束镜到达半波片,经半波片调制后射入至空间光调制器,再经空间光调制器调制后射回至分束镜,再经分束镜反射至反射镜,然后经光阑滤除高阶衍射杂光后射到双面45°反光镜的另一侧反光镜面,经双面45°反光镜再次反射后经过偏振分束镜射至XY扫描振镜上,然后经XY扫描振镜反射至扫描透镜,再经过管镜后至1/4波片的偏振面翻转90°被物镜聚焦后打在待测样品上;聚焦后的光斑再从待测样品上反射回物镜,经过1/4波片,偏振面再次翻转90°后依次经过管镜、扫面透镜和XY扫描振镜,射在偏振分束镜上,再经偏振分束镜反射至收集透镜并聚焦到光纤小孔,最后被光电倍增管收集;共焦显微镜模式像差矫正方法,其特征是,包括以下步骤:
步骤a,建立Zernike多项式;将像差分解为Zernike多项式4-11阶组合;Zernike多项式4-11阶多项式形式如下:
其中i代表Zernike像差阶数,作为判定循环是否继续的条件,r为极坐标系径向坐标,θ为极坐标系角向坐标;
步骤c,计算五幅图像的灰度方差函数;
步骤d,s=0:质心法计算像差系数ai,s=1:线性计算像差系数ai’;
步骤e,使i=i+1,若i<12,返回步骤b,若i>12则进行步骤f;
步骤f,若s=0,则使i=4,s=1返回步骤b;若s=1则进行步骤g;
步骤g,计算最佳相差系数,进行共焦像差矫正。
优选的:所述步骤b中将Zernike多项式i阶系数为像差分别转换为bmp图像加载至附加于共焦显微系统的空间光调制器,s=0时,空间光调制器初始值为系数0;s=1时,空间光调制器初始值为清零后的系数a。在空间光调制器5次加载时分别利用共焦系统对样品成像,得到5幅图像作为计算图像灰度方差的数据基础。
优选的:所述步骤c中的图像灰度方差函数为:
其中∑pixels代表对m×n图像像素点值累加求和,<I>为求图像平均值,M为图像灰度方差函数值,m为图像行数,n为图像列数,x为图像中像素点横坐标,y为图像中像素点纵坐标,对步骤b中所得5幅图像分别进行计算,得到5个图像灰度方差值,作为拟合图像像差二项式曲线的数据基础。
6、优选的:所述步骤d中,s=0时,利用5幅图像的灰度方差函数值,做质心法形式如下:
M1-5和x1-5分别对应第一幅至第五幅图像的灰度方差函数值和相应偏执系数,ai为系统像差系数ai负值;对于s=1时,已知若共焦显微镜不含相差,则加入正负相同幅值偏置,其图像灰度方差函数值之差为零,且其差的值与共焦显微镜中所含像差系数成线性比例,利用像差偏置x=±b所对应的两幅图像的灰度方差函数值之差计算共焦显微镜像差系数ai’。
优选的:所述步骤f中,为减少各阶Zernike像差间的串扰,在第一次循环结束后,即i=12时,得到共焦系统的4-11阶像差系数,比较得到的像差系数,其中系数a<0.7rad项像差模式清零,系数a≧0.7rad项像差模式预先矫正,即在空间光调制器中加入清零后的系数a像差的bmp图像为初始值,以此基础加重新进行4-11阶偏置,进入第二次循环:i=4,s=1。
优选的:所述步骤g中,将重新计算得到最佳像差系数a”,其中a”=a’+a,将Φ(r,θ)=∑-ai″Zernikei(r,θ)转换为bmp图像加载至空间光调制器,达到矫正共焦系统像差的效果。
本发明与现有产品相比具有以下效果:本发明根据像差理论,将像差分Zernike多项式4-11阶组合,依靠空间光调制器依次偏置不同幅值不同Zernike阶像差,依靠图像灰度方差函数计算像差系数,利用空间光调制器将其加载至共焦显微系统中,进行像差矫正,可在不增加样品处理操作和装置最简化的情况下有效克服像差对其成像的影响,达到成像质量提升17%,分辨率和图像幅值提升的效果。
附图说明
图1是本发明所述的共焦显微镜模式像差矫正方法的流程图;
图2是共焦系统的结构示意图。
图中:1-激光器、2-双面45°反光镜、3-分束镜、4-半波片、5-空间光调制器、6-反射镜、7-光阑、8-偏振分束镜、9-XY扫描振镜、10-扫描透镜、11-管镜、12-1/4波片、13-物镜、14-待测样品、15-收集透镜、16-光纤小孔、17-光电倍增管。
具体实施方式
下面根据附图详细阐述本发明优选的实施方式。
如图2所示,本发明所述的共焦显微镜模式像差矫正方法所依托的共焦系统,其组成包括激光器1、双面45°反光镜2、分束镜3、半波片4、空间光调制器5、反射镜6、光阑7、偏振分束镜8、XY扫描振镜9、扫描透镜10、管镜11、物镜12、1/4波片13、收集透镜15、光纤小孔16和光电倍增管17,激光器1发出的激光经双面45°反光镜2的一侧反光镜面反射后途经分束镜3到达半波片4,经半波片4调制后射入至空间光调制器5,再经空间光调制器5调制后射回至分束镜3,再经分束镜3反射至反射镜6,然后经光阑7滤除高阶衍射杂光后射到双面45°反光镜2的另一侧反光镜面,经双面45°反光镜2再次反射后经过偏振分束镜8射至XY扫描振镜9上,然后经XY扫描振镜9反射至扫描透镜10,再经过管镜11后经过1/4波片12偏振面翻转90°后被物镜13聚焦至的打在待测样品14上;聚焦后的光斑再从待测样品14上反射回物镜13,经过1/4波片12,偏振面再次翻转90°后依次经过、管镜11、扫面透镜10和XY扫描振镜9,射在偏振分束镜8上,再经偏振分束镜8反射至收集透镜15并聚焦到光纤小孔16,最后被光电倍增管17收集;
所述的激光器1、偏振分束镜8、XY扫描振镜9、扫描透镜10、管镜11、物镜12和1/4波片13组成照明装置,所述的照明装置按照照明光传播方向依次为:激光器1、偏振分束镜8、XY扫描振镜9、扫描透镜10、管镜11、1/4波片12和物镜13。
所述的偏振分束镜8、XY扫描振镜9、扫描透镜10、管镜11、1/4波片12、物镜13、收集透镜15、光纤小孔16和光电倍增管17组成探测装置,所述的探测装置按照照明光传播方向依次为:物镜13、1/4波片12、管镜11、扫描透镜10、XY扫描振镜9、偏振分束镜8、收集透镜15、光纤小孔16和光电倍增管17。
所述的双面45°反光镜2、分束镜3、半波片4、空间光调制器5、反射镜6、光阑7组成共焦系统,所述的共焦系统按照照明光传播方向依次为:双面45°反射镜2、分束镜3、半波片4、空间光调制器5、反射镜6和光阑7。
如图1所示,共焦显微镜模式像差矫正方法,其特征是,包括以下步骤:
步骤a,建立Zernike多项式;将像差分解为Zernike多项式4-11阶组合;Zernike多项式4-11阶多项式形式如下:
其中i代表Zernike像差阶数,作为判定循环是否继续的条件,r、θ分别为极坐标系径向、角向坐标;
步骤c,计算五幅图像的灰度方差函数;
步骤d,s=0:质心法计算像差系数ai,s=1:线性计算像差系数ai’;
步骤e,使i=i+1,若i<12,返回步骤b,若i>12则进行步骤f;
步骤f,若s=0,则使i=4,s=1返回步骤b;若s=1则进行步骤g;
步骤g,计算最佳相差系数,进行共焦像差矫正。
进一步:所述步骤b中将Zernike多项式i阶系数为像差分别转换为bmp图像加载至附加于共焦显微系统的空间光调制器,s=0时,空间光调制器初始值为系数0;s=1时,空间光调制器初始值为清零后的系数a。在空间光调制器5次加载时分别利用共焦系统对样品成像,得到5幅图像作为计算图像灰度方差的数据基础。
进一步:所述步骤c中的图像灰度方差函数为:
其中∑pixels代表对m×n图像像素点值累加求和,<I>为求图像平均值,M为图像灰度方差函数值,m为图像行数,n为图像列数,x为图像中像素点横坐标,y为图像中像素点纵坐标,对步骤b中所得5幅图像分别进行计算,得到5个图像灰度方差值,作为拟合图像像差二项式曲线的数据基础。
7、进一步:所述步骤d中,s=0时,利用5幅图像的灰度方差函数值,做质心法形式如下:
M1-5和x1-5分别对应第一幅至第五幅图像的灰度方差函数值和相应偏执系数,ai为系统像差系数;对于s=1时,已知若共焦显微镜不含相差,则加入正负相同幅值偏置,其图像灰度方差函数值之差为零,且其差的值与共焦显微镜中所含像差系数成线性比例,利用像差偏置x=±b所对应的两幅图像的灰度方差函数值之差计算共焦显微镜像差系数ai’。
进一步:所述步骤f中,为减少各阶Zernike像差间的串扰,在第一次循环结束后,即i=12时,得到共焦系统的4-11阶像差系数,比较得到的像差系数,其中系数ai<0.7rad项像差模式清零,系数ai≧0.7rad项像差模式预先矫正,即在空间光调制器中加入清零后的系数a像差的bmp图像为初始值,以此基础加重新进行4-11阶偏置,进入第二次循环:i=4,s=1。
进一步:所述步骤g中,将重新计算得到最佳像差系数a”,其中a”=a’+a,将Φ(r,θ)=∑-ai″Zernikei(r,θ)转换为bmp图像加载至空间光调制器,达到矫正共焦系统像差的效果。
本实施方式只是对本专利的示例性说明,并不限定它的保护范围,本领域技术人员还可以对其局部进行改变,只要没有超出本专利的精神实质,都在本专利的保护范围内。
Claims (5)
1.共焦显微镜模式像差矫正方法,其所依托的共焦系统,其组成包括激光器(1)、双面45°反光镜(2)、分束镜(3)、半波片(4)、空间光调制器(5)、反射镜(6)、光阑(7)、偏振分束镜(8)、XY扫描振镜(9)、扫描透镜(10)、管镜(11)、物镜(12)、1/4波片(13)、收集透镜(15)、光纤小孔(16)和光电倍增管(17),激光器(1)发出的激光经双面45°反光镜(2)的一侧反光镜面反射后途经分束镜(3)到达半波片(4),经半波片(4)调制后射入至空间光调制器(5),再经空间光调制器(5)调制后射回至分束镜(3),再经分束镜(3)反射至反射镜(6),然后经光阑(7)滤除高阶衍射杂光后射到双面45°反光镜(2)的另一侧反光镜面,经双面45°反光镜(2)再次反射后经过偏振分束镜(8)射至XY扫描振镜(9)上,然后经XY扫描振镜(9)反射至扫描透镜(10),经过管镜(11)后经过1/4波片(12)偏振面翻转90°后被物镜(13)聚焦至的打在待测样品(14)上;聚焦后的光斑再从待测样品(14)上反射回物镜(13),经过1/4波片(12),偏振面再次翻转90°后依次经过、管镜(11)、扫描透镜(10)和XY扫描振镜(9),射在偏振分束镜(8)上,再经偏振分束镜(8)反射至收集透镜(15)并聚焦到光纤小孔(16),最后被光电倍增管(17)收集;共焦显微镜模式像差矫正方法,其特征是,包括以下步骤:
步骤a,建立Zernike多项式;将像差分解为Zernike多项式4-11阶组合;Zernike多项式4-11阶多项式形式如下:
其中i代表Zernike像差阶数,作为判定循环是否继续的条件,r为极坐标系径向坐标,θ为极坐标系角向坐标;
步骤c,计算五幅图像的灰度方差函数;
步骤d,s=0:质心法计算像差系数ai,s=1:线性计算像差系数ai’;
步骤e,使i=i+1,若i<12,返回步骤b,若i>12则进行步骤f;
步骤f,若s=0,则使i=4,s=1返回步骤b;若s=1则进行步骤g;为减少各阶Zernike像差间的串扰,在第一次循环结束后,即i=12时,得到共焦系统的4-11阶像差系数,比较得到的像差系数,其中系数a<0.7rad项像差模式清零,系数a≧0.7rad项像差模式预先矫正,即在空间光调制器中加入清零后的系数a像差的bmp图像为初始值,以此基础加重新进行4-11阶偏置,进入第二次循环:i=4,s=1;
步骤g,计算最佳相差系数,进行共焦像差矫正;
所述步骤d中,s=0时,利用5幅图像的灰度方差函数值,做质心法形式如下:
M1-5和x1-5分别对应第一幅至第五幅图像的灰度方差函数值和相应偏执系数,ai为系统像差系数;对于s=1时,已知若共焦显微镜不含相差,则加入正负相同幅值偏置,其图像灰度方差函数值之差为零,且其差的值与共焦显微镜中所含像差系数成线性比例,利用像差偏置x=±b所对应的两幅图像的灰度方差函数值之差计算共焦显微镜像差系数ai’。
5.根据权利要求4所述的共焦显微镜模式像差矫正方法,其特征在于:所述步骤g中,将重新计算得到最佳像差系数a”,其中a”=a’+a,
将Φ(r,θ)=∑-ai”Zernikei(r,θ)转换为bmp图像加载至空间光调制器,达到矫正共焦系统像差的效果。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711242335.4A CN108088653B (zh) | 2017-11-30 | 2017-11-30 | 共焦显微镜模式像差矫正方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711242335.4A CN108088653B (zh) | 2017-11-30 | 2017-11-30 | 共焦显微镜模式像差矫正方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108088653A CN108088653A (zh) | 2018-05-29 |
CN108088653B true CN108088653B (zh) | 2020-11-20 |
Family
ID=62173529
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711242335.4A Active CN108088653B (zh) | 2017-11-30 | 2017-11-30 | 共焦显微镜模式像差矫正方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108088653B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109540946A (zh) * | 2018-10-25 | 2019-03-29 | 华东师范大学 | 基于人工智能算法的透射电子显微镜样品漂移矫正方法 |
CN109445082B (zh) * | 2018-11-30 | 2020-05-01 | 南京智博医疗器械有限公司 | 一种共焦扫面显微镜光源调制方法 |
CN109813533B (zh) * | 2019-02-25 | 2020-08-11 | 深圳市安思疆科技有限公司 | 一种批量测试doe衍射效率和均匀性的方法及装置 |
CN113741123B (zh) * | 2020-05-28 | 2022-12-30 | 成都理想境界科技有限公司 | 一种投影矫正方法及光纤扫描成像设备 |
CN117589081A (zh) * | 2023-08-29 | 2024-02-23 | 哈尔滨工业大学 | 一种原位共焦显微三维形貌畸变矫正方法 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103226238B (zh) * | 2013-05-14 | 2014-12-17 | 哈尔滨工业大学 | 基于反射式液晶空间光调制器的共焦横向扫描方法 |
CN103424861B (zh) * | 2013-08-19 | 2015-09-16 | 北京信息科技大学 | 基于柱偏振涡旋光束的超分辨共焦显微成像方法及装置 |
CN106767400B (zh) * | 2016-11-23 | 2019-05-10 | 哈尔滨工业大学 | 基于空间光调制器的结构探测共焦显微成像方法及装置 |
-
2017
- 2017-11-30 CN CN201711242335.4A patent/CN108088653B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN108088653A (zh) | 2018-05-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108088653B (zh) | 共焦显微镜模式像差矫正方法 | |
JP4377404B2 (ja) | 画像向上機能を備えたカメラ | |
CN107843982B (zh) | 基于实时相位差异技术的无波前探测自适应光学系统 | |
Bowman et al. | An SLM-based Shack–Hartmann wavefront sensor for aberration correction in optical tweezers | |
US8830351B2 (en) | Image processing method and image processing apparatus for image restoration to reduce a detected color shift | |
US10247933B2 (en) | Image capturing device and method for image capturing | |
CN108646407B (zh) | 一种无波前感知自适应光学成像系统及其成像方法 | |
KR20090004428A (ko) | 광학 설계 방법 및 시스템과 광학 수차를 갖는 광학 요소를이용한 촬상 소자 | |
CN1768346A (zh) | 用于最小化成像系统中的光程差效应的系统和方法 | |
JP2009188676A (ja) | 撮像装置、撮像方法、およびプログラム | |
CN107077722A (zh) | 图像捕获装置及图像捕获方法 | |
CN107991235B (zh) | 共焦显微镜模式像差矫正装置 | |
JP5587264B2 (ja) | 調整方法、調整装置、及び、撮像装置の製造方法 | |
WO2020258434A1 (zh) | 一种基于tie的相位成像方法、装置及可读存储介质 | |
CN105203213A (zh) | 一种计算复合波前传感自适应光学系统复原电压的方法 | |
CN108037074A (zh) | 共焦显微镜模式像差矫正方法 | |
Zepp et al. | Optimization of the holographic wavefront sensor for open-loop adaptive optics under realistic turbulence. Part I: simulations | |
CN207882564U (zh) | 基于实时相位差异技术的无波前探测自适应光学系统 | |
CN108037075A (zh) | 共焦显微镜模式像差矫正方法 | |
JP5224976B2 (ja) | 画像補正装置、および、画像補正方法、および、プログラム、および、記録媒体 | |
CN108089318A (zh) | 共焦显微镜模式像差矫正方法 | |
JP3584285B2 (ja) | 歪画像補正方法及び装置 | |
CN114460735A (zh) | 基于空间光调制器波前校正的显微成像装置与方法 | |
CN115131201A (zh) | 基于叠层重建的远场漫反射合成孔径超分辨成像方法 | |
CN113192146B (zh) | 一种基于深度学习消除栅线投影相移技术中离焦误差的方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |