球面波初相位系数的确定方法、全息图生成方法及设备

技术领域

本申请涉及计算全息技术领域，更具体地，涉及一种球面波初相位系数的确定方法、基于球面波的全息图生成方法、电子设备及可读存储介质。

背景技术

现阶段，大部分AR/VR显示设备中采用基于双目视差的“伪3D”显示技术，长时间佩戴这种显示设备会使用户产生诸如模糊、眩晕等视觉疲劳现象，严重时甚至可能对人的视觉系统造成不可逆转的伤害。为了解决上述问题，能够真实再现具有深度信息的3D场景的计算全息技术，成为了近年来的研究热点。

计算全息技术通过全息图生成算法，例如，Gerchberg—Saxton算法(GS算法)，将三维图像(RGBD图像)信息编码到仅包含相位信息的相位型全息图或仅包含振幅信息的振幅型全息图中。相比于振幅型全息图，相位型全息图凭借极高的衍射效率成为了现阶段最为常用的全息图类型。接着，将相位型全息图加载于空间光调制器(spatial lightmodulator，SLM)上，通过特定的光学系统即可在指定的空间中生成三维图像(即，再现像)。如何抑制再现像中的散斑与高频伪影，提高再现像观感是目前研究的热点。

发明内容

本申请提供了一种可至少部分解决相关技术中存在的上述问题的球面波初相位系数的确定方法、全息图生成方法及设备。

本申请一方面提供一种球面波初相位系数的确定方法，包括：设定球面波的初相位系数的循环计算规则；按照所设定的循环计算规则对球面波的初相位系数进行循环计算；其中，在每一次计算过程中包括：基于本次计算过程中的初相位系数形成样本图像的全息图；获取全息图的再现像以及再现像的频谱；确定频谱中的主峰位置；以及在循环计算终止后，将主峰位置与中心频谱之间的距离最大时对应的计算过程中的初相位系数作为球面波的初相位系数。

在一些实施方式中，基于本次计算过程中的初相位系数形成样本图像的全息图包括：基于本次计算过程中的初相位系数形成本次计算过程中球面波的初相位；为样本图像附加所形成的初相位构建样本图像的物平面光场分布；以及根据所构建出的物平面光场分布计算出样本图像的全息图。

在一些实施方式中，根据所构建出的物平面光场分布计算出样本图像的全息图包括：根据所构建出的物平面光场分布，采用GS算法计算出样本图像的相位型全息图。

在一些实施方式中，获取全息图的再现像包括：获取光学成像系统拍摄全息图形成的图像。

在一些实施方式中，确定频谱中的主峰位置之后，在每一次计算过程中还包括：确定频谱的主峰位置与中心频谱之间的距离；以及建立距离与本次计算过程中的初相位系数之间的对应关系。

在一些实施方式中，将主峰位置与中心频谱之间的距离最大时对应的计算过程中的初相位系数作为球面波的初相位系数包括：根据所建立的对应关系，确定出主峰位置与中心频谱的距离最大时对应的计算过程中的初相位系数；以及将所确定出的初相位系数作为球面波的初相位系数。

在一些实施方式中，样本图像的灰度值为255。

在一些实施方式中，设定球面波的初相位系数的循环计算规则包括：以等差数列的形式设定球面波的初相位系数在循环过程中的的计算公式；以及设定最大循环计算次数作为循环计算终止的判断条件。

一种基于球面波的全息图生成方法，包括：采用如如上所述的球面波初相位系数的确定方法确定出球面波的初相位系数；根据所确定出的初相位系数形成球面波的初相位；为输入图像附加球面波的初相位构建输入图像的物平面光场分布；以及根据所构建出的物平面光场分布计算出输入图像的全息图。

一种基于球面波的全息图生成系统，包括：

初相位系数确定模块，用于采用如上所述的球面波初相位系数的确定方法确定出球面波的初相位系数；初相位形成模块，用于根据所确定出的初相位系数形成所述球面波的初相位；物平面光场分布构建模块，用于为输入图像附加所述球面波的初相位构建所述输入图像的物平面光场分布；以及全息图生成模块，用于根据所构建出的物平面光场分布生成所述输入图像的全息图。

