微结构光纤高分辨率三维折射率测试方法与流程

(一)技术领域

本发明涉及的是一种微结构光纤高分辨率三维折射率测试方法，可用于微结构光纤的高分辨率三维显微成像，可广泛应用于光学透明物体的无损、无标记、非接触式的三维层析成像等，属于显微成像技术领域。

(二)

背景技术：

微结构光纤的三维折射率分布，是其一种重要的固有属性，对于光学透明的微结构光纤来说，三维折射率分布，可以反映样品的微结构、密度等信息，因此要实现无损、无标记、非接触式的三维层析成像。

在当代的生命科学研究中，通常采用荧光成像对待测样品进行标记。但是在进行标记的过程中，会对待测样品本身产生一定的影响，从而会影响最终的研究结果。而数字全息层析成像技术是一种新型的无损、无标记、非接触式的新型成像技术，可以重建得到微结构光纤的三维折射率分布信息，是近年来的一个研究热点。

数字全息层析成像技术结合了数字全息显微成像技术和计算机断层扫描技术，是近年来提出的新技术。近年来，虽然提出了多种应用数字全息层析成像技术的成像方法，但是大部分的思路都是结合马赫曾德尔干涉光路进行数字全息记录。

马赫曾德尔干涉光路成像所用的器件更多，系统稳定性要求更高，操作成为复杂。基于马赫曾德尔干涉光路的方法，对器件的要求更高，而且光路更为复杂且难以调试，因此迫切需要一种新的成像方法，其所用器件更少，光路更为简单，系统稳定性更高，操作更为简便，而且对微结构光纤测量的灵敏度更高，成像分辨率更高。

本发明提出了一种微结构光纤高分辨率三维折射率测试方法，本发明采用f-p腔干涉光路，其光路简单，所用的器件更少，因此系统的稳定性更高。从原理上来说，光束在f-p腔多次反射，并多次经过微结构光纤，光程差依次累积，因此测量的灵敏度更高。f-p腔干涉形成的全息图，精细度相比马赫曾德尔干涉更高，全息图的记录更为精确。

专利cn201310082100.9公开了一种数字全息成像在线重构显示系统及方法，其特点是采用了马赫曾德尔干涉数字全息记录光路，相比本发明提出的微结构光纤高分辨率三维折射率测试方法，本发明提出的新方法，测量灵敏度更高。

专利cn201610911993.7公开了一种双波长相位显微成像系统和方法、以及对应相位恢复方法，其特点是采用马赫曾德尔干涉光路来实现双波长同轴相移干涉显微系统，与本发明有着本质性的区别，相比本发明提出度新方法，也没法本发明的测量灵敏度高。

专利cn201710518263.5公开了一种三折射率层析显微成像系统及其方法，其虽然可以还原出样本的三维折射率信息，但是在成像光路上，与本发明有着本质区别，且所需器件更为复杂。

专利cn201710904860.1公开了一种光学相干断层扫描成像系统，该成像系统采用了马赫曾德尔干涉光路，其特点是采用了光纤来简化系统，降低成本，但是相比f-p腔的光路结构，仍然较为复杂。

专利cn201810145657.5公开了一种高分辨率数字全息衍射层析成像，其特点是采用马赫曾德尔干涉光路结构，利用合成孔径方法得到n幅合成高分辨率全息图，进而获得被测样品的高分辨率三维折射率再现。相对来说，结构更为复杂，与本发明专利有着本质区别。

专利cn201910136421.x公开了一种超分辨率数字全息成像系统和成像方法，其成像系统的特点是在传统的马赫曾德尔干涉光路前加入一块透射式空间光调制器，对光源进行调制。与采用f-p腔的光路结构的本发明相比，有着本质区别。

