CN114167707B - 一种离轴数字全息参考光估计方法及系统 - Google Patents

一种离轴数字全息参考光估计方法及系统 Download PDF

Info

Publication number
CN114167707B
CN114167707B CN202111490469.4A CN202111490469A CN114167707B CN 114167707 B CN114167707 B CN 114167707B CN 202111490469 A CN202111490469 A CN 202111490469A CN 114167707 B CN114167707 B CN 114167707B
Authority
CN
China
Prior art keywords
index
matrix
spectrum
reference light
frequency maximum
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
Application number
CN202111490469.4A
Other languages
English (en)
Other versions
CN114167707A (zh
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.)
Laser Fusion Research Center China Academy of Engineering Physics
Original Assignee
Laser Fusion Research Center China Academy of Engineering Physics
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
Application filed by Laser Fusion Research Center China Academy of Engineering Physics filed Critical Laser Fusion Research Center China Academy of Engineering Physics
Priority to CN202111490469.4A priority Critical patent/CN114167707B/zh
Publication of CN114167707A publication Critical patent/CN114167707A/zh
Application granted granted Critical
Publication of CN114167707B publication Critical patent/CN114167707B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • 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
    • G03H1/04—Processes or apparatus for producing holograms
    • G03H1/08—Synthesising holograms, i.e. holograms synthesized from objects or objects from holograms
    • G03H1/0866—Digital holographic imaging, i.e. synthesizing holobjects from holograms
    • 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
    • G03H1/04—Processes or apparatus for producing holograms
    • G03H1/0443—Digital holography, i.e. recording holograms with digital recording means
    • 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
    • G03H1/04—Processes or apparatus for producing holograms
    • G03H1/08—Synthesising holograms, i.e. holograms synthesized from objects or objects from holograms
    • G03H1/0808—Methods of numerical synthesis, e.g. coherent ray tracing [CRT], diffraction specific
    • 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
    • G03H1/04—Processes or apparatus for producing holograms
    • G03H1/08—Synthesising holograms, i.e. holograms synthesized from objects or objects from holograms
    • G03H1/0808—Methods of numerical synthesis, e.g. coherent ray tracing [CRT], diffraction specific
    • G03H2001/0833—Look up table
    • 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
    • G03H1/04—Processes or apparatus for producing holograms
    • G03H1/08—Synthesising holograms, i.e. holograms synthesized from objects or objects from holograms
    • G03H1/0866—Digital holographic imaging, i.e. synthesizing holobjects from holograms
    • G03H2001/0883—Reconstruction aspect, e.g. numerical focusing

Abstract

本发明涉及一种离轴数字全息参考光估计方法及系统，属于数字全息成像领域，对全息图像进行离散傅里叶变换，并将零频分量转移到全息图像的频谱中心，得到+1级像频谱中频率最大值的矩阵索引；并确定傅里叶变换上采样频谱图中局部的索引区域；从中寻找索引区域中频率最大值的位置，并确定索引区域中频率最大值位置的行索引和列索引；根据索引区域中频率最大值位置的行索引和列索引，计算傅里叶变换上采样频谱图中+1级像频谱中频率最大值位置的行索引和列索引；根据+1级像频谱中频率最大值位置的行索引和列索引，计算参考光的角度余弦信息；根据参考光的角度余弦信息，计算估计的参考光的光场分布信息，能够提高参考光的估计精度和重建质量。

Description

