CN102625921A - 用于恢复波场的相位的方法和设备 - Google Patents

用于恢复波场的相位的方法和设备 Download PDF

Info

Publication number
CN102625921A
CN102625921A CN2010800410434A CN201080041043A CN102625921A CN 102625921 A CN102625921 A CN 102625921A CN 2010800410434 A CN2010800410434 A CN 2010800410434A CN 201080041043 A CN201080041043 A CN 201080041043A CN 102625921 A CN102625921 A CN 102625921A
Authority
CN
China
Prior art keywords
wave field
plane
described method
function
wave
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.)
Pending
Application number
CN2010800410434A
Other languages
English (en)
Inventor
张福才
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
University of Sheffield
Original Assignee
University of Sheffield
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Priority claimed from GB0916167A external-priority patent/GB0916167D0/en
Priority claimed from GB0917930A external-priority patent/GB2474442A/en
Application filed by University of Sheffield filed Critical University of Sheffield
Publication of CN102625921A publication Critical patent/CN102625921A/zh
Pending legal-status Critical Current

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/0075—Optical systems or apparatus not provided for by any of the groups G02B1/00 - G02B26/00, G02B30/00 with means for altering, e.g. increasing, the depth of field or depth of focus
    • 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/42—Diffraction optics, i.e. systems including a diffractive element being designed for providing a diffractive effect
    • G02B27/44—Grating systems; Zone plate systems
    • G—PHYSICS
    • G03—PHOTOGRAPHY; CINEMATOGRAPHY; ANALOGOUS TECHNIQUES USING WAVES OTHER THAN OPTICAL WAVES; ELECTROGRAPHY; HOLOGRAPHY
    • G03H—HOLOGRAPHIC PROCESSES OR APPARATUS
    • G03H1/00—Holographic processes or apparatus using light, infrared or ultraviolet waves for obtaining holograms or for obtaining an image from them; Details peculiar thereto
    • H—ELECTRICITY
    • H01—ELECTRIC ELEMENTS
    • H01J—ELECTRIC DISCHARGE TUBES OR DISCHARGE LAMPS
    • H01J25/00—Transit-time tubes, e.g. klystrons, travelling-wave tubes, magnetrons
    • G—PHYSICS
    • G02—OPTICS
    • G02B—OPTICAL ELEMENTS, SYSTEMS OR APPARATUS
    • G02B26/00—Optical devices or arrangements for the control of light using movable or deformable optical elements
    • G02B26/08—Optical devices or arrangements for the control of light using movable or deformable optical elements for controlling the direction of light
    • G02B26/0816—Optical devices or arrangements for the control of light using movable or deformable optical elements for controlling the direction of light by means of one or more reflecting elements
    • G02B26/0833—Optical devices or arrangements for the control of light using movable or deformable optical elements for controlling the direction of light by means of one or more reflecting elements the reflecting element being a micromechanical device, e.g. a MEMS mirror, DMD

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Optics & Photonics (AREA)
  • Investigating Or Analysing Materials By Optical Means (AREA)
  • Radar Systems Or Details Thereof (AREA)
  • Holo Graphy (AREA)
  • Polarising Elements (AREA)
  • Diffracting Gratings Or Hologram Optical Elements (AREA)

Abstract

一种恢复波场的相位的方法，包括下列步骤：提供初始平面处的波场的估计φ 0 ；以及经由波场变换装置在进入平面与检测器平面之间来回传播波场，进入平面是具有将波场限定到其内的面积的平面，其中在进入平面应用支持约束，而在检测器平面应用幅值约束，波场变换装置设置成将波场变换函数应用于波场，波场变换函数的特征在于与透镜函数的有限偏离。

Description

