CN110133849B - 一种点扩散函数重建方法 - Google Patents
一种点扩散函数重建方法 Download PDFInfo
- Publication number
- CN110133849B CN110133849B CN201910427576.9A CN201910427576A CN110133849B CN 110133849 B CN110133849 B CN 110133849B CN 201910427576 A CN201910427576 A CN 201910427576A CN 110133849 B CN110133849 B CN 110133849B
- Authority
- CN
- China
- Prior art keywords
- amplitude
- spread function
- point spread
- point
- phase
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G02—OPTICS
- G02B—OPTICAL ELEMENTS, SYSTEMS OR APPARATUS
- G02B27/00—Optical systems or apparatus not provided for by any of the groups G02B1/00 - G02B26/00, G02B30/00
- G02B27/0012—Optical design, e.g. procedures, algorithms, optimisation routines
Abstract
本发明公开了一种点扩散函数重建方法,包括:(1)输入目标点扩散函数;(2)针对输入的目标点扩散函数进行采样;(3)将采样得到的每个点看作一个偶极子进行初始相位和初始振幅估计;(4)计算得到偶极子在瞳面上的辐射场;(5)通过辐射场得到瞳面上的振幅和相位分布;(6)通过得到的振幅和相位生成点扩散函数;(7)将点扩散函数与目标点扩散函数比较,若满足要求,则输出得到的振幅和相位分布;若不满足要求,则更新初始相位和初始振幅,返回步骤(4)。本发明直接由目标点扩散函数得到瞳面上的相位及振幅分布,大大提高了计算效率。并且,由于偶极子辐射的光是向随机方向传播的,该方法可以扩展到多个具有多个光瞳面的显微***。
Description
技术领域
本发明设计光学工程领域,具体是涉及一种点扩散函数重建方法。
背景技术
聚焦光斑(Tightly Focused Spots)广泛应用于光学存储、光学操控、光学微细加工、全息显示器和远场光学显微镜。对于聚焦光斑的不同的应用场景,人们需要对聚焦光斑进行相应的操控。在超分辨显微成像领域,特别是在受激辐射光淬灭显微镜(STEDMicroscopy)中,对于聚焦光斑的操控尤为重要。如果有一个好的聚焦光斑,那么显微镜的分辨率将会更高。
聚焦光斑又被称为点扩散函数(Point Spread Function,简称PSF),是描述光学***对点光源解析能力的函数。因为点光源在经过任何光学***后都会由于衍射而形成一个扩大的像点,这个像点就是点扩散函数。
在给定输入参数的情况下,利用矢量衍射公式得到***的点扩散函数是容易的,但是还无法对其进行操控。为了得到人们所需要的点扩散函数,科学家们研究出了点扩散函数的操控技术(也称为点扩散函数重建技术),即,给定一个目标点扩散函数,通过一定的算法,得到产生该点扩散函数所需的输入参数,这些输入参数构成的函数被称为瞳函数,这些输入参数一般指入射光在瞳面上的振幅和相位。
近些年来,点扩散函数的操控得到了很大的发展。比如游尚挺等人2015年在Optics Letters上面发表的题为《Iterative phase-retrieval method for generatingstereo array of polarization-controlled focal spots》的文章提出了一种基于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个偶极子向瞳面上对应坐标为(θ,φ)的点辐射光波的波矢);且满足:
若M为无限大,上式则可以转化为积分的形式。
步骤(5)中,由步骤(4)中得到的辐射场可得到瞳面上的振幅和相位分布。
瞳面上点的振幅A和相位α为:
其中,K=2π/λ为波数,是以镜头(比如物镜)的理想焦点为原点建立的柱坐标系点坐标;θ是光束的孔径角;是瞳面上的方位角;是入射光在瞳面上的振幅值的分布,即步骤(5)所得的振幅;是入射光的琼斯矩阵;是入射光在瞳面上的相位分布,即步骤(5)所得的相位;是物方到像方的坐标转置矩阵,为已知参数。
步骤(7)中,将步骤(6)中生成的点扩散函数与目标点扩散函数比较,若结果不理想,则改变初始相位通过迭代,再计算瞳面上的振幅和相位分布,直到获得的瞳面上的振幅和相位能够产生符合预期的点扩散函数。步骤(7)中,具体比较时,计算当前得到的点扩散函数与目标点扩散函数各点强度值的均方差,若均方差在允许范围内,则判定为满足要求;反之判定为不满足要求。
该方法的原理是:通过计算结果点扩散函数各点强度值(即振幅的平方)与目标点扩散函数各点强度值的均方差,若均方差在允许范围内,则结果符合预期;若均方差不在允许范围内,则开始迭代,具体的迭代方法为,将计算所得到的振幅函数替换为一个高斯分布或者均匀分布的振幅函数,然后向前计算点扩散函数,并保留相位分布作为初始估计,然后再重新计算瞳函数,这样迭代,直到点扩散函数达到设计目标。
步骤(7)中,当前得到的点扩散函数不满足要求,更新初始相位和初始振幅时:将得到的振幅函数替换为一个高斯分布或者均匀分布的振幅函数,完成初始振幅的更新;将当前得到的相位分布作为初始相位,完成初始相位的更新。
作为优选,若输入的目标点扩散函数为图片,则应将其转换为对应的矩阵,再对其进行稀疏采样。
本发明的方法可以利用计算机完成,即通过计算机编程进行迭代运算,利用本发明的方法得到最佳的点扩散函数参数(入射光在瞳面上的振幅和相位等)后,可以利用下述装置对得到的点扩散函数参数进行验证:
基于上述方法,本发明还提出了一种实现点扩散函数操控的装置。该装置包括:
激光光源、光调制模块、透镜、接收模块。
本发明中,正向入射光为相干准直线偏光,经光调制模块,然后通过透镜出射,到达接收模块。
本发明中,所述光调制模块位于透镜的前焦面上,接收模块位于透镜的后焦面上。
本发明中,所述光调制模块的作用是对入射光进行振幅、相位和偏振调制;优选为空间光调制器。
本发明中,所述透镜为双胶合消色差凸透镜;优选焦距为100mm。若用透镜组替代该透镜,则需要在其前面加一4f***,光调制模块位于4f***第一透镜的后焦面上,4f***的作用是将光调制模块成像到透镜组的入瞳面。
本发明中,所述接收模块用于检测生成的点扩散函数;优选为电荷耦合器件。
本发明的重建方法将目标点扩散函数看作位于不同位置的多个具有不同初始相位和振幅的偶极子的组合,通过对偶极子进行稀疏采样,再计算这些采样得到的偶极子在瞳面上的辐射场,在通过瞳面上的辐射场计算出其振幅和相位分布,再通过避免迭代,或者用少次数迭代,从而快速获得重构点扩散函数所需的瞳函数。上述方法可以避免多次迭代需要的大量时间,并且更加直接地由目标点扩散函数得到瞳面上的相位及振幅分布,因此大大提高了计算效率。并且,由于偶极子辐射的光是向随机方向传播的,该方法可以扩展到多个具有多个光瞳面的显微***。
相对于现有技术,本发明具有以下有益的技术效果:
(1)计算速率高;
(2)结构简单;
(3)可以扩展到具有多个瞳面的光学***。
附图说明
图1为本发明提供的一种用于操纵点扩散函数的相位反演方法流程图;
图2和图3为应用本发明提供的方法得到的试验效果图(实施例);
图4为应用本发明的光路结构示意图(实施例)。
图中:
激光光源1、光调制模块2、透镜3、接收模块4。
具体实施方式
以下结合附图和实施例来详细说明本发明的技术方案,但本发明并不仅限于此。
本发明是在Matlab R2018b语言环境下实现的。首先输入目标点扩散函数,再对目标点扩散函数进行稀疏采样,获得采样点的位置数据;然后对每个采样点分别给予初始振幅和初始相位;再限制瞳面的大小,将每个偶极子带入到辐射场计算公式中,算出这些偶极子在瞳面上的辐射场分布;再通过辐射场分布,得到瞳面上的振幅和相位分布;然后利用瞳面上的振幅和相位分布,通过矢量衍射理论评估所得点扩散函数的质量,若不符合预期,那么通过迭代来使生成的点扩散函数符合预期。
本发明是一种用于操纵点扩散函数的相位反演方法,该方法的流程见图1所示,该方法包括以下步骤:
我们将任意一个目标点扩散函数看作位于不同位置的多个具有不同初始相位和振幅的偶极子的组合。
将点扩散函数的每一个点(像素)看作一个偶极子,将每个点的强度值看作对应偶极子的振幅。
由于在进行图片的拍摄和计算机模拟时,图片即为离散的点阵,所以目标点扩散函数可以被看作位于不同位置的多个具有不同初始相位和振幅的偶极子的组合。
由于一个点扩散函数可以用一张图片或者一个矩阵来描述,所以可以把点扩散函数图片的每一个像素或者矩阵的每一个元素看作一个偶极子。并且其像素值的大小可以看作其振幅的大小。
第一步,输入目标点扩散函数,对其进行稀疏采样,获取采样点的位置数据。
(1)在Matlab R2018b语言环境下输入目标点扩散函数,若输入的目标点扩散函数为图片,则应将其转换为对应的矩阵,再对其进行稀疏采样,得到对应偶极子的坐标信息和强度信息。
很多情况下,采样率不需要太高,举例说明,对于某一个目标点扩散函数,虽然可以用一个1000*1000的矩阵来精确描述一个点扩散函数的强度分布,然而由于光的衍射性质有时候会增加两个偶极子之间的串扰,所以偶极子的分布越密集,计算难度越大,从而有可能导致无法获得一个符合预期的点扩散函数。所以我们可以采用稀疏采样,采集适当数量的偶极子数据,避免串扰的同时,也可以提高计算效率。比如我们可以利用一个位置间隔为0.6λ的3*3的偶极子点阵来计算上述点扩散函数。
第二步,对采集得到的每个偶极子进行初始的振幅和相位估计。
作为一种实施方式,为了简化计算,可以将采样得到的偶极子的初始振幅都设为1,且假设方向相同,初始相位设为π。
第三步,带入辐射场计算公式,得出这些偶极子在瞳面上的辐射场分布。
计算这些偶极子的辐射场,这些偶极子的远场辐射可以分解为平面波在不同方向上的传播的叠加(为波矢),θ是入射光光束的孔径角,设定θ的最大值θmax为镜头的汇聚角(对于点扩散函数来说是发散角),NA表示镜头(可以是物镜)的数值孔径,n是镜头的物方折射率,是瞳面上的方位角。
其中,M为偶极子的总数,AN和ψN分别是第N个偶极子的初始振幅和相位,为第N个偶极子向瞳面上对应坐标为的点辐射光波的波矢,θ和代表瞳面上点的位置坐标,是第N个偶极子的位置矢量,是第N个偶极子的位置坐标。
若M为无限大,上式则可以转化为积分的形式。
第四步,由步骤三中得到的辐射场可得到瞳面上的振幅和相位分布。瞳面上点的振幅A和相位α为:
第五步,将瞳面上的振幅和相位分布带入到矢量衍射理论公式中,算出一个点扩散函数。
在物镜具有高数值孔径的条件下,根据矢量衍射理论,物镜焦点附近的电场矢量可以表示为:
其中,K=2π/λ为波数,是以物镜的理想焦点为原点建立的柱坐标系,θ是光束的孔径角,是瞳面上的方位角。是步骤四所得的入射光在瞳面上的振幅值的分布,是入射光的偏振矩阵或琼斯矩阵,是步骤四所得的入射光在瞳面上的相位分布,是物方到像方的坐标转置矩阵。
再利用公式:
即可得到整个点扩散函数,即上式中的I。
第七步,将步骤六中生成的点扩散函数与目标点扩散函数比较,若结果不理想,则改变初始相位通过迭代,再计算瞳面上的振幅和相位分布,直到获得的瞳面上的振幅和相位能够产生符合预期的点扩散函数。
该方法的原理是:通过计算结果点扩散函数与目标点扩散函数各点强度值的均方差,若均方差在允许范围内,则结果符合预期;若均方差不在允许范围内,则开始迭代,具体的迭代方法为,将计算所得到的振幅函数替换为一个高斯分布或者均匀分布的振幅函数,然后向前计算点扩散函数,并保留相位分布作为初始估计,然后再重新计算瞳函数,这样迭代,直到点扩散函数达到设计目标。
实施例1
如图2所示为本发明的试验效果图YALE。
图中:目标点扩散函数、振幅、相位、结果点扩散函数。
本发明中,将目标点扩散函数YALE输入,对其进行稀疏采样,获得相应偶极子的位置,再对每个偶极子进行初始振幅和相位估计,经过相位反演,得到光瞳函数的振幅和相位,再将振幅和相位带入到矢量衍射理论公式中,得出结果点扩散函数YALE。
实施例2
如图3所示为本发明的试验效果图斗牛犬。
图中:目标点扩散函数、振幅、相位、结果点扩散函数。
本发明中,将目标点扩散函数斗牛犬输入,对其进行稀疏采样,获得相应偶极子的位置,再对每个偶极子进行初始振幅和相位估计,经过相位反演,得到光瞳函数的振幅和相位,再将振幅和相位带入到矢量衍射理论公式中,得出结果点扩散函数斗牛犬。
实施例3
如图4所示为本发明的点扩散函数操控装置图。
图中:激光光源1、光调制模块2、透镜3、接收模块4。
本发明中,正向入射光为相干准直线偏光,经光调制模块,然后通过透镜出射,到达接收模块。
本发明中,所述光调制模块位于透镜的前焦面上,接收模块位于透镜的后焦面上。
本发明中,所述光调制模块的作用是对入射光进行振幅、相位和偏振调制;优选为空间光调制器。
本发明中,所述透镜为双胶合消色差凸透镜;优选焦距为100mm。若用透镜组替代该透镜,则需要在其前面加一4f***,光调制模块位于4f***第一透镜的后焦面上,4f***的作用是将光调制模块成像到透镜组的入瞳面。
本发明中,所述接收模块用于检测生成的点扩散函数;优选为电荷耦合器件。
Claims (10)
1.一种点扩散函数重建方法,其特征在于,包括如下步骤:
(1)输入目标点扩散函数,读取每个点的坐标信息和强度信息;
(2)针对输入的目标点扩散函数进行采样;
(3)将采样得到的每个点看作一个偶极子进行初始相位和初始振幅估计;
(4)计算得到偶极子在镜头瞳面上的辐射场;
(5)通过辐射场得到镜头瞳面上的振幅和相位分布;
(6)通过得到的振幅和相位,经镜头聚焦生成点扩散函数;
(7)将步骤(6)中生成的点扩散函数与目标点扩散函数比较,若满足要求,则输出得到的振幅和相位分布;若不满足要求,则更新初始相位和初始振幅,返回步骤(4)。
2.根据权利要求1所述的点扩散函数重建方法,其特征在于,所述步骤(2)中,针对输入的目标点扩散函数进行采样时采用稀疏采样。
3.根据权利要求2所述的点扩散函数重建方法,其特征在于,以位置间隔为0.5~0.7λ进行稀疏采样,λ为入射光的波长。
4.根据权利要求1所述的点扩散函数重建方法,其特征在于,步骤(3)中,在进行初始振幅估计时,假设每个偶极子的振幅相同。
8.根据权利要求1所述的点扩散函数重建方法,其特征在于,步骤(7)中,计算当前得到的点扩散函数与目标点扩散函数各点强度值的均方差,若均方差在允许范围内,则判定为满足要求;反之判定为不满足要求。
9.根据权利要求1所述的点扩散函数重建方法,其特征在于,步骤(7)中,当前得到的点扩散函数不满足要求,更新初始相位和初始振幅时:将得到的振幅函数替换为一个高斯分布或者均匀分布的振幅函数,完成初始振幅的更新;将当前得到的相位分布作为初始相位,完成初始相位的更新。
10.根据权利要求1所述的点扩散函数重建方法,其特征在于,若输入的目标点扩散函数为图片,则应将其转换为对应的矩阵,再对其进行稀疏采样。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910427576.9A CN110133849B (zh) | 2019-05-22 | 2019-05-22 | 一种点扩散函数重建方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910427576.9A CN110133849B (zh) | 2019-05-22 | 2019-05-22 | 一种点扩散函数重建方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110133849A CN110133849A (zh) | 2019-08-16 |
CN110133849B true CN110133849B (zh) | 2020-06-09 |
Family
ID=67572085
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910427576.9A Active CN110133849B (zh) | 2019-05-22 | 2019-05-22 | 一种点扩散函数重建方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110133849B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112620931B (zh) * | 2020-12-09 | 2022-08-05 | 青岛理工大学 | 一种光束整形的激光加工装置及方法 |
CN112945835A (zh) * | 2021-01-22 | 2021-06-11 | 西安交通大学 | 基于深度变化点扩散函数的细胞精确三维形态恢复方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102540474A (zh) * | 2012-01-11 | 2012-07-04 | 哈尔滨工业大学 | 一种实现边缘陡峭且光强波动低的平顶光束整形控制方法及其整形装置 |
CN103645562A (zh) * | 2013-12-20 | 2014-03-19 | 中南大学 | 一种同时调制振幅和相位的光束整形方法 |
WO2014074138A1 (en) * | 2012-11-09 | 2014-05-15 | Nikon Corporation | Globally dominant point spread function estimation |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
AU2003217081A1 (en) * | 2002-03-17 | 2003-09-29 | Gareth Paul Bell | Optimising point spread function of spatial filter |
WO2007116365A2 (en) * | 2006-04-07 | 2007-10-18 | Ecole Polytechnique Federale De Lausanne (Epfl) | Method and apparatus to measure and compute the amplitude point spread function and associated parameters of a coherent optical imaging system |
US9202265B2 (en) * | 2012-04-30 | 2015-12-01 | Nikon Corporation | Point spread function cost function with non-uniform weights |
CN104777612B (zh) * | 2015-04-28 | 2017-01-25 | 中南大学 | 一种对光束的振幅和相位同时整形的迭代方法 |
-
2019
- 2019-05-22 CN CN201910427576.9A patent/CN110133849B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102540474A (zh) * | 2012-01-11 | 2012-07-04 | 哈尔滨工业大学 | 一种实现边缘陡峭且光强波动低的平顶光束整形控制方法及其整形装置 |
WO2014074138A1 (en) * | 2012-11-09 | 2014-05-15 | Nikon Corporation | Globally dominant point spread function estimation |
CN103645562A (zh) * | 2013-12-20 | 2014-03-19 | 中南大学 | 一种同时调制振幅和相位的光束整形方法 |
Non-Patent Citations (1)
Title |
---|
Point-spread function optimization in isoSTED nanoscopy;Xiang Hao;《Optics letter》;20150729;全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN110133849A (zh) | 2019-08-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Ang et al. | Quantum limit for two-dimensional resolution of two incoherent optical point sources | |
Hou et al. | Deep-learning-based phase control method for tiled aperture coherent beam combining systems | |
Chen et al. | Self-splitting properties of a Hermite-Gaussian correlated Schell-model beam | |
Zhang et al. | Phase retrieval based on wave-front relay and modulation | |
CN110133849B (zh) | 一种点扩散函数重建方法 | |
CN105379253A (zh) | 景深3d成像空间光调制器显微镜 | |
WO2022258075A1 (zh) | 一种基于dmd生成多参量可调光场的方法、装置及*** | |
Underwood et al. | Wave optics approach for incoherent imaging simulation through distributed turbulence | |
Stuart et al. | Fast algorithms for generating binary holograms | |
Foreman et al. | Inversion of the Debye-Wolf diffraction integral using an eigenfunction representation of the electric fields in the focal region | |
Volostnikov et al. | Diffractive elements based on spiral beams as devices for determining the depth of bedding of radiation sources | |
Ferdman et al. | Diffractive optical system design by cascaded propagation | |
Katkovnik et al. | Computational wavelength resolution for in-line lensless holography: phase-coded diffraction patterns and wavefront group-sparsity | |
CN112099228A (zh) | 一种基于rcwa的偏振体光栅衍射光线追迹仿真***及方法 | |
Wang et al. | A phase distribution design method for phased arrays multibeam independently generating and 3-D scanning | |
Chen et al. | Superresolution microscopy imaging based on full-wave modeling and image reconstruction | |
CN113589560B (zh) | 基于角向二值化相位的轨道角动量光梳生成方法与*** | |
Khonina et al. | Defined distribution forming in the near diffraction zone based on expansion of finite propagation operator eigenfunctions | |
US6343307B1 (en) | Synthesis of light beams | |
CN111562707B (zh) | 一种用于超晶格非线性成像过程的迂回编码设计方法 | |
JP7274181B2 (ja) | 円柱光学系を用いた楕円ビーム設計 | |
JPH08110749A (ja) | 位相格子の設計システム | |
Shostka et al. | Assessment of spatial optical trapping with digital holographic sensing for biological micro-particles studying | |
Voelz et al. | Wave and ray optics simulations of short exposure incoherent imaging in atmospheric turbulence | |
Li et al. | Comprehensive holographic parallel beam modulation inside material based on automatic differentiation |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |