局域轨道函数的三维重构方法及装置

1.本技术涉及三维成像技术领域，特别涉及一种局域轨道函数的三维重构方法及装置。


背景技术：

2.通常情况下样品的电子显微镜照片中蕴含着样品丰富的、难以直观获取的信息，应用三维重构算法从样品一系列电子显微镜投影照片中获取样品的三维结构信息对在根本的层面上理解材料成分和性能的关系有很大帮助。
3.近年来，随着数据采集方法、迭代三维重建算法和后处理方法的发展，aet(atomic electrical tomography，原子尺度电子层析成像)已成为三维和四维原子尺度结构表征的有力工具，它提供了在原子水平上关联材料结构和性质的能力，相关技术中的aet能够以亚埃精度确定三维原子坐标和元素种类，并揭示其在动态过程中的原子尺度时间演化。
4.然而，相关技术中，aet算法从重构出的三维密度矩阵中通过寻峰来获取原子的三维坐标，对硬件的需求较高，且在寻峰过程中，需要进行繁琐的人为干预，难以避免因人为干预造成的误差，同时，在直接对原子的三维坐标进行重构的同时，无法对采集图像时存在的样品漂移、样品台的机械倾转误差进行矫正，准确度较差，有待改善。


技术实现要素：

5.本技术提供一种局域轨道函数的三维重构方法及装置，以解决相关技术中只能从重构出的三维密度矩阵中获取原子的三维坐标，且无法对误差进行校正，导致重构的过程对硬件需求较高，且重构的三维坐标精度较差的技术问题。
6.本技术第一方面实施例提供一种局域轨道函数的三维重构方法，包括以下步骤：采集多个倾转角度下的样品的图像数据；基于所述图像数据，在实空间等间隔进行撒点并利用线性累加，得到所述多个倾转角度的每个倾转角度下的计算图像；以及根据每个倾转角度下的计算图像计算损失函数，并获取所述损失函数关于待优化参数的梯度，且根据所述梯度优化所述待优化参数，筛选出满足预设条件的原子，重新计算新的损失函数，直至满足收敛条件，在所述实空间重构局域轨道函数的中心的三维空间坐标和局域轨道函数的形状，得到三维重构结果。
7.可选地，在本技术的一个实施例中，采集所述多个倾转角度下的样品的图像数据，包括：获取所述多个倾转角度下的样品的初始图像数据；对所述初始图像数据进行对中合轴和降噪处理，并由处理后的图像归一化，得到所述图像数据。
8.可选地，在本技术的一个实施例中，所述损失函数的计算公式为：
[0009][0010]
其中，w代表损失函数，m代表倾转角度的总数，j代表倾转角度的序号，i代表局域轨道函数序号，p代表图像的边长，u代表图像中每个像素的横坐标，v代表图像中每个像素的纵坐标，fj(u,v)代表第j个角度下计算得到的图像，bj(u,v)代表第j个角度下实验所得的
图像。
[0011]
可选地，在本技术的一个实施例中，所述筛选出满足预设条件的原子，包括：在每一步迭代更新完参数后，检测所有局域轨道函数的参数在检测到任一局域轨道函数的参数小于由所述所有局域轨道函数参数得到的阈值时，删除所述任一局域轨道函数的同时，通过建立二叉树得到局域轨道函数中心距离小于预设像素的局域轨道函数对索引，删除局域轨道函数对中的任意一个局域轨道函数；每一步迭代更新完参数并删除完所述局域轨道函数后，将每个局域轨道函数的参数均以预设概率减小为预设倍数，且设置保护时间，使得在所述保护时间内，不允许执行筛选操作和删除操作。
[0012]
可选地，在本技术的一个实施例中，所述待优化参数包括每个局域轨道函数中心的三维空间坐标、描述其形状的参数、每个转角对应的三个欧拉角、每个角度下样品的漂移、样品台的机械倾转偏差中的至少一项。
[0013]
本技术第二方面实施例提供一种局域轨道函数的三维重构装置，包括：采集模块，用于采集多个倾转角度下的样品的图像数据；累加模块，用于基于所述图像数据，在实空间等间隔进行撒点并利用线性累加，得到所述多个倾转角度的每个倾转角度下的计算图像；以及重构模块，用于根据每个倾转角度下的计算图像计算损失函数，并获取所述损失函数关于待优化参数的梯度，且根据所述梯度优化所述待优化参数，筛选出满足预设条件的原子，重新计算新的损失函数，直至满足收敛条件，在所述实空间重构局域轨道函数的中心的三维空间坐标和局域轨道函数的形状，得到三维重构结果。
[0014]
可选地，在本技术的一个实施例中，所述采集模块包括：获取单元，用于获取所述多个倾转角度下的样品的初始图像数据；降噪单元，用于对所述初始图像数据进行对中合轴和降噪处理，并由处理后的图像归一化，得到所述图像数据。
[0015]
可选地，在本技术的一个实施例中，所述损失函数的计算公式为：
[0016][0017]
其中，w代表损失函数，m代表倾转角度的总数，j代表倾转角度的序号，i代表局域轨道函数序号，p代表图像的边长，u代表图像中每个像素的横坐标，v代表图像中每个像素的纵坐标，fj(u,v)代表第j个角度下计算得到的图像，bj(u,v)代表第j个角度下实验所得的图像。
[0018]
可选地，在本技术的一个实施例中，所述重构模块包括：检测单元，用于在每一步迭代更新完参数后，检测所有局域轨道函数的参数删除单元，用于在检测到任一局域轨道函数的参数小于由所述所有局域轨道函数参数得到的阈值时，删除所述任一局域轨道函数的同时，通过建立二叉树得到局域轨道函数中心距离小于预设像素的局域轨道函数对索引，删除局域轨道函数对中的任意一个局域轨道函数；保护单元，用于在每一步迭代更新完参数并删除完所述局域轨道函数后，将每个局域轨道函数的参数均以预设概率减小为预设倍数，且设置保护时间，使得在所述保护时间内，不允许执行筛选操作和删除操作。
[0019]
可选地，在本技术的一个实施例中，所述待优化参数包括每个局域轨道函数中心的三维空间坐标、描述其形状的参数、每个转角对应的三个欧拉角、每个角度下样品的漂移、样品台的机械倾转偏差中的至少一项。
[0020]
本技术第三方面实施例提供一种电子设备，包括：存储器、处理器及存储在所述存
储器上并可在所述处理器上运行的计算机程序，所述处理器执行所述程序，以实现如上述实施例所述的局域轨道函数的三维重构方法。
[0021]
本技术第四方面实施例提供一种计算机可读存储介质，所述计算机可读存储介质存储计算机指令，所述计算机指令用于使所述计算机执行如上述实施例所述的局域轨道函数的三维重构方法。
[0022]
本技术实施例基于采集的多个倾转角度下的样品的图像数据，获得多个倾转角度下的计算图像，进而计算出计算图像的损失函数，求取优化参数的梯度，从而获得优化参数，在经过反复筛选和计算后，直至损失函数满足收敛条件，进而获得三维重构结果，可以精简三维坐标的重构过程，不仅降低了对硬件的要求，还能减少繁琐的人为干预，节约劳动力成本，同时在迭代过程中，还可以矫正样品漂移和样品台的机械倾转误差，进而提高三维坐标重构的准确性。由此，解决了相关技术中只能从重构出的三维密度矩阵中获取原子的三维坐标，且无法对误差进行校正，导致重构的过程对硬件需求较高，且重构的三维坐标精度较差的技术问题。
[0023]
本技术附加的方面和优点将在下面的描述中部分给出，部分将从下面的描述中变得明显，或通过本技术的实践了解到。
附图说明
[0024]
本技术上述的和/或附加的方面和优点从下面结合附图对实施例的描述中将变得明显和容易理解，其中：
[0025]
图1为根据本技术实施例提供的一种局域轨道函数的三维重构方法的流程图；
[0026]
图2为根据本技术一个实施例的局域轨道函数的三维重构方法的流程图；
[0027]
图3为根据本技术一个实施例的待重构的10000个原子构成的小颗粒在25
°
、0
°
和-25
°
的模拟示意图；
[0028]
图4为根据本技术一个实施例的局域轨道函数的三维重构方法的初始输入的原子坐标构成的散点示意图；
[0029]
图5为根据本技术一个实施例的局域轨道函数的三维重构方法的初始输入的局域轨道函数在倾转角度分别为25
°
、0
°
和-25
°
的计算示意图；
[0030]
图6为根据本技术一个实施例的局域轨道函数的三维重构方法的迭代过程中损失函数的值的折线示意图；
[0031]
图7为根据本技术一个实施例的局域轨道函数的三维重构方法的收敛后获得的原子模型示意图；
[0032]
图8为根据本技术一个实施例的局域轨道函数的三维重构方法的收敛时计算图像与实验图像的差值；
[0033]
图9为根据本技术一个实施例的局域轨道函数的三维重构方法的收敛计算所得原子坐标与真实坐标之间的距离直方图；
[0034]
图10为根据本技术实施例提供的一种局域轨道函数的三维重构装置的结构示意图；
[0035]
图11为根据本技术实施例提供的电子设备的结构示意图。
具体实施方式
[0036]
下面详细描述本技术的实施例，所述实施例的示例在附图中示出，其中自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。下面通过参考附图描述的实施例是示例性的，旨在用于解释本技术，而不能理解为对本技术的限制。
[0037]
下面参考附图描述本技术实施例的局域轨道函数的三维重构方法及装置。针对上述背景技术中心提到的相关技术只能从重构出的三维密度矩阵中获取原子的三维坐标，且无法对误差进行校正，导致重构的过程对硬件需求较高，且重构的三维坐标精度较差的技术问题，本技术提供了一种局域轨道函数的三维重构方法，在该方法中，基于采集的多个倾转角度下的样品的图像数据，获得多个倾转角度下的计算图像，进而计算出计算图像的损失函数，求取优化参数的梯度，从而获得优化参数，在经过反复筛选和计算后，直至损失函数满足收敛条件，进而获得三维重构结果，可以精简三维坐标的重构过程，不仅降低了对硬件的要求，还能减少繁琐的人为干预，节约劳动力成本，同时在迭代过程中，还可以矫正样品漂移和样品台的机械倾转误差，进而提高三维坐标重构的准确性。由此，解决了相关技术中只能从重构出的三维密度矩阵中获取原子的三维坐标，且无法对误差进行校正，导致重构的过程对硬件需求较高，且重构的三维坐标精度较差的技术问题。
[0038]
具体而言，图1为本技术实施例所提供的一种局域轨道函数的三维重构方法的流程示意图。
[0039]
如图1所示，该局域轨道函数的三维重构方法包括以下步骤：
[0040]
在步骤s101中，采集多个倾转角度下的样品的图像数据。
[0041]
在实际执行过程中，本技术实施例可以拍摄样品在一系列不同转角下的haadf(high-angle annular dark-field imaging，高角环形暗场像)图像，进而获得样品多个倾转角度下的样品的图像数据，便于后续对图像数据进行处理，并进行三维坐标重构。
[0042]
可选地，在本技术的一个实施例中，采集多个倾转角度下的样品的图像数据，包括：获取多个倾转角度下的样品的初始图像数据；对初始图像数据进行对中合轴和降噪处理，并由处理后的图像归一化，得到图像数据。
[0043]
具体地，本技术实施例通过拍摄获得多个倾转角度下的样品的初始图像数据后，为保证后续三维坐标重构的准确性，可以对初始图像数据进行对中和轴和降噪处理，在便于后续计算的同时，避免噪点对三维坐标重构的影响，再通过线性归一化，得到可用于后续计算的图像数据。
[0044]
其中，线性归一化的具体过程如下：
[0045][0046]
其中，m代表倾转角度的总数，j代表倾转角度的序号，u代表图像中每个像素的横坐标，v代表图像中每个像素的纵坐标，bj(u,v)代表第j个角度下实验得到的图像。
[0047]
在步骤s102中，基于图像数据，在实空间等间隔进行撒点并利用线性累加，得到多个倾转角度的每个倾转角度下的计算图像。
[0048]
作为一种可能实现的方式，本技术实施例首先可以基于上述步骤获得的图像数据，在实空间以一定规律设置离散化、有限数目的局域轨道函数作为迭代过程的初始输入，其中，每个局域轨道函数的强度h可以设置为1e-5
，宽度可以设置为1.4，各角度下的样品漂
移量(uj,vj)和欧拉转角的偏移量皆可以设置为0。
[0049]
具体地，本技术实施例可以将不同角度下的局域轨道函数中心的三维坐标表示为：
[0050][0051]
其中，ψj、θj、为第j个角度对应的分布绕x轴、y轴、z轴的欧拉角，(xi,yi,zi)为第i个局域轨道函数在未旋转(三个欧拉角皆为0)时的中心三维位置坐标，(u
ij
,v
ij
,w
ij
)为第i个局域轨道函数在第j个角度下的中心三维位置坐标。
[0052]
进一步地，本技术实施例可以计算各角度下局域轨道中心的坐标：
[0053][0054]
进而，本技术实施例可以根据上述坐标，将局域轨道函数表示为三维高斯函数，计算每个局域轨道函数在各个角度下的值：
[0055][0056]
其中，d
ij
表示第i个局域轨道函数在第j个角度下在实空间(u,v,w)位置处的值，hi和bi为待优化参数，分别代表第i个局域轨道函数的强度与宽度，(u,v,w)为实空间三维位置坐标，(u
ij
,v
ij
,w
ij
)为第i个局域轨道函数在第j个角度下的中心三维位置坐标，(uj,vj,wj)为第j个角度下局域轨道函数中心相对于真实位置(u
ij
,v
ij
,w
ij
)的三维方向漂移。
[0057]
再通过线性累加，本技术实施例可以得到计算图像：
[0058][0059]
其中，n代表局域轨道函数总个数，d
ij
表示第i个局域轨道函数在第j个角度下在实空间(u,v,w)位置处的值。
[0060]
在步骤s103中，根据每个倾转角度下的计算图像计算损失函数，并获取损失函数关于待优化参数的梯度，且根据梯度优化待优化参数，筛选出满足预设条件的原子，重新计算新的损失函数，直至满足收敛条件，在实空间重构局域轨道函数的中心的三维空间坐标和局域轨道函数的形状，得到三维重构结果。
[0061]
在实际执行过程中，本技术实施例可以利用梯度的优化算法求解参数，并通过多次计算和筛选，在实空间重构局域轨道函数的中心的三维空间坐标和局域轨道函数的形状，得到三维重构结果。
[0062]
具体地，本技术实施例可以利用上述步骤中计算获得的计算图像，计算损失函数，进而获取损失函数关于待优化参数的梯度，利用梯度的优化算法，优化待优化参数，进而筛选出满足预设条件的原子，并重复计算损失函数，直至满足收敛条件，从而在实空间重构局域轨道函数的中心的三维空间坐标和局域轨道函数的形状，得到三维重构结果。
[0063]
本技术实施例可以直接获得样品原子的三维坐标，跳过了先得到三维密度矩阵再寻峰得到三维坐标的过程，不仅降低了对硬件的要求，还去除了寻峰过程中的繁琐人为干预，节约了劳动成本，同时在迭代过程中，还可以矫正样品漂移和样品台的机械倾转误差，进而提高三维坐标重构的准确性。
[0064]
需要注意的是，预设条件可以由本领域技术人员根据实际情况进行设置，在此不做具体限制。
[0065]
可选地，在本技术的一个实施例中，损失函数的计算公式为：
[0066][0067]
其中，w代表损失函数，m代表倾转角度的总数，j代表倾转角度的序号，i代表局域轨道函数序号，p代表图像的边长，u代表图像中每个像素的横坐标，v代表图像中每个像素的纵坐标，fj(u,v)代表第j个角度下计算得到的图像，bj(u,v)代表第j个角度下实验所得的图像。
[0068]
具体地，本技术实施例可以通过上述步骤中的计算图像计算损失函数，并将损失函数写为关于计算图像和实验图像的函数：
[0069][0070]
其中，w代表损失函数，m代表倾转角度的总数，j代表倾转角度的序号，i代表局域轨道函数序号，p代表图像的边长，u代表图像中每个像素的横坐标，v代表图像中每个像素的纵坐标，fj(u,v)代表第j个角度下计算得到的图像，bj(u,v)代表第j个角度下实验所得的图像，mj(u,v)为第j个角度下计算图像与实验图像之差。
[0071]
可选地，在本技术的一个实施例中，待优化参数包括每个局域轨道函数中心的三维空间坐标、描述其形状的参数，每个转角对应的三个欧拉角，每个角度下样品的漂移、样品台的机械倾转偏差中的至少一项。
[0072]
在获取损失函数后，本技术实施例可以求出损失函数关于局域轨道函数中心三维坐标(xi,yi,zi)、强度hi、宽度bi、每个角度下样品漂移量(uj,uj)、样品台角度偏差量ψj、θj、等参数的梯度。
[0073]
进一步地，本技术实施例可以利用计算获得的梯度对目标参数进行更新：
[0074][0075][0076][0077][0078][0079]
[0080][0081][0082][0083][0084]
其中，和是各个参数的学习率。
[0085]
需要注意的是，梯度的求取可以利用具有自动求导功能的软件库实现，也可以通过解析表达式实现。
[0086]
可选地，在本技术的一个实施例中，筛选出满足预设条件的原子，包括：在每一步迭代更新完参数后，检测所有局域轨道函数的参数在检测到任一局域轨道函数的参数小于由所有局域轨道函数参数得到的阈值时，删除任一局域轨道函数的同时，通过建立二叉树得到局域轨道函数中心距离小于预设像素的局域轨道函数对索引，删除局域轨道函数对中的任意一个局域轨道函数；每一步迭代更新完参数并删除完局域轨道函数后，将每个局域轨道函数的参数均以预设概率减小为预设倍数，且设置保护时间，使得在保护时间内，不允许执行筛选操作和删除操作。
[0087]
可以理解的是，筛选出满足预设条件的原子包括删除冗余局域轨道函数和筛选局域轨道函数。
[0088]
其中，删除冗余局域轨道函数的方法为：在每一步迭代更新完参数后，检查所有局域轨道函数的参数h，若某局域轨道函数的参数h小于小于由所有局域轨道函数参数得到的阈值时，则删除该局域轨道函数，同时通过建立二叉树的方法，得到局域轨道函数中心距离预设像素的局域轨道函数对索引，删除局域轨道函数对中的任意一个局域轨道函数。
[0089]
需要注意的是，阈值可以由本领域技术人员根据实际情况进行设置，也可以设置为参考值，如所有局域轨道函数参数h最大值的0.01倍；预设像素可以由本领域技术人员根据实际情况进行设置，也可以设置为参考值，如2像素。
[0090]
筛选局域轨道函数的方法为：在每一步迭代更新完参数并删除完局域轨道函数后，每个局域轨道函数的参数h预设概率减小为预设倍数，且设置保护时间，在保护时间内局域轨道函数不可被筛选和删除。
[0091]
需要注意的是，预设概率和预设倍数可以由本领域技术人员根据实际情况进行设置，也可以设置为参考值，如每个局域轨道函数的参数h皆以0.02的概率减小为0.1倍；保护事件可以由本领域技术人员根据实际情况进行设置，也可以设置为参考值，如设置保护时间t＝50。
[0092]
下面结合图2至图9所示，以一个具体实施例对本技术实施例的局域轨道函数的三维重构方法的工作原理进行详细阐述。
[0093]
如图2所示，以重构10000个原子组成的小颗粒为例，本技术实施例包括以下步骤：
[0094]
步骤s201：采集图像数据。如图3所示，本技术实施例需要重构10000个原子组成的小颗粒，通过在角度为
±
25
°
、
±
20
°
、
±
15
°
、
±5°
、0
°
时倾转，获得其模拟图像，图3中第1、2、
3列分别对应25
°
、0
°
和-25
°
的图像。
[0095]
步骤s202：计算图像。本技术实施例可以基于图像数据，在实空间等间隔进行撒点并利用线性累加，得到多个倾转角度的每个倾转角度下的计算图像。如图4所示，初始撒点集合为等间隔原子构成的圆柱形，并参考图5所示，利用公式得到迭代初始的各个角度下的计算图像，其中，第1、2、3列分别对应25
°
、0
°
和-25
°
的图像。
[0096]
步骤s203：计算损失函数，并对目标参数进行迭代更新。本技术实施例可以计算损失函数及目标参数的梯度，并利用以下公式对目标参数进行迭代更新：
[0097][0098][0099][0100][0101][0102][0103][0104][0105][0106][0107]
其中，和是各个参数的学习率。
[0108]
步骤s204：删除、筛选局域轨道函数。进一步地，本技术实施例可以删除强度h小于最强局域轨道函数强度hmax*0.01的局域轨道函数，同时以0.02的概率重置该局域轨道函数，使局域轨道函数强度h减小为原来的0.1倍，设置保护周期t＝50。
[0109]
通过循环计算损失函数，并更新各目标参数，删除和筛选局域轨道函数的操作，直至局域轨道函数数目收敛，实现三维坐标重构。
[0110]
其中，损失函数的值随迭代过程的折线图如图6所示；最终收敛时得到的原子模型的示意图如图7所示；计算图像与实验图像的差值参考图如图8所示，其中，第1列为计算图，第2列为实验图，第3列为二者差值，第1行为倾转角度为25
°
下的数据，第2行为倾转角度为20
°
下的数据；计算所得原子坐标与真实坐标距离直方图如图9所示。
[0111]
根据本技术实施例提出的局域轨道函数的三维重构方法，基于采集的多个倾转角
度下的样品的图像数据，获得多个倾转角度下的计算图像，进而计算出计算图像的损失函数，求取优化参数的梯度，从而获得优化参数，在经过反复筛选和计算后，直至损失函数满足收敛条件，进而获得三维重构结果，可以精简三维坐标的重构过程，不仅降低了对硬件的要求，还能减少繁琐的人为干预，节约劳动力成本，同时在迭代过程中，还可以矫正样品漂移和样品台的机械倾转误差，进而提高三维坐标重构的准确性。由此，解决了相关技术中只能从重构出的三维密度矩阵中获取原子的三维坐标，且无法对误差进行校正，导致重构的过程对硬件需求较高，且重构的三维坐标精度较差的技术问题。
[0112]
其次参照附图描述根据本技术实施例提出的局域轨道函数的三维重构装置。
[0113]
图10是本技术实施例的局域轨道函数的三维重构装置的方框示意图。
[0114]
如图10所示，该局域轨道函数的三维重构装置10包括：采集模块100、累加模块200和重构模块300。
[0115]
具体地，采集模块100，用于采集多个倾转角度下的样品的图像数据。
[0116]
累加模块200，用于基于图像数据，在实空间等间隔进行撒点并利用线性累加，得到多个倾转角度的每个倾转角度下的计算图像。
[0117]
重构模块300，用于根据每个倾转角度下的计算图像计算损失函数，并获取损失函数关于待优化参数的梯度，且根据梯度优化待优化参数，筛选出满足预设条件的原子，重新计算新的损失函数，直至满足收敛条件，在实空间重构局域轨道函数的中心的三维空间坐标和局域轨道函数的形状，得到三维重构结果。
[0118]
可选地，在本技术的一个实施例中，采集模块100包括：获取单元和降噪单元。
[0119]
其中，获取单元，用于获取多个倾转角度下的样品的初始图像数据。
[0120]
降噪单元，用于对初始图像数据进行对中合轴和降噪处理，并由处理后的图像归一化，得到图像数据。
[0121]
可选地，在本技术的一个实施例中，损失函数的计算公式为：
[0122][0123]
其中，w代表损失函数，m代表倾转角度的总数，j代表倾转角度的序号，i代表局域轨道函数序号，p代表图像的边长，u代表图像中每个像素的横坐标，v代表图像中每个像素的纵坐标，fj(u,v)代表第j个角度下计算得到的图像，bj(u,v)代表第j个角度下实验所得的图像。
[0124]
可选地，在本技术的一个实施例中，重构模块300包括：检测单元和保护单元。
[0125]
其中，检测单元，用于在每一步迭代更新完参数后，检测所有局域轨道函数的参数删除单元，用于在检测到任一局域轨道函数的参数小于由所有局域轨道函数参数得到的阈值时，删除任一局域轨道函数的同时，通过建立二叉树得到局域轨道函数中心距离小于预设像素的局域轨道函数对索引，删除局域轨道函数对中的任意一个局域轨道函数。
[0126]
保护单元，用于在每一步迭代更新完参数并删除完局域轨道函数后，将每个局域轨道函数的参数均以预设概率减小为预设倍数，且设置保护时间，使得在保护时间内，不允许执行筛选操作和删除操作。
[0127]
可选地，在本技术的一个实施例中，待优化参数包括每个局域轨道函数中心的三维空间坐标、描述其形状的参数，每个转角对应的三个欧拉角，每个角度下样品的漂移、样
品台的机械倾转偏差中的至少一项。
[0128]
需要说明的是，前述对局域轨道函数的三维重构方法实施例的解释说明也适用于该实施例的局域轨道函数的三维重构装置，此处不再赘述。
[0129]
根据本技术实施例提出的局域轨道函数的三维重构装置，基于采集的多个倾转角度下的样品的图像数据，获得多个倾转角度下的计算图像，进而计算出计算图像的损失函数，求取优化参数的梯度，从而获得优化参数，在经过反复筛选和计算后，直至损失函数满足收敛条件，进而获得三维重构结果，可以精简三维坐标的重构过程，不仅降低了对硬件的要求，还能减少繁琐的人为干预，节约劳动力成本，同时在迭代过程中，还可以矫正样品漂移和样品台的机械倾转误差，进而提高三维坐标重构的准确性。由此，解决了相关技术中只能从重构出的三维密度矩阵中获取原子的三维坐标，且无法对误差进行校正，导致重构的过程对硬件需求较高，且重构的三维坐标精度较差的技术问题。
[0130]
图11为本技术实施例提供的电子设备的结构示意图。该电子设备可以包括：
[0131]
存储器1101、处理器1102及存储在存储器1101上并可在处理器1102上运行的计算机程序。
[0132]
处理器1102执行程序时实现上述实施例中提供的局域轨道函数的三维重构方法。
[0133]
进一步地，电子设备还包括：
[0134]
通信接口1103，用于存储器1101和处理器1102之间的通信。
[0135]
存储器1101，用于存放可在处理器1102上运行的计算机程序。
[0136]
存储器1101可能包含高速ram存储器，也可能还包括非易失性存储器(non-volatile memory)，例如至少一个磁盘存储器。
[0137]
如果存储器1101、处理器1102和通信接口1103独立实现，则通信接口1103、存储器1101和处理器1102可以通过总线相互连接并完成相互间的通信。总线可以是工业标准体系结构(industry standard architecture，简称为isa)总线、外部设备互连(peripheral component，简称为pci)总线或扩展工业标准体系结构(extended industry standard architecture，简称为eisa)总线等。总线可以分为地址总线、数据总线、控制总线等。为便于表示，图11中仅用一条粗线表示，但并不表示仅有一根总线或一种类型的总线。
[0138]
可选地，在具体实现上，如果存储器1101、处理器1102及通信接口1103，集成在一块芯片上实现，则存储器1101、处理器1102及通信接口1103可以通过内部接口完成相互间的通信。
[0139]
处理器1102可能是一个中央处理器(central processing unit，简称为cpu)，或者是特定集成电路(application specific integrated circuit，简称为asic)，或者是被配置成实施本技术实施例的一个或多个集成电路。
[0140]
本实施例还提供一种计算机可读存储介质，其上存储有计算机程序，该程序被处理器执行时实现如上的局域轨道函数的三维重构方法。
[0141]
在本说明书的描述中，参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本技术的至少一个实施例或示例中。在本说明书中，对上述术语的示意性表述不必须针对的是相同的实施例或示例。而且，描述的具体特征、结构、材料或者特点可以在任一个或n个实施例或示例中以合适的方式结合。此外，在不相互矛盾的情况下，本领域的技
术人员可以将本说明书中描述的不同实施例或示例以及不同实施例或示例的特征进行结合和组合。
[0142]
此外，术语“第一”、“第二”仅用于描述目的，而不能理解为指示或暗示相对重要性或者隐含指明所指示的技术特征的数量。由此，限定有“第一”、“第二”的特征可以明示或者隐含地包括至少一个该特征。在本技术的描述中，“n个”的含义是至少两个，例如两个，三个等，除非另有明确具体的限定。
[0143]
流程图中或在此以其他方式描述的任何过程或方法描述可以被理解为，表示包括一个或更n个用于实现定制逻辑功能或过程的步骤的可执行指令的代码的模块、片段或部分，并且本技术的优选实施方式的范围包括另外的实现，其中可以不按所示出或讨论的顺序，包括根据所涉及的功能按基本同时的方式或按相反的顺序，来执行功能，这应被本技术的实施例所属技术领域的技术人员所理解。
[0144]
在流程图中表示或在此以其他方式描述的逻辑和/或步骤，例如，可以被认为是用于实现逻辑功能的可执行指令的定序列表，可以具体实现在任何计算机可读介质中，以供指令执行系统、装置或设备(如基于计算机的系统、包括处理器的系统或其他可以从指令执行系统、装置或设备取指令并执行指令的系统)使用，或结合这些指令执行系统、装置或设备而使用。就本说明书而言，"计算机可读介质"可以是任何可以包含、存储、通信、传播或传输程序以供指令执行系统、装置或设备或结合这些指令执行系统、装置或设备而使用的装置。计算机可读介质的更具体的示例(非穷尽性列表)包括以下：具有一个或n个布线的电连接部(电子装置)，便携式计算机盘盒(磁装置)，随机存取存储器(ram)，只读存储器(rom)，可擦除可编辑只读存储器(eprom或闪速存储器)，光纤装置，以及便携式光盘只读存储器(cdrom)。另外，计算机可读介质甚至可以是可在其上打印所述程序的纸或其他合适的介质，因为可以例如通过对纸或其他介质进行光学扫描，接着进行编辑、解译或必要时以其他合适方式进行处理来以电子方式获得所述程序，然后将其存储在计算机存储器中。
[0145]
应当理解，本技术的各部分可以用硬件、软件、固件或它们的组合来实现。在上述实施方式中，n个步骤或方法可以用存储在存储器中且由合适的指令执行系统执行的软件或固件来实现。如，如果用硬件来实现和在另一实施方式中一样，可用本领域公知的下列技术中的任一项或他们的组合来实现：具有用于对数据信号实现逻辑功能的逻辑门电路的离散逻辑电路，具有合适的组合逻辑门电路的专用集成电路，可编程门阵列(pga)，现场可编程门阵列(fpga)等。
[0146]
本技术领域的普通技术人员可以理解实现上述实施例方法携带的全部或部分步骤是可以通过程序来指令相关的硬件完成，所述的程序可以存储于一种计算机可读存储介质中，该程序在执行时，包括方法实施例的步骤之一或其组合。
[0147]
此外，在本技术各个实施例中的各功能单元可以集成在一个处理模块中，也可以是各个单元单独物理存在，也可以两个或两个以上单元集成在一个模块中。上述集成的模块既可以采用硬件的形式实现，也可以采用软件功能模块的形式实现。所述集成的模块如果以软件功能模块的形式实现并作为独立的产品销售或使用时，也可以存储在一个计算机可读取存储介质中。
[0148]
上述提到的存储介质可以是只读存储器，磁盘或光盘等。尽管上面已经示出和描述了本技术的实施例，可以理解的是，上述实施例是示例性的，不能理解为对本技术的限
制，本领域的普通技术人员在本技术的范围内可以对上述实施例进行变化、修改、替换和变型。