一种离轴数字全息参考光估计方法及系统
技术领域
本发明涉及数字全息成像领域，特别是涉及一种离轴数字全息参考光估计方法及系统。
背景技术
离轴数字全息由于引入了倾斜平面参考光，使得零级像、+1级像和-1级像能够在频域中分离，其图像重建过程中需要准确得到平面参考光倾角以消除倾斜相位畸变。传统的获取离轴全息平面参考光信息的方法是将全息图变换到频域，通过定位+1级像或-1级像频谱的最大值点来直接提取参考光倾角。然而，由于像素尺寸以及工作波长的综合影响，这种方法对平面参考光倾斜角度估计不准确，从而提取的参考光信息会存在较大误差，进而使得重建质量下降。特别是在太赫兹波段，其工作波长与像素尺寸的比值较可见光波段大数10倍，导致在太赫兹波段以传统方法估计的参考光倾角误差范围是可见光波段下的数10倍。因此，目前亟需一种离轴数字全息参考光估计方法及系统，以解决现有技术中对平面参考光倾斜角度的估计精度低，重建质量差的问题。
发明内容
本发明的目的是提供一种离轴数字全息参考光估计方法及系统，以提高对平面参考光倾斜角度的估计精度，从而提升离轴数字全息的重建质量。
为实现上述目的，本发明提供了如下方案：
一方面，本发明提出了一种离轴数字全息参考光估计方法，包括：
对全息图像进行离散傅里叶变换，并将离散傅里叶变换后的全息图像的零频分量转移到全息图像的频谱中心，对所述频谱中心两侧±1级像频谱中频率最大值的位置坐标进行定位，得到+1级像频谱中频率最大值的矩阵索引；所述±1级像频谱包括+1级像频谱和-1级像频谱，所述+1级像频谱表示物光波与参考光波共轭的乘积项的频谱，所述-1级像频谱表示物光波共轭与参考光波的乘积项的频谱；
根据所述矩阵索引的行索引fr、列索引fc以及全息图像的离散傅里叶变换的上采样率f，确定全息图像的傅里叶变换上采样频谱图fHU中局部的索引区域fH_in；所述索引区域fH_in为由行索引范围(fr-2)*f+1到fr*f，列索引范围(fc-2)*f+1到fc*f构成的大小为2f*2f的索引区域；
从所述索引区域fH_in中寻找频率最大值的位置，得到频率最大值位置的行索引fr_in和列索引fc_in；
根据所述索引区域中频率最大值位置的行索引fr_in和列索引fc_in，计算所述傅里叶变换上采样频谱图fHU中+1级像频谱中频率最大值位置的行索引fr_u和列索引fc_u；
根据所述+1级像频谱中频率最大值位置的行索引fr_u和列索引fc_u，计算参考光的角度余弦信息；
根据所述参考光的角度余弦信息，计算得到估计的参考光的光场分布信息。
另一方面，本发明还提出了一种离轴数字全息参考光估计系统，包括：
+1级像频谱频率最大值矩阵索引获取模块，用于对全息图像进行离散傅里叶变换，并将离散傅里叶变换后的全息图像的零频分量转移到全息图像的频谱中心，对所述频谱中心两侧±1级像频率最大值的位置坐标进行定位，得到+1级像频谱中频率最大值的矩阵索引；所述±1级像频谱包括+1级像频谱和-1级像频谱，所述+1级像频谱表示物光波与参考光波共轭的乘积项的频谱，所述-1级像频谱表示物光波共轭与参考光波的乘积项的频谱；
索引区域确定模块，用于根据所述矩阵索引的行索引fr、列索引fc以及全息图像的离散傅里叶变换的上采样率f，确定全息图像的傅里叶变换上采样频谱图fHU中局部的索引区域fH_in；所述索引区域fH_in为由行索引范围(fr-2)*f+1到fr*f，列索引范围(fc-2)*f+1到fc*f构成的大小为2f*2f的索引区域；
索引区域频率最大值位置确定模块，用于从所述索引区域fH_in中寻找频率最大值的位置，得到频率最大值位置的行索引fr_in和列索引fc_in；
+1级像频谱频率最大值位置精确获取模块，用于根据所述索引区域中频率最大值位置的行索引fr_in和列索引fc_in，计算所述傅里叶变换上采样频谱图fHU中+1级像频谱中频率最大值位置的行索引fr_u和列索引fc_u；
参考光的角度余弦信息计算模块，用于根据所述+1级像频谱中频率最大值位置的行索引fr_u和列索引fc_u，计算参考光的角度余弦信息；
参考光光场分布估计模块，用于根据所述参考光的角度余弦信息，计算得到估计的参考光的光场分布信息。
根据本发明提供的具体实施例，本发明公开了以下技术效果：
本发明提出了一种离轴数字全息参考光估计方法，通过对全息图像进行离散傅里叶变换，由于离散傅里叶变换后零频分量位于图像的左上角，因此将零频分量转移到全息图像的频谱中心，通过对频谱中心两侧±1级像频谱中频率最大值位置坐标进行定位，从而得到+1级像频谱中频率最大值的矩阵索引，实现了对+1级像频谱中频率最大值位置的粗定位。然后基于傅里叶变换矩阵乘法的思想，对+1级像频谱中心进行局部上采样，确定全息图像的傅里叶变换上采样频谱图中由行索引范围和列索引范围构成的局部索引区域；通过从该索引区域中寻找索引区域中频率最大值的位置的方法，确定出+1级像频谱中频率最大值位置的行索引和列索引，从而通过对全息图像进行局部上采样傅里叶变换，实现了对频谱中心处的+1级像频谱中频率最大值位置的精定位，能够准确、可靠地确定+1级像频谱中频率最大值位置。最终根据确定出的+1级像频谱中频率最大值位置的行索引和列索引，计算出参考光的角度余弦信息，并通过角度余弦信息计算的方法估计出参考光的光场分布信息，从而实现了高效率、高精度的平面参考光的估计，能够提高对平面参考光倾斜角度的估计精度，进而提升了离轴数字全息的重建质量。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。以下附图并未刻意按实际尺寸等比例缩放绘制，重点在于示出本发明的主旨。
图1为本发明实施例1提供的离轴数字全息参考光估计方法的流程图；
图2为本发明实施例1提供的振幅型USAF1951样品的振幅分布图；
图3为本发明实施例1提供的采用传统方法进行参考光估计后得到的振幅重建结果的示意图；
图4为本发明实施例1提供的采用本发明方法进行参考光估计后得到的振幅重建结果的示意图；
图5为本发明实施例2提供的离轴数字全息参考光估计系统的结构框图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动的前提下所获得的所有其他实施例，都属于本发明保护的范围。
如本发明和权利要求书中所示，除非上下文明确提示例外情形，“一”、“一个”、“一种”和/或“该”等词并非特指单数，也可包括复数。一般说来，术语“包括”与“包含”仅提示包括已明确标识的步骤和元素，而这些步骤和元素不构成一个排它性的罗列，方法或者设备也可能包含其他的步骤或元素。
虽然本发明对根据本发明的实施例的系统中的某些模块做出了各种引用，然而，任何数量的不同模块可以被使用并运行在用户终端和/或服务器上。所述模块仅是说明性的，并且所述系统和方法的不同方面可以使用不同模块。
本发明中使用了流程图用来说明根据本发明的实施例的系统所执行的操作。应当理解的是，前面或下面操作不一定按照顺序来精确地执行。相反，根据需要，可以按照倒序或同时处理各种步骤。同时，也可以将其他操作添加到这些过程中，或从这些过程移除某一步或数步操作。
本发明的目的是提供一种离轴数字全息参考光估计方法及系统，以提高对平面参考光倾斜角度的估计精度，从而提升离轴数字全息的重建质量。
基于此技术目的，本发明首先通过对+1级像频谱中频率最大值位置进行粗定位，然后基于傅里叶变换矩阵乘法的思想，对+1级像频谱中心区域进行局部上采样，以对+1级像频谱中频率最大值位置的精定位，从而实现了高效率、高精度的平面参考光的估计，能够提高对平面参考光倾斜角度的估计精度，从而提升离轴数字全息的重建质量。
为使本发明的上述目的、特征和优点能够更加明显易懂，下面结合附图和具体实施方式对本发明作进一步详细的说明。
实施例1
如图1所示，本实施例提供了一种离轴数字全息参考光估计方法，该方法具体步骤如下：
步骤S1、对全息图像进行离散傅里叶变换，并将离散傅里叶变换后的全息图像的零频分量转移到全息图像的频谱中心，对所述频谱中心两侧±1级像频谱中频率最大值的位置坐标进行定位，得到+1级像频谱中频率最大值的矩阵索引。
本实施例中，首先对全息图像进行离散傅里叶变换，由于离散傅里叶变换后零频分量位于全息图像的左上角位置，因此，本实施例将零频分量转移到全息图像的频谱中心，通过对频谱中心两侧±1级像频谱的频率最大值位置坐标进行定位，粗略地确定出以频谱中心即零频中心位置为中心，+1级像频谱中频率最大值的矩阵索引，实现了对+1级像频谱中频率最大值位置的粗定位。
其中，±1级像频谱包括+1级像频谱和-1级像频谱，所述+1级像频谱表示物光波与参考光波共轭的乘积项的频谱，所述-1级像频谱表示物光波共轭与参考光波的乘积项的频谱。本实施例中，在整个全息图像离散傅里叶变换后的频谱中心的两侧分别为+1级像频谱和-1级像频谱，在对频谱中心两侧±1级像频谱的频率最大值位置坐标进行定位后，选取±1级像频谱中的+1级像频谱进行后续的分析处理，根据定位后得到的+1级像频谱中频率最大值的位置坐标建立矩阵索引，设定该矩阵索引中的行索引fr和列索引fc。
本实施例中，在对全息图像进行离散傅里叶变换后，可以采用Matlab软件中的FFTShift函数，对离散傅里叶变换后的全息图像的零频分量进行FFT Shift变换，从而将离散傅里叶变换后的全息图像的零频分量转移到全息图像的频谱中心位置。
步骤S2、根据所述矩阵索引的行索引fr、列索引fc以及全息图像的离散傅里叶变换的上采样率f，确定全息图像的傅里叶变换上采样频谱图fHU中局部的索引区域fH_in；所述索引区域fH_in为由行索引范围(fr-2)*f+1到fr*f，列索引范围(fc-2)*f+1到fc*f构成的大小为2f*2f的索引区域。
本发明在对+1级像频谱中频率最大值位置进行粗定位后，基于傅里叶变换矩阵乘法的思想，设置离散傅里叶变换的上采样率f，在基于全息图像的傅里叶变换上采样频谱图fHU中构建出了多个矩阵，这些矩阵结合起来形成了索引区域fH_in。具体包括：
步骤S2.1、利用公式分别构建大小为2f*2f的矩阵EXP_HU和大小为M*N的矩阵EXP_H，矩阵中第i行第j列的元素分别表示为：
EXP_HU(i,j)＝exp(iπ((fr-2)*f+i+(fc-2)*f+j))；
Figure BDA0003399105030000061
其中，M表示全息图像H的矩阵总行数，N表示全息图像H的矩阵总列数；i表示矩阵中第i行，j表示矩阵中第j列。
步骤S2.2、利用公式构建大小为2f*M的矩阵RMatrixK，矩阵中第i行第j列的元素表示为：
Figure BDA0003399105030000062
步骤S2.3、利用公式构建大小为N*2f的矩阵CMatrixK，矩阵中第i行第j列的元素表示为：
Figure BDA0003399105030000063
步骤S2.4、根据矩阵EXP_HU、矩阵EXP_H、矩阵RMatrixK和矩阵CMatrixK，计算得到全息图像的傅里叶变换上采样频谱图fHU中由行索引范围(fr-2)*f+1到fr*f，列索引范围(fc-2)*f+1到fc*f构成的大小为2f*2f的索引区域fH_in，表示为：
fH_in＝EXP_HU⊙(RMatrixK*(H⊙EXP_H)*CMatrix)；
其中，H表示全息图像，⊙为矩阵对应元素相乘符号。
步骤S3、从所述索引区域fH_in中寻找频率最大值的位置，得到频率最大值位置的行索引fr_in和列索引fc_in。本实施例中，从所述索引区域fH_in中寻找频率最大值的位置，也就是从索引区域fH_in中确定+1级像频谱中频率最大值所在的矩阵索引。
步骤S4、根据所述索引区域中频率最大值位置的行索引fr_in和列索引fc_in，计算全息图像的傅里叶变换上采样频谱图fHU中+1级像频谱中频率最大值位置的行索引fr_u和列索引fc_u。具体包括：
步骤S4.1、根据索引区域中频率最大值位置的fr_in，计算+1级像频谱中频率最大值位置的行索引fr_u，表示为：
fr_u＝(fr-2)*f+fr_in；
步骤S4.2、根据索引区域中频率最大值位置的列索引fc_in，计算+1级像频谱中频率最大值位置的列索引fc_u，表示为：
fc_u＝(fc-2)*f+fc_in。
步骤S5、根据所述+1级像频谱中频率最大值位置的行索引fr_u和列索引fc_u，计算参考光的角度余弦信息。具体包括：
步骤S5.1、根据所述+1级像频谱中频率最大值位置的行索引fr_u和列索引fc_u，计算零频分量转移到频谱中心的位移量。
本实施例中，计算零频分量转移到频谱中心的位移的公式为：
Figure BDA0003399105030000081
其中，M表示全息图像H的矩阵总行数，N表示全息图像H的矩阵总列数；fx_u和fy_u分别表示零频分量转移到频谱中心在x轴上产生的位移以及在y轴上产生的位移。
步骤S5.2、根据所述零频分量转移到频谱中心的位移量，计算参考光的角度余弦信息。
本实施例中，计算参考光的角度余弦信息的公式为：
Figure BDA0003399105030000082
其中，cosα表示参考光与x轴之间角度的余弦值，cosβ表示参考光与y轴之间角度的余弦值，λ表示波长，Δp表示探测器像素尺寸。
步骤S6、根据所述参考光的角度余弦信息，计算得到估计的参考光的光场分布信息。
本实施例中，计算估计的参考光的光场分布信息的公式为：
Figure BDA0003399105030000083
其中，ReferenceWave(x,y)表示估计的参考光光场分布坐标，i为复数，i2＝-1。
本发明通过对全息图像进行离散傅里叶变换，将零频分量转移到全息图像的频谱中心，通过对频谱中心两侧的频率最大值位置坐标进行定位，从而得到+1级像频谱中频率最大值的矩阵索引，实现了对+1级像频谱中频率最大值位置的粗定位。然后基于傅里叶变换矩阵乘法的思想，在频谱中心局部进行上采样傅里叶变换，实现了对频谱中心处的+1级像频谱中频率最大值位置的精定位，能够准确、可靠地确定+1级像频谱中频率最大值位置，最终根据确定出的+1级像频谱中频率最大值位置的行索引和列索引，计算出参考光的角度余弦信息，并通过角度余弦信息计算的方法估计出参考光的光场分布信息，从而实现了高效率、高精度的平面参考光的估计，能够提高对平面参考光倾斜角度的估计精度，进而提升了离轴数字全息的重建质量。
为了进一步说明本发明相对于现有技术的效果，现列举一实例具体说明：
以振幅型USAF1951样品为例在太赫兹波段开展仿真计算，样品振幅分布如图2所示。离轴全息仿真过程中，设定波长为118.8μm，探测器像素尺寸为17μm，物体离探测器距离10mm，参考光与x，y轴角度均为60°。采用传统参考光信息提取方法得到参考光角度为59.668°，角谱重建结果如图3所示。利用本发明参考光估计方法，当上采样率设为1000时，得到参考光角度为60.0063°，角谱重建结果如图4所示。对比图3和图4可直观看出，图4中的图像重建质量比图3好很多。因此，相比传统估计方法，本发明方法能够高精度地估计出参考光的光场分布信息，能够显著提高重建质量。
容易理解的是，本实施例设置的上述样品型号、波长、探测器像素尺寸、物体离探测器距离、参考光与x，y轴角度以及上采样率等具体数值，仅仅是列举出的一组数值，这些数值并不是固定的、唯一的，还可以使用其他的数值，可根据实际情况自行设定。
实施例2
如图5所示，本实施例提供了一种离轴数字全息参考光估计系统，该系统各个模块的功能与实施例1方法各个步骤相同且一一对应，该系统具体包括：
+1级像频谱频率最大值矩阵索引获取模块M1，用于对全息图像进行离散傅里叶变换，并将离散傅里叶变换后的全息图像的零频分量转移到全息图像的频谱中心，对所述频谱中心两侧±1级像频谱中频率最大值的位置坐标进行定位，得到+1级像频谱中频率最大值的矩阵索引；所述±1级像频谱包括+1级像频谱和-1级像频谱，所述+1级像频谱表示物光波与参考光波共轭的乘积项的频谱，所述-1级像频谱表示物光波共轭与参考光波的乘积项的频谱；
索引区域确定模块M2，用于根据所述矩阵索引的行索引fr、列索引fc以及所述离散傅里叶变换的上采样率f，确定全息图像的傅里叶变换上采样频谱图fHU中局部的索引区域fH_in；所述索引区域fH_in为由行索引范围(fr-2)*f+1到fr*f，列索引范围(fc-2)*f+1到fc*f构成的大小为2f*2f的索引区域；
索引区域频率最大值位置确定模块M3，用于从所述索引区域fH_in中寻找频率最大值的位置，得到频率最大值位置的行索引fr_in和列索引fc_in；
+1级像频谱频率最大值位置精确获取模块M4，用于根据所述索引区域中频率最大值位置的行索引fr_in和列索引fc_in，计算全息图像的傅里叶变换上采样频谱图fHU中+1级像频谱中频率最大值位置的行索引fr_u和列索引fc_u；
参考光的角度余弦信息计算模块M5，用于根据所述+1级像频谱中频率最大值位置的行索引fr_u和列索引fc_u，计算参考光的角度余弦信息；
参考光光场分布估计模块M6，用于根据所述参考光的角度余弦信息，计算得到估计的参考光的光场分布信息。
本实施例中，所述索引区域确定模块M2包括：
矩阵构建单元，用于分别构建大小为2f*2f的矩阵EXP_HU和大小为M*N的矩阵EXP_H，矩阵中第i行第j列的元素分别表示为：
EXP_HU(i,j)＝exp(iπ((fr-2)*f+i+(fc-2)*f+j))；
Figure BDA0003399105030000101
其中，M表示全息图像H的矩阵总行数，N表示全息图像H的矩阵总列数；i表示矩阵中第i行，j表示矩阵中第j列；
所述矩阵构建单元，还用于构建大小为2f*M的矩阵RMatrixK，矩阵中第i行第j列的元素表示为：
Figure BDA0003399105030000102
所述矩阵构建单元，还用于构建大小为N*2f的矩阵CMatrixK，矩阵中第i行第j列的元素表示为：
Figure BDA0003399105030000111
索引区域获取单元，用于根据矩阵EXP_HU、矩阵EXP_H、矩阵RMatrixK和矩阵CMatrixK，计算得到索引区域fH_in：
fH_in＝EXP_HU⊙(RMatrixK*(H⊙EXP_H)*CMatrix)；
其中，H表示全息图像，⊙为矩阵对应元素相乘符号。
本实施例中，所述索引区域频率最大值位置确定模块M3包括：
+1级像频谱中频率最大值位置行索引计算单元，用于根据索引区域中频率最大值位置的fr_in，计算+1级像频谱中频率最大值位置的行索引fr_u，表示为：
fr_u＝(fr-2)*f+fr_in；
+1级像频谱频率最大值位置列索引计算单元，用于根据索引区域中频率最大值位置的列索引fc_in，计算+1级像频谱中频率最大值位置的列索引fc_u，表示为：
fc_u＝(fc-2)*f+fc_in。
除非另有定义，这里使用的所有术语(包括技术和科学术语)具有与本发明所属领域的普通技术人员共同理解的相同含义。还应当理解，诸如在通常字典里定义的那些术语应当被解释为具有与它们在相关技术的上下文中的含义相一致的含义，而不应用理想化或极度形式化的意义来解释，除非这里明确地这样定义。
上面是对本发明的说明，而不应被认为是对其的限制。尽管描述了本发明的若干示例性实施例，但本领域技术人员将容易地理解，在不背离本发明的新颖教学和优点的前提下可以对示例性实施例进行许多修改。因此，所有这些修改都意图包含在权利要求书所限定的本发明范围内。应当理解，上面是对本发明的说明，而不应被认为是限于所公开的特定实施例，并且对所公开的实施例以及其他实施例的修改意图包含在所附权利要求书的范围内。本发明由权利要求书及其等效物限定。