用于恢复波场的相位的方法和设备
技术领域
本发明涉及用于恢复波场相位的方法及对应设备。在一些实施例中，提供一种基于通过从物体所散射的辐射所形成的衍射图案的强度测量来构成物体的图像的方法。
背景技术
人们认识到，物体的图像可从该物体所散射的波场的相位和强度的测量来构成。但是，图像检测器通常不能测量波场的相位，而是仅提供强度的测量。所谓的‘相位问题’，即确定波场的相位的问题，一直是更为关注的主题。
对相位问题的解决方案，通常涉及基于衍射图案的强度测量的波场迭代计算。
Sayre提出了单独从有限远物体的衍射强度，来复原从该有限远物体所衍射的波场相位的可能性(D.Sayre，Acta Crystallographica 5，843(1952))。利用有限支持作为物体约束的算法由Fienup(J.R.Fienup，Optics Letters 3，27(1978))于1978年提出。Fienup和其它研究人员广泛地研究了这些算法的依赖于边界形状、对称性和锐度的解的唯一性和收敛性，参见例如R.Barakat和G.Newsam的Journal of Mathematical Physics 25，3190(1984)，R.H.T.Bates，Optics(Jena)61，247(1982)以及J.R.Fienup和C.C.Wackerman的Journal of the Optical Society of America A 3，1897(1986)。
给定紧支撑和特定物体波场分布，已经表明使用模拟数据的重构是可能的。但是，尽管这些理论有进步，也很少有令人信服质量的实验结果被证实，直到Miao以及合作人员于1999年的黄金图案的x射线重构(J.Miao，P.Charalambous，J.Kirz和D.Sayre，Nature(London)400，342(1999))。
后来的实验将这种技术应用于包括单纯金颗粒的样本(J.Miao，T.Ishikawa，B.Johnson，E.H.Anderson，B.Lai和K.O.Hodgson，Phys. Rev. Lett. 89，088303(2002))、复合酵母细胞(D.Shapiro等人，Proc.Natl.Sci. U.S.A.，102，15343(2005))和纳米晶体内部的应变场的断层扫描映射(M.A.Pfeifer，G.J.Williams，I.A.Vartanyants，R.Harder和I.K.Robinson，Nature 442，63(2006))。
同时，还通过各种方式来改进主要基于混合输入输出(HIO)算法的算法。可通过收缩包围算法来避免紧支撑要求，收缩包围算法能够在迭代过程中连同物体一起来细化支持的估计(S.Marchesini等人，Phys.Rev.B 68，140101(2003))。
从迄今为止所证明的单个衍射图案的相位恢复，仍然局限于小的隔离标本。这对它在材料和生物科学中的广泛使用，呈现出了根本的限制。
这种技术中所要求的有限出射波还能够通过采用有限探头照射扩展物体来提供。但是，这种情况下的相位恢复因平滑边界和非负性损失而面临许多困难(J.R.Fienup，J.Opt.Soc.Am.A 4，118(1987)，J.M.Rodenburg和H.M.L.Faulkner，Appl.Opt. 85，4795(2004))。
一般来说，这些算法固有地伴随着平移和厄米对称模糊性，如J.R.Fienup和C.C.Wackerman在J.Opt.Soc.Am. A 3，1897(1986)中所论述的那样。这些‘平凡’解的竞争可引起缓慢收敛。在某些情况下，近来已经表明，如果照射被弯曲并且还是准确已知的，则扩展物体的重构是可能的(B.Abbey等人，Nat.Phys.4，394(2008))。
当前相干衍射成像技术(CDI)中的另一个障碍是对检测器的动态范围和噪声性能所施加的极严格要求。为了覆盖典型衍射图案的全范围，要求数量级为220的动态范围的检测器。这完全超出诸如电荷耦合器件(CCD)之类的常用检测器的能力。束阑必须用于阻挡中心射束，但是这引起所谓的缺失数据问题，缺失数据必须保持较低或者使用通过其它手段所测量的数据来修正。
当前CDI中的另一个基本障碍是收集高角衍射数据的困难。在x射线或电子辐射的情况下，大多数实际感兴趣样本较弱地衍射。因此，只有入射束能量的极小部分将被衍射到高角区域中；衍射到高角(或高阶)区域以便得到更高分辨率图像。甚至对于没有噪声的理想检测器，仍然要求较长时间段使检测器得到充分数量的计数。一些样本无法耐受所需时间长度的辐射。
通过使用更亮和高成本辐射源，能够将数据采集时间减小一些，但是它将引起甚至是增大的样本损坏。
图1(a)示出了已知的实验布置，其中来自源10的辐射束（也可描述为波场）由设置在物平面中的物体20来散射。散射波场30然后入射到设置成测量波场强度的检测器40上。
图1(b)示出迭代计算波场30的相位和幅度的已知方法。支持约束应用于物平面中的波场的估计，其然后传播到检测器平面，在那应用幅值约束(检测器平面处的波的幅值由检测器40对波场强度的测量来确定)。
US6369932、US6545790和US6906839公开用于复原波前的相位信息的系统和方法。这些文献公开了采用准直辐射来照射材料的标本，并且使辐射经过具有预定阻挡图案的光阑或者一个或多个滤波器。记录经过光阑或滤波器的辐射强度分布。这个过程重复进行至少五次，其中每次使用不同的相应光阑或滤波器。
发明内容
在本发明的第一方面，提供一种恢复波场的相位的方法，包括下列步骤：提供初始平面的波场的估计φ 0 ；以及传播来自初始平面的波场，并且在进入平面与检测器平面之间来回传播，该进入平面是一种具有将波场限定到其内的面积的平面，其中在进入平面应用支持约束，而在检测器平面应用幅值约束，并且波场变换装置设置在进入平面与检测器平面之间的波场的路径中，波场变换装置设置成向波场应有波场变换函数，波场变换函数的特征在于与透镜函数的有限偏离。
‘有限’含有与透镜函数的大或小的偏离。换言之，波场变换函数不是理想透镜函数。波场变换函数也不是简单的自由空间传播函数。优选地，与透镜函数的偏离较大。优选地，与透镜函数的偏离足够大，以在该方法规定的迭代次数之内实现波场的相位恢复。
透镜函数表示一种装置的函数，其能够将从一个点传出的波场转换成表现为从不同的点传出或者收敛到不同的点的另一个波场。具体来说，透镜函数是能够将球面波(或者其低阶近似)转换为另一个球面波(或者其低阶近似)的函数。这里提到的术语‘点’还应当被理解为在衍射受限意义上涵盖‘光点’。
为了避免疑问，透镜函数没有包括自由空间传播函数。
有利地，波场变换函数与透镜函数显著偏离。
要理解，词语‘来回’表示波场在进入平面与检测器之间传播，根据实现算法的编程人员或者该方法用户的选择而在第一方向或第二方向开始。
此外，要理解，初始平面可处于进入平面与检测器平面处或之间。作为选择，初始平面可在进入平面的上游。
进入平面是其中已知的将波场通量限定到某个面积的平面。
要理解，本发明的实施例的优点在于，波场的相位可从与波场的强度对应的单个数据集来恢复。换言之，在本发明的一些实施例，不需要得到与检测器所进行的波场强度的不同相应记录对应的多个数据集。这种‘一次通过’的特征具有的优点在于，它使得能够执行图像中强度实时变化的动态事件的相位恢复。例如，可执行从就地实验所得到的图像的相位恢复，以及在发现样本在波场的辐射下作为时间的函数而退化的情况下。
本发明相对于US6369932、US6545790和US6906839具有的优点在于，相位恢复可以仅使用与波场强度的分布对应的单个数据集来执行。本发明还与US6369932、US6545790和US6906839不同，因为本发明要求应用支持约束。
先前，支持约束的应用要求小隔离样本。这种限制通过本发明的实施例来克服，其中采用通常具有强调制性质的波场变换装置。当波场从波场变换装置后向传播到其中应用了支持约束的进入平面时，波场估计的不正确分量被传播到支持的外部，使得支持约束在细化波场的估计中具有增加的有效性。
优选地，在波场沿第一方向经过进入平面与检测器平面之间时波场变换函数被应用于波场，而在波场沿与第一方向相反的第二方向经过进入平面与检测器平面之间时波场变换函数的逆被应用于波场。
该方法还包括通过在进入平面与检测器平面之间来回重复地传播波场，来迭代地计算波场的相位。
在一些实施例中，波场变换装置设置成，使得经过该装置的射线光程长度使产生于其相干加强的强度与理想透镜的强度充分不同，以便在规定迭代次数之内使相位恢复能够执行到所需空间分辨率。
波场变换函数的特征可在于将傅立叶变换应用于波场，此后将波场与调制函数相乘，随后将另一个傅立叶变换应用于波场，调制函数是一种具有与透镜函数有限偏离的函数。
表述‘特征在于’表示波场变换函数可直接采取所述步骤的形式或等效的形式，而无需显式执行所述步骤。
本发明的实施例的优点在于，能够得到高角信号强度，比通过当前CDI技术可得到的那些高角信号强度要高若干数量级，而无需更强辐射源并且不会增加样本损坏。
如上所述，对于x射线或电子辐射，大多数样本较弱地散射，并且因此较低数量的光子或电子散射到高角。期望的是，较大数量的光子或电子散射到高角区域中。
现有技术解决方案涉及使用能够产生高达1033光子/脉冲的诸如自由电子激光器之类的亮辐射源。自由电子激光器会耗费数千万或数亿美元来构建。但是，大量能量将以定向射束的形式而丢失，从而导致样本损坏增加。
本发明实施例的优点在于，可对于具有软边缘和大相位变化的波场或者采用有限探头所照射的扩展物体来恢复相位。
本发明的实施例的优点还在于，相位可从强度的单个记录来恢复，例如衍射图案强度的单个记录。此外，该方法允许从所记录的衍射数据的相位恢复，而无需具有大的动态范围的检测器。
相反，具有数量级为210-214动态范围的检测器通常是足够的。
该方法的一个关键方面是采用具有已知透射或反射函数(或‘传递函数’)的波场变换装置(WTD)。WTD可以是已知的乘性调制函数(一般为复合调制函数)装置，或者已知脉冲响应函数(傅立叶域中的传递函数)系统。可采用多个WTD，以便形成波场变换系统。
本发明的实施例适合与具有下面波长范围和类型的辐射的波场一起使用，如可见光、红外光、紫外光、太赫兹频率的辐射、x射线辐射、电子辐射、中子辐射、以及任何其它适当辐射。在一些实施例中还可采用声波场。
支持约束S优选地按照下式来应用：
其中，φ n 和φ ’ n 分别是第n次迭代的进入波场的当前和更新估计，φ ’ n 在第一次迭代时设置成等于φ n 。
优选地，其中在假定待测量波场具有有效值的像素处S取单位值，而在其它情况下取零值。
可应用幅值约束，以便按照方程式
Figure 426826DEST_PATH_IMAGE002
来确定检测器平面φ D n+1 处的波场的估计，其中A n+1 和? n+1 分别表示检测器处的波场的幅度和相位。
优选地按照方程式来应用幅值约束，其中φ ’D n+1 是应用幅值约束之后在检测器处的波场，I是记录的衍射强度，以及P(I)是强度I的函数。
优选地，P(I)采取形式P(I)=I γ ，其中γ是常数。
该方法优选地包括将γ设置为大致在从大约0.5至2的范围中的值的步骤。
该方法可包括下列步骤：以γ的第一值执行n 1 次迭代，随后以γ的第二值执行n 2 次迭代。
γ的第一值可大于第二值。
γ的第二值可以为0.5。
该方法可包括将β选择为具有在从大约0.4至大约0.8的范围中的值的步骤。
该方法可包括将β选择为具有值0.62的步骤。
该方法可包括选择γ、β和S的值，以使信号误差比SER在大约100次迭代之后具有小于或基本上等于10-5的值的步骤。
该方法的前面可加上提供初始平面处波场的初始估计的步骤。
初始平面可设置在从下述中选择的一个位置处，与进入平面一致的位置、与检测器平面一致的位置、以及进入平面与检测器平面之间的位置。
波场变换装置可设置成呈现从下述之中选取的一个，入射波场的线性响应和非线性响应。
波场变换装置可设置成具有复合传输，该复合传输是一种既呈现损失又呈现相位延迟的传输。
波场变换装置可包括从相位板、一维光栅、二维光栅、晶体板和空间光调制器之中选取的至少一个。
波场变换装置可包括多个交叉耦合的光纤，设置成将从波场变换装置的入口平面入射的光传送到波场变换装置的出口平面。
交叉耦合表示来自一个光纤的光能够耦合到一个或多个其它光纤。
该装置可设置成在装置的入口和出口平面的相应位置之间传送光，使得波场变换函数的特征在于，由于从装置的入口和出口平面处相应光纤端部的相应位置之间的对应性和相应光纤的长度之中选取的至少一个，而与透镜函数的有限偏离。
因此，在一些实施例中，光纤全部具有大致相同的长度，但是设置成通过在波场从入口平面传送到出口平面时有效地交换WTD入口平面处的波场的像素位置来‘加扰’相位和幅度。可选地，光纤还可具有不同长度，由此将相移引入由给定光纤所传送的波场。
进一步备选地，光纤可具有不同长度，但是入口与出口平面之间的像素的位置可以基本上未改变。
波场变换装置可设置为从入射辐射的透射与入射辐射的反射之中选取的一个。
波场变换装置可包括多个像素元件。
该方法可包括相对于检测器和/或进入平面来调整波场变换装置的取向的步骤。
该方法还可包括提供多个波形变换装置的步骤。
该方法可包括提供按照级联配置的多个波形变换装置的步骤，由此波场设置成经由多个装置的每一个而通过进入平面与检测器平面之间。
波场变换函数可以是从离散的和连续的函数之中选取的一个。
该方法可包括提供波形变换装置的步骤，其中波形变换装置包括从像差透镜和具有不可忽视像差的复合透镜系统之中选取的至少一个。
波形变换函数可以是从可逆运算符和非乘性运算符之中选取的一个。
由检测器所记录的数据可设置成对应于从Fraunhofer衍射图案、Fresnel衍射图案和像差图像之中选取的一个。
支持约束可以是从1D区域的长度、2D面积的边界和3D体积之中选取的一个。
支持约束可应用于多个空间分隔的区域。
波场可包括多个空间分隔的区域。
波场可包括从3D波场、2D波场和1D信号之中选取的一个。
波场变换装置可设置成散射从电磁辐射、可见光子、x射线光子、电子、中子和质子之中选取的至少一个。
优选地，波场包括从太赫兹频率辐射、红外辐射、可见光辐射、深紫外线辐射、软x射线辐射和硬x射线辐射之中选取的电磁辐射。
波场可设置成包括基本上相干的辐射。 
波场可设置成基本上由相干辐射组成。 
波场可以是由物体所散射的波场。
该方法可包括计算物体的一个或多个平面中的波场的相位的步骤。
该方法可包括计算波场路径的所需位置处的波场相位和幅度。
在本发明的第二方面，提供一种恢复波场的相位的方法，包括下列步骤：提供设置成将波场变换函数应用于波场的波场变换装置，其中波场变换函数的特征在于是具有与透镜函数的有限偏离的函数；经由波场变换装置将波场从进入平面传递到检测器平面；以及通过检测器来记录检测器平面处的波场的强度，该方法还包括下列步骤：经由波场变换装置在进入平面与检测器平面之间以虚拟方式来回传播波场，其中在进入平面应用支持约束，以及在检测器平面应用幅值约束，幅值约束对应于检测器所记录的强度。
入射到波场变换装置上的波场可以是由物体散射的波场。
在进入平面与检测器平面之间以虚拟方式来回传播波场的步骤之前，可加上提供初始平面处的波场的估计φ 0 的步骤。
在本发明的第三方面，提供用于恢复波场的相位的设备，包括：设置成将波场变换函数应用于波场的波场变换装置，其中波场变换函数的特征在于是具有与透镜函数有限偏离的函数；响应波场的强度的检测器；以及计算机系统，该设备设置成允许波场经由波场变换装置从进入平面传播到检测器，进入平面是一种将波场限定到有限面积的平面，该计算机系统设置成在进入平面与检测器之间来回传播虚拟波场，该计算机系统设置成将支持约束应用于进入平面处的波场以及将幅值约束应用于检测器处的波场，该系统设置成当波场在进入平面与检测器之间传播时将波场变换函数应用于波场，由此恢复波场的所需位置处的波场相位。
计算机系统可设置成通过在进入平面与检测器平面之间来回重复地传播波场，来迭代地计算波场的相位。
计算机系统可设置成在波场沿第一方向经过进入平面与检测器之间时将波场变换函数应用于波场，而在波场沿与第一方向相反的第二方向经过进入平面与检测器之间时应用波场变换函数的逆。
波场变换函数的特征可在于与透镜函数的有限偏离。
在进入平面与检测器之间来回传播虚拟波场的步骤之前，优选地可加入从初始平面传播波场的步骤。
计算机系统可设置成提供初始平面处的波场的估计，并且随后从初始平面并且在进入平面与检测器之间传播波场。
作为替代或补充，计算机系统可设置成提示用户输入初始平面处的波场的估计，并且随后从初始平面以及在进入平面与检测器之间传播波场。
要理解，波场的初始估计可以是波场的源与检测器之间的基本上任何位置处的波场的初始估计。算法则开始于对应这样位置的阶段，即在其上进行波场的初始估计的位置。
因此，如果初始估计是入射平面与WTD平面之间的位置处的估计，则算法可设置成将波场传播到WTD的平面，并且然后在按照图4的流程图继续进行之前应用波场变换函数。备选地，算法可设置成将波场传播到入射平面，并且在按照图4的流程图继续进行之前应用支持约束。
附图说明
现在将参照附图来描述本发明的实施例，附图包括：
图1示出(a) 用于测量散射波场强度的组件的现有技术布置，以及(b) 确定散射波场的幅度和相位的现有技术方法；
图2示出传输几何结构中按照本发明实施例的设备组件的布置；
图3示出其中进入平面、调制器平面和检测器平面是(a) 平行的以及(b) 非平行的反射几何结构中按照本发明一个实施例的设备组件的布置；
图4示出按照本发明一个实施例的、用于确定波场的幅度和相位的算法的步骤；
图5示出按照本发明一个实施例的设备的布置；
图6是按照本发明实施例的波场变换装置的设计的示意图，示出一种装置的一部分表面的 (a) 平面图和(b) 放大的透视图；图6(c)示出结合多个光纤的装置的透视图，而图6(d)和图6(e)分别示出沿从入口平面到出口平面的方向所查看的入口平面和出口平面；
图7示出用于提供一个示例中使用的波场的(a) 幅度和(b) 相位的值的图像，以便证明使用按照本发明的方法的相位恢复；
图8示出作为支持松度Θ的三个值的迭代次数函数的信号误差比(SER)的图表；
图9示出在单叶子植物样本的重构过程期间所记录的数据，示出(a) 采用放大视图插图的由检测器所记录的衍射图案；(b) 其中在迭代计算这个平面处的幅度和相位的过程期间应用支持约束的平面(‘约束应用平面’)处的幅度图；以及样本平面处的(c) 幅度图和(d) 相位图；
图10示出(a) 沿与波场传播到波场变换装置的方向平行的方向所查看的具有2D周期相位结构的波场变换装置的平面图，以及(b) 作为按照本发明一个实施例使用这个波场变换装置来确定波场的幅度和相位的过程期间的迭代次数函数的信号误差比SER的图表；
图11示出(a) 沿与波场传播到波场变换装置的方向平行的方向所查看的具有1D周期相位结构的波前调制器的平面图，以及(b) 作为按照本发明一个实施例使用这个波场变换装置来确定波场的幅度和相位的过程期间的迭代次数函数的信号误差比SER的图表；以及
图12示出(a)用作波场变换装置的成像系统幅度传递函数的相位图；(b) 作为按照本发明一个实施例使用这个波场变换装置来确定幅度和相位的过程期间的迭代次数函数的SER的图表；(c) 物体平面处的波场幅度的重构图；以及(d) 物体平面处的波场相位的重构图。
图13示出用于证明本发明实施例的操作的1D测试信号的(a) 幅度和(b) 相位。
图14示出按照本发明一个实施例的、在用于确定波场幅度和相位的算法的单次迭代之后图13的信号的(a) 幅度和(b) 相位。
图15示出按照本发明一个实施例的、在用于确定波场幅度和相位的算法的多次迭代之后图13的信号的(a) 幅度和(b) 相位。
具体实施方式
在本发明的一个实施例中，提供具有如图2(a)所示所设置的组件的设备100。设备100具有设置成采用辐射来照射物体120的照射源110。由物体120所散射的辐射设置成入射到波场变换装置(WTD)130，并且透过该波场变换装置，WTD 130又可称作波前调制装置。在图2(a)的实例的情况下，WTD 130采取相位板的形式。
由WTD 130所散射的辐射设置成入射到检测器140上。辐射在本文中将描述为一种波场，其特征在于在空间中任何给定位置处的幅度值和相位值。
图2(a)和(b)的配置可称作透射操作模式，因为辐射透过WTD 130，并且入射到位于 WTD 130的相对侧上的检测器140。
图2(a)和(b)中还示出入射平面125。入射平面125是其中应用支持约束的平面。在一些实施例中，入射平面125可称作进入平面125。
入射平面125可在物体120的上游或下游。在一些实施例中，入射平面125设置成处于来自照射源110的辐射束的跨越点处或附近，由此限制入射平面处的照射面积。
在一些备选实施例中，入射平面选择成与物体120的平面一致。入射平面125的其它多个位置也是有用的。
包括多个空间光调制器(SLM)的其它WTD也是有用的。
图2(b)示出与图2(a)相似的布置，其中WTD 130示出为具有相互间隔开的入口平面130A和出口平面130B。
图3(a)和(b)示出以反射操作模式设置的设备。在图3(a)的特定布置中，入射平面225、WTD平面230和检测器平面240各自基本上相互平行。在图3(b)的特定布置中，入射平面225、WTD平面230和检测器平面完全没有相互平行。
要理解，在图3(a)和(b)的布置中，来自源210的辐射由物体220散射到WTD 230，并且散射到设置成检测由WTD 230所‘反射’的光线的检测器240。
WTD 230可以是已知结构的强调制相位分析器，并且可在常规CDI设置中插入样本下游的位置。在这种情况下，设备的布置在某些方面可与Zhang等人(参见以上文献)在多图像重构算法中所公开的相似。
使用WTD的效果是双重的。首先，由检测器所产生的强度样本的相依性得到增强，因为由所照射物体的更大数量的点所散射的波场部分入射到其中发生部分干扰的检测器上。
增强的相依性对相位逆问题提供超定机制。这与其中超定主要通过不相关的多个记录来提供的Zhang等人的公开相反。
其次，衍射图案散布到比原始物体更宽的倒易空间的体积中，从而提供降低的动态范围的更均匀强度分布。与现有技术方法相比，这又具有降低检测器所需动态范围的优点。
在描述图2(a)的布置中，规定至少三个平面的位置。它们是WTD 130所在的平面、检测器140所在的平面以及“入射平面”125。
在图2(a)和(b)的布置中，入射平面125、WTD 130的平面以及检测器140的平面基本上平行，并且沿波场的传播方向分别间隔开距离d 1 和d 2 。在一些备选实施例中，入射平面125、WTD 130的平面以及检测器140的平面不是相互平行的。
在一些实施例中，WTD的位置根据WTD的入口平面130A和WTD的出口平面130B来确定，参见例如图2(b)。
按照本发明方法的要求在于，在入射平面125，波场的范围是有限的，但是波场的边界无需呈现强度的突然变化，即边界可以是软的。
为了实现平面之间的波传播，可采用Fresnel算法。通过适当选择d 1 和d 2 的值(图2) ，能够易于实现Fresnel近似条件。
如果设备被配置成用于远场条件，如同x射线和电子波衍射的情况中那样，则傅立叶变换可用作射束传播器。
在Fresnel算法中，在不同平面的取样间隔是相关的。作为举例，如果检测器具有N×N像素，其各为正方形并且具有边Δx D ，则在WTD的取样间隔Δx M 和入射平面Δx为：
Figure 878667DEST_PATH_IMAGE004
其中λ是所采用的辐射的波长。
如果入射平面125定义为与物体120一致，则Δx也是物体120图像的可实现空间分辨率，它是使用从检测器140所得到的幅度信息以及按照本方法所确定的相位信息所重构的。
在上述推导方程式1和2中，没有做与在三个平面处的波场如何快速改变相关的假设。在平滑波场的情况下，粗略取样间隔可以是足够的，并且可有利地使用其它已知射束传播算法，如角频谱方法。例如，它们可允许更大的视场。
在一些实施例中，限制WTD 130的传输面积，以便所产生衍射波能够由检测器140充分地取样。在图2的实施例中，WTD 130的边长度设置为具有长度NΔx，即刚好足够大以满足检测器140平面中的波场的尼奎斯特取样要求。
相位复原方法采取图4示意性地所示的算法400的形式。在一些实施例中，该方法开始于作为进入平面处波场的进入波场的估计。进入波场可写作φ 0 (pΔx,qΔx)，其中p和q是离散空间坐标。
该方法如下所述进行。首先，应用410支持约束，以便得到波场的另一个估计：
Figure 477139DEST_PATH_IMAGE001
其中，φ n 和φ ' n 分别是第n次迭代中的入射波场的当前和更新估计。为了清晰和简明起见，省略了空间坐标。这个公式又可称作‘更新公式’。
φ ' n  在第一次迭代设置为等于φ n 。S表示支持约束，并且对于假定待测量波场具有有效值的像素取单位值，而在其它情况下取零值。参数β能够被调整成以改变反馈强度，并且取大约0.4至大约0.8的范围中的值。值β＝0.62一直用于在这里所呈现的模拟和实验重构，除非另加说明。
其次，将波场传播420到WTD 130的平面。
随后，将波场与WTD的复合传输相乘430，以便确定遇到WTD之后的波前的预计形式。
然后，将波场传播440到检测器140的平面，从而产生
Figure DEST_PATH_IMAGE005
；其中A n+1 和φ n+1 分别表示幅度和相位。
由检测器140所测量的辐射的强度(作为幅值的平方)是已知的，并且相应地，按照该方法的下一个步骤是应用450幅值约束：
Figure 960685DEST_PATH_IMAGE006
其中I是记录的衍射强度。参数γ能够在从大约0.5至大约2的范围中调整。可发现，收敛与γ的值极大地相关；大值引起逐次迭代的解的大变化，并且能够快速确定波场的轮廓，但是具有较差质量。
在一些实施例中，首先的n 1 次迭代以γ的大值来执行，最后的n 2 次迭代以γ的较小值来执行，如γ=0.5。因此，总迭代次数为n 1 +n 2 。
然后，将波场后向传播460到WTD，并且通过将WTD处的波场除以470 WTD的传输函数，去除WTD的影响。
随后，将波场后向传播480到进入平面，从而产生进入场的更新估计φ ' n+1 。
上述方法步骤按照迭代方式重复进行，直到连续估计之间的改进变为充分小，或者直到已经执行给定迭代次数。
如上所述，要理解的是，波场的初始估计可以是波场的源与检测器之间的基本上任何位置处的波场的初始估计。算法则开始于对应这样位置的阶段，即在其上进行波场的初始估计的位置。
因此，如果初始估计是入射/进入平面与WTD平面之间的位置处的估计，则算法可设置成将波场传播到WTD的平面，并然后在按照图4的流程图继续进行之前应用波场变换函数。备选地，算法可设置成将波场传播到入射平面，并且在按照图4的流程图继续进行之前应用支持约束。
更新公式，即上述方程式(3)，与HIO算法中使用的不同。该公式选择为符合上述方程式(4)中的参数γ的引入。
在步骤5，即450，参数γ可逐步改变。在一些实施例中，可以发现，可通过随迭代进行而将γ逐渐减小到值0.5，来得到更好的总收敛速率。在一些实施例中，参数γ没有改变。在一些实施例中，参数γ没有在每一次迭代时改变；而是在规定次数处改变参数γ，例如在预定迭代数之后，如交替迭代。
其它形式的幅值约束也是有用的。
上述方法400在计算装置中通过算法来实现，并以表示波场强度和相位的人工构成的数据集来运行。
该算法采用各种波场来测试，包括具有较硬边界的波场、具有较软边界的波场、具有基本上平坦的相位变化的波场、以及具有越过波场范围[-π, π]的相位大变化的波场。
图5示出其中聚光透镜305用于把来自源的波场聚焦到样本320的布置的示例。由样本320所散射的辐射经过采取薄板形式的WTD 330，并且随后入射到对波场的强度敏感的检测器340上。
图6示出按照本发明一个实施例的WTD的示例。图6(a)的暗像素化图案示出WTD的部分，其设置成降低由那些部分所透射的波场的幅度。备选地，图6(a)的暗像素化图案示出WTD的部分，其设置成延迟由那些部分所透射的波场的相位。图6(b)示出相位板表面的一部分的透视图，示出相位板厚度的离散变化，厚度变化的阶跃变化处于相邻像素之间的边界处。
图6(c)示出按照本发明一个实施例的另一个WTD 630的示例。WTD 630具有相互间隔开并且具有延伸于它们之间的多个光纤632的入口平面630A和出口平面630B。入射到入口平面630上的波场由光纤632传送到出口平面630B。光纤632设置成使得在沿相同方向（例如沿入口平面630A到出口平面630B的方向）查看平面630A、630B时，从出口平面出现的波场经过变换函数将入口平面630A的指定坐标(X,Y)处的至少一个光纤632映射到出口平面630B的不同指定坐标(X+a,X+b)。
因此，再现在出口平面630B的波场实际上是入口平面630A处的波场的‘加扰’版本。
图6(d)是WTD 630的入口平面630A的示意图。分别在位置(X,Y)和(X ' ,Y ' )示出第一和第二光纤632A、632B的每个的第一自由端。
图6(e)是如沿从入口平面630A到出口平面630B的方向所查看的WTD 630的出口平面630B的示意图。在位置(X+a,Y+b)示出第一光纤632A的第二自由端，其中a和b是非零值，以及在位置 (X ' +a ' ,Y ' +b ' )示出第二光纤632B的第二自由端，其中a ' 和b ' 是非零值。
因此，入射到入口平面630A上的波场出现在出口平面630B，其中与没有引入‘加扰’的波场相比具有强度和相位的空间重新布置，例如a、b、a ' 和b ' 全部为零的情况。
其它布置也是有用的。
在一些实施例中，计算装置可设置成提供从进入平面到检测器或者进入平面的上游的任何规定平面处的幅度值和相位值。例如，计算装置可设置成提供样本中任何所需位置处的幅度值和相位值，由此提供样本内部体积的图像或者任何其它所需图像，例如样本的整个体积的传输图像。
示例1
呈现从具有较软边缘和较强相位变化的波场测量所得到的结果。
常规CDI方法解决这种情况有困难。
波场使用图7(a)和(b)所示的图像来生成。将图7(a)的图像像素的强度值与爱里斑的对应像素的强度值相乘，以便定义波场的幅度。结合爱里斑对图像提供软边界。幅度缩放到范围[0,1]。
波场相位的对应值使用缩放到[0,2π]范围相位的图7(b)的图像来定义。重叠在图7(b)图像上的虚线圆表示应用于图7(a)的爱里斑的第一零值轮廓的对应位置。
用于算法的设置参数为：λ=635 nm；d 1 =9.7 mm；d 2 =47.7 mm，以及Δx D =7.4μm。
对于这个示例，WTD选择为带有设计图案的相位板，其具有在相位延迟上的基本随机的空间变化，该板的各位置具有0或π的相位延迟。
计算衍射图案，该图案具有量化为212个等级的256×256样本。
实际上，可能难以定位波的真实边界。在这种情况下，可引入说明支持约束的松度的因子：
Figure 2010800410434100002DEST_PATH_IMAGE007
其中，B和D分别是波场范围和支持的线性维数，如图7(a)所示。
在本示例中，真实进入波是已知的，并因此算法的收敛能够直接使用信号误差比来测量
Figure 455252DEST_PATH_IMAGE008
其中，求和针对所有取样指数。
SER是通常使用的归一化RMS误差量度的倒数。但是，由于在本算法中引入了参数γ，衍射平面中的误差量度不适合在这里使用。对于首先的n 1 次迭代，所计算的RMS非常大并且是无意义的；相反，SER提供小的值。作为迭代次数函数的SER的测量值，还能够从作为迭代次数函数的SER图表的斜率来提供关于算法的收敛行为的信息。
图8针对Θ的三个值示出作为迭代次数函数的SER的图表。从图表中可以理解，上述算法甚至当存在所提供支持的大量不定因素时，也迅速收敛。在Θ=1.4的情况下，B和D分别取值116和138。
对于本示例，支持比围绕所有边的实际边界要宽11个像素。随着松度Θ的增加，最终SER减小，并且收敛斜率减慢。当相继的SER的相对变化小于10-5时，迭代终止。
使用n1的不同值，以便得到三个曲线。作为一般规则，当Θ较大时，n1的大值是优选的。一旦SER>100，重构幅度与原始值凭肉眼是不可区分的。此外，相位的计算值与原始值相差恒定的偏差，并因此在这里没有示出重构图像。
还执行了具有对称幅度和相位分布的波场的模拟。也得到类似收敛性能。
由于该方法不要求波场的明确定义的边界，所以其中应用支持约束的入射平面不一定必须位于物体平面。
在一些实施例中，发现即使所使用的支持大小(单位为像素)保持相同，使用波场具有最小范围的平面也给出最快收敛和最佳图像质量。
在一些实施例中，发现进入平面处定义为FF=(B/L) 2 的填充因子是确定收敛中的最关键参数，其中L是进入平面处的视场的线性维数。
在上述模拟中，使用具有随机相位图和像素大小Δx M 的WTD。实际上，期望WTD的特征大小较大，以便促进WTD的易于制造。如果辐射源的相干性不是限制，WTD像素大小能够通过改变距离d 1 和d 2 来自由选择。在一些实施例中，WTD的特征大小能够选择为比Δx M 大许多。
执行了实验，其中执行了WTD的像素的合成。对于测试性进入波场，图2的幅度和相位图被重新调整大小，以给出填充因子1/9。对于4×和6×Δx M 值的板像素大小，得到100的SER值所需迭代次数分别为44和223。
对于具有弯曲照射的弱相位物体，有可能在WTD中使用甚至更大的像素尺寸。考虑8 keV的可能的x射线实验，其中预期分辨率Δx=20 nm，d 2 =8m，N=256，以及Δx D =24μm。按照方程式(1)和(2)计算距离d 1 为3.3 mm和Δx M =200 nm。
具有大约1μm特征大小的WTD是采用当前制造技术易于实现的，并且能够在这种配置中使用。
为了算法的更好收敛或者为了整个系统的更好能量效率，在WTD的传输波形设计中，还存在大的灵活性。
如果在制造之后能够得到其准确调制函数或传递函数，则在算法中能够考虑WTD中的任何制造误差。
术语‘调制函数’应理解为指的是乘性波场变换函数(透射或反射)，而术语‘传递函数’应理解为指的是例如作为透镜的特性的卷积波场变换函数。
本文中提到的‘变换函数’应理解为包括提到的调制函数或传递函数。
在采用作为相移板的WTD的一些实施例中，调制函数能够直接从WTD的表面形貌的测量来计算。在一些实施例中，表面形貌测量可使用共焦显微镜或表面轮廓仪来进行。
图10(a)是具有二维周期相位结构的WTD的示意图。图10(b)示出作为按照本发明一个实施例的、使用具有这种二维周期相位结构的WTD算法的迭代次数函数的信号误差比(SER)的图表。
图11(a)是具有一维周期相位结构的WTD的示意图。图11(b)示出作为按照本发明一个实施例的、使用具有这种一维周期相位结构的WTD算法的迭代次数函数的信号误差比(SER)的对应图表。从图10(b)和图11(b)能够看到，与采用具有一维周期相位结构的WTD相比，采用具有二维周期相位结构的WTD要求更低迭代次数。
可使用具有周期调制函数的一维或二维WTD的事实，在其中采用x射线或电子辐射的应用中是重要的，因为结晶材料板可用作WTD。在一些实施例中，其它的材料板也是有用的，包括单晶和多晶材料板。
还可使用反射WTD，如C.Kohler、F.Zhang和W.Osten所述论述的(Applied Optics 48，4003(2009))。这对于要求使用这样一种波长的应用是特别重要的，对于该波长其反射分量比折射分量更易于可用，如较短波长。
图12(a)是具有像差并且用作WTD的成像系统的幅度传递函数的相位图。图12(b)示出作为按照本发明一个实施例的、使用图12(a)所示WTD算法的迭代次数函数的SER的对应图表。图12(c)和图12(d)分别示出物体平面处的(c) 幅度和(d) 相位的重构图。
示例2
来自635 nm激光二极管的光束由焦距为50 mm的透镜来准直和会聚，以便提供如图5所示的照射探头。将WTD放在射束的跨越点后面大约18.45 mm的距离。WTD从蚀刻成具有变化厚度的硅石玻璃来形成，以便传递所需相位延迟。
WTD形成为具有1100×1100像素，按照图6所示的实施例每个的形状为正方形并且为16μm宽。每个像素设置有针孔，像素阵列由此提供相位图。针孔具有6μm的孔尺寸，以使因像素之间的过渡边缘而引起的伪像为最小。
应理解，不需要针孔，以及没有包括针孔的布置也是有用的。
将具有各边为7.4μm的正方形像素的CCD照相机放在WTD下游70 mm，以便记录衍射图案。
微观单叶子植物标本用作测试样本。测试样本放在WTD上游19.88 mm的位置。
图9(a)示出记录的衍射图案，其示出因使用WTD而引起的均匀分布的、完全展开的斑纹图案。图案的一部分经过放大，并且示为插图。
如上所述，按照本发明实施例的WTD设置成散射入射波场，使得衍射图案中心射束的图像强度因散射到衍射图案的较高角区域而降低，由此显著增强暗场信号。
这对于辐射敏感样本会是特别有利的。例如，某些样本暴露于x射线或电子辐射能够引起巨大的辐射损坏。因此，现有技术要求使用光阑以便记录高角衍射数据，这可通过本发明的一些实施例来克服。
本发明的实施例还允许采用动态范围降低的检测器。这再次是暗场信号强度的增强和中心射束强度降低的结果。
在本示例中，衍射图案的中心376×376样本用于重构样本的图像。数量376按照方程式1来计算，以便满足所需缩放关系。
图9(b)示出其中在n1=30、n2=20次迭代之后，应用支持约束的入射平面中幅度的重构。入射平面选择为探头的跨越点平面。所使用的支持的边界由图9(b)的虚线表示。
从图9(b)中应理解，支持实际上比较宽松(具有较高‘松度’)，因为支持约束的面积比其上入射波场具有大信号强度的面积大许多。
在这个迭代过程中没有应用支持细化算法，如收缩包围算法。但是，采用支持细化算法可引起算法的甚至更迅速收敛。
图9(c)示出样本平面处的波场的幅度的图，以及图9(d)示出这个位置处对应的波场的相位图。
示例3
众所周知的是，一维(1D)信号的相位恢复比较高维情况(2D或3D)要难许多，因为相位问题本身变得更可能欠确定。恢复一维信号的相位具有许多应用，例如超短脉冲的形状确定以及在大地测量等等中。这个示例证明，按照本发明实施例的方法通过数值实验还能够同样适用于1D信号。
测试了不同种类的信号。在这里，选择如图13(a)和(b)所示的具有强变化相位和软边缘的信号。对于这种信号，现有方法面临严重困难。调制器的调制函数仅在相位方面具有变化，它均匀地分布在0与2π的范围之内。强度图通过Fresnel射束传播算法来生成。相位恢复的过程开始于猜测波，其幅度和相位如图14(a)和(b)所示。幅度是具有滚下边缘的调制高斯脉冲；相位是截取的正弦波。图14(a)中的矩形表示应用支持约束的区域1400。图15示出120次迭代之后的重构幅度和相位。能够理解，应用的支持约束1400实际上大于实际信号程度。对支持紧密度的不敏感性证明，这种方法具有优于其它基于Fienup算法方法的极大优点。
作为举例来选择图13所示的信号。还测试了其它信号，包括具有强随机相位的信号–这是实际上会遇到的最困难情况。得到了类似收敛性能。
本发明的实施例提供了用于测量波场相位的新方法。这种技术适合于具有弱或强相位变化的复值场。该技术克服了使用单个衍射图案测量的当前CDI方法的隔离样本要求。与当前CDI方法相比，它还极大地增强收集高角衍射数据的能力。宽松支持足以提供快速收敛。该方法涉及较简单的实验布置，并且对外部振动是不敏感的，并因此可易于适用于现场应用，显示了优于基于二射束干扰技术方法的优点，如干涉测量和轴外全息摄影。本发明的实施例适合于实时应用以及在短时标上发生的现象的研究。
该方法与一种极有希望的解决方案相兼容，该解决方案使用脉冲激光器解决样本损坏问题，参见例如H.N.Chapman等人，Nat.Phys. 2，839(2006)。
要理解，WTD可采取许多不同形式中的一个或多个。
例如，WTD可以是板，并且波前变换函数可采取乘性透射/反射函数的形式。
WTD还可以是系统。波前变换函数可采取脉冲响应函数(或者傅立叶域中的传递函数)的形式。
WTD的其它形式也是有用的，包括具有线性或非线性响应的装置。
在简单乘性装置的情况下，WTD可以是相位板，例如具有复合传输(即，具有损失和相位延迟)的相位板。
在一些实施例中，WTD是一维或二维光栅。WTD例如可以是设置成散射辐射的晶体板，诸如x射线辐射、电子、中子、质子或者任何其它适当辐射。
WTD可包括空间光调制器。WTD可以是反射和/或透射的。
WTD可按照相对于检测器的倾斜取向来设置。
在一些实施例中，WTD具有像素化结构。WTD可设置成为与一个或多个其它WTD结合的级联。
与WTD关联的调制函数可以是离散的。备选地，调制函数可以是连续的。
在复合系统(通过卷积所述)的情况下，WTD可包括像差透镜和/或具有一定量像差的复合透镜系统。
WTD可以是可逆运算符。WTD可以是非乘性的。
要理解，由检测器所记录的数据可对应于从Fraunhofer衍射图案、Fresnel衍射图案和像差图像之中选取的一个。
模数约束可涉及强度的一般非线性函数。
本发明的实施例对于大范围应用提供对波场的一般相位测量问题的解决方案。除了将CDI转变为供材料科学和生物医学科学中使用的例行技术的可能性之外，该方法除了其它应用之外还在计量学和波场感测方面得到应用。
在一个实施例中，波场变换装置设置成与检测器邻接。在一些实施例中，物体设置成与波场变换装置邻接。在一些实施例中，物体、波场变换装置和检测器各设置成相互邻接。因此，提到波场的传播，意指的是考虑WTD函数的虚拟波场。
在本说明书的描述和权利要求书中，词语“包括”和“包含”及其变化，如“包含”，通篇表示“包括但不限于”，并且不是要(并且不是)排除其它部分、添加物、组件、整体或步骤。
在本说明书的描述和权利要求书中，单数通篇包含复数，除非上下文另有要求。具体来说，在使用不定冠词的情况下，本说明书应被理解为考虑复数以及单数，除非上下文另有要求。
结合本发明的具体方面、实施例或示例所述的特征、整体、特性、化合物、化学部分或编组，要被理解为可适用于本文所述的任何其它方面、实施例或示例，除非与其不相容。 