一种电子设备，包括：处理器；以及存储器，与处理器通讯连接，其中，存储器存储有可被处理器执行的程序，当程序被处理器执行时，处理器能够执行根据如上所述的球面波初相位系数的确定方法或者根据如上所述的基于球面波的全息图生成方法。

一种可读存储介质，可读存储介质上存储有计算机程序，计算机程序被处理器执行时实现根据如上所述的球面波初相位系数的确定方法或者根据如上所述的基于球面波的全息图生成方法。

本申请至少一个实施方式提供的球面波初相位系数的确定方法，通过将基于样本图像的全息图的再现像频谱中的主峰位置与中心频谱的距离最大时对应的计算过程中的初相位系数作为球面波的初相位系数，使得所确定的球面波初相位系数为考虑了全息图再现像的散斑尺寸与周期特性后的数值，从而使得为输入图像附加采用该初相位系数的球面波初相位计算全息图时，可抑制输入图像的全息图的再现像中的散斑与高频伪影噪声现象，提高再现像的观感。

附图说明

通过阅读参照以下附图所作的对非限制性实施例的详细描述，本申请的其它特征、目的和优点将会变得更明显。其中：

图1A是根据相关技术为输入图像附加随机初相位时，所得全息图的再现像的示意图；

图1B是图1A中虚线框内的局部示意图；

图2A是根据相关技术为输入图像附加传统球面波初相位时，所得全息图的再现像的示意图；

图2B是图2A中虚线框内的局部示意图；

图3A是根据本申请示例性实施方式提供的球面波初相位系数的确定方法的流程示意图；

图3B是根据本申请示例性实施方式提供的基于球面波的全息图生成方法的流程示意图；

图4A是根据本申请实施方式的步骤S120的流程示意图

图4B是根据本申请实施方式的步骤S121的流程示意图；

图5A是根据本申请实施方式的样板图像的示意图；

图5B是根据图5A的样板图像的全息图的再现像示意图；

图5C是根据图5B中再现像的频谱示意图；

图5D是根据图5C中再现像中心水平方向上的单边谱示意图；

图6是根据本申请实施方式的步骤S130的流程示意图；

图7A和图7B是在不同的主峰位置与中心频谱的距离时选取的初相位系数形成的全息图的再现像示意图；

图8A和图8B分别是采用本申请提出的基于球面波的全息图生成方法生成的全息图对应的再现像及其局部放大图；

图9A和图9B分别是基于传统球面波初相位的全息图生成算法生成的全息图对应的再现像及其局部放大图；

图10A和图10B分别是基于随机初相位的全息图生成算法生成的全息图对应的再现像及其局部放大图；

图11是本申请的示例性实施方式的基于球面波的全息图生成系统示意图；

图12是本申请的示例性实施方式的电子设备的结构示意图；以及

图13是本申请的示例性实施方式的电子设备的结构示意图。

具体实施方式

为了更好地理解本申请，将参考附图对本申请的各个方面做出更详细的说明。应理解，这些详细说明只是对本申请的示例性实施方式的描述，而非以任何方式限制本申请的范围。在说明书全文中，相同的附图标号指代相同的元件。表述“和/或”包括相关联的所列项目中的一个或多个的任何和全部组合。

应注意，在本说明书中，第一、第二、第三等的表述仅用于将一个特征与另一个特征区域分开来，而不表示对特征的任何限制，尤其不表示任何的先后顺序。因此，在不背离本申请的教导的情况下，本申请中讨论的第一部分也可被称作第二部分，反之亦然。

在附图中，为了便于说明，已稍微调整了部件的厚度、尺寸和形状。附图仅为示例而并非严格按比例绘制。如在本文中使用的，用语“大致”、“大约”以及类似的用语用作表近似的用语，而不用作表程度的用语，并且旨在说明将由本领域普通技术人员认识到的、测量值或计算值中的固有偏差。

