CN106895917B - 用于运动伪影校正的空间光谱编码并行oct系统和方法 - Google Patents
用于运动伪影校正的空间光谱编码并行oct系统和方法 Download PDFInfo
- Publication number
- CN106895917B CN106895917B CN201710077740.9A CN201710077740A CN106895917B CN 106895917 B CN106895917 B CN 106895917B CN 201710077740 A CN201710077740 A CN 201710077740A CN 106895917 B CN106895917 B CN 106895917B
- Authority
- CN
- China
- Prior art keywords
- spatial spectral
- light
- spectrum
- light beam
- parallel oct
- 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
- 230000003595 spectral effect Effects 0.000 title claims abstract description 108
- 230000033001 locomotion Effects 0.000 title claims abstract description 57
- 238000000034 method Methods 0.000 title claims abstract description 53
- 238000012937 correction Methods 0.000 title claims abstract description 35
- 238000001228 spectrum Methods 0.000 claims abstract description 57
- 230000008569 process Effects 0.000 claims abstract description 28
- 238000005305 interferometry Methods 0.000 claims abstract description 4
- 239000000523 sample Substances 0.000 claims description 50
- 230000028161 membrane depolarization Effects 0.000 claims description 31
- 239000006185 dispersion Substances 0.000 claims description 24
- 238000013519 translation Methods 0.000 claims description 16
- 238000000926 separation method Methods 0.000 claims description 14
- 230000007935 neutral effect Effects 0.000 claims description 11
- 238000012545 processing Methods 0.000 claims description 9
- 230000005540 biological transmission Effects 0.000 claims description 8
- 238000004364 calculation method Methods 0.000 claims description 8
- 238000001514 detection method Methods 0.000 claims description 8
- 238000005286 illumination Methods 0.000 claims description 8
- 230000009466 transformation Effects 0.000 claims description 8
- 150000001875 compounds Chemical class 0.000 claims description 7
- 230000003287 optical effect Effects 0.000 claims description 7
- 238000009738 saturating Methods 0.000 claims description 4
- 230000015572 biosynthetic process Effects 0.000 claims description 3
- 238000013459 approach Methods 0.000 claims description 2
- 230000009885 systemic effect Effects 0.000 abstract description 2
- 238000005516 engineering process Methods 0.000 description 6
- 238000005070 sampling Methods 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- 230000007246 mechanism Effects 0.000 description 2
- 230000000149 penetrating effect Effects 0.000 description 2
- 230000010287 polarization Effects 0.000 description 2
- 230000035945 sensitivity Effects 0.000 description 2
- 230000002411 adverse Effects 0.000 description 1
- 230000004075 alteration Effects 0.000 description 1
- 239000000571 coke Substances 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 230000004069 differentiation Effects 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 238000010191 image analysis Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 230000000241 respiratory effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
- G01J3/00—Spectrometry; Spectrophotometry; Monochromators; Measuring colours
- G01J3/28—Investigating the spectrum
- G01J3/2823—Imaging spectrometer
Landscapes
- Physics & Mathematics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- General Physics & Mathematics (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
本发明公开了一种用于运动伪影校正的空间光谱编码并行OCT系统和方法。本发明在宽带光源和并行OCT系统之间加入了空间光谱编码模块，该模块沿光束扫描方向对宽带光源出射的宽带光谱进行空间光谱编码；本发明利用相邻两步扫描过程中，两个矩形照明区域的重叠区域所对应的数据进行基于互相关算法的运动伪影校正，来计算和补偿相邻次扫描过程中样品的随机运动量。利用相邻多步扫描过程所获得的多幅干涉光谱中对同一位置的不同编码光谱，拼接出该位置的完整光谱，完成空间光谱解码，从而恢复并行OCT系统的理论横向分辨率和轴向分辨率。本发明能够在保证原有并行OCT系统分辨率的同时，提供高精度、高准确度的样品随机运动量校正。
Description
技术领域
本发明属于光学相干测量领域，具体涉及一种用于运动伪影校正的空间光谱编码并行OCT系统和方法。
技术背景
OCT探测光束扫描过程中，样品存在的随机运动(如活体生物的呼吸运动，电动平移台的随机抖动等)将不可避免的造成OCT图像出现运动伪影，进而对后续基于图像的定量分析产生不利的影响。因而，如何抑制该随机运动的影响一直是OCT技术中的关键问题。
目前，OCT技术中主要有两类方法来抑制样品随机运动量：1.提高OCT的成像速度；2.检定并补偿该随机运动量。提高OCT成像速度最有效的方法是采用全场OCT成像技术，该技术采用扫频宽带光源和面阵CCD，从而避免了光束扫描机构的使用，因而能够实现MHz的OCT成像速度，极大的抑制了成像周期内样品的随机运动量对OCT图像质量的影响。然而，该全场OCT成像技术存在着探测单元之间信号互相串扰的固有问题，从而导致图像分辨率和信号灵敏度的降低。因此，为了兼顾OCT成像速度和信号灵敏度，采用线性扫描的并行OCT技术将更适合应用于抑制样品的随机运动量。并行OCT由于只是同时获得单次B-Scan的信号，因此对相邻B-Scan仍然需要检定并补偿样品随机运动量。目前检定并补偿样品随机运动量的方法主要分为两种：1.OCT探测光束扫描过程中实时检定和补偿，2.后续算法处理过程中算法检定和补偿。方法1中需要引入额外的激光光束实时追踪样品的随机运动量，并通过复杂的回馈装置，通过改进扫描机构的扫描模式来补偿该随机运动量。该方法增加了系统的复杂性和系统成本。方法2通过对相邻过采样数据之间的互相关计算，来求取OCT探测光束扫描过程中样品的随机运动量，并进行补偿。然而该方法需要更高密度的数据采样点，并且需要假定样品表面相对平滑。
发明内容
本发明针对现有技术的不足，提出一种用于运动伪影校正的空间光谱编码并行OCT系统和方法。
一种用于运动伪影校正的空间光谱编码并行OCT系统，该系统包括宽带光源、空间光谱编码模块以及并行OCT模块。从宽带光源出射的宽带光经过空间光谱编码模块后，宽带光沿着整个系统光束扫描方向(x轴方向)形成空间光谱编码光，空间光谱编码光随后进入并行OCT模块，聚焦照射在待测样品上形成空间光谱编码的矩形照明区域，经由待测样品反射的信号光重新进入并行OCT模块，与并行OCT模块中的参考光发生干涉，干涉形成的干涉光由并行OCT模块中的探测单元进行光谱分解和并行探测，所测数据最终传入并行OCT模块的计算机中进行后续的算法处理。
所述空间光谱编码模块包括：色散棱镜、第一消色差透镜、刀口棱镜、第一平面反射镜、第二平面反射镜、第一消偏振分光棱镜、第二消色差透镜和第三平面反射镜。色散棱镜对进入空间光谱编码模块的宽带光沿着整个系统光束扫描方向(x轴方向)进行空间线性色散，色散后的宽带光通过第一消色差透镜和刀口棱镜后发生左半边和右半边光谱的对称分离，左右分离后的宽带光分别经由第一平面反射镜、第二平面反射镜反射，再入射到第一消偏振分光棱镜后重新汇集，再经由第二消色差透镜准直，此时左右半边光谱在空间上完全重合，色散方向相反，从而形成空间光谱编码光。最后，空间光谱编码光经由第三平面反射镜后出射空间光谱编码模块。
所述并行OCT模块包括：第一柱面透镜、第二消偏振分光棱镜、中性密度衰减片、第二柱面透镜、第三消色差透镜、第四平面反射镜、第四消色差透镜、电动平移台、组合透镜、狭缝、第五消色差透镜、透射光栅、第六消色差透镜、面阵CMOS、高速图像采集卡、电动平移台控制器以及计算机。从空间光谱编码模块出射的空间光谱编码光进入并行OCT模块，首先经由第一柱面透镜后变为柱面波光束出射，该柱面波光束随后由第二消偏振分光棱镜分成第一光束和第二光束。第一光束经由中性密度衰减片、第二柱面透镜以及第三消色差透镜后聚焦于第四平面反射镜上，该光束随后再逆向经由第三消色差透镜、第二柱面透镜和中性密度衰减片回到第二消偏振分光棱镜。第二光束经由第四消色差透镜聚焦于放置于电动平移台的待测样品上，从待测样品返回的光束再逆向经由第四消色差透镜回到第二消偏振分光棱镜。返回的第一光束和第二光束在第二消偏振分光棱镜发生干涉，随后经由组合透镜、狭缝、第五消色差透镜、透射光栅、第六消色差透镜后由面阵CMOS进行光谱分解和并行探测。探测所得数据经由高速图像采集卡传入计算机中进行后续的算法处理。
一种用于运动伪影校正的空间光谱编码并行OCT方法，该方法包括以下步骤：
步骤Ⅰ、利用空间光谱编码模块对宽带光源出射的宽带光进行空间光谱编码，使得最终经过并行OCT模块照射到待测样品上时，形成空间光谱编码的矩形照明区域；
步骤Ⅱ、利用相邻两步扫描过程中，两个矩形照明区域的重叠区域所对应的数据进行基于互相关算法的运动伪影校正，来计算和补偿相邻次扫描过程中样品的随机运动量；
步骤Ⅲ、待校正过运动伪影后，利用相邻多步扫描过程所获得的数据完成空间光谱解码，从而恢复并行OCT系统的理论横向分辨率和轴向分辨率。
所述步骤Ⅰ包括：在空间光谱编码模块中，利用色散棱镜对入射的宽带光沿着整个系统光束扫描方向(x轴方向)进行空间线性色散，然后利用第一消色差透镜和刀口棱镜对色散后的宽带光做左半边和右半边光谱的对称分离，再利用第一平面反射镜、第二平面反射镜和第一消偏振分光棱镜使分离的光谱重新汇集，最后利用第二消色差透镜将汇集后的光准直出射，此时左右半边光谱在空间上完全重合，色散方向相反，从而形成空间光谱编码光。空间光谱编码光经过并行OCT模块照射到待测样品上，形成空间光谱编码的矩形照明区域，在x轴方向上，不同的位置对应于不同的光谱段，即x坐标被编码于光谱坐标中。
所述步骤Ⅱ包括：设定适宜的扫描步长使相邻两步光束扫描的矩形照明区域有一定比例的重叠，对相邻两步的重叠区域所对应的数据做基于互相关算法的运动伪影校正，以前一步的数据为基准，计算得到后一步的数据所需要补偿的错移量。所述的基于互相关算法的运动伪影校正由像素级错移校正和相位级错移校正两个步骤所组成。首先对相邻两步重叠区域所对应的数据分别做快速傅里叶变换，对变换后的数据做互相关算法处理，通过结果的模值计算得到样品在后一步光束扫描中的像素级错移量，并对该步数据进行该像素级错移量补偿。然后把前一步数据和校正过像素级错移量的后一步数据分别分成两个部分，分别做快速傅里叶变换，再对变换后的数据分别对应地做互相关算法处理，对处理后的两个部分的结果求和，通过求和结果的相位值计算得到样品在后一步光束扫描中的相位级错移量，并对该步数据进行该相位级错移量补偿。
所述步骤Ⅲ包括：待校正过运动伪影后，利用相邻多步扫描过程所获得的多幅干涉光谱中对应同一位置的不同编码光谱，拼接出该位置的完整光谱，完成空间光谱解码，从而恢复并行OCT系统的理论横向分辨率和轴向分辨率。
与背景技术相比，本发明具有的有益效果是：
1.相比于全场OCT成像技术，本发明采用了空间光谱编码模块和并行OCT模块相结合的系统架构，因而同时兼顾了较高的OCT成像速度和较高的信号信噪比。
2.相比于OCT探测光束扫描过程中实时检定和补偿的方法，本发明由于不需要额外的激光光束以及复杂的回馈装置用于校正样品的随机运动量，因此系统将显得更为简洁，从而有效的降低了系统成本。
3.相比于技术背景中所使用相邻过采样数据做互相关计算的方法，本发明采用了空间光谱编码的方法，相邻2步扫描的重叠区域编码于不同的光谱段，因此避免了过采样，有效地减少了数据采集密度。
4.相比于技术背景中所使用相邻过采样数据做互相关计算的方法，本发明采用了空间光谱编码的方法，可以保证用于互相关计算的数据是严格对应于同一位置的，因而能够适用于非平滑的样品。
附图说明
图1是本发明的空间光谱编码并行OCT系统示意图；
图2是本发明的空间光谱编码和扫描示意图；
图3是本发明的运动伪影矫正和空间光谱解码示意图。
图1中：1、宽带光源，2、色散棱镜，3、第一消色差透镜，4、刀口棱镜，5、第一平面反射镜，6、第二平面反射镜，7、第一消偏振分光棱镜，8、第二消色差透镜，9、第三平面反射镜，10、第一柱面透镜，11、第二消偏振分光棱镜，12、中性密度衰减片，13、第二柱面透镜，14、第三消色差透镜，15、第四平面反射镜，16、第四消色差透镜，17、待测样品，18、电动平移台，19、组合透镜，20、狭缝，21、第五消色差透镜，22、透射光栅，23、第六消色差透镜，24、面阵CMOS，25、高速图像采集卡，26、计算机，27、电动平移台控制器，28、空间光谱编码模块，29、并行OCT模块。
具体实施方式
本发明的目的在于提供一种用于运动伪影校正的空间光谱编码并行OCT系统和方法。本发明相较于传统的并行OCT系统，在宽带光源和并行OCT系统之间加入了空间光谱编码模块，该模块沿光束扫描方向对宽带光源出射的宽带光谱进行空间光谱编码，从而实现空间光谱编码的矩形区域照明。然后，利用相邻两步扫描过程中，两个矩形照明区域的重叠区域所对应的数据进行基于互相关算法的运动伪影校正，来计算和补偿相邻次扫描过程中样品的随机运动量。利用相邻多步扫描过程所获得的多幅干涉光谱中对同一位置的不同编码光谱，拼接出该位置的完整光谱，完成空间光谱解码，从而恢复并行OCT系统的理论横向分辨率和轴向分辨率。因而，本发明能够在保证原有并行OCT系统分辨率的同时，提供高精度、高准确度的样品随机运动量校正。
本发明的目的是通过如下技术方案实现的：
一种用于运动伪影校正的空间光谱编码并行OCT系统：
该系统包括宽带光源、空间光谱编码模块以及并行OCT模块。从宽带光源出射的宽带光经过空间光谱编码模块后，宽带光沿着整个系统光束扫描方向(x轴方向)形成空间光谱编码光，空间光谱编码光随后进入并行OCT模块，聚焦照射在待测样品上形成空间光谱编码的矩形照明区域，经由待测样品反射的信号光重新进入并行OCT模块，与并行OCT模块中的参考光发生干涉，干涉形成的干涉光由并行OCT模块中的探测单元进行光谱分解和并行探测，所测数据最终传入并行OCT模块的计算机中进行后续的算法处理。
所述空间光谱编码模块包括：色散棱镜、第一消色差透镜、刀口棱镜、第一平面反射镜、第二平面反射镜、第一消偏振分光棱镜、第二消色差透镜和第三平面反射镜。色散棱镜对进入空间光谱编码模块的宽带光沿着整个系统光束扫描方向(x轴方向)进行空间线性色散，色散后的宽带光通过第一消色差透镜和刀口棱镜后发生左半边和右半边光谱的对称分离，左右分离后的宽带光分别经由第一平面反射镜、第二平面反射镜和第一消偏振分光棱镜后重新汇集，再经由第二消色差透镜准直，此时左右半边光谱在空间上完全重合，色散方向相反，从而形成空间光谱编码光。最后，空间光谱编码光经由第三平面反射镜后出射空间光谱编码模块。
所述并行OCT模块包括：第一柱面透镜、第二消偏振分光棱镜、中性密度衰减片、第二柱面透镜、第三消色差透镜、第四平面反射镜、第四消色差透镜、电动平移台、组合透镜、狭缝、第五消色差透镜、透射光栅、第六消色差透镜、面阵CMOS、高速图像采集卡、电动平移台控制器以及计算机。从空间光谱编码模块出射的空间光谱编码光进入并行OCT模块，首先经由第一柱面透镜后变为柱面波光束出射，该柱面波光束随后由第二消偏振分光棱镜分成第一光束和第二光束。第一光束经由中性密度衰减片、第二柱面透镜以及第三消色差透镜后聚焦于第四平面反射镜上，该光束随后再逆向经由第三消色差透镜、第二柱面透镜和中性密度衰减片回到第二消偏振分光棱镜。第二光束经由第四消色差透镜聚焦于放置于电动平移台的待测样品上，从待测样品返回的光束再逆向经由第四消色差透镜回到第二消偏振分光棱镜。返回的第一光束和第二光束在第二消偏振分光棱镜发生干涉，随后经由组合透镜、狭缝、第五消色差透镜、透射光栅、第六消色差透镜后由面阵CMOS进行光谱分解和并行探测。探测所得数据经由高速图像采集卡传入计算机中进行后续的算法处理。
一种用于运动伪影校正的空间光谱编码并行OCT方法：
1)在空间光谱编码模块中，利用色散棱镜对入射的宽带光沿着整个系统光束扫描方向(x轴方向)进行空间线性色散，然后利用第一消色差透镜和刀口棱镜对色散后的宽带光做左半边和右半边光谱的对称分离，再利用第一平面反射镜、第二平面反射镜和第一消偏振分光棱镜使分离的光谱重新汇集，最后利用第二消色差透镜将汇集后的光准直出射，此时左右半边光谱在空间上完全重合，色散方向相反，从而形成空间光谱编码光。空间光谱编码光经过并行OCT模块照射到待测样品上，形成空间光谱编码的矩形照明区域，在x轴方向上，不同的位置对应于不同的光谱段，即x坐标被编码于光谱坐标中。
2)设定适宜的扫描步长使相邻两步光束扫描的矩形照明区域有一定比例的重叠，对相邻两步的重叠区域所对应的数据做基于互相关算法的运动伪影校正，以前1步的数据为基准，计算得到后一步的数据所需要补偿的错移量。所述的基于互相关算法的运动伪影校正由像素级错移校正和相位级错移校正两个步骤所组成。首先对相邻两步重叠区域所对应的数据分别做快速傅里叶变换，对变换后的数据做互相关算法处理，通过结果的模值计算得到样品在后一步光束扫描中的像素级错移量，并对该步数据进行该像素级错移量补偿。然后把前一步数据和校正过像素级错移量的后一步数据分别分成两个部分，分别做快速傅里叶变换，对变换后的数据分别对应地做互相关算法处理，对处理后的两个部分的结果求和，通过求和结果的相位值计算得到样品在后一步光束扫描中的相位级错移量，并对该步数据进行该相位级错移量补偿。
3)待校正过运动伪影后，利用相邻多步扫描过程所获得的多幅干涉光谱中对应同一位置的不同编码光谱，拼接出该位置的完整光谱，完成空间光谱解码，从而恢复并行OCT系统的理论横向分辨率和轴向分辨率。
下面结合附图和实施示例对本发明作进一步的说明：
如图1所示，一种用于运动伪影校正的空间光谱编码并行OCT系统，包括宽带光源1、色散棱镜2、第一消色差透镜3、刀口棱镜4、第一平面反射镜5、第二平面反射镜6、第一消偏振分光棱镜7、第二消色差透镜8、第三平面反射镜9、第一柱面透镜10、第二消偏振分光棱镜11、中性密度衰减片12、第二柱面透镜13、第三消色差透镜14、第四平面反射镜15、第四消色差透镜16、待测样品17、电动平移台18、组合透镜19、狭缝20、第五消色差透镜21、透射光栅22、第六消色差透镜23、面阵CMOS24、高速图像采集卡25、计算机26、电动平移台控制器27、空间光谱编码模块28、并行OCT模块29。
从宽带光源1发出的宽带光，进入空间光谱编码模块28。在空间光谱编码模块28中，宽带光源1发出的宽带光经过色散棱镜2后沿着整个系统光束扫描方向(x轴方向)发生空间线性色散，色散后的宽带光通过第一消色差透镜3和刀口棱镜4后发生左半边(k-m～k0)和右半边(k0～km)光谱的对称分离，如图1中空心箭头和实心箭头所示。左右分离后的宽带光分别经由第一平面反射镜5、第二平面反射镜6和第一消偏振分光棱镜7后重新汇集，再经由第二消色差透镜8准直，此时左右半边光谱在空间上完全重合，色散方向相反，从而形成空间光谱编码光。最后，空间光谱编码光经由第三平面反射镜9后出射空间光谱编码模块28。
从空间光谱编码模块28出射的宽带光进入并行OCT模块29，首先经由第一柱面透镜10后变为柱面波光束出射，该柱面波光束随后由第二消偏振分光棱镜11分成第一光束和第二光束。第一光束经由中性密度衰减片12、第二柱面透镜13以及第三消色差透镜14后聚焦于第四平面反射镜15上，该光束随后再逆向经由第三消色差透镜14、第二柱面透镜13和中性密度衰减片12回到第二消偏振分光棱镜11。第二光束经由第四消色差透镜16聚焦于放置于电动平移台18的待测样品17上，从待测样品17返回的光束再逆向经由第四消色差透镜16回到第二消偏振分光棱镜11。返回的第一光束和第二光束在第二消偏振分光棱镜11发生干涉，随后经由组合透镜19、狭缝20、第五消色差透镜21、透射光栅22、第六消色差透镜后23由面阵CMOS24进行光谱分解和并行探测。探测所得数据经由高速图像采集卡25传入计算机26中进行后续的算法处理。
如图2所示，为本发明的空间光谱编码和扫描示意图。图中所示y轴对应于传统并行OCT系统的并行探测方向。由于空间光谱编码模块对宽带光沿着整个系统光束扫描方向(x方向)发生空间光谱编码，产生左右半边光谱在空间上完全重合但色散方向相反的空间光谱编码光，因此沿x轴方向，待测样品不同位置将由不同的光谱段照射，但照射任一位置的两个光谱段的平均波数又是相同的，即实现了对待测样品的空间光谱编码。干涉光在探测单元中通过光谱分解形成的干涉光谱由图1中所示面阵CMOS24探测。通过以上所述空间光谱编码的方式，在一个干涉光谱中同时记录了照明区域内x位置的信息，并以光谱段的不同作为区分。为了简化说明，以扫描步长为照明区域宽度的三分之一为例，图中只列出了k-3、k-2、k-1、k0、k1、k2、k3，7个光谱坐标。如图2中第n步扫描时，照明的区域在x轴上的坐标分别为xn-1、xn、xn+1，对应的光谱段分别为(k-3～k-2、k2～k3)、(k-2～k-1、k1～k2)以及(k-1～k0、k0～k1)，平均波数均为k0。相同的x和k的对应关系也体现在面阵CMOS24所探测到干涉光谱中。图1中计算机26发出控制信号，经由电动平移台控制器27后，传送给电动平移台18，从而实现图2中x轴扫描步进n到x轴扫描步进n+1的移动。第n步和第n+1步扫描所照明的区域在x轴上的坐标分别为xn-1、xn、xn+1和xn、xn+1、xn+2，重叠区域为xn、xn+1。
如图3所示，为运动伪影矫正和空间光谱解码示意图。其包括步骤Ⅱ和步骤Ⅲ，在步骤Ⅱ中，利用窗口函数H’和H”分别选取图2中两步扫描所获得的数据In和In+1中，对应于重叠区域的部分数据，并分别对其进行快速傅里叶变换处理对快速傅里叶变换后的数据做互相关计算①，通过结果的模值计算得到样品相邻次光束扫描过程中的像素级错移量对数据In+1补偿该像素级错移量②。然后把In和校正过像素级错移量的In+1分成两个部分并分别做傅里叶变换，得到结果和和和分别做互相关算法处理，再对处理所得的结果进行求和运算③，从而消除与深度位置有关的相位量，进而求取相位级错移量并对其In+1进行该相位级错移量补偿④，从而得到校正过样品随机运动量的数据I'n+1。在步骤Ⅲ中，对校正过运动伪影的相邻3步扫描的数据In-1、In和In+1，分别提取出编码了位置xn的离散光谱数据(k-1～k1)、(k-2～k-1、k1～k2)以及(k-3～k-2、k2～k3)，从而拼接出相应光束扫描位置xn处的完整光谱(k-3～k3)，进而恢复并行OCT系统的理论横向分辨率和轴向分辨率。
Claims (7)
1.用于运动伪影校正的空间光谱编码并行OCT系统，其特征在于：该系统包括宽带光源、空间光谱编码模块以及并行OCT模块；从宽带光源出射的宽带光经过空间光谱编码模块后，宽带光沿着整个系统光束扫描方向形成空间光谱编码光，空间光谱编码光随后进入并行OCT模块，聚焦照射在待测样品上形成空间光谱编码的矩形照明区域，经由待测样品反射的信号光重新进入并行OCT模块，与并行OCT模块中的参考光发生干涉，干涉形成的干涉光由并行OCT模块中的探测单元进行光谱分解和并行探测，所测数据最终传入并行OCT模块的计算机中进行后续的算法处理。
2.如权利要求1所述的用于运动伪影校正的空间光谱编码并行OCT系统，其特征在于：所述空间光谱编码模块包括：色散棱镜、第一消色差透镜、刀口棱镜、第一平面反射镜、第二平面反射镜、第一消偏振分光棱镜、第二消色差透镜和第三平面反射镜；色散棱镜对进入空间光谱编码模块的宽带光沿着整个系统光束扫描方向进行空间线性色散，色散后的宽带光通过第一消色差透镜和刀口棱镜后发生左半边和右半边光谱的对称分离，左右分离后的宽带光分别经由第一平面反射镜、第二平面反射镜反射，再入射到第一消偏振分光棱镜后重新汇集，再经由第二消色差透镜准直，此时左右半边光谱在空间上完全重合，色散方向相反，从而形成空间光谱编码光；最后，空间光谱编码光经由第三平面反射镜后出射空间光谱编码模块。
3.如权利要求1所述的用于运动伪影校正的空间光谱编码并行OCT系统，其特征在于：所述并行OCT模块包括：第一柱面透镜、第二消偏振分光棱镜、中性密度衰减片、第二柱面透镜、第三消色差透镜、第四平面反射镜、第四消色差透镜、电动平移台、组合透镜、狭缝、第五消色差透镜、透射光栅、第六消色差透镜、面阵CMOS、高速图像采集卡、电动平移台控制器以及计算机；从空间光谱编码模块出射的空间光谱编码光进入并行OCT模块，首先经由第一柱面透镜后变为柱面波光束出射，该柱面波光束随后由第二消偏振分光棱镜分成第一光束和第二光束；第一光束经由中性密度衰减片、第二柱面透镜以及第三消色差透镜后聚焦于第四平面反射镜上，该光束随后再逆向经由第三消色差透镜、第二柱面透镜和中性密度衰减片回到第二消偏振分光棱镜；第二光束经由第四消色差透镜聚焦于放置于电动平移台的待测样品上，从待测样品返回的光束再逆向经由第四消色差透镜回到第二消偏振分光棱镜；返回的第一光束和第二光束在第二消偏振分光棱镜发生干涉，随后经由组合透镜、狭缝、第五消色差透镜、透射光栅、第六消色差透镜后由面阵CMOS进行光谱分解和并行探测；探测所得数据经由高速图像采集卡传入计算机中进行后续的算法处理。
4.根据权利要求1所述的用于运动伪影校正的空间光谱编码并行OCT系统的方法，其特征在于：该方法包括以下步骤：
步骤Ⅰ、利用空间光谱编码模块对宽带光源出射的宽带光进行空间光谱编码，使得最终经过并行OCT模块照射到待测样品上时，形成空间光谱编码的矩形照明区域；
步骤Ⅱ、利用相邻两步扫描过程中，两个矩形照明区域的重叠区域所对应的数据进行基于互相关算法的运动伪影校正，来计算和补偿相邻次扫描过程中样品的随机运动量；
步骤Ⅲ、待校正过运动伪影后，利用相邻多步扫描过程所获得的数据完成空间光谱解码，从而恢复并行OCT系统的理论横向分辨率和轴向分辨率。
5.根据权利要求4所述的用于运动伪影校正的空间光谱编码并行OCT系统的方法，其特征在于：所述步骤Ⅰ包括：在空间光谱编码模块中，利用色散棱镜对入射的宽带光沿着整个系统光束扫描方向进行空间线性色散，然后利用第一消色差透镜和刀口棱镜对色散后的宽带光做左半边和右半边光谱的对称分离，再利用第一平面反射镜、第二平面反射镜和第一消偏振分光棱镜使分离的光谱重新汇集，最后利用第二消色差透镜将汇集后的光准直出射，此时左右半边光谱在空间上完全重合，色散方向相反，从而形成空间光谱编码光；空间光谱编码光经过并行OCT模块照射到待测样品上，形成空间光谱编码的矩形照明区域，在x轴方向上，不同的位置对应于不同的光谱段，即x坐标被编码于光谱坐标中。
6.根据权利要求4所述的用于运动伪影校正的空间光谱编码并行OCT系统的方法，其特征在于：所述步骤Ⅱ包括：设定适宜的扫描步长使相邻两步光束扫描的矩形照明区域有一定比例的重叠，对相邻两步的重叠区域所对应的数据做基于互相关算法的运动伪影校正，以前一步的数据为基准，计算得到后一步的数据所需要补偿的错移量；所述的基于互相关算法的运动伪影校正由像素级错移校正和相位级错移校正两个步骤所组成；首先对相邻两步重叠区域所对应的数据分别做快速傅里叶变换，对变换后的数据做互相关算法处理，通过结果的模值计算得到样品在后一步光束扫描中的像素级错移量，并对该步数据进行该像素级错移量补偿；然后把前一步数据和校正过像素级错移量的后一步数据分别分成两个部分，分别做快速傅里叶变换，对变换后的数据分别对应地做互相关算法处理，对处理后的两个部分的结果求和，通过求和结果的相位值计算得到样品在后一步光束扫描中的相位级错移量，并对该步数据进行该相位级错移量补偿。
7.根据权利要求4所述的用于运动伪影校正的空间光谱编码并行OCT系统的方法，其特征在于：所述步骤Ⅲ包括：待校正过运动伪影后，利用相邻多步扫描过程所获得的多幅干涉光谱中对应同一位置的不同编码光谱，拼接出该位置的完整光谱，完成空间光谱解码，从而恢复并行OCT系统的理论横向分辨率和轴向分辨率。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710077740.9A CN106895917B (zh) | 2017-02-14 | 2017-02-14 | 用于运动伪影校正的空间光谱编码并行oct系统和方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710077740.9A CN106895917B (zh) | 2017-02-14 | 2017-02-14 | 用于运动伪影校正的空间光谱编码并行oct系统和方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106895917A CN106895917A (zh) | 2017-06-27 |
| CN106895917B true CN106895917B (zh) | 2018-09-28 |
Family
ID=59198494
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710077740.9A Active CN106895917B (zh) | 2017-02-14 | 2017-02-14 | 用于运动伪影校正的空间光谱编码并行oct系统和方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106895917B (zh) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112326585B (zh) * | 2020-10-28 | 2021-10-22 | 南京农业大学 | 一种用于快速检测草莓白粉病的红外光谱装置及检测方法 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6122241A (en) * | 1995-01-31 | 2000-09-19 | Ricoh Company, Ltd. | Optical pick-up device for processing central and peripheral optical signal components |
| CN101797146A (zh) * | 2010-01-29 | 2010-08-11 | 浙江大学 | 无扫描光谱编码的内窥成像的方法和系统 |
| CN102628799A (zh) * | 2012-04-25 | 2012-08-08 | 中国科学院上海光学精密机械研究所 | 无需深度扫描的时域光学相干层析成像方法与系统 |
| CN206496841U (zh) * | 2017-02-14 | 2017-09-15 | 浙江大学 | 用于运动伪影校正的空间光谱编码并行oct系统 |
-
2017
- 2017-02-14 CN CN201710077740.9A patent/CN106895917B/zh active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6122241A (en) * | 1995-01-31 | 2000-09-19 | Ricoh Company, Ltd. | Optical pick-up device for processing central and peripheral optical signal components |
| CN101797146A (zh) * | 2010-01-29 | 2010-08-11 | 浙江大学 | 无扫描光谱编码的内窥成像的方法和系统 |
| CN102628799A (zh) * | 2012-04-25 | 2012-08-08 | 中国科学院上海光学精密机械研究所 | 无需深度扫描的时域光学相干层析成像方法与系统 |
| CN206496841U (zh) * | 2017-02-14 | 2017-09-15 | 浙江大学 | 用于运动伪影校正的空间光谱编码并行oct系统 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106895917A (zh) | 2017-06-27 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| DE102015207431B4 (de) | Erfassungsvorrichtung und -verfahren einer Oberflächenform verwendend eine Beugungswellenfront einer Aperturblendenstitchmessung | |
| DE69722876T2 (de) | Vorrichtung und Verfahren zur optischen Profilmessung | |
| KR101336048B1 (ko) | 광 단층촬영의 촬상방법 및 그 장치 | |
| US9194811B1 (en) | Apparatus and methods for improving defect detection sensitivity | |
| CN104964982B (zh) | 基于oct复信号的玻璃表面真伪缺陷识别方法及系统 | |
| CN102494623A (zh) | 镜头中光学表面中心间距的非接触式测量方法及测量装置 | |
| CN107543508A (zh) | 光学系统及使用该系统的物体表面三维形貌侦测方法 | |
| CN104698068B (zh) | 高空间分辨激光双轴差动共焦光谱‑质谱显微成像方法与装置 | |
| CN102439395A (zh) | 形状测定装置、观察装置及影像处理方法 | |
| CN109186759B (zh) | 一种光栅光谱仪像质测量方法和装置 | |
| CN106768280A (zh) | 一种基于多波长无透镜傅里叶变换数字全息的振动检测装置 | |
| DE19832175B4 (de) | Verfahren und Anordnungen zur transversalen optischen Kohärenztomographie | |
| CN105785566A (zh) | 一种利用空间光调制器改善光声成像有限视角的方法 | |
| CN109916331A (zh) | 一种基于复合光栅的结构光微纳结构三维检测方法 | |
| CN116183568B (zh) | 一种三维结构光照明超分辨显微成像的高保真重构的方法和装置 | |
| RU2006140854A (ru) | Способ контроля шероховатости деревянной поверхности | |
| CN111369484B (zh) | 钢轨廓形检测方法及装置 | |
| CN106895917B (zh) | 用于运动伪影校正的空间光谱编码并行oct系统和方法 | |
| CN106872469A (zh) | 一种基于共角干涉的层析相位显微方法与装置 | |
| CN206496841U (zh) | 用于运动伪影校正的空间光谱编码并行oct系统 | |
| CN105510229A (zh) | 一种超分辨虚拟结构光照明成像装置及其成像方法 | |
| CN109932162B (zh) | 一种基于白光配准的腔模参数检测装置及检测方法 | |
| CN101477325A (zh) | 多彩离轴数字全像系统及其成像方法 | |
| EP1931970B1 (de) | Vorrichtung und verfahren zur tiefenaufgelösten optischen erfassung einer probe | |
| CN105259668A (zh) | 基于黑色支撑的叠层成像技术 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |