叠层成像中物体片层透射函数的恢复方法、装置和存储介质制造方法及图纸

本发明专利技术提供一种叠层成像中物体片层透射函数的恢复方法、装置和存储介质，方法包括：获得被测物体不同扫描位置的多张实验衍射像；通过迭代操作由实验衍射像恢复物体各片层的透射函数，各迭代操作包括：在当前迭代中确定各片层待修正透射函数的振幅和相位，并估计入射波函数；基于实验衍射像、传播函数以及各片层的透射函数的振幅和相位，利用成像公式获得被测物体在不同扫描位置的模拟衍射像；基于模拟衍射像和实验衍射像的差值，计算衍射像的误差函数；基于误差函数分别对各片层透射函数的振幅与相位求梯度公式，计算各片层的透射函数的振幅和相位的修正量；基于修正量来修正各片层透射函数的振幅和相位，获得当前迭代后的模拟衍射像。衍射像。衍射像。

【技术实现步骤摘要】
叠层成像中物体片层透射函数的恢复方法、装置和存储介质


[0001]本专利技术电子显微镜成像
，尤其涉及一种叠层成像中物体片层透射函数和入射波的恢复方法、装置和存储介质。

技术介绍

[0002]叠层成像(Ptychography)是一种迭代相位恢复算法，用于X射线及电子显微像的透射成像技术中，可以通过算法，恢复待测物的波函数以及入射光的函数形式。
[0003]在量子力学中，量子系统的量子态可以用波函数(wave function)来描述，波函数是一种复数函数。入射电子波函数、物体的透射函数、出射波函数，都是这样的复数函数。其中，入射电子波函数可写作：
[0004]Φ(R)＝∫Aper(k)exp[
‑
iχ(k)]exp(
‑
ikR)dk；
[0005]其中，Aper(k)为光阑函数；k表示频率域空间的矢量，R是空间域的矢量，二者可通过傅里叶变换相互转换；χ(k)代表聚光镜的像差函数，包括焦距及像散等各种成像因素。近平行的入射束透过光阑后，经过聚光镜会照射到样品表面，入射时的光斑通常形成较小的光斑，与焦距、聚光镜的像散有关。
[0006]叠层成像技术是通过对待测样品以小于入射光直径的步长进行扫描，采集穿透过样品后的衍射像，利用迭代计算可以重建出照明光和样品波函数分布，是一种理论分辨率为衍射极限、利用算法恢复物体透过函数的无透镜相位成像技术。近些年随着相关研究的蓬勃发展，人们对叠层成像算法的不断优化，在可见光、X射线和电子束等领域已被广泛应用于相位成像、波前恢复和光学计量。
[0007]基于叠层成像的针对物体的透射函数和入射波的恢复算法，包括有：仅考虑了一个待测对象的PIE(Ptychographic Iterative Engine，层叠衍射迭代引擎)、ePIE(extending Ptychographic Iterative Engine，扩展的层叠衍射迭代引擎)、DM(Difference Map，差异图)、WDD(wigner distribution deconvolution，Wigner分布反卷积)和LSQ
‑
ML叠层成像(Linear least
‑
square maximum
‑
likelihood phychography，线性最小二乘最大似然叠层成像)等算法，以及考虑了物体为三维多片层结构的3ePIE、3DM(3D Difference Map)和3ML(3D maximum likelihood，三维最大似然)算法。这些算法都是从记录的衍射像(或称衍射图)来恢复物体透射函数及入射波的方法。最早期的PIE算法，不能对入射波进行修正，并且无法对多层物体的透射函数进行恢复；而早期的WDD算法是以去卷积方式进行恢复。
[0008]PIE、ePIE、DM、LSQ
‑
MAL、3ePIE、3DM以及3ML算法在恢复目标函数的时候，都是使用迭代恢复算法。其中，PIE算法仅考虑了一个待测对象，ePIE和3ePIE算法，从已求得(或简单估计)的物体透射函数和入射波函数出发计算像面上的波函数，并以实验衍射像的振幅的开方替换该像面上的波函数的振幅部分，再逆向向物体和入射光方向传播，从而修正之前求得的物体和入射光的波函数。ePIE和3ePIE算法需要的迭代次数多，收敛慢导致计算速度慢。相比于ePIE算法，虽然3ePIE算法可以对多层物体的透射函数进行恢复，但由于需要迭
代的次数仍然过多，使得整个算法运行速度依然较慢，导致在应用上受到一定限制。
[0009]用于物体透射函数和入射波恢复的DM、3DM、3ML以及LSQ
‑
MAL算法是以模拟衍射像与实验衍射像的偏差来更新迭代物体与入射波函数的方法。图1为基于模拟衍射像与实验衍射像之间偏差的迭代恢复方法的示意性原理图，如图1所示，这些恢复方法中，首先对物体的目标函数O(R)(或O(k))和入射光的目标函数P(R)(或P(k))在空间域R(或频率域k)空间设置个猜测的初值，之后利用已知的成像公式，得到每个扫描点对应的衍射像I(k)(或也可以表示为图像空间的衍射像I(R))，与获得并输入的实验衍射像进行比较，从而得到模拟衍射像和实验衍射像两者的差ΔI(R)或ΔI(q)(此处，q表示频率域)，再基于该差值计算物体波函数修正量ΔO(R)和入射波修正量ΔP(R)。
[0010]在现有技术中，DM和3DM算法是使用Fienup的混合输入输出算法(Fienup
’
s hybrid input
‑
output algorithm)计算物体透射函数的修正。
[0011]LSQ
‑
MAL算法中，其目标函数是梯度形式的最大化概率密度函数，相当于是取衍射像的开根号的差为最小化的目标函数，且该算法中只考虑了样品只有一个物体波函数的推导，因此LSQ
‑
MAL算法的目标函数恢复效率和效果都不够理想。
[0012]3ML算法是对3ePIE算法的延伸，其将多层切片技术用于解释样品内的传播效应，使成像超越景深限制。多层切片技术允许重建多个目标切片和入射光照，以及未知物体厚度的检索。数值研究表明，较小的扫描步长会显著增加景深，使用多层切片方法可进一步扩展景深。3ML算法中，最小化的目标函数相当于是取实验衍射像与模拟衍射像的开根号的差之和；优化的目标是每层波函数，是复数形式，包括了实部和虚部两个部分。现有的3ML算法虽然已经可以对三维物体进行恢复，但是该算法在推导过程中是对O(物体)和P(照明光)来求偏导，仍需要较多的迭代次数，其第一次的恢复结果无法产生正确的波函数相位，从而影响图像恢复的速度。
[0013]如何提供一种能够进一步降低迭代次数从而能够对三维物体的透射函数进行快速恢复的方法，是一个有待解决的问题。

