一种点扩散函数重建方法技术

本发明专利技术公开了一种点扩散函数重建方法，包括：(1)输入目标点扩散函数；(2)针对输入的目标点扩散函数进行采样；(3)将采样得到的每个点看作一个偶极子进行初始相位和初始振幅估计；(4)计算得到偶极子在瞳面上的辐射场；(5)通过辐射场得到瞳面上的振幅和相位分布；(6)通过得到的振幅和相位生成点扩散函数；(7)将点扩散函数与目标点扩散函数比较，若满足要求，则输出得到的振幅和相位分布；若不满足要求，则更新初始相位和初始振幅，返回步骤(4)。本发明专利技术直接由目标点扩散函数得到瞳面上的相位及振幅分布，大大提高了计算效率。并且，由于偶极子辐射的光是向随机方向传播的，该方法可以扩展到多个具有多个光瞳面的显微系统。

A Point Spread Function Reconstruction Method

【技术实现步骤摘要】
一种点扩散函数重建方法
本专利技术设计光学工程领域，具体是涉及一种点扩散函数重建方法。
技术介绍
聚焦光斑(TightlyFocusedSpots)广泛应用于光学存储、光学操控、光学微细加工、全息显示器和远场光学显微镜。对于聚焦光斑的不同的应用场景，人们需要对聚焦光斑进行相应的操控。在超分辨显微成像领域，特别是在受激辐射光淬灭显微镜(STEDMicroscopy)中，对于聚焦光斑的操控尤为重要。如果有一个好的聚焦光斑，那么显微镜的分辨率将会更高。聚焦光斑又被称为点扩散函数(PointSpreadFunction，简称PSF)，是描述光学系统对点光源解析能力的函数。因为点光源在经过任何光学系统后都会由于衍射而形成一个扩大的像点，这个像点就是点扩散函数。在给定输入参数的情况下，利用矢量衍射公式得到系统的点扩散函数是容易的，但是还无法对其进行操控。为了得到人们所需要的点扩散函数，科学家们研究出了点扩散函数的操控技术(也称为点扩散函数重建技术)，即，给定一个目标点扩散函数，通过一定的算法，得到产生该点扩散函数所需的输入参数，这些输入参数构成的函数被称为瞳函数，这些输入参数一般指入射光在瞳面上的振幅和相位。近些年来，点扩散函数的操控得到了很大的发展。比如游尚挺等人2015年在OpticsLetters上面发表的题为《Iterativephase-retrievalmethodforgeneratingstereoarrayofpolarization-controlledfocalspots》的文章提出了一种基于G-S算法的相位反演方法，该方法虽然能够实现相位反演，但是其经过了大量的迭代计算，因此所需的计算时间极长，并且只适用于单个瞳函数系统，很难扩展到需要同时构建多个瞳函数的情况，比如4Pi显微系统。
技术实现思路
针对上述问题，本专利技术提供了一种点扩散函数重建方法。在保证重建点扩散函数的质量的前提下，通过避免迭代，或者用少次数迭代，从而快速获得重构点扩散函数所需的瞳函数，加快了运算速率。本专利技术采用将任意一个目标点扩散函数看作位于不同位置的多个具有不同初始相位和振幅的偶极子的组合。将目标点扩散函数的每一个点(像素)看作一个偶极子，将每个点的强度值看作对应偶极子的振幅，得到一种计算效率高的点扩散函数重建方法。具体如下：一种点扩散函数重建方法，包括如下步骤：(1)输入目标点扩散函数；(2)针对输入的目标点扩散函数进行采样，读取每个点的坐标信息和强度信息(主要是指点的像素值)；(3)将采样得到的每个点看作一个偶极子进行初始相位和初始振幅估计；(4)计算得到偶极子在瞳面上的辐射场；(5)通过辐射场得到瞳面上的振幅和相位分布；(6)通过得到的振幅和相位生成点扩散函数；(7)将步骤(6)中生成的点扩散函数与目标点扩散函数比较，若满足要求，则输出得到的振幅和相位分布；若不满足要求，则更新初始相位和初始振幅，返回步骤(4)～(7)。步骤(1)中，输入的目标点扩散函数，包括目标点扩散函数内各点的坐标信息，强度信息。光的衍射性质会增加两个偶极子之间的串扰，所以很多情况下，采样率不需要太高，举例说明，针对某一目标点扩散函数，虽然可以用一个1000*1000的矩阵来精确描述一个点扩散函数的强度分布，但是实际上，为了进一步降低上述串扰带来的不利影响，只需要一个位置间隔为0.6λ的3*3的偶极子点阵来计算即可，否则，由于串扰的存在，偶极子的分布越密集，计算难度越大，从而导致无法获得一个符合预期的点扩散函数。λ为入射光的波长。作为优选，所述步骤(2)中，针对输入的目标点扩散函数进行采样时采用稀疏采样。作为优选，步骤(2)中，以位置间隔为0.5～0.7λ进行稀疏采样，λ为入射光的波长。作为优选，步骤(3)中，在进行初始振幅估计时，假设每个偶极子的振幅相同。这里的相同，包括数值相同，方向也相同。作为进一步优选，初始相位可以设置为π，初始振幅可以设置为1。本专利技术中，步骤(4)中，计算这些偶极子的辐射场，这些偶极子的远场辐射可以分解为平面波在不同方向上的传播的叠加为位置坐标为θ，的波矢，即第N个偶极子向瞳面上对应坐标为(θ,φ)的点辐射光波的波矢)；且满足：所述辐射场利用复振幅进行表示，则方向上的复振幅可以表示为：θ是入射光光束的孔径角；是瞳面上的方位角；；设定θ的最大值θmax为镜头的汇聚角(对于点扩散函数来说是发散角)，NA表示镜头的数值孔径，n是镜头的像方折射率；M为偶极子的总数；AN和ψN分别是第N个偶极子的初始振幅和初始相位；是第N个偶极子的位置坐标；是第N个偶极子的位置矢量。若M为无限大，上式则可以转化为积分的形式。步骤(5)中，由步骤(4)中得到的辐射场可得到瞳面上的振幅和相位分布。瞳面上点的振幅A和相位α为：为瞳面上坐标为的点对应的复振幅；||代表对取模；angle[]表示对取幅角。由于是相位反演，所以相位项带了一个负号。步骤(6)中，通过得到的振幅和相位生成点扩散函数。点扩散函数通过如下公式计算：其中为镜头焦点附近的电场矢量，“||”为求模运算符；在镜头(比如物镜)具有高数值孔径的条件下，根据矢量衍射理论，物镜焦点附近的电场矢量可以表示为：其中，K＝2π/λ为波数，是以镜头(比如物镜)的理想焦点为原点建立的柱坐标系点坐标；θ是光束的孔径角；是瞳面上的方位角；是入射光在瞳面上的振幅值的分布，即步骤(5)所得的振幅；是入射光的琼斯矩阵；是入射光在瞳面上的相位分布，即步骤(5)所得的相位；是物方到像方的坐标转置矩阵，为已知参数。本专利技术通过即可得到整个点扩散函数。步骤(7)中，将步骤(6)中生成的点扩散函数与目标点扩散函数比较，若结果不理想，则改变初始相位通过迭代，再计算瞳面上的振幅和相位分布，直到获得的瞳面上的振幅和相位能够产生符合预期的点扩散函数。步骤(7)中，具体比较时，计算当前得到的点扩散函数与目标点扩散函数各点强度值的均方差，若均方差在允许范围内，则判定为满足要求；反之判定为不满足要求。该方法的原理是：通过计算结果点扩散函数各点强度值(即振幅的平方)与目标点扩散函数各点强度值的均方差，若均方差在允许范围内，则结果符合预期；若均方差不在允许范围内，则开始迭代，具体的迭代方法为，将计算所得到的振幅函数替换为一个高斯分布或者均匀分布的振幅函数，然后向前计算点扩散函数，并保留相位分布作为初始估计，然后再重新计算瞳函数，这样迭代，直到点扩散函数达到设计目标。步骤(7)中，当前得到的点扩散函数不满足要求，更新初始相位和初始振幅时：将得到的振幅函数替换为一个高斯分布或者均匀分布的振幅函数，完成初始振幅的更新；将当前得到的相位分布作为初始相位，完成初始相位的更新。作为优选，若输入的目标点扩散函数为图片，则应将其转换为对应的矩阵，再对其进行稀疏采样。本专利技术的方法可以利用计算机完成，即通过计算机编程进行迭代运算，利用本专利技术的方法得到最佳的点扩散函数参数(入射光在瞳面上的振幅和相位等)后，可以利用下述装置对得到的点扩散函数参数进行验证：基于上述方法，本专利技术还提出了一种实现点扩散函数操控的装置。该装置包括：激光光源、光调制模块、透镜、接收模块。本专利技术中，正向入射光为相干准直线偏光，经光调制模块，然后通过透镜出射，到达接收模块。本专利技术中，所述光调制本文档来自技高网...

【技术保护点】
1.一种点扩散函数重建方法，其特征在于，包括如下步骤：(1)输入目标点扩散函数，读取每个点的坐标信息和强度信息；(2)针对输入的目标点扩散函数进行采样；(3)将采样得到的每个点看作一个偶极子进行初始相位和初始振幅估计；(4)计算得到偶极子在瞳面上的辐射场；(5)通过辐射场得到瞳面上的振幅和相位分布；(6)通过得到的振幅和相位生成点扩散函数；(7)将步骤(6)中生成的点扩散函数与目标点扩散函数比较，若满足要求，则输出得到的振幅和相位分布；若不满足要求，则更新初始相位和初始振幅，返回步骤(4)。

【技术特征摘要】
1.一种点扩散函数重建方法，其特征在于，包括如下步骤：(1)输入目标点扩散函数，读取每个点的坐标信息和强度信息；(2)针对输入的目标点扩散函数进行采样；(3)将采样得到的每个点看作一个偶极子进行初始相位和初始振幅估计；(4)计算得到偶极子在瞳面上的辐射场；(5)通过辐射场得到瞳面上的振幅和相位分布；(6)通过得到的振幅和相位生成点扩散函数；(7)将步骤(6)中生成的点扩散函数与目标点扩散函数比较，若满足要求，则输出得到的振幅和相位分布；若不满足要求，则更新初始相位和初始振幅，返回步骤(4)。2.根据权利要求1所述的点扩散函数重建方法，其特征在于，所述步骤(2)中，针对输入的目标点扩散函数进行采样时采用稀疏采样。3.根据权利要求2所述的点扩散函数重建方法，其特征在于，以位置间隔为0.5～0.7λ进行稀疏采样，λ为入射光的波长。4.根据权利要求1所述的点扩散函数重建方法，其特征在于，步骤(3)中，在进行初始振幅估计时，假设每个偶极子的振幅相同。5.根据权利要求1所述的点扩散函数重建方法，其特征在于，步骤(4)中，所述辐射场利用复振幅进行表示：θ是入射光光束的孔径角；是瞳面上的方位角；且θ、满足：为第N个偶极子向瞳面上对应坐标为的点辐射光波的波矢；设定θ的最大值为θmax，即镜头的汇聚角，NA表示镜头的数值孔径，n是镜头的像方折射率；M为偶极子的总数；AN和ψN分别是第N个偶极子的初...

【专利技术属性】
技术研发人员：郝翔，刘鑫，匡翠方，刘旭，
申请(专利权)人：浙江大学，
类型：发明
国别省市：浙江,33