Claims (10)

1.一种离轴数字全息参考光估计方法，其特征在于，包括：
对全息图像进行离散傅里叶变换，并将离散傅里叶变换后的全息图像的零频分量转移到全息图像的频谱中心，对所述频谱中心两侧±1级像频谱中频率最大值的位置坐标进行定位，得到+1级像频谱中频率最大值的矩阵索引；所述±1级像频谱包括+1级像频谱和-1级像频谱，所述+1级像频谱表示物光波与参考光波共轭的乘积项的频谱，所述-1级像频谱表示物光波共轭与参考光波的乘积项的频谱；
根据所述矩阵索引的行索引fr、列索引fc以及所述离散傅里叶变换的上采样率f，确定全息图像的傅里叶变换上采样频谱图fHU中的局部的索引区域fH_in；所述索引区域fH_in为由行索引范围(fr-2)*f+1到fr*f，列索引范围(fc-2)*f+1到fc*f构成的大小为2f*2f的索引区域；
从所述索引区域fH_in中寻找频率最大值的位置，得到频率最大值位置的行索引fr_in和列索引fc_in；
根据所述索引区域中频率最大值位置的行索引fr_in和列索引fc_in，计算全息图像的傅里叶变换上采样频谱图fHU中+1级像频谱中频率最大值位置的行索引fr_u和列索引fc_u；
根据所述+1级像频谱中频率最大值位置的行索引fr_u和列索引fc_u，计算参考光的角度余弦信息；
根据所述参考光的角度余弦信息，计算得到估计的参考光的光场分布信息。
2.根据权利要求1所述的离轴数字全息参考光估计方法，其特征在于，采用Matlab软件中的FFT Shift函数，对离散傅里叶变换后的全息图像的零频分量进行FFT Shift变换，以将所述零频分量转移到全息图像的频谱中心位置。
3.根据权利要求1所述的离轴数字全息参考光估计方法，其特征在于，所述根据所述矩阵索引的行索引fr、列索引fc以及所述离散傅里叶变换的上采样率f，确定全息图像的傅里叶变换上采样频谱图fHU中的局部的索引区域fH_in，具体包括：
分别构建大小为2f*2f的矩阵EXP_HU和大小为M*N的矩阵EXP_H，矩阵中第i行第j列的元素分别表示为：
EXP_HU(i,j)＝exp(iπ((fr-2)*f+i+(fc-2)*f+j))；
Figure FDA0003399105020000021
其中，M表示全息图像H的矩阵总行数，N表示全息图像H的矩阵总列数；i表示矩阵中第i行，j表示矩阵中第j列；
构建大小为2f*M的矩阵RMatrixK，矩阵中第i行第j列的元素表示为：
Figure FDA0003399105020000022
构建大小为N*2f的矩阵CMatrixK，矩阵中第i行第j列的元素表示为：
Figure FDA0003399105020000023
根据矩阵EXP_HU、矩阵EXP_H、矩阵RMatrixK和矩阵CMatrixK，计算得到索引区域fH_in：
fH_in＝EXP_HU⊙(RMatrixK*(H⊙EXP_H)*CMatrix)；
其中，H表示全息图像，⊙为矩阵对应元素相乘符号。
4.根据权利要求1所述的离轴数字全息参考光估计方法，其特征在于，所述根据所述索引区域最大值位置的行索引fr_in和列索引fc_in，计算全息图像的傅里叶变换上采样频谱图fHU中+1级像频谱中频率最大值位置的行索引fr_u和列索引fc_u，具体包括：
根据索引区域中频率最大值位置的fr_in，计算+1级像频谱中频率最大值位置的行索引fr_u，表示为：
fr_u＝(fr-2)*f+fr_in；
根据索引区域中频率最大值位置的列索引fc_in，计算+1级像频谱中频率最大值位置的列索引fc_u，表示为：
fc_u＝(fc-2)*f+fc_in。
5.根据权利要求1所述的离轴数字全息参考光估计方法，其特征在于，所述根据所述+1级像频谱中频率最大值位置的行索引fr_u和列索引fc_u，计算参考光的角度余弦信息，具体包括：
根据所述+1级像频谱中频率最大值位置的行索引fr_u和列索引fc_u，计算+1级像频谱中频率最大值位置与频谱中心的距离；
根据所述距离，计算参考光的角度余弦信息。
6.根据权利要求5所述的离轴数字全息参考光估计方法，其特征在于，计算+1级像频谱中频率最大值位置与频谱中心的距离的公式为：
Figure FDA0003399105020000031
其中，M表示全息图像H的矩阵总行数，N表示全息图像H的矩阵总列数；fx_u和fy_u分别表示+1级像频谱中频率最大值位置与频谱中心在x轴上距离以及在y轴上的距离；
计算参考光的角度余弦信息的公式为：
Figure FDA0003399105020000032
其中，cosα表示参考光与x轴之间角度的余弦值，cosβ表示参考光与y轴之间角度的余弦值，λ表示波长，Δp表示探测器像素尺寸。
7.根据权利要求6所述的离轴数字全息参考光估计方法，其特征在于，计算估计的参考光的光场分布信息的公式为：
Figure FDA0003399105020000041
其中，ReferenceWave(x,y)表示估计的参考光的光场分布坐标，i为复数，i2＝-1。
8.一种离轴数字全息参考光估计系统，其特征在于，包括：
+1级像频谱频率最大值矩阵索引获取模块，用于对全息图像进行离散傅里叶变换，并将离散傅里叶变换后的全息图像的零频分量转移到全息图像的频谱中心，对所述频谱中心两侧±1级像频谱中频率最大值的位置坐标进行定位，得到+1级像频谱中频率最大值的矩阵索引；所述±1级像频谱包括+1级像频谱和-1级像频谱，所述+1级像频谱表示物光波与参考光波共轭的乘积项的频谱，所述-1级像频谱表示物光波共轭与参考光波的乘积项的频谱；
索引区域确定模块，用于根据所述矩阵索引的行索引fr、列索引fc以及所述离散傅里叶变换的上采样率f，确定全息图像的傅里叶变换上采样频谱图fHU中的局部的索引区域fH_in；所述索引区域fH_in为由行索引范围(fr-2)*f+1到fr*f，列索引范围(fc-2)*f+1到fc*f构成的大小为2f*2f的索引区域；
索引区域频率最大值位置确定模块，用于从所述索引区域fH_in中寻找频率最大值的位置，得到频率最大值位置的行索引fr_in和列索引fc_in；
+1级像频谱频率最大值位置精确获取模块，用于根据所述索引区域中频率最大值位置的行索引fr_in和列索引fc_in，计算全息图像的傅里叶变换上采样频谱图fHU中+1级像频谱中频率最大值位置的行索引fr_u和列索引fc_u；
参考光的角度余弦信息计算模块，用于根据所述+1级像频谱中频率最大值位置的行索引fr_u和列索引fc_u，计算参考光的角度余弦信息；
参考光光场分布估计模块，用于根据所述参考光的角度余弦信息，计算得到估计的参考光的光场分布信息。
9.根据权利要求8所述的离轴数字全息参考光估计系统，其特征在于，所述索引区域确定模块包括：
矩阵构建单元，用于分别构建大小为2f*2f的矩阵EXP_HU和大小为M*N的矩阵EXP_H，矩阵中第i行第j列的元素分别表示为：
EXP_HU(i,j)＝exp(iπ((fr-2)*f+i+(fc-2)*f+j))；
Figure FDA0003399105020000051
其中，M表示全息图像H的矩阵总行数，N表示全息图像H的矩阵总列数；i表示矩阵中第i行，j表示矩阵中第j列；
所述矩阵构建单元，还用于构建大小为2f*M的矩阵RMatrixK，矩阵中第i行第j列的元素表示为：
Figure FDA0003399105020000052
所述矩阵构建单元，还用于构建大小为N*2f的矩阵CMatrixK，矩阵中第i行第j列的元素表示为：
Figure FDA0003399105020000053
索引区域获取单元，用于根据矩阵EXP_HU、矩阵EXP_H、矩阵RMatrixK和矩阵CMatrixK，计算得到索引区域fH_in：
fH_in＝EXP_HU⊙(RMatrixK*(H⊙EXP_H)*CMatrix)；
其中，H表示全息图像，⊙为矩阵对应元素相乘符号。
10.根据权利要求8所述的离轴数字全息参考光估计系统，其特征在于，所述索引区域频率最大值位置确定模块包括：
+1级像频谱中频率最大值位置行索引计算单元，用于根据索引区域中频率最大值位置的fr_in，计算+1级像频谱中频率最大值位置的行索引fr_u，表示为：
fr_u＝(fr-2)*f+fr_in；
+1级像频谱频率最大值位置列索引计算单元，用于根据索引区域中频率最大值位置的列索引fc_in，计算+1级像频谱中频率最大值位置的列索引fc_u，表示为：
fc_u＝(fc-2)*f+fc_in。
CN202111490469.4A 2021-12-08 2021-12-08 一种离轴数字全息参考光估计方法及系统 Active CN114167707B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111490469.4A CN114167707B (zh) 2021-12-08 2021-12-08 一种离轴数字全息参考光估计方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111490469.4A CN114167707B (zh) 2021-12-08 2021-12-08 一种离轴数字全息参考光估计方法及系统

Publications (2)

Publication Number Publication Date
CN114167707A CN114167707A (zh) 2022-03-11
CN114167707B true CN114167707B (zh) 2023-05-16

Family

ID=80484343

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111490469.4A Active CN114167707B (zh) 2021-12-08 2021-12-08 一种离轴数字全息参考光估计方法及系统

Country Status (1)

Country Link
CN (1) CN114167707B (zh)

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4082431A (en) * 1975-04-22 1978-04-04 Minnesota Mining And Manufacturing Company Image processing system using incoherent radiation and spatial filter hologram
CN102317873A (zh) * 2008-12-10 2012-01-11 赫勒拉德有限责任公司 用于彩色运动全息术的系统和方法
CN103235495A (zh) * 2013-04-07 2013-08-07 哈尔滨工程大学 一种基于高斯核函数离轴数字全息零级像的抑制方法
CN104407507A (zh) * 2014-09-18 2015-03-11 河北工程大学 一种基于Hilbert变换的数字全息高精度位相重建方法
CN109116709A (zh) * 2018-08-14 2019-01-01 中国石油大学(华东) 一种弱离轴相移数字全息相移值抽取与物光恢复方法
CN110095968A (zh) * 2019-04-24 2019-08-06 清华大学深圳研究生院 一种离轴数字全息图像重建装置和方法及显微成像系统
CN111273533A (zh) * 2019-09-27 2020-06-12 中国工程物理研究院激光聚变研究中心 一种同轴数字全息自动聚焦方法及系统
EP3677969A1 (en) * 2017-08-30 2020-07-08 University of Hyogo Holographic imaging device and data processing method therefor

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4082431A (en) * 1975-04-22 1978-04-04 Minnesota Mining And Manufacturing Company Image processing system using incoherent radiation and spatial filter hologram
CN102317873A (zh) * 2008-12-10 2012-01-11 赫勒拉德有限责任公司 用于彩色运动全息术的系统和方法
CN103235495A (zh) * 2013-04-07 2013-08-07 哈尔滨工程大学 一种基于高斯核函数离轴数字全息零级像的抑制方法
CN104407507A (zh) * 2014-09-18 2015-03-11 河北工程大学 一种基于Hilbert变换的数字全息高精度位相重建方法
EP3677969A1 (en) * 2017-08-30 2020-07-08 University of Hyogo Holographic imaging device and data processing method therefor
CN109116709A (zh) * 2018-08-14 2019-01-01 中国石油大学(华东) 一种弱离轴相移数字全息相移值抽取与物光恢复方法
CN110095968A (zh) * 2019-04-24 2019-08-06 清华大学深圳研究生院 一种离轴数字全息图像重建装置和方法及显微成像系统
CN111273533A (zh) * 2019-09-27 2020-06-12 中国工程物理研究院激光聚变研究中心 一种同轴数字全息自动聚焦方法及系统