Claims (53)

1. 一种恢复波场的相位的方法，包括下列步骤：
提供初始平面处的波场的估计φ 0 ；以及
经由波场变换装置(130，230，330，630)在进入平面(125，225)与检测器平面(140，240，340)之间来回传播所述波场，所述进入平面是具有将所述波场限定到其内的面积的平面，
其中，在所述进入平面应用(410)支持约束，而在所述检测器平面应用幅值约束(450)，所述波场变换装置设置成将波场变换函数(430，470)应用于所述波场，其中所述波场变换函数的特征在于与透镜函数的有限偏离。
2. 如权利要求1所述的方法，其中，在所述波场沿第一方向经过所述进入平面与所述检测器平面之间时，所述波场变换函数被应用(430)于所述波场；而在所述波场沿与所述第一方向相反的第二方向经过所述进入平面与所述检测器平面之间时，所述装置将所述波场变换函数的逆应用(470)于所述波场。
3. 如权利要求1或2所述的方法，包括步骤：通过在所述进入平面与所述检测器平面之间来回重复传播所述波场，来迭代地计算所述波场的相位。
4. 如以上权利要求中的任一项所述的方法，其中，所述波场变换函数的特征在于：将傅立叶变换应用于所述波场，此后将所述波场与调制函数相乘，随后将另一个傅立叶变换应用于所述波场，所述调制函数是具有与透镜函数有限偏离的函数。
5. 如以上权利要求中的任一项所述的方法，其中，所述支持约束S按照下述方程式来应用：
其中φ n 和φ ’ n 分别是第n次迭代的进入波场的当前和更新估计，φ ’ n 在第一次迭代被设置成等于φ n 。
6. 如权利要求5所述的方法，其中，S在假定待测量的所述波场具有有效值的像素处取单位值，而在其它情况下取零值。
7. 如以上权利要求中的任一项所述的方法，应用所述幅值约束，以便按照方程式
Figure DEST_PATH_IMAGE004
来确定所述检测器平面φ D n+1 处的所述波场的估计，其中A n+1 和? n+1 分别表示所述检测器处的所述波场的幅度和相位。
8. 如权利要求7所述的方法，其中，按照等方程式
Figure DEST_PATH_IMAGE006
来应用所述幅值约束，其中φ 'D n+1 是应用所述幅值约束之后在所述检测器处的所述波场，I是记录的衍射强度，以及P(I)是强度I的函数。
9. 如权利要求8所述的方法，其中，P(I)采取形式P(I)=I γ ，其中γ是常数。
10. 如权利要求9所述的方法，包括步骤：将γ设置为大致在从大约0.5至2的范围中的值。
11. 如权利要求9或10所述的方法，包括下列步骤：以γ的第一值执行n 1 次迭代，随后以γ的第二值执行n 2 次迭代。
12. 如权利要求11所述的方法，其中，γ的所述第一值大于所述第二值。
13. 如权利要求11或12所述的方法，其中，γ的所述第二值为0.5。
14. 如权利要求5或者权利要求6至从属到权利要求6的权利要求13中的任一项所述的方法，包括步骤：将β选择成具有从大约0.4至大约0.8的范围中的值。
15. 如权利要求14所述的方法，包括步骤：将β选择成具有值0.62。
16. 如从属至权利要求5的权利要求9或者通过从属于权利要求9而从属于权利要求5的权利要求10至15中的任一项所述的方法，包括将γ、β和S的值选择成使信号误差比SER在大约100次迭代之后具有小于或者基本上等于10-5的值。
17. 如以上权利要求中的任一项所述的方法，在该方法前面可添加提供初始平面处的所述波场的初始估计的步骤。
18. 如权利要求17所述的方法，其中，所述初始平面设置在从下述位置中选择的一个之处，与所述进入平面一致、与所述检测器平面一致、以及在所述进入平面与所述检测器平面之间。
19. 如以上权利要求中的任一项所述的方法，其中，所述波场变换装置设置成呈现从对入射波场的线性响应和非线性响应之中选取的一个。
20. 如以上权利要求中的任一项所述的方法，其中，所述波场变换装置设置成具有复合传输，该复合传输是既呈现损失和又呈现相位延迟的传输。
21. 如以上权利要求中的任一项所述的方法，其中，所述波场变换装置包括从相位板、一维光栅、二维光栅、晶体板和空间光调制器之中选取的至少一个。
22. 如以上权利要求中的任一项所述的方法，其中，所述波场变换装置包括多个交叉耦合的光纤，该交叉耦合的光纤设置成将从所述波场变换装置的入口平面入射的光传送到所述波场变换装置的出口平面。
23. 如权利要求22所述的方法，其中，所述装置设置成在所述装置的入口和出口平面的相应位置之间传送光，使得所述波场变换函数的特征在于，由于从所述装置的入口和出口平面处的相应光纤端部的相应位置之间的对应性和相应光纤长度之中选取的至少一个，而与透镜函数的有限偏离。
24. 如以上权利要求中的任一项所述的方法，其中，所述波场变换装置设置为从入射辐射的透射和入射辐射的反射之中选取的一个。
25. 如以上权利要求中的任一项所述的方法，其中，所述波场变换装置包括多个像素元件。
26. 如以上权利要求中的任一项所述的方法，包括步骤：相对于所述检测器和/或所述进入平面来调整所述波场变换装置的取向。
27. 如以上权利要求中的任一项所述的方法，包括步骤:提供多个波形变换装置。
28. 如以上权利要求中的任一项所述的方法，包括步骤:提供按照级联配置的多个波形变换装置，由此波场设置成经由所述多个装置的每一个而通过进入平面与检测器平面之间。
29. 如以上权利要求中的任一项所述的方法，其中，所述波场变换函数是从离散的和连续的之中选取的一个。
30. 如以上权利要求中的任一项所述的方法，包括步骤:提供波形变换装置，其包括从像差透镜和具有不可忽视像差的复合透镜系统之中选取的至少一个。
31. 如以上权利要求中的任一项所述的方法，其中，所述波形变换函数是从可逆运算符和非乘性运算符之中选取的一个。
32. 如以上权利要求中的任一项所述的方法，其中，由所述检测器记录的数据设置成对应于从Fraunhofer衍射图案、Fresnel衍射图案和像差图像之中选取的一个。
33. 如以上权利要求中的任一项所述的方法，其中，所述支持约束是从1D区域的长度、2D面积的边界和3D体积之中选取的一个。
34. 如以上权利要求中的任一项所述的方法，其中，将所述支持约束应用于多个空间分隔区域。
35. 如以上权利要求中的任一项所述的方法，其中，所述波场是从3D波场、2D波场和1D信号之中选取的一个。
36. 如以上权利要求中的任一项所述的方法，其中，所述波场变换装置设置成散射从电磁辐射、可见光子、x射线光子、电子、中子和质子之中选取的至少一个。
37. 如以上权利要求中的任一项所述的方法，其中，所述波场包括从太赫兹频率辐射、红外辐射、可见光辐射、深紫外线辐射、软X射线辐射和硬X射线辐射之中选取的电磁辐射。
38. 如以上权利要求中的任一项所述的方法，其中，所述波场设置成包括基本上相干的辐射。
39. 如以上权利要求中的任一项所述的方法，其中，所述波场设置成基本上由相干辐射组成。
40. 如以上权利要求中的任一项所述的方法，其中，所述波场是由物体所散射的波场。
41. 如权利要求40所述的方法，包括步骤:计算物体的一个或多个平面中的所述波场的相位。
42. 如以上权利要求中的任一项所述的方法，包括在所述波场路径的所需位置处计算所述波场的相位和幅度。
43. 一种恢复波场的相位的方法，包括下列步骤：
提供设置成将波场变换函数(430，470)应用于波场的波场变换装置(130，230，330，630)，其中所述波场变换函数的特征在于是具有与透镜函数的有限偏离的函数；
经由所述波场变换装置将波场从进入平面(125，225)传递到检测器平面(140，240，340)，并且通过检测器(340)来记录所述检测器平面处的所述波场的强度，
所述方法还包括下列步骤：
经由所述波场变换装置在所述进入平面与所述检测器平面之间以虚拟方式来回传播所述波场，其中在所述进入平面应用(410)支持约束，以及在所述检测器平面应用(450)幅值约束，所述幅值约束对应于所述检测器记录的强度。
44. 如权利要求43所述的方法，其中，入射到所述波场变换装置的所述波场是由物体(120，220，320)散射的波场。
45. 如权利要求43或44中的任一项所述的方法，其中，在所述进入平面与所述检测器平面之间以虚拟方式来回传播所述波场的步骤之前，添加提供初始平面处的所述波场的估计φ 0 的步骤。
46. 一种用于恢复波场的相位的设备，包括：
设置成将波场变换函数(430，470)应用于波场的波场变换装置(130，230，330，630)，其中所述波场变换函数的特征在于是具有与透镜函数的有限偏离的函数；
响应所述波场强度的检测器(140，240，340)；以及
计算机系统，
所述设备设置成允许波场经由所述波场变换装置从进入平面(125，225)传播到所述检测器，所述进入平面是一种将所述波场限定到有限面积的平面，
所述计算机系统设置成在所述进入平面与所述检测器之间来回传播虚拟波场，
所述计算机系统设置成将支持约束应用(410)于所述进入平面处的波场并且将幅值约束应用(450)于所述检测器处的波场，所述系统设置成在所述进入平面与所述检测器之间传播所述波场时将所述波场变换函数应用于所述波场，由此恢复所述波场的所需位置处的所述波场的相位。
47. 如权利要求46所述的设备，其中，所述计算机系统设置成通过在所述进入平面与所述检测器平面之间来回重复地传播所述波场，来迭代地计算所述波场的相位。
48. 如权利要求46或47所述的设备，其中，所述系统设置成：在所述波场沿第一方向经过所述进入平面与所述检测器之间时，将所述波场变换函数应用于所述波场；而在所述波场沿与所述第一方向相反的第二方向经过所述进入平面与所述检测器之间时，应用所述波场变换函数的逆。
49. 如权利要求46至48中的任一项所述的设备，其中，在所述进入平面与所述检测器之间来回传播虚拟波场的步骤之前，添加从初始平面传播所述波场的步骤。
50. 如权利要求49所述的设备，其中，所述计算机系统设置成提供所述初始平面处的所述波场的估计，并且随后从所述初始平面并在所述进入平面与所述检测器之间传播所述波场。
51. 如权利要求49或50中的任一项所述的设备，其中，所述计算机系统设置成提示用户输入所述初始平面处的所述波场的估计，并且随后从所述初始平面并在所述进入平面与所述检测器之间传播所述波场。
52. 一种计算机程序，包括使计算机执行如权利要求1至45中的任一项所述的方法的程序指令。
53. 一种其上具有计算机程序代码部件的计算机程序产品，当所述程序被加载时使所述计算机按照权利要求1至45中任一项所述的方法来恢复波场的相位。
CN2010800410434A 2009-09-15 2010-09-10 用于恢复波场的相位的方法和设备 Pending CN102625921A (zh)