还应理解的是，诸如“包括”、“包括有”、“具有”、“包含”和/或“包含有”等表述在本说明书中是开放性而非封闭性的表述，其表示存在所陈述的特征、元件和/或部件，但不排除一个或多个其它特征、元件、部件和/或它们的组合的存在。此外，当诸如“...中的至少一个”的表述出现在所列特征的列表之后时，其修饰整列特征，而非仅仅修饰列表中的单独元件。此外，当描述本申请的实施方式时，使用“可”表示“本申请的一个或多个实施方式”。并且，用语“示例性的”旨在指代示例或举例说明。

除非另外限定，否则本文中使用的所有措辞(包括工程术语和科技术语)均具有与本申请所属领域普通技术人员的通常理解相同的含义。还应理解的是，除非本申请中有明确的说明，否则在常用词典中定义的词语应被解释为具有与它们在相关技术的上下文中的含义一致的含义，而不应以理想化或过于形式化的意义解释。

需要说明的是，在不冲突的情况下，本申请中的实施例及实施例中的特征可以相互组合。下面将参考附图并结合实施例来详细说明本申请。

如前所述，相比于振幅型全息图，相位型全息图凭借极高的衍射效率成为了现阶段最为常用的全息图类型。

在现有的全息图生成算法中，为输入图像附加初相位(initial phase)能够有效地将输入图像中的低频分量编码到全息图中，进而使再现像尽可能地接近输入图像。

其中，随机初相位是最为常见的方法，以输入图像为二维图像为例，其具体表达形式如下：

式中，

然而，受限于激光光源的高相干特性，随机初相位会在再现像中引入无规则排序的颗粒状散斑，如图1A及其局部放大图1B所示。经研究，采用球面波初相位可减小初相位中相邻像素点间的差值，从而可实现弱化上述散斑的目的。

同样，以二维图像为例，球面波初相位的表达形式如下：

式中，

可以看出，上述初相位定义方式能够使相位沿中心向四周有序变化，从而弱化散斑。

球面波的初相位系数a决定了全息图再现像的纹理特征，然而传统的球面波初相位定义方式却忽略了这一点，仅简单地将系数a的取值方式定义为π/M，M为全息图的像素尺寸。

当为输入图像附加传统的球面波初相位时，所得全息图的再现像及其局部放大图分别如图2A和图2B所示。可以看出，传统的球面波初相位虽然弱化了全息图再现像中的散斑，却出现了众多的高频伪影，同样影响观感。综上所述，随机初相位系数会在全息图再现像中引入散斑噪声，而传统的球面波初相位系数则会在全息图再现像中引入高频伪影。

基于此，本申请的示例性实施方式提供一种球面波初相位系数的确定方法，包括：设定球面波的初相位系数的循环计算规则；按照所设定的循环计算规则对球面波的初相位系数进行循环计算；其中，在每一次计算过程中包括：基于本次计算过程中的初相位系数形成样本图像的全息图；获取全息图的再现像以及再现像的频谱；确定频谱中的主峰位置；以及在循环计算终止后，将主峰位置与中心频谱的距离最大时对应的计算过程中的初相位系数作为球面波的初相位系数。本申请的上述方案中，在循环计算球面波的初相位系数过程中，通过将基于样本图像的全息图的再现像频谱中的主峰位置与中心频谱的距离最大时对应的计算过程中的初相位系数作为球面波的初相位系数，使得所确定的球面波初相位系数为考虑了全息图再现像的散斑尺寸与周期特性后的数值，从而使得为输入图像附加采用该初相位系数的球面波初相位计算全息图时，可抑制输入图像的全息图的再现像中的散斑与高频伪影噪声现象，提高再现像的观感。

在下文中，将结合附图更详细地描述本方案的具体示例。

如图3A所示，本申请第一实施方式提供了一种球面波初相位系数的确定方法1000。如图3A所示，在球面波初相位系数的确定方法1000中包括如下步骤：

S110、设定球面波的初相位系数的循环计算规则；