本发明公开了一种微结构光纤高分辨率三维折射率测试方法，本发明可用于微结构光纤的高分辨率显微成像，可广泛应用于光学透明物体的无损、无标记、非接触式的层析成像等领域。该微结构光纤高分辨率三维折射率测试方法采用基于f-p腔的数字全息记录光路，运用数字全息层析成像技术，对f-p腔光路记录的数字全息图，进行层析重建，还原出微结构光纤的三维折射率分布。与在先技术相比，由于采用了基于f-p腔的数字全息记录光路，光束在f-p腔中多次反射经过微结构光纤，测量的灵敏度更高。该微结构光纤高分辨率三维折射率测试方法具有结构简单、灵敏度高、分辨率更高的优点。

(三)

技术实现要素：

本发明的目的在于提供一种结构简单、灵敏度高、分辨率更高的微结构光纤高分辨率三维折射率测试方法。

本发明的目的是这样实现的：

该微结构光纤高分辨率三维折射率测试方法包括基于f-p腔的数字全息图记录、数值重建、解包裹、误差处理、三维相位分布重建和三维折射率转换。记录入射f-p腔光波的复振幅a0，及未放入微结构光纤样品时的数字全息图，之后在f-p腔中放入微结构光纤，在θ＝0时，记录含有微结构光纤样品信息的最优的数字全息图，对记录的数字全息图进行数值重建，得到透射光的复振幅分布ut，代入已知的数值，可以得到该角度时的相位分布δi，对相位分布进行解包裹处理，含有微结构光纤样品信息的相位分布δi减去未放入微结构光纤样品时的相位分布δ0，即得到只含有微结构光纤样品信息的相位分布δ。旋转微结构光纤样品，记录不同角度时的最优的数字全息图，依次得到该角度时的只含有微结构光纤样品信息的相位分布。对所有角度的相位分布做iradon变换，即可重建微结构光纤的三维相位分布δ(x,y,z)，由推导公式可以转换得到微结构光纤的高分辨率三维折射率分布n(x,y,z)。

本发明提出的微结构光纤高分辨率三维折射率测试方法，适用于包括基于f-p腔的数字全息记录光路11、微结构光纤12、控制模块13和计算显示模块14的测量系统，如图2a所示。所述基于f-p腔的数字全息记录光路11包括光源，扩束器，f-p腔，显微物镜，图像采集器等。所述微结构光纤12为具有各种微结构的光纤。所述控制模块13由计算机、仪器控制单元和仪器控制接口组成，对图像采集器、旋转控制平台等进行控制和操作，完成包含微结构光纤信息的数字全息图记录。所述计算显示模块14，对记录的数字全息图进行程序处理，并在线显示微结构光纤12的三维折射率分布信息。

所述基于f-p腔的数字全息记录光路11包括光源，扩束器，f-p腔，显微物镜，图像采集器等。优选的，光源采用波长为532nm的激光光源，并选用相应波长的扩束器，选用腔长固定的f-p标准具，放大倍数为20倍的显微物镜，像素尺寸为3.45μm×3.45μm的chargecoupleddevice(ccd)等器件。

所述微结构光纤12位于基于f-p腔的数字全息记录光路11中的f-p腔中，光束在f-p腔多次反射，并多次经过微结构光纤，光程差依次累积，因此测量的灵敏度更高。

所述微结构光纤高分辨率三维折射率测试方法中，f-p腔的腔长大于微结构光纤12的直径，使得微结构光纤12能够放入f-p腔中。将微结构光纤12去除涂覆层，并填充折射率匹配液。

所述微结构光纤高分辨率三维折射率测试方法中，f-p腔的界面反射率在0～1之间。

所述微结构光纤高分辨率三维折射率测试方法中，f-p腔具有不同长度的腔长，且可以具有不同的形状和尺寸。f-p腔的腔长可以是可调的，当f-p腔的腔长固定时，即是f-p标准具。

所述的微结构光纤12为具有各种微结构的光纤，是具有任意三维折射率分布的普通单模光纤，特种光纤等。微结构光纤12为微小物体，光束经过微结构光纤12时，会产生柱透镜效应，因此需要选择并填充折射率与微结构光纤12最外层折射率相当的折射率匹配液。微结构光纤12可以是具有不同尺寸和结构分布的普通单模光纤，特种光纤等。

