基于FS-ORAP和相位补偿的非迭代三维全息图生成方法及装置

基于fs-orap和相位补偿的非迭代三维全息图生成方法及装置
技术领域
1.本发明涉及图像技术领域，具体涉及一种基于fs-orap和相位补偿的非迭代三维全息图生成方法及装置。


背景技术：

2.三维(3d)全息图可以重建三维场景的整个光波场；因此，能够提供人眼可以感知的所有深度线索。随着空间光调制器(slm)和计算技术的发展，可以利用计算机生成全息图(cgh)对三维场景进行动态重建。cgh是一个广泛的研究领域，它包含将再现光场所需要的信息编码为全息图的各种计算方法。
3.3d cgh常用迭代方法生成纯相位全息图，三维顺序gs算法和三维全局gs算法是两种基于gerchberg-saxton(g-s)算法的3d cgh方法，前者倾向于使用最后一个平面作为约束，因此无法选择重建平面。此外，上述两种方法都是针对具体目标进行迭代的，当需要生成不同对象的相位全息图时，需要重新迭代生成，非常耗时。2018年，alejandro velez zea等人提出了一种优化随机相位(optimized random phase，orap)方法，用于非迭代生成纯相位全息图。orap具有与目标窗口相同的支撑，但目标支撑的大小和位置不能改变。否则，需要生成新的orap。
4.针对传统gs方法耗时和orap方法灵活性差的问题，本文提出基于单个满支撑优化随机相位和相位补偿的非迭代三维计算机生成全息图。


技术实现要素：