S120，按照所设定的循环计算规则对球面波的初相位系数进行循环计算；其中，在每一次计算过程中包括基于本次计算过程中的初相位系数形成样本图像的全息图；获取全息图的再现像以及再现像的频谱；确定频谱中的主峰位置；以及

S130，在循环计算终止后，将主峰位置与中心频谱的距离最大时对应的计算过程中的初相位系数作为球面波的初相位系数。

应理解的是，在球面波初相位系数的确定方法1000中所示的步骤不是排它性的，还可以在所示步骤中的任何步骤之前、之后或之间执行其它步骤。此外，上述步骤中的一些步骤可以是同时地执行的或者可以是按照不同于图3A所示的顺序执行的。下面结合图3A、图4A至图7B进一步描述上述的步骤S110至S130。

S110

如前式(2)中所示出的球面波初相位

式中，a为球面波的初相位系数，x和y为相应的像素坐标。初相位系数a的大小直接决定了全息图再现像的散斑特征。

可选的，在步骤S110中以等差数列的形式设定球面波的初相位系数的循环计算公式如下：

a

式中，a

在第1次计算过程中，a

在第2次计算过程中，a

在第3次计算过程中，a

依次类推，在第n次计算过程中，a

可选的，预设的初始值a

步骤S110中还确定出初相位系数的循环计算的终止条件。

在一些实施方式中，循环计算的终止条件为n＞N，其中N为采样个数。N可以为大于1的任意正整数，N越大，搜索范围越大，结果越准确。示例性的，N为10至10000之间的任一数值。可选的，N为50至5000之间的任一数值。

因此，本申请的球面波初相位系数a的取值范围为a∈[a

此处需要说明的是，根据预设的初始值以及终止条件的不同，初相位系数a的选取范围可以从0至无穷大。

此外，虽然本申请以等差数列的形式介绍相关技术方案，但是对于以等比数列或者变步长的形式确定的一组球面波的初相位系数同样适用于本申请的技术方案，对此本申请不做限制。

S120

在步骤S110中确定球面波的初相位系数的循环计算公式以及终止条件之后，在步骤S120中对球面波的初相位系数进行循环计算。

如图4A所示，步骤S120中的每一次计算过程均依次包括以下步骤：

S121、基于本次计算过程中的初相位系数形成样本图像的全息图；

S122、获取全息图的再现像以及再现像的频谱；以及

S123、确定频谱中的主峰位置。

通过上述步骤确定出基于样本图像的全息图的再现像频谱中的主峰位置，以便在后续步骤中将主峰位置与中心频谱的距离最大时对应的计算过程中的初相位系数作为球面波的初相位系数。

在一些实施方式中，如图4B所示，步骤S121中包括：

S1211、基于本次计算过程中的初相位系数形成本次计算过程中球面波的初相位；

S1212、为样本图像附加所形成的初相位构建样本图像的物平面光场分布；以及

S1213、根据所构建出的物平面光场分布计算出样本图像的全息图。

具体而言，计算开始后，当n＝1时为第1次计算过程，在第1次计算过程中的初相位系数a

在步骤S1211中基于本次计算过程中的初相位系数a

可选的，步骤S1213中根据步骤S1212中所构建出的物平面光场分布，采用GS算法计算出样本图像的相位型全息图

此处需要说明的是，根据物平面光场分布，采用GS算法计算出样本图像的相位型全息图

在步骤S121中形成样本图像的全息图后，在步骤S122中获取全息图的再现像以及再现像的频谱。

在一些实施方式中，步骤S122中获取全息图的再现像包括接收光学成像系统拍摄全息图形成的图像。

具体而言，将步骤S1213中形成的全息图

本申请的上述方案引入了相机实拍再现像的策略，使得全息图的再现像能够包含因光学系统像差等因素引入的噪声成分，保障了确定球面波初相位系数的准确性。这是由于在相关技术中，一般采用仿真的方式获得全息图的再现像，而再现像的仿真结果无法准确反应其散斑特征，同时，仿真时采用的光传播模型与实际光学系统有一定的差距，因而本申请采用相机实拍再现像的方法解决上述问题。

继续在步骤S122中获取再现像的频谱，以样本图像为二维图像为例，样本图像全息图的再现像的二维频谱参照图5C所示。需要说明的是，虽然本申请此处以二维图像为例进行介绍，但是本申请涉及的所有方法对三维图像是同样适用的。

在步骤S123中确定上述频谱的主峰位置。可选的，考虑到频谱的偶函数特性，仅需对中心水平方向上的单边谱进行分析，如图5C中虚线框内所示。需要说明的是，图5C中示出的是对频谱取log计算后的结果。图5D为全息再现像中心水平方向上的单边谱。如图5D所示，频谱中的主峰如图5D中“*”号位置所示，本申请此处的主峰的定义方式为：1)频谱中的极大值点，即该点的强度大于左右点的强度；2)所有极大值点中强度最大的点。