所述控制模块13由计算机、仪器控制单元和仪器控制接口组成，对图像采集器、旋转控制平台等进行控制和操作，完成包含微结构光纤信息的数字全息图记录。搭建完成基于f-p腔的数字全息记录光路11后，利用控制模块13，控制旋转控制平台带着微结构光纤12转动一周，并控制ccd采集包含微结构光纤12信息的数字全息图且存储图。

所述计算显示模块14，对存储的数字全息图进行程序处理，并在线显示微结构光纤12的三维折射率分布信息。根据本发明提供的显微成像新方法，对ccd采集的数字全息图进行处理，可以重建得到微结构光纤12的高分辨率三维折射率分布。

基于f-p腔的数字全息记录光路11中，光束在f-p腔中多次反射、透射，如图3b所示，多次经过微结构光纤，记录微结构光纤的三维信息，光程差依次累积，最终透过f-p腔的光束的复振幅为：

其中，ut为透射光的复振幅，a0是入射f-p腔的光波33复振幅，r为f-p腔的两平行的平面玻璃板内侧的表面反射率，如图3a中的31和32，δ为微结构光纤的相位分布。

则f-p腔多光束干涉时，通过数字全息获得的相位分布为：

其中，n是腔内介质的折射率，d是f-p腔的厚度，λ为光源的波长，φ为入射f-p腔的光波33的入射角34。

光束沿着传播方向，经过微结构光纤12内部的每一点的折射率的累积，就是通过数字全息图获得的相位分布，当微结构光纤12内部及与微结构光纤12周围的环境介质之间的折射率差较小时，光程差是沿光束路径方向折射率的累积，则相位分布与微结构光纤12的折射率分布的关系为：

其中，n(x,y,z)为微结构光纤12内部的折射率分布，z轴是光束传播的方向，n0为微结构光纤12周围的环境介质的折射率。

对于常见的基于马赫曾德尔干涉的数字全息记录光路，光源经分光棱镜分成两束，一束光经过微结构光纤作为物光波，另一束光作为参考光。两束光又经分束棱镜合束，在ccd记录平面相遇干涉形成全息图。光束只经过一次微结构光纤，记录微结构光纤的三维信息，光程差累积一次。对ccd记录的数字全息图，进行数值重建，再现物光波的复振幅分布为umz。同时得到物光波的相位分布为：

其中，umz表示再现的物光波的复振幅分布，re表示取复振幅的实部，im表示取虚部。

基于马赫曾德尔干涉的数字全息记录光路，光束沿着传播方向，经过微结构光纤，相位差在微结构光纤内部沿光束传播方向累积，光束只经过一次微结构光纤，因此相位分布与微结构光纤的折射率分布的关系为：

其中，n(x,y,z)为微结构光纤内部的折射率分布，z轴是光束传播的方向，n0为微结构光纤周围的环境介质的折射率。

对比公式(3)和公式(5)，可以发现，本发明提出的基于f-p腔的微结构光纤高分辨率三维折射率测试方法，光束在f-p腔中多次反射、透射，多次经过微结构光纤，记录微结构光纤的三维信息，光程差依次累积，测量灵敏度更高，因此分辨率更高，从原理上优于常见的基于马赫曾德尔干涉的折射率测试方法。

因此在基于f-p腔的数字全息记录光路11，首先记录f-p腔中入射光的复振幅，记录未放入微结构光纤时数字全息图，在f-p标准具中放入微结构光纤。记录该角度含有微结构光纤12信息的透射光的最优数字全息图，对记录的数字全息图进行数值重建得到其复振幅分布，代入f-p腔的反射率r，根据公式(1)即可获得不含有微结构光纤12信息的相位分布和该角度时含有微结构光纤12信息的相位分布。含有微结构光纤12信息的相位分布减去不含有微结构光纤12信息的相位分布，得到该角度只含有微结构光纤12信息的相位分布。控制模块13，控制旋转控制平台带着微结构光纤12转动一周，同时控制ccd采集包含微结构光纤12信息的数字全息图。利用计算显示模块14，由公式(1)，即可依次得到每一角度的只含有微结构光纤12信息的相位分布，依次将各个截面上所有角度的微结构光纤12的相位分布，进行iradon变换，即可重建得到微结构光纤的三维相位分布δ(x,y,z)。根据公式(3)，可转换得到微结构光纤12的高分辨率三维折射率分布，通过计算显示模块14，即可在线显示微结构光纤12的高分辨率三维折射率分布n(x,y,z)。