技术实现思路

[0014]鉴于此，本专利技术实施例提供了一种叠层成像中物体片层透射函数和入射波的恢复方法和装置，以消除或改善现有技术中存在的一个或更多个缺陷，实现对三维物体的透射函数的快速恢复，同时可以实现对入射波的快速恢复。
[0015]本专利技术的一个方面提供了一种叠层成像中物体波和入射波的恢复方法，该方法包括以下步骤：
[0016]获得被测物体不同扫描位置的多张实验衍射像；
[0017]通过多次迭代操作由实验衍射像恢复物体各片层的透射函数，各次迭代操作包括如下步骤：
[0018]在当前迭代中确定被测物体各片层待修正透射函数的振幅和相位，并估计入射波函数；
[0019]基于实验衍射像、传播函数以及被测物体各片层的透射函数的振幅和相位，利用成像公式获得被测物体在不同扫描位置的模拟衍射像；
[0020]基于模拟衍射像和实验衍射像的差值，计算衍射像的误差函数；
[0021]基于衍射像的误差函数分别对各片层透射函数的振幅与相位求梯度公式，利用求得的梯度公式计算各片层的透射函数的振幅和相位的修正本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种叠层成像中物体片层透射函数的恢复方法，其特征在于，该方法包括以下步骤：获得被测物体不同扫描位置的多张实验衍射像；通过多次迭代操作由实验衍射像恢复物体各片层的透射函数，各次迭代操作包括如下步骤：在当前迭代中确定被测物体各片层待修正透射函数的振幅和相位，并估计入射波函数；基于实验衍射像、传播函数以及被测物体各片层的透射函数的振幅和相位，利用成像公式获得被测物体在不同扫描位置的模拟衍射像；基于模拟衍射像和实验衍射像的差值，计算衍射像的误差函数；基于衍射像的误差函数分别对各片层透射函数的振幅与相位求梯度公式，利用求得的梯度公式计算各片层的透射函数的振幅和相位的修正量；以及基于获得的修正量来修正各片层透射函数的振幅和相位，从而获得当前迭代优化后的模拟衍射像。2.根据权利要求1所述的方法，其特征在于，在初次迭代中估计的入射波函数表示为：Φ0(R)＝∫A
aper
(k)exp[
‑
i
χ
(k)]exp(
‑
ikR)dk；其中，A
aper
为聚光镜的光阑函数，χ为聚光镜的像差函数，k表示频率域空间的矢量，R是空间域空间的矢量；利用成像公式得出的不同扫描位置的模拟衍射像，模拟衍射像表示为I
N
(k)或I
N
(R)：I
N
(k)＝φ
N
(k)[φ
N
(k)]
*
；其中，I
N
(k)是频率域的模拟衍射像，为实数形式的矩阵；φ
N
(k)为频率域的第N个片层的出射波函数；I
N
(R)为空间域的模拟衍射像，其为通过反傅里叶变换从I
N
(k)计算得到的复数形式的矩阵；Φ
N
(R)和为空间域的第N片层的出射波函数及其共轭形式；R和R
′
均为空间域的位置矢量；Φ
N
(R)表示为：其中，P
n
(R)表示第n层的真空传播函数；t
n
(R)表示第n层的透射函数，n＝1,2,3,
…
,N；第n层透射函数t
n
(R)为复数矩阵。3.根据权利要求2所述的方法，其特征在于，所述透射函数符合如下公式：t
n
(R)＝A
n
(R)exp[iV
n
(R)]；其中，A
n
(R)为第n片层透射函数的振幅部分，V
n
(R)为第n片层透射函数的相位部分。4.根据权利要求1所述的方法，其特征在于，所述方法还包括：基于衍射像的误差函数对入射波函数求梯度公式，利用求得的梯度公式计算入射波函数的修正量，并基于计算的入射波函数的修正量来修正估计的入射波函数。5.根据权利要求4所述的方法，其特征在于，所述基于获得的修正量来修正各片层透射函数的振幅和相位包括：基于如下公式修正
各片层透射函数的振幅和相位：各片层透射函数的振幅和相位：所述基于计算的入射波函数的修正量来修正估计的入射波函数包括：基于如下公式来修正估计的入射波函数：其中，和分别表示整个样品区域第n片层在第j+1次迭代和第j次迭代中透射函数的相位部分；表示整个样品区域第n片层在第j次迭代的透射函数的相位修正量；和分别表示整个样品区域第n片层在第j+1次迭代和第j次迭代中透射函数的振幅部分；表示整个样品区域第n片层在第j次迭代的透射函数的振幅修正量；和分别表示第j+1次和第j次迭代中的入射波函数；表示第j次迭代的入射波函数的修正量；以及分别为第j次迭代中透射函数的相位、振幅以及入射波函数的常数修正项。6.根据权利要求5所述的方法，其特征在于，第N
‑
n+1片层的透射函数的相位和振幅的修正量分别满足以下公式：ΔV
N
‑
n+1
(R)＝4Im{t
N
‑
n+1
(R)Φ
N
‑
n
(R)W
n
(R)}；ΔA
N
‑
n+1
(R)＝
‑
4Re{exp[iV
N...

【专利技术属性】
技术研发人员：林芳，高博文，张晓涵，林上港，谭穗妍，
申请(专利权)人：华南农业大学，
类型：发明
国别省市：