Non-Patent Citations (6)

* Cited by examiner, † Cited by third party
Title
Phase aberration compensation by spectrum centering in digital holographic microscopy;Huakun Cui等;《Optics Communications 》;第4152-4155页 *
Phase basedmethodforlocationofthecentersofsidebandsinspatial frequencydomaininoff-axisdigitalholographicmicrocopy;Yun Liu等;《OpticsandLasersinEngineering》;第115-124页 *
Spatial filtering for zero-order and twin-image elimination in digital off-axis holography;Etienne Cuche等;《APPLIED OPTICS》;第39卷(第23期);第4070-4075页 *
一种有效解决离轴数字全息相图倾斜畸变的数字 参考平面方法;李芳等;《物理学报》;第67卷(第9期);第094202-1至094202-9页 *
数字全息技术及散斑干涉技术在形变测量领域的发展及应用;闫浩等;《红外与激光工程》;第48卷(第6期);第0603010-1至0603010-13页 *
连续太赫兹波同轴数字全息相衬成像;王大勇等;《中国激光》;第41卷(第8期);第0809003-1至0809003-6页 *

Also Published As

Publication number Publication date
CN114167707A (zh) 2022-03-11

Similar Documents

Publication Publication Date Title
Wang et al. Quaternion polar complex exponential transform for invariant color image description
CN109901281B (zh) 一种穆勒显微偏振像差校准方法及装置
CN105157561A (zh) 基于降采样相位畸变主成分分析的快速相位畸变补偿方法
Jahid et al. Fast algorithm of 3D discrete image orthogonal moments computation based on 3D cuboid
CN114167707B (zh) 一种离轴数字全息参考光估计方法及系统
CN110675318B (zh) 一种基于主结构分离的稀疏表示图像超分辨率重建方法
CN114333983A (zh) 一种基于深度学习的磁共振波谱智能量化分析方法
Deng et al. Multiple frame splicing and degradation learning for hyperspectral imagery super-resolution
Auclair et al. Separation of dust emission from the cosmic infrared background in Herschel observations with wavelet phase harmonics
Hibbard et al. Computed three‐dimensional reconstruction of median‐eminence capillary modules: image alignment and correlation
CN111882532A (zh) 下肢x光影像中关键点提取方法
CN104574315B (zh) 基于光强传输矩阵的光学系统成像复原方法
WO2023142205A1 (zh) 一种InSAR时序相位的优化方法及装置
Sciammarella et al. Generalization of the Poincare sphere to process 2D displacement signals
McGaughey et al. Generating two-dimensional fractional Brownian motion using the fractional Gaussian process (FGp) algorithm
Armitage-Caplan et al. PReBeaM for Planck: A Polarized Regularized Beam Deconvolution Map-Making Method
Matei Design procedure for 2D Gaussian anisotropic filters
Hu et al. Single MR image super-resolution via mixed self-similarity attention network
CN113554555A (zh) 相位畸变补偿方法、装置、计算机设备和存储介质
Elfagrich Fringe pattern demodulation using the fan two-dimensional continuous wavelet transform conjugate to Shannon entropy
CN109856942A (zh) 一种基于不规则孔径正交多项式拟合的数字全息显微检测高阶像差校正方法
CN114237000B (zh) 一种离轴数字全息优化重建方法及系统
CN111121615B (zh) 一种基于希尔伯特变换的相移干涉条纹图批量选取方法
Matei Analytical Design Technique for the 2D IIR Differentiator
Liu et al. Testing the LSST Difference Image Analysis Pipeline Using Synthetic Source Injection Analysis

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