本发明所述的，微结构光纤高分辨率三维折射率测试方法，主要包括以下步骤：

第一步：记录f-p腔中入射光的复振幅分布，记录f-p腔中未放入微结构光纤12的数字全息图。在f-p腔中放入微结构光纤12。

第二步：控制模块13控制旋转控制平台，使微结构光纤12转动，记录该角度时放入微结构光纤12后的最优数字全息图，并截取一定大小。

第三步：对记录的数字全息图进行数值重建，根据公式(1)得到相位分布，并对得到的相位分布进行解包裹处理。

第四步：将放入微结构光纤解包裹得到的相位分布减去未放入微结构光纤解包裹得到的相位分布，得到该角度只含有微结构光纤样品相位信息的相位分布。

第五步：控制模块13控制旋转控制平台，使微结构光纤12转动一周，依次重复第二步～第四步，即可得到微结构光纤12不同角度时各个截面上的相位分布。

第六步：通过计算显示模块14，依次将各个截面上所有角度的微结构光纤12的相位分布，进行iradon变换，即可重建得到微结构光纤12的三维相位分布δ(x,y,z)。

第七步：通过计算显示模块14，根据公式(3)，可转换得到微结构光纤12的高分辨率三维折射率分布，通过计算显示模块14，即可在线显示微结构光纤12的高分辨率三维折射率分布n(x,y,z)。

本发明提供的微结构光纤高分辨率三维折射率测试方法，包括基于f-p腔的数字全息图记录、数值重建、解包裹、误差处理、三维相位分布重建和三维折射率转换。与在先技术相比，由于采用了基于f-p腔的数字全息记录光路，光束在f-p腔中多次反射经过微结构光纤，测量的灵敏度更高。该微结构光纤高分辨率三维折射率测试方法具有结构简单、灵敏度高、分辨率更高的优点。

(四)附图说明

图1是一种微结构光纤高分辨率三维折射率测试方法的示意图。包括基于f-p腔的数字全息图记录、数值重建、解包裹、误差处理、三维相位分布重建和三维折射率转换。

图2a是本发明提出的一种微结构光纤高分辨率三维折射率测试方法适用的测量系统结构图。包括基于f-p腔的数字全息记录光路11、微结构光纤12、控制模块13和计算显示模块14。图2b是本发明实施例中基于f-p腔的一种数字全息记录光路的示意图。该实施例中，一种基于f-p腔的数字全息记录光路包括光源21、衰减器22、扩束器23、f-p标准具24、微结构光纤25、显微物镜26、ccd27和计算机28。

图3a是本发明中f-p腔的示意图。f-p腔由两块平行的平面玻璃板组成，在两块平面玻璃板的向里的两侧31和32处镀有反射率为r的介质膜。图3b是光束入射f-p腔，在f-p腔中多次反射、透射的示意图。33是入射光波，34是入射角，35、37、39是透射光的示意图，36、38是反射光的示意图。

图4是本发明实施例中f-p腔与微结构光纤的示意图。旋转控制平台41控制微结构光纤25旋转，微结构光纤25插入填充有折射率与微结构光纤25最外层折射率相当的折射率匹配液的比色皿42中，光束传播方向如43所示。

图5是本发明微结构光纤高分辨率三维折射率测试方法的步骤流程图。

(五)具体实施方式

下面结合具体的实施例来进一步阐述本发明。

实施例1：