5.为了解决现有技术中的不足，本发明提供了一种基于fs-orap和相位补偿的非迭代三维全息图生成方法及装置。
6.在本发明的第一方面，提供了一种基于fs-orap和相位补偿的非迭代三维全息图生成方法，该方法包括：
7.对待处理的三维目标沿深度分层，生成多个二维图像；
8.根据所述多个二维图像的振幅与事先生成的满支撑优化随机相位，生成各层纯相位全息图；
9.对所述各层纯相位全息图分别添加补偿相位，并将添加补偿相位后得到的各层全息图的复振幅叠加，生成复合相位全息图；
10.对所述复合相位全息图进行菲涅尔变换重建所述三维目标。
11.如上所述的方面和任一可能的实现方式，进一步提供一种实现方式，所述满支撑优化随机相位的生成，包括：
12.预先在重建平面生成和slm平面相同大小的满支撑单位振幅；
13.根据所述满支撑单位振幅和初始化随机相位生成满支撑优化随机相位。
14.如上所述的方面和任一可能的实现方式，进一步提供一种实现方式，所述根据所
述满支撑单位振幅和初始化随机相位生成满支撑优化随机相位，包括：
15.将所述满支撑单位振幅与所述初始化随机相位相乘，得到所述重建平面上的复振幅；
16.将所述重建平面上的复振幅进行傅里叶逆变换，传播到slm平面；
17.在所述slm平面和所述重建平面之间执行若干次标准的g-s循环，得到所述满支撑优化随机相位。
18.如上所述的方面和任一可能的实现方式，进一步提供一种实现方式，根据所述多个二维图像的振幅与事先生成的满支撑优化随机相位，生成各层纯相位全息图，包括：
19.将所述多个二维图像的振幅分别与所述满支撑优化随机相位相乘，并将得到的乘积进行傅里叶逆变换，生成每层的纯相位全息图。
20.如上所述的方面和任一可能的实现方式，进一步提供一种实现方式，所述对所述各层纯相位全息图分别添加补偿相位，包括：
21.将各层纯相位全息图相位分别与不同焦距的菲涅耳相位掩膜叠加，得到新的相位分布。
22.如上所述的方面和任一可能的实现方式，进一步提供一种实现方式，所述对所述复合相位全息图进行菲涅尔变换重建所述三维目标物体，包括：
23.对所述复合相位全息图添加恢复相位，并根据添加恢复相位后的复合相位全息图重建所述三维目标。
24.在本发明的第二方面，提供了一种基于fs-orap和相位补偿的非迭代三维全息图生成装置，该装置包括：
25.划分模块，用于对待处理的三维目标沿深度分层，生成多个二维图像；
26.计算模块，用于根据所述多个二维图像的振幅与事先生成的满支撑优化随机相位，生成各层纯相位全息图；
27.补偿相位模块，用于对所述各层纯相位全息图分别添加补偿相位，并将添加补偿相位后得到的各层全息图的复振幅叠加，生成复合相位全息图；
28.重建模块，用于对所述复合相位全息图进行菲涅尔变换重建所述三维目标。
29.和现有技术相比，本发明的优点为：
30.(1)本发明实现了使用单幅fs-orap非迭代快速生成三维目标的全息图，利用fs-orap，可以生成具有任意位置和大小的目标的幅值相位全息图，克服了原有orap方法在目标位置和大小上的限制，无需迭代，极大的加快了三维目标相位全息图的生成。
31.(2)本发明采用菲涅尔相位掩膜进行相位补偿，可以自由选择重建平面来重建三维目标物体。
32.应当理解，发明内容部分中所描述的内容并非旨在限定本发明的关键或重要特征，亦非用于限制本发明的范围。本发明的其它特征将通过以下的描述变得容易理解。
附图说明
33.结合附图并参考以下详细说明，本发明各实施例的上述和其他特征、优点及方面将变得更加明显。附图用于更好地理解本方案，不构成对本发明的限定在附图中，相同或相似的附图标记表示相同或相似的元素，其中：
34.图1示出了根据本发明的基于fs-orap和相位补偿的非迭代三维全息图生成方法的流程图；
35.图2示出了根据本发明的三维全息图基于层的方法生成示意图；
36.图3示出了根据本发明的fs-orap方法生成纯相位全息图流程图；
37.图4示出了根据本发明的各层全息图生成模块示意图；
38.图5示出了根据本发明的补偿相位模块示意图；
39.图6示出了根据本发明的全息投影系统原理示意图；
40.图7示出了根据本发明的三维目标全息图重建过程示意图；
41.图8示出了根据本发明的sfs-orap-pc-3d整体流程图；
42.图9示出了根据本发明的基于fs-orap和相位补偿的非迭代三维全息图生成装置的框图；
43.图10示出了用来实现本发明实施例的基于fs-orap和相位补偿的非迭代三维全息图生成方法的电子设备的框图。
具体实施方式
44.为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的全部其他实施例，都属于本发明保护的范围。
45.另外，本文中术语“和/或”，仅仅是一种描述关联对象的关联关系，表示可以存在三种关系，例如，a和/或b，可以表示：单独存在a，同时存在a和b，单独存在b这三种情况。另外，本文中字符“/”，一般表示前后关联对象是一种“或”的关系。
46.本发明中，为了改善传统gs方法耗时和orap方法灵活性差的问题，提出了一种基于fs-orap和相位补偿的非迭代三维全息图生成方法，直接应用于三维全息显示全息图的生成，突破了原始orap方法在空间域中目标振幅固定支撑的大小和位置的限制。该方法使用单个fs-orap生成各层的相位全息图，无需迭代，极大的加快了三维目标相位全息图的生成。
47.图1示出了根据基于fs-orap和相位补偿的非迭代三维全息图生成方法100的流程图，如图1所示：
48.s101，对待处理的三维目标物体沿深度分层，生成多个二维图像。
49.在一些实施例中，全息投影显示过程可分为波前记录和波前重建两个步骤。波前记录是物体波前和参考波前的干涉记录，即利用光的干涉原理，将物体发射的特定光波以干涉条纹的形式记录下来，使物光波前的全部信息全部记录在介质中，所记录的干涉条纹图像称为全息图；当全息图被再现的光波照射时，由于衍射原理能再现原始物光波，从而形成与原物体逼真的三维图像，全息投影系统原理如图2。
50.具体地，如图2所示，通过激光器laser发射一束激光，经beam expander扩散后，通过lens将扩散后的激光转换为平面波，通过平面波照射slm平面的纯相位全息图，在image plane上观察重建图像。
51.在一些实施例中，可以采用基于层的方法生成三维目标的全息图，首先需要根据
深度信息，将三维物体分割成多个二维层，再生成各个层的相位全息图，最后通过复用将全息图进行组合，图3为三维全息图基于层的方法生成示意图，如图3所示：
52.将三维目标分层，得到多个二维层，各个二维层到slm平面的距离分别为z1、z2……zn
，在slm平面生成各个层的相位全息图，最后通过复用将全息图进行组合重建。
53.s102，根据所述多个二维图像的振幅与事先生成的满支撑随机相位，生成各层纯相位全息图。
54.在一些实施例中，alejandro velez zea等人提出的优化随机相位(orap)方法只能生成特定大小和指定位置目标的全息图，当生成不同尺寸和位置目标的全息图，则需要重新迭代生成orap，因此，当该方法用于三维目标的全息图生成时，需要针对每层振幅的大小和尺寸迭代生成orap，生成各层的相位全息图，这增加了三维目标的全息图的生成时间，并且灵活性低。
55.为解决orap仅支持固定大小位置目标全息图的生成，本发明使用满支撑优化随机相位方法(fs-orap)。除具有orap的优点外，fs-orap在需要生成不同尺寸或位置目标的全息图时只需生成一次fs-orap，使用该fs-orap即可进行目标重建。fs-orap是满支撑的，支持生成任意大小和位置的目标的全息图，极大的节省了生成多幅不同大小和位置目标的全息图的时间。
56.在一些实施例中，fs-orap的生成，可以首先生成一个与slm平面大小相同的满支撑单位振幅窗口作为重构平面，使用初始化随机相位乘以此重构平面，对乘积进行逆傅里叶变换(ift)，然后进行标准g-s算法循环的多次迭代，将得到的复振幅的振幅替换为目标窗口振幅。
57.在slm平面中，使用与slm区域大小相同的均匀振幅，在重建平面中，使用先前创建的目标窗口振幅，迭代保证重建的窗口振幅和目标窗口振幅之间有足够的近似，迭代完成后，生成fs-orap，并将分层后的二维图像振幅与fs-orap相乘，对得到的乘积做傅里叶逆变换，生成每层的纯相位全息图，通过fs-orap方法生成纯相位全息图流程如图4。
58.具体地，本发明中各层全息图生成模块如图5所示，首先将带有“某某大学”字样的三维图像分割成分别带有“某”、“某”、“大”、“学”字样的二维图像layer1-layer4，并将上述二维图像与事先生成的fs-orap相乘，并对乘积做傅里叶逆变换，得到各层二维图像的振幅全息图ah1-ah4以及各层二维图像的纯相位信息图ph1-ph4。
59.s103，对所述各层纯相位全息图分别添加补偿相位，并将添加补偿相位后得到的各层全息图的复振幅叠加，生成复合相位全息图。
60.在一些实施例中，由于在步骤s102中得到的纯相位全息图无法体现各层二维图像的深度信息，无法描述各层与slm平面的距离也无法选择重建平面，因此在重建时不能还原三维目标，因此，对各层的相位全息图添加补偿相位，以增添深度信息。
61.本发明将计算得到的全息图相位与菲涅耳相位掩膜叠加后加载到slm平面上，不同焦距fi(i＝1，2，3，4)的菲涅耳相位掩膜，分别与相位全息图进行叠加，得到新的相位分布为：
[0062][0063]
其中λ为波长，u,v为坐标。
[0064]
进一步地，将添加补偿相位后的各层全息图的复振幅通过叠加组合在一起，生成三维目标的全息图，并将该过程模块化，相位补偿模块具体如图6所示，对各层的纯相位全息图ph1-ph4分别与菲涅耳相位掩膜fpm1-fpm4叠加，得到各复合相位全息图cph1-cph4。
[0065]
在一些实施例中，补偿相位为正焦距菲涅尔相位掩膜。
[0066]
s104，对所述复合相位全息图进行菲涅尔变换重建所述三维目标物体。
[0067]
在一些实施例中，将复合相位全息图转变为复振幅之后进行叠加，然后进行重建。对全息图复振幅做菲涅尔变换即可重建原三维目标，重建过程具体如图7。
[0068]
具体地，如图7所示，将复合相位全息图复振幅ca_cph做不同距离的菲涅尔变换，重建三维目标。
[0069]
根据本公开的另一实施例对本公开三维全息图生成流程进行详细说明，三维全息图流程如图8所示：
[0070]
首先将带有“某某大学”字样的三维目标分成带有“某”、“某”、“大”、“学”的二维平面，并将所述二维平面的振幅与事先生成的fs-opar相乘，并对乘积做傅里叶逆变换，得到各层纯相位全息图ph1-ph4，并对各层纯相位全息图ph1-ph4添加补偿相位，得到复合相位全息图ca-cph，根据得到的复合相位全息图，进行菲涅尔变换重建三维目标。
[0071]
根据本发明，实现了以下技术效果：
[0072]
结合满支撑优化随机相位(fs-orap)方法和三维分层思想，首先基于单幅fs-orap非迭代快速生成三维目标的全息图，即只需要生成一个fs-orap便可实现各层相位全息图的生成，极大地提高了效率；其次利用fs-orap，可以生成具有任意位置和大小的目标的幅值相位全息图，克服了原有orap方法在目标位置和大小上的限制。最后，采用菲涅尔相位掩膜进行相位补偿，可以自由选择重建平面。
[0073]
需要说明的是，对于前述的各方法实施例，为了简单描述，故将其都表述为一系列的动作组合，但是本领域技术人员应该知悉，本发明并不受所描述的动作顺序的限制，因为依据本发明，某些步骤可以采用其他顺序或者同时进行。其次，本领域技术人员也应该知悉，说明书中所描述的实施例均属于可选实施例，所涉及的动作和模块并不一定是本发明所必须的。
[0074]
以上是关于方法实施例的介绍，以下通过装置实施例，对本发明所述方案进行进一步说明。
[0075]
图9示出了根据本发明的基于fs-orap和相位补偿的非迭代三维全息图生成装置900的方框图。如图9所示，装置900包括：
[0076]
划分模块910，用于对待处理的三维目标物体沿深度分层，生成多个二维图像；
[0077]
计算模块920，用于根据所述多个二维图像的振幅与事先生成的满支撑随机相位，生成各层纯相位全息图；
[0078]
补偿相位模块930，用于对所述各层纯相位全息图分别添加补偿相位，并将添加补偿相位后得到的各层全息图的复振幅叠加，生成复合相位全息图；
[0079]
重建模块940，用于对所述复合相位全息图进行不同距离菲涅尔变换重建所述三维目标物体。所属领域的技术人员可以清楚地了解到，为描述的方便和简洁，所述描述的模块的具体工作过程，可以参考前述方法实施例中的对应过程，在此不再赘述。
[0080]
本发明的技术方案中，所涉及的用户个人信息的获取，存储和应用等，均符合相关
法律法规的规定，且不违背公序良俗。
[0081]
根据本发明，本发明还提供了一种电子设备、一种可读存储介质和一种计算机程序产品。
[0082]
图10示出了可以用来实施本发明的电子设备110的示意性框图。电子设备旨在表示各种形式的数字计算机，诸如，膝上型计算机、台式计算机、工作台、个人数字助理、服务器、刀片式服务器、大型计算机、和其它适合的计算机。电子设备还可以表示各种形式的移动装置，诸如，个人数字处理、蜂窝电话、智能电话、可穿戴设备和其它类似的计算装置。本文所示的部件、它们的连接和关系、以及它们的功能仅仅作为示例，并且不意在限制本文中描述的和/或者要求的本发明的实现。
[0083]
电子设备110包括计算单元111，其可以根据存储在rom112中的计算机程序或者从存储单元118加载到ram113中的计算机程序，来执行各种适当的动作和处理。在ram113中，还可存储电子设备110操作所需的各种程序和数据。计算单元111、rom112以及ram113通过总线114彼此相连。i/o接口115也连接至总线114。
[0084]
电子设备110中的多个部件连接至i/o接口115，包括：输入单元116，例如键盘、鼠标等；输出单元117，例如各种类型的显示器、扬声器等；存储单元118，例如磁盘、光盘等；以及通信单元119，例如网卡、调制解调器、无线通信收发机等。通信单元119允许电子设备110通过诸如因特网的计算机网络和/或各种电信网络与其他设备交换信息/数据。
[0085]
计算单元111可以是各种具有处理和计算能力的通用和/或专用处理组件。计算单元111的一些示例包括但不限于中央处理单元(cpu)、图形处理单元(gpu)、各种专用的人工智能(ai)计算芯片、各种运行机器学习模型算法的计算单元、数字信号处理器(dsp)、以及任何适当的处理器、控制器、微控制器等。计算单元111执行上文所描述的各个方法和处理，例如方法100。例如，在一些实施例中，方法100可被实现为计算机软件程序，其被有形地包含于机器可读介质，例如存储单元118。在一些实施例中，计算机程序的部分或者全部可以经由rom112和/或通信单元119而被载入和/或安装到电子设备110上。当计算机程序加载到ram 113并由计算单元111执行时，可以执行上文描述的方法100的一个或多个步骤。备选地，在其他实施例中，计算单元111可以通过其他任何适当的方式(例如，借助于固件)而被配置为执行方法100。
[0086]
本文中以上描述的系统和技术的各种实施方式可以在数字电子电路系统、集成电路系统、现场可编程门阵列(fpga)、专用集成电路(asic)、专用标准产品(assp)、芯片上系统的系统(soc)、负载可编程逻辑设备(cpld)、计算机硬件、固件、软件、和/或它们的组合中实现。这些各种实施方式可以包括：实施在一个或者多个计算机程序中，该一个或者多个计算机程序可在包括至少一个可编程处理器的可编程系统上执行和/或解释，该可编程处理器可以是专用或者通用可编程处理器，可以从存储系统、至少一个输入装置、和至少一个输出装置接收数据和指令，并且将数据和指令传输至该存储系统、该至少一个输入装置、和该至少一个输出装置。
[0087]
用于实施本发明的方法的程序代码可以采用一个或多个编程语言的任何组合来编写。这些程序代码可以提供给通用计算机、专用计算机或其他可编程数据处理装置的处理器或控制器，使得程序代码当由处理器或控制器执行时使流程图和/或框图中所规定的功能/操作被实施。程序代码可以完全在机器上执行、部分地在机器上执行，作为独立软件
包部分地在机器上执行且部分地在远程机器上执行或完全在远程机器或服务器上执行。
[0088]
在本发明的上下文中，机器可读介质可以是有形的介质，其可以包含或存储以供指令执行系统、装置或设备使用或与指令执行系统、装置或设备结合地使用的程序。机器可读介质可以是机器可读信号介质或机器可读储存介质。机器可读介质可以包括但不限于电子的、磁性的、光学的、电磁的、红外的、或半导体系统、装置或设备，或者上述内容的任何合适组合。机器可读存储介质的更具体示例会包括基于一个或多个线的电气连接、便携式计算机盘、硬盘、随机存取存储器(ram)、只读存储器(rom)、可擦除可编程只读存储器(eprom或快闪存储器)、光纤、便捷式紧凑盘只读存储器(cd-rom)、光学储存设备、磁储存设备、或上述内容的任何合适组合。
[0089]
为了提供与用户的交互，可以在计算机上实施此处描述的系统和技术，该计算机具有：用于向用户显示信息的显示装置；以及键盘和指向装置(例如，鼠标或者轨迹球)，用户可以通过该键盘和该指向装置来将输入提供给计算机。其它种类的装置还可以用于提供与用户的交互；例如，提供给用户的反馈可以是任何形式的传感反馈(例如，视觉反馈、听觉反馈、或者触觉反馈)；并且可以用任何形式(包括声输入、语音输入或者、触觉输入)来接收来自用户的输入。
[0090]
可以将此处描述的系统和技术实施在包括后台部件的计算系统(例如，作为数据服务器)、或者包括中间件部件的计算系统(例如，应用服务器)、或者包括前端部件的计算系统(例如，具有图形用户界面或者网络浏览器的用户计算机，用户可以通过该图形用户界面或者该网络浏览器来与此处描述的系统和技术的实施方式交互)、或者包括这种后台部件、中间件部件、或者前端部件的任何组合的计算系统中。可以通过任何形式或者介质的数字数据通信(例如，通信网络)来将系统的部件相互连接。通信网络的示例包括：局域网(lan)、广域网(wan)和互联网。
[0091]
计算机系统可以包括客户端和服务器。客户端和服务器一般远离彼此并且通常通过通信网络进行交互。通过在相应的计算机上运行并且彼此具有客户端-服务器关系的计算机程序来产生客户端和服务器的关系。服务器可以是云服务器，也可以为分布式系统的服务器，或者是结合了区块链的服务器。
[0092]
应该理解，可以使用上面所示的各种形式的流程，重新排序、增加或删除步骤。例如，本发明中记载的各步骤可以并行地执行也可以顺序地执行也可以不同的次序执行，只要能够实现本发明的技术方案所期望的结果，本文在此不进行限制。
[0093]
上述具体实施方式，并不构成对本发明保护范围的限制。本领域技术人员应该明白的是，根据设计要求和其他因素，可以进行各种修改、组合、子组合和替代。任何在本发明的精神和原则之内所作的修改、等同替换和改进等，均应包含在本发明保护范围之内。