Applications Claiming Priority (5)

Application Number Priority Date Filing Date Title
GB0916167.0 2009-09-15
GB0916167A GB0916167D0 (en) 2009-09-15 2009-09-15 Method and apparatus for retrieving a phase of a wavefield
GB0917930A GB2474442A (en) 2009-10-13 2009-10-13 Retrieving a phase of a wavefield
GB0917930.0 2009-10-13
PCT/GB2010/051516 WO2011033287A1 (en) 2009-09-15 2010-09-10 Method and apparatus for retrieving a phase of a wavefield

Publications (1)

Publication Number Publication Date
CN102625921A true CN102625921A (zh) 2012-08-01

Family

ID=43084462

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2010800410434A Pending CN102625921A (zh) 2009-09-15 2010-09-10 用于恢复波场的相位的方法和设备

Country Status (8)

Country Link
US (1) US20120179425A1 (zh)
EP (1) EP2478407A1 (zh)
JP (1) JP2013504780A (zh)
CN (1) CN102625921A (zh)
AU (1) AU2010297015A1 (zh)
CA (1) CA2774082A1 (zh)
GB (1) GB2486369A (zh)
WO (1) WO2011033287A1 (zh)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105548080A (zh) * 2016-01-15 2016-05-04 北京工业大学 一种连续太赫兹波空间扫描相干衍射成像系统及方法
CN106338823A (zh) * 2016-10-27 2017-01-18 中国科学院光电技术研究所 一种基于混合焦距菲涅耳波带片的相位反演方法
CN108027475A (zh) * 2015-06-01 2018-05-11 光学技术注册协会莱布尼兹研究所 与多模光纤有关的改进
CN112236725A (zh) * 2018-08-23 2021-01-15 杜尔利塔斯有限公司 全息图计算的方法
CN115047636A (zh) * 2022-06-20 2022-09-13 北京无线电计量测试研究所 一种高斯型太赫兹宽带空间滤波器