图2b给出了一种基于f-p腔的数字全息记录光路的实施例。该光路包括光源21、衰减器22、扩束器23、f-p标准具24、微结构光纤25、显微物镜26、ccd27和计算机28。光束由光源21出发，经过衰减器22，光束能量减弱，光束经过扩束器23，光束直径扩大，光束经过放有微结构光纤25的f-p标准具24时，在f-p标准具24多次反射，多次经过微结构光纤25，经过f-p标准具24的每一束透射光在通过显微物镜26后，在ccd27上干涉叠加，由ccd27记录干涉的数字全息图，并保存在计算机28中。

在本实施例中，优选的，光源21采用波长为532nm的激光光源，并选用相应波长的扩束器23，选用可调的衰减器22，选用腔长为60mm的f-p标准具24，放大倍数为20倍的显微物镜26，像素尺寸为3.45μm×3.45μm的ccd27等器件。

在本实施例中，优选的，f-p标准具24如图3a所示，选用内侧表面镀有反射率为0.8的介质膜的两块平行的平面玻璃板，组成腔长为60mm的f-p标准具24。优选的，微结构光纤25为双芯光纤，并将去除涂覆层的双芯光纤，放入填充有折射率匹配液的比色皿42中，将比色皿42放入f-p标准具24腔中，如图4所示，调整数字全息记录光路，使光束43的光轴经过微结构光纤25。

在本实施例中，优选的，旋转控制平台41选用精密的旋转电机，通过控制模块，精密的控制微结构光纤25绕光纤轴向方向旋转一周。比色皿42上带有陶瓷插芯，可以固定微结构光纤25绕轴转动，填充的折射率匹配液折射率与光纤包层的折射率相同。

在本实施例中，所述基于f-p腔的数字全息记录光路包括光源21、衰减器22、扩束器23、f-p标准具24、显微物镜26、ccd27等。所述微结构光纤25为双芯光纤。所述控制模块由计算机28、仪器控制单元和仪器控制接口组成，对ccd27、控制微结构光纤25转动的旋转控制平台41等进行控制和操作，完成包含微结构光纤25信息的数字全息图记录。所述计算显示模块，对记录的数字全息图按本发明提出的显微成像方法进行程序处理，并在线显示微结构光纤的三维折射率信息。

本实施例中，运用微结构光纤高分辨率三维折射率测试方法，对微结构光纤25的数字全息图进行处理，并在线显示微结构光纤25的三维折射率信息，如图5所示，主要包括以下步骤：

第一步51：记录光束在f-p标准具24中入射光的复振幅分布，并记录f-p标准具24中未放入微结构光纤25时的数字全息图。在f-p标准具24中放入微结构光纤25。

第二步52：控制模块控制精密的旋转电机旋转，使微结构光纤25转动，记录该角度时放入微结构光纤25后的最优数字全息图，并截取一定大小的最优数字全息图。

第三步53：对记录的数字全息图进行数值重建。

第四步54：根据公式(1)得到微结构光纤25的相位分布，并对得到的相位分布进行解包裹处理。

第五步55：将放入微结构光纤25解包裹得到的相位分布减去未放入微结构光纤25解包裹得到的相位分布，得到只含有微结构光纤25相位信息的相位分布。

第六步56：控制模块控制旋转控制平台，使微结构光纤25转动一周，依次重复第二步～第五步，即可得到微结构光纤25不同角度时各个截面上的相位分布。

第七步57：通过计算显示模块，依次将各个截面上所有角度的微结构光纤25的相位分布，进行iradon变换，即可重建得到微结构光纤25的三维相位分布δ(x,y,z)。

第八步58：通过计算显示模块，根据公式(3)，可转换得到微结构光纤25的高分辨率三维折射率分布，通过计算显示模块，即可在线显示微结构光纤25的高分辨率三维折射率分布n(x,y,z)。

以上所述，仅为本发明的优选实施例，但本发明的保护范围并不局限于此。任何本领域的技术人员根据本发明的精神和范围，对本发明进行各种改动和变化，均应包含在本发明权利要求保护范围内。