技术特征：
1.一种基于fs-orap和相位补偿的非迭代三维全息图生成方法，其特征在于，该方法包括：对待处理的三维目标沿深度分层，生成多个二维图像；根据所述多个二维图像的振幅与事先生成的满支撑优化随机相位，生成各层纯相位全息图；对所述各层纯相位全息图分别添加补偿相位，并将添加补偿相位后得到的各层全息图的复振幅叠加，生成复合相位全息图；对所述复合相位全息图进行菲涅尔变换重建所述三维目标。2.根据权利要求1所述的方法，其特征在于，所述满支撑优化随机相位的生成，包括：预先在重建平面生成和slm平面相同大小的满支撑单位振幅；根据所述满支撑单位振幅和初始化随机相位生成满支撑优化随机相位。3.根据权利要求2所述的方法，其特征在于，所述根据所述满支撑单位振幅和初始化随机相位生成满支撑优化随机相位，包括：将所述满支撑单位振幅与所述初始化随机相位相乘，得到所述重建平面上的复振幅；将所述重建平面上的复振幅进行傅里叶逆变换，传播到slm平面；在所述slm平面和所述重建平面之间执行若干次标准的g-s循环，得到所述满支撑优化随机相位。4.根据权利要求1所述的方法，其特征在于，所述根据所述多个二维图像的振幅与事先生成的满支撑优化随机相位，生成各层纯相位全息图，包括：将所述多个二维图像的振幅分别与所述满支撑优化随机相位相乘，并将得到的乘积进行傅里叶逆变换，生成每层的纯相位全息图。5.根据权利要求1所述的方法，其特征在于，所述对所述各层纯相位全息图分别添加补偿相位，包括：将各层纯相位全息图相位分别与不同焦距的菲涅耳相位掩膜叠加，得到新的相位分布。6.根据权利要求1所述的方法，其特征在于，所述对所述复合相位全息图进行菲涅尔变换重建所述三维目标物体，包括：对所述复合相位全息图添加恢复相位，并根据添加恢复相位后的复合相位全息图重建所述三维目标。7.一种基于fs-orap和相位补偿的非迭代三维全息图生成装置，其特征在于，该装置包括：划分模块，用于对待处理的三维目标沿深度分层，生成多个二维图像；计算模块，用于根据所述多个二维图像的振幅与事先生成的满支撑优化随机相位，生成各层纯相位全息图；补偿相位模块，用于对所述各层纯相位全息图分别添加补偿相位，并将添加补偿相位后得到的各层全息图的复振幅叠加，生成复合相位全息图；重建模块，用于对所述复合相位全息图进行菲涅尔变换重建所述三维目标。

技术总结
本发明提供了一种基于FS-ORAP和相位补偿的非迭代三维全息图生成方法及装置。所述方法包括：对待处理的三维目标沿深度分层，生成多个二维图像；根据所述多个二维图像的振幅与事先生成的满支撑优化随机相位，生成各层纯相位全息图；对所述各层纯相位全息图分别添加补偿相位，并将添加补偿相位后得到的各层全息图的复振幅叠加，生成复合相位全息图；对所述复合相位全息图进行菲涅尔变换重建所述三维目标。以此方式，突破了原始ORAP方法在空间域中目标振幅固定支撑的大小和位置的限制，且无需迭代，极大的加快了三维目标相位全息图的生成。极大的加快了三维目标相位全息图的生成。极大的加快了三维目标相位全息图的生成。


技术研发人员：张成 韩鹏 史继森 周昊 吴峰 章权兵 程鸿 沈川
受保护的技术使用者：安徽大学
技术研发日：2023.06.20
技术公布日：2023/9/20