Families Citing this family (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB201020516D0 (en) 2010-12-03 2011-01-19 Univ Sheffield Improvements in providing image data
GB201107053D0 (en) 2011-04-27 2011-06-08 Univ Sheffield Improvements in providing image data
GB201112119D0 (en) 2011-07-14 2011-08-31 Univ Sheffield Method and apparatus for position determination
GB201201140D0 (en) * 2012-01-24 2012-03-07 Phase Focus Ltd Method and apparatus for determining object characteristics
GB201207800D0 (en) 2012-05-03 2012-06-13 Phase Focus Ltd Improvements in providing image data
JP6283875B2 (ja) * 2013-09-05 2018-02-28 キヤノンメディカルシステムズ株式会社 医用画像処理装置、ｘ線診断装置およびｘ線コンピュータ断層撮影装置
CN104215339B (zh) * 2014-09-14 2017-02-15 中国科学院光电技术研究所 一种基于连续远场的波前复原系统及方法
CN104266769B (zh) * 2014-10-23 2017-03-22 北京理工大学 一种相位恢复方法
US9791390B2 (en) * 2015-01-22 2017-10-17 EDAX, Incorporated Devices and systems for spatial averaging of electron backscatter diffraction patterns
CN105675151B (zh) * 2016-01-06 2018-10-23 上海大学 一种基于光强传输方程相位恢复应用装置
US9911180B2 (en) * 2016-04-04 2018-03-06 Raytheon Company Computational imaging with uncalibrated pupil phase
CN111474188B (zh) * 2019-01-24 2021-07-27 中国科学院上海光学精密机械研究所 基于动态调制的单次曝光波前重建和相位成像装置和方法
CN113093499B (zh) * 2021-04-15 2022-03-01 中国地质大学（北京） 一种离散孔径内推的太赫兹数字全息成像方法和系统

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6545790B2 (en) * 1999-11-08 2003-04-08 Ralph W. Gerchberg System and method for recovering phase information of a wave front

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CA2390546A1 (en) 1999-11-08 2001-05-17 Wavefront Analysis Inc. System and method for recovering phase information of a wave front

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6545790B2 (en) * 1999-11-08 2003-04-08 Ralph W. Gerchberg System and method for recovering phase information of a wave front

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
C. KOHLER ET AL.: "Characterization of a spatial light modulator and its application in phase retrieval", 《APPLIED OPTICS》, vol. 48, no. 20, 10 July 2009 (2009-07-10), pages 4003 - 4008, XP001524393, DOI: doi:10.1364/AO.48.004003 *
FUCAI ZHANG ET AL.: "Phase retrieval of arbitrary complex-valued fields through aperture-plane modulation", 《PHYSICAL REVIEW A》, vol. 75, no. 4, 4 April 2007 (2007-04-04), XP002611492, DOI: doi:10.1103/PhysRevA.75.043805 *
G. J. WILLIAMS ET AL.: "Fresnel Coherent Diffractive Imaging", 《PHYSICAL REVIEW LETTERS》, vol. 97, no. 2, 14 July 2006 (2006-07-14), XP002611494, DOI: doi:10.1103/PhysRevLett.97.025506 *
I. JOHNSON ET AL.: "Coherent Diffractive Imaging Using Phase Front Modifications", 《PHYSICAL REVIEW LETTERS》, vol. 100, no. 15, 18 April 2008 (2008-04-18), XP002611493, DOI: doi:10.1103/PHYSREVLETT.100.155503 *

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108027475A (zh) * 2015-06-01 2018-05-11 光学技术注册协会莱布尼兹研究所 与多模光纤有关的改进
CN105548080A (zh) * 2016-01-15 2016-05-04 北京工业大学 一种连续太赫兹波空间扫描相干衍射成像系统及方法
CN105548080B (zh) * 2016-01-15 2018-07-31 北京工业大学 一种连续太赫兹波空间扫描相干衍射成像系统及方法
CN106338823A (zh) * 2016-10-27 2017-01-18 中国科学院光电技术研究所 一种基于混合焦距菲涅耳波带片的相位反演方法
CN106338823B (zh) * 2016-10-27 2020-10-13 中国科学院光电技术研究所 一种基于混合焦距菲涅耳波带片的相位反演方法
CN112236725A (zh) * 2018-08-23 2021-01-15 杜尔利塔斯有限公司 全息图计算的方法
CN115047636A (zh) * 2022-06-20 2022-09-13 北京无线电计量测试研究所 一种高斯型太赫兹宽带空间滤波器

Also Published As

Publication number Publication date
GB2486369A (en) 2012-06-13
AU2010297015A1 (en) 2012-04-05
WO2011033287A1 (en) 2011-03-24
CA2774082A1 (en) 2011-03-24
JP2013504780A (ja) 2013-02-07
EP2478407A1 (en) 2012-07-25
US20120179425A1 (en) 2012-07-12
GB201204434D0 (en) 2012-04-25

Similar Documents

Publication Publication Date Title
CN102625921A (zh) 用于恢复波场的相位的方法和设备
US6545790B2 (en) System and method for recovering phase information of a wave front
JP4997334B2 (ja) 超解像を達成するための新規なデジタル方法を有する光学顕微鏡
US5384573A (en) Image synthesis using time sequential holography
US11293806B2 (en) Wavefront sensor and method of using it
Lu et al. Noniterative spatially partially coherent diffractive imaging using pinhole array mask
US6369932B1 (en) System and method for recovering phase information of a wave front
CN108332866A (zh) 基于振幅调制的波前在线快速重建装置和重建方法
CN103559686A (zh) 一种基于多面图像信息的同轴相衬成像相位恢复方法及系统
CN103323938B (zh) 一种获取立体显微图像的方法
CN103559698A (zh) 一种基于混合迭代的同轴相衬成像相位恢复方法及系统
Hallada et al. Fresnel zone plate light field spectral imaging
Spencer et al. Digital holography wave-front sensing in the presence of strong atmospheric turbulence and thermal blooming
Pedersen et al. Numerical study of Bragg CDI on thick polycrystalline specimens
Fevola et al. A Monte Carlo ray-tracing simulation of coherent X-ray diffractive imaging
Banet et al. Speckle decorrelation effects on motion-compensated, multi-wavelength 3D digital holography: theory and simulations
GB2474442A (en) Retrieving a phase of a wavefield
Eriksson Development of a Software Tool for Mid-Spatial Frequency Analysis
Tappin et al. Simulating diffraction effects in heliospheric imagers
Wei et al. Parameter retrieval of small particles in dark-field Fourier ptychography and a rectangle in real-space ptychography
Zhang et al. A numerical method of improved bandwidth adaptability for simulating diffraction integral of a spatial light modulator
Wynne et al. Enhanced optical alignment of a digital micro mirror device through Bayesian adaptive exploration
Monnier et al. GPU-based simulation of optical propagation through turbulence for active and passive imaging
Kazakov et al. Using the statistical properties of phase fluctuations of the light field in a Shack-Hartmann wave front sensor
Wei et al. Parameter retrieval methods in ptychography

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C02 Deemed withdrawal of patent application after publication (patent law 2001)
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20120801