在一些实施方式中，继续参照图4A所示，步骤S123中确定频谱中的主峰位置之后，在每一次计算过程中还包括：

S124、确定频谱的主峰位置与中心频谱之间的距离；以及

S125、建立上述距离与本次计算过程中的初相位系数之间的对应关系。

具体而言，在第1次计算过程中，确定主峰后，即可在步骤S124中计算出主峰相对于中心频谱的距离dis

至此，第1次计算过程结束。

然后以循环计算的终止条件判断循环计算是否结束，如果n<＝N，则n＝n+1，并返回步骤S120中继续进行下一次循环计算。如果n>N，则结束循环计算。

示例性的，当循环计算结束后，可获得多组(系数a，距离dis)的数据对。例如：

第1次循环计算所获得的结果：(系数a

第2次循环计算所获得的结果：(系数a

第3次循环计算所获得的结果：(系数a

第n次循环计算所获得的结果：(系数a

……

第N次循环计算所获得的结果：(系数a

S130

步骤S130中，将主峰位置与中心频谱之间的距离最大时对应的计算过程中的初相位系数作为球面波的初相位系数。

在一些实施方式中，如图6所示，步骤S130包括：

S131、根据所建立的对应关系，确定出主峰位置与中心频谱的距离最大时对应的循环计算过程中的初相位系数；以及

S132、将所确定出的初相位系数作为球面波的初相位系数。

频谱上的主峰位置反应了全息图再现像上散斑的尺寸与周期特性，主峰位置距离中心频谱越远，散斑尺寸越小且周期越小，因而使全息图的再现像看起来更为细腻。

根据大量的实验结果，对于像素尺寸为1080*1080的全息图来说，初相位系数a的选取范围大致在150至200之间较为合适。此外，基于本申请上述方法确定的球面波的初相位系数a通常大于100，而基于现有相关技术确定的球面波初相位系数通常小于1，因此，也可将本申请中的球面波初相位称为大系数球面波初相位。

图7A和图7B分别示出了在距离dis＝150lp/mm和dis＝50lp/mm时，样板图像的全息图对应的再现像。不难看出，图7A中散斑的周期与尺寸要明显小于图7B中散斑的周期与尺寸，证明了本申请上述球面波系数选择方法的可行性。

本申请的上述方案具有如下特点：

1)确定球面波初相位系数时，将纯白图像作为输入，能够使再现像清楚地呈现散斑特征而无其余干扰因素，简化了再现像频谱分析的难度；

2)通过分析再现像的频谱特征，来判断散斑尺寸与周期特性，进而指导球面波初相位系数的选择，从而抑制了全息图的再现像中的散斑与高频伪影噪声；

3)通过相机实拍的方式能够记录光学系统像差等因素对再现像质量的影响，使得在选择球面波初相位系数时能够充分考虑实际光学系统的误差，保障了选取球面波初相位系数的准确性。

本申请示例性第二实施方式还提供一种基于球面波的全息图生成方法2000，如图3B所示，包括以下步骤：

S210、采用球面波初相位系数的确定方法确定出球面波的初相位系数，并根据所确定出的初相位系数形成球面波的初相位；

S220、为输入图像附加球面波的初相位构建输入图像的物平面光场分布；以及

S230、根据所构建出的物平面光场分布计算出输入图像的全息图。

具体而言，根据上述实施方式提供的如图3A示出的球面波初相位系数的确定方法1000所确定的球面波初相位系数a，构建大系数球面波初相位

式中的u

然后，采用GS算法根据物平面光场分布u计算出输入图像对应的全息图

图8A和图8B分别示出了采用本申请提出的基于球面波的全息图生成方法生成的全息图对应的再现像及其局部放大图。

为对比说明本申请的与相关技术的区别，图9A和图9B以及图10A和图10B则分别为基于传统球面波初相位的全息图生成算法以及基于随机初相位的全息图生成算法生成的全息图对应的再现像及其局部放大图。

从图8A和图8B、图9A和图9B以及图10A和图10B中可以看出，本申请生成的全息图对应的再现像整体观感明显优于传统算法，其散斑呈均匀周期性分布，散斑尺寸也更小。

本申请示例性第三实施方式还提供一种基于球面波的全息图生成系统3000，如图11所示，包括：

初相位系数确定模块310，用于采用如上所述的球面波初相位系数的确定方法1000确定出球面波的初相位系数；

初相位形成模块320，用于根据所确定出的初相位系数形成所述球面波的初相位；

物平面光场分布构建模块330，用于为输入图像附加所述球面波的初相位构建所述输入图像的物平面光场分布；以及

全息图生成模块340，用于根据所构建出的物平面光场分布生成所述输入图像的全息图。

本申请提供的基于球面波的全息图生成系统3000用于实施如图3B所示出的基于球面波的全息图生成方法2000，因此，基于球面波的全息图生成系统3000的具体实施过程可参照上述方案，本申请此处不再赘述。

此外，基于上述球面波初相位系数的确定方法或者基于球面波的全息图生成方法，本申请的实施方式还提供一种电子设备，例如，服务器、云端服务器等。

图12示出了本申请的第一示例性实施方式的电子设备的结构示意图。

如图12所示，电子设备包括：至少一个处理器701；以及与至少一个处理器701通信连接的存储器702；其中，存储器存储有可被至少一个处理器701执行的指令，指令被至少一个处理器701执行，以使至少一个处理器701能够执行上述实施方式提及的球面波初相位系数的确定方法或者基于球面波的全息图生成方法。其中，电子设备旨在表示各种形式的数字计算机，诸如，膝上型计算机、台式计算机、工作台、个人数字助理、服务器、刀片式服务器、大型计算机、和其它适合的计算机。电子设备还可以表示各种形式的移动装置，诸如，个人数字处理、蜂窝电话、智能电话和其它类似的计算装置。本文所示的部件、它们的连接和关系、以及它们的功能仅仅作为示例，并且不意在限制本文中描述的和/或者要求的本申请的实现。

图13示出了本申请的第二示例性实施方式的电子设备的结构示意图。

如图13所示，电子设备还可例如包括：I/O接口703、输入单元704、输出单元705、通信单元706、只读存储器(ROM)707和随机存取存储器(RAM)708。具体地，处理器701可根据存储在ROM707中的计算机程序或者从存储器702加载到RAM708中的计算机程序，来执行各种适当的动作和处理。在RAM708中，还可存储电子设备操作所需的各种程序和数据。处理器701、ROM707以及RAM708通过总线709彼此相连。I/O接口(输入/输出接口)703也连接至总线709。

电子设备中的多个部件连接至I/O接口703，包括：输入单元704，例如键盘、鼠标等；输出单元705，例如各种类型的显示器、扬声器等；存储器702，例如磁盘、光盘等；以及通信单元706，例如网卡、调制解调器、无线通信收发机等。通信单元706允许电子设备通过诸如因特网的计算机网络和/或各种电信网络与其他设备交换信息/数据。

处理器701可以是各种具有处理和计算能力的通用和/或专用处理组件。处理器701的一些示例包括但不限于中央处理单元(CPU)、图形处理单元(GPU)、各种专用的人工智能(AI)计算芯片、各种运行机器学习模型算法的处理器、数字信号处理器(DSP)、以及任何适当的处理器、控制器、微控制器等。处理器701执行上文所描述的各个方法和处理，例如球面波初相位系数的确定方法或者基于球面波的全息图生成方法。例如，在一些实施例中，球面波初相位系数的确定方法或者基于球面波的全息图生成方法可被实现为计算机软件程序，其被有形地包含于计算机可读存储介质，例如存储器702。在一些实施例中，计算机程序的部分或者全部可以经由ROM707和/或通信单元706而被载入和/或安装到电子设备上。当计算机程序加载到RAM708并由处理器701执行时，可以执行上文描述的球面波初相位系数的确定方法或者基于球面波的全息图生成方法的一个或多个步骤。备选地，在其他实施例中，处理器701可以通过其他任何适当的方式(例如，借助于固件)而被配置为执行球面波初相位系数的确定方法或者基于球面波的全息图生成方法。

本文中以上描述的系统和技术的各种实施方式可以在数字电子电路系统、集成电路系统、场可编程门阵列(FPGA)、专用集成电路(ASIC)、专用标准产品(ASSP)、芯片上系统的系统(SOC)、复杂可编程逻辑设备(CPLD)、计算机硬件、固件、软件、和/或它们的组合中实现。这些各种实施方式可以包括：实施在一个或者多个计算机程序中，该一个或者多个计算机程序可在包括至少一个可编程处理器的可编程系统上执行和/或解释，该可编程处理器可以是专用或者通用可编程处理器，可以从存储系统、至少一个输入装置、和至少一个输出装置接收数据和指令，并且将数据和指令传输至该存储系统、该至少一个输入装置、和该至少一个输出装置。

于实施本申请的方法的程序代码可以采用一个或多个编程语言的任何组合来编写。上述程序代码可以封装成计算机程序产品。这些程序代码或计算机程序产品可以提供给通用计算机、专用计算机或其他可编程数据处理装置的处理器或控制器，使得程序代码当由处理器701执行时使流程图和/或框图中所规定的功能/操作被实施。程序代码可以完全在计算机上执行、部分地在计算机上执行，作为独立软件包部分地在计算机上执行且部分地在远程计算机上执行或完全在远程计算机或服务器上执行。

电子设备的具体描述及有益效果，可以参考上述球面波初相位系数的确定方法或者基于球面波的全息图生成方法的描述，不再赘述。

此外，这里需要指出的是：本申请另一方面还提供了一种计算机可读存储介质，且所述计算机可读存储介质中存储有前文提及的球面波初相位系数的确定方法或者基于球面波的全息图生成方法所执行的计算机程序，且所述计算机程序包括程序指令，当所述处理器执行所述程序指令时，能够执行前文所述球面波初相位系数的确定方法或者基于球面波的全息图生成方法的描述，因此，这里将不再进行赘述。另外，对采用相同方法的有益效果描述，也不再进行赘述。对于本申请所涉及的计算机可读存储介质实施例中未披露的技术细节，请参照本申请方法实施例的描述。

本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程，是可以通过计算机程序来指令相关的硬件来完成，程序可存储于一计算机可读存储介质中，该程序在执行时，可包括如上述各方法的实施例的流程。其中，存储介质可为磁碟、光盘、只读存储记忆体(Read-Only Memory，ROM)或随机存储记忆体(Random Access Memory，RAM)等。

以上描述仅为本申请的实施方式以及对所运用技术原理的说明。本领域技术人员应当理解，本申请中所涉及的保护范围，并不限于上述技术特征的特定组合而成的技术方案，同时也应涵盖在不脱离技术构思的情况下，由上述技术特征或其等同特征进行任意组合而形成的其它技术方案。例如上述特征与本申请中公开的(但不限于)具有类似功能的技术特征进行互相替换而形成的技术方案。