CN106530223B - 基于频域调制的快速Fourier鬼成像方法及系统 - Google Patents

基于频域调制的快速Fourier鬼成像方法及系统 Download PDF

Info

Publication number
CN106530223B
CN106530223B CN201611070198.6A CN201611070198A CN106530223B CN 106530223 B CN106530223 B CN 106530223B CN 201611070198 A CN201611070198 A CN 201611070198A CN 106530223 B CN106530223 B CN 106530223B
Authority
CN
China
Prior art keywords
frequency domain
binary
fourier
domain modulation
fringes
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
CN201611070198.6A
Other languages
English (en)
Other versions
CN106530223A (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.)
Tsinghua University
Original Assignee
Tsinghua University
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 Tsinghua University filed Critical Tsinghua University
Priority to CN201611070198.6A priority Critical patent/CN106530223B/zh
Publication of CN106530223A publication Critical patent/CN106530223A/zh
Application granted granted Critical
Publication of CN106530223B publication Critical patent/CN106530223B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T3/00Geometric image transformations in the plane of the image
    • G06T3/18Image warping, e.g. rearranging pixels individually

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Image Analysis (AREA)
  • Image Processing (AREA)

Abstract

本发明公开了一种基于频域调制的Fourier鬼成像方法及系统,其中,方法包括:通过频域调制得到二值条纹,以根据二值条纹生成灰度条纹;通过灰度条纹采集图像频谱的幅值和相位;根据图像频谱的幅值和相位恢复出频谱;通过Fourier反变换得到成像图像。该成像方法可以通过频域调制的手段提升Fourier鬼成像的速度,从而快速实现高精度鬼成像,提高成像的实用性。

Description

基于频域调制的快速Fourier鬼成像方法及系统
技术领域
本发明涉及计算摄像学技术领域,特别涉及一种基于频域调制的快速Fourier鬼成像方法及系统。
背景技术
相关技术,Fourier鬼成像通过直接采集图像频谱的幅值和相位恢复出频谱,然后再通过Fourier反变换得到图像。Fourier鬼成像不同于传统单像素成像的地方在于,该工作中使用4步相位平移正弦信号的调制图案来完成对成像场景的傅里叶频谱的采集,而非传统意义上采用随机图案来完成成像过程。由于同一正弦信号在不同相移下在具有相同的总强度,4种不同的相移相互之间比较即可去除频率远远高于或者远远低于调制频率的环境光噪声。但是,由于仪器生成正弦信号远远慢于生成二值信号,导致Fourier鬼成像速度较慢。
发明内容
本发明旨在至少在一定程度上解决相关技术中的技术问题之一。
为此,本发明的一个目的在于提出一种基于频域调制的Fourier鬼成像方法,该方法可以提升Fourier鬼成像的速度,从而快速实现高精度鬼成像。
本发明的另一个目的在于提出一种基于频域调制的Fourier鬼成像系统。
为达到上述目的,本发明一方面实施例提出了一种基于频域调制的Fourier鬼成像方法,包括以下步骤:通过频域调制得到二值条纹,以根据所述二值条纹生成灰度条纹;通过所述灰度条纹采集图像频谱的幅值和相位;根据所述图像频谱的幅值和相位恢复出频谱;通过Fourier反变换得到成像图像。
本发明实施例的基于频域调制的Fourier鬼成像方法,可以通过频域调制的手段提升Fourier鬼成像的速度,从而快速实现高精度鬼成像,通过频域调制的方法使得在Fourier鬼成像中采用的正弦条纹能由二值条纹代替,提高成像的实用性和可靠性。
另外,根据本发明上述实施例的基于频域调制的Fourier鬼成像方法还可以具有以下附加的技术特征:
进一步地,在本发明的一个实施例中,滤除所述二值条纹中高次频率点,以得到对应于所述二值条纹的正弦条纹。
进一步地,在本发明的一个实施例中,所述对应于所述二值条纹的正弦条纹的获取步骤包括:获取二值图像;通过在频域截取多个像素为半径的圆得到所述正弦条纹。
进一步地,在本发明的一个实施例中,根据垂直方向的灰度值判断所述正弦条纹是否满足预设条件。
进一步地,在本发明的一个实施例中,通过空间光调制器DMD频域调制所述二值条纹以得到所述正弦条纹。
为达到上述目的,本发明另一方面实施例提出了一种基于频域调制的Fourier鬼成像系统,包括:调制模块,用于通过频域调制得到二值条纹,以根据所述二值条纹生成灰度条纹;采集模块,用于通过所述灰度条纹采集图像频谱的幅值和相位;恢复模块,用于根据所述图像频谱的幅值和相位恢复出频谱;成像模块,用于通过Fourier反变换得到成像图像。
本发明实施例的基于频域调制的Fourier鬼成像系统,可以通过频域调制的手段提升Fourier鬼成像的速度,从而快速实现高精度鬼成像,通过频域调制的方法使得在Fourier鬼成像中采用的正弦条纹能由二值条纹代替,提高成像的实用性和可靠性。
另外,根据本发明上述实施例的基于频域调制的Fourier鬼成像系统还可以具有以下附加的技术特征:
进一步地,在本发明的一个实施例中,还包括:滤除模块,用于滤除所述二值条纹中高次频率点,以得到对应于所述二值条纹的正弦条纹。
进一步地,在本发明的一个实施例中,所述滤除模块还用于获取二值图像,并且通过在频域截取多个像素为半径的圆得到所述正弦条纹。
进一步地,在本发明的一个实施例中,所述滤除模块还用于根据垂直方向的灰度值判断所述正弦条纹是否满足预设条件。
进一步地,在本发明的一个实施例中,通过空间光调制器DMD频域调制所述二值条纹以得到所述正弦条纹。
本发明附加的方面和优点将在下面的描述中部分给出,部分将从下面的描述中变得明显,或通过本发明的实践了解到。
附图说明
本发明上述的和/或附加的方面和优点从下面结合附图对实施例的描述中将变得明显和容易理解,其中:
图1为根据本发明实施例的基于频域调制的Fourier鬼成像方法的流程图;
图2为根据本发明一个实施例的正弦条纹及其频域的示意图;
图3为根据本发明一个实施例的二值条纹及其频域的示意图;
图4为根据本发明一个实施例的频域选通带以及选通后反变换得到的条纹示意图;
图5为根据本发明一个实施例的选通后反变换得到的条纹延垂直于条纹方向的灰度值示意图;
图6为根据本发明一个实施例的快速鬼成像设计光路示意图;
图7为根据本发明一个实施例的样本图片及其Fourier变换的示意图;
图8为根据本发明一个实施例的重建图片及其Fourier变换的示意图;
图9为本发明实施例的基于频域调制的Fourier鬼成像系统的结构示意图。
具体实施方式
下面详细描述本发明的实施例,所述实施例的示例在附图中示出,其中自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。下面通过参考附图描述的实施例是示例性的,旨在用于解释本发明,而不能理解为对本发明的限制。
下面参照附图描述根据本发明实施例提出的基于频域调制的Fourier鬼成像方法及系统,首先将参照附图描述根据本发明实施例提出的基于频域调制的Fourier鬼成像方法。
图1是本发明实施例的基于频域调制的Fourier鬼成像方法的流程图。
如图1所示,该基于频域调制的Fourier鬼成像方法包括以下步骤:
在步骤S101中,通过频域调制得到二值条纹,以根据二值条纹生成灰度条纹。
其中,在本发明的一个实施例中,滤除二值条纹中高次频率点,以得到对应于二值条纹的正弦条纹。
可以理解的是,如图2和图3所示,观察到二维平面上的正弦条纹,其Fourier平面上只有三个点,分别是中心的、代表均值的点,以及对称于中心的代表频率的两个点。同样级次的二值条纹,其频谱图相比正弦条纹而言多了很多点,这些点是高频分量,对应着高频的正弦条纹。滤除高次频率的点,就可以得到相应的正弦条纹。
在步骤S102中,通过灰度条纹采集图像频谱的幅值和相位。
在步骤S103中,根据图像频谱的幅值和相位恢复出频谱。
在步骤S104中,通过Fourier反变换得到成像图像。
也就是说,在本发明的实施例中,可以得到由二值条纹生成的灰度条纹,随后进行Fourier鬼成像的一般步骤,从而得到成像图像。
其中,在本发明的一个实施例中,对应于二值条纹的正弦条纹的获取步骤包括:获取二值图像;通过在频域截取多个像素为半径的圆得到正弦条纹。
进一步地,在本发明的一个实施例中,根据垂直方向的灰度值判断正弦条纹是否满足预设条件。
需要说明的是,预设条件可以根据实际情况进行设置,在此不作具体限制。
进一步地,在本发明的一个实施例中,通过空间光调制器DMD频域调制二值条纹以得到正弦条纹。
可以理解的是,通过Matlab仿真,发现通过在频域上截取一个圆形区域,将低频部分选通,高频部分滤除,可以达到很好的正弦化效果。如图4和图5所示,有4条亮纹的二值图像,通过在频域截取5个像素为半径的圆,即可得到同级次的正弦条纹。通过测定其垂直方向的灰度值,说明得到的条纹是理想的正弦条纹。
由此可得,通过数字微镜阵列频域调制二值条纹从而可以得到正弦条纹。
举例而言,如图6所示,数字微镜阵列1负责产生二值条纹,经过透镜1后在其焦面上形成Fourier变换,随后由数字微镜阵列2截取低频部分,再经过透镜2做Fourier反变换后得到正弦条纹。随后正弦条纹打在样本上,最后由单像素接收,完成整个采集过程。
进一步地,下面由上述光路设计全系统的仿真。选择64x64的样本图样,然后利用上述得到的由二值条纹生成的灰度条纹,随后进行Fourier鬼成像的一般步骤,可得到如图7和图8所示的重建图像,由此可知,本发明实施例的方法虽然会降低一些分辨率,但是确实是有效的。
根据本发明实施例的基于频域调制的Fourier鬼成像方法,可以通过频域调制的手段提升Fourier鬼成像的速度,从而快速实现高精度鬼成像,通过频域调制的方法使得在Fourier鬼成像中采用的正弦条纹能由二值条纹代替,提高成像的实用性和可靠性。
其次参照附图描述根据本发明实施例提出的基于频域调制的Fourier鬼成像系统。
图9是本发明实施例的基于频域调制的Fourier鬼成像系统的结构示意图。
如图9所示,该基于频域调制的Fourier鬼成像系统10包括:调制模块100、采集模块200、恢复模块300和成像模块400。
其中,调制模块100用于通过频域调制得到二值条纹,以根据二值条纹生成灰度条纹。采集模块200用于通过灰度条纹采集图像频谱的幅值和相位。恢复模块300用于根据图像频谱的幅值和相位恢复出频谱。成像模块400用于通过Fourier反变换得到成像图像。本发明实施例的成像系统10可以通过频域调制的手段提升Fourier鬼成像的速度,从而快速实现高精度鬼成像,提高成像的实用性。
进一步地,在本发明的一个实施例中,本发明实施例的成像系统还包括:滤除模块。其中,滤除模块用于滤除二值条纹中高次频率点,以得到对应于二值条纹的正弦条纹。
进一步地,在本发明的一个实施例中,滤除模块还用于获取二值图像,并且通过在频域截取多个像素为半径的圆得到正弦条纹。
进一步地,在本发明的一个实施例中,滤除模块还用于根据垂直方向的灰度值判断正弦条纹是否满足预设条件。
进一步地,在本发明的一个实施例中,通过空间光调制器DMD频域调制二值条纹以得到正弦条纹。
需要说明的是,前述对基于频域调制的Fourier鬼成像方法实施例的解释说明也适用于该实施例的基于频域调制的Fourier鬼成像系统,此处不再赘述。
根据本发明实施例的基于频域调制的Fourier鬼成像系统,可以通过频域调制的手段提升Fourier鬼成像的速度,从而快速实现高精度鬼成像,通过频域调制的方法使得在Fourier鬼成像中采用的正弦条纹能由二值条纹代替,提高成像的实用性和可靠性。
在本发明的描述中,需要理解的是,术语“中心”、“纵向”、“横向”、“长度”、“宽度”、“厚度”、“上”、“下”、“前”、“后”、“左”、“右”、“竖直”、“水平”、“顶”、“底”“内”、“外”、“顺时针”、“逆时针”、“轴向”、“径向”、“周向”等指示的方位或位置关系为基于附图所示的方位或位置关系,仅是为了便于描述本发明和简化描述,而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作,因此不能理解为对本发明的限制。
此外,术语“第一”、“第二”仅用于描述目的,而不能理解为指示或暗示相对重要性或者隐含指明所指示的技术特征的数量。由此,限定有“第一”、“第二”的特征可以明示或者隐含地包括至少一个该特征。在本发明的描述中,“多个”的含义是至少两个,例如两个,三个等,除非另有明确具体的限定。
在本发明中,除非另有明确的规定和限定,术语“安装”、“相连”、“连接”、“固定”等术语应做广义理解,例如,可以是固定连接,也可以是可拆卸连接,或成一体;可以是机械连接,也可以是电连接;可以是直接相连,也可以通过中间媒介间接相连,可以是两个元件内部的连通或两个元件的相互作用关系,除非另有明确的限定。对于本领域的普通技术人员而言,可以根据具体情况理解上述术语在本发明中的具体含义。
在本发明中,除非另有明确的规定和限定,第一特征在第二特征“上”或“下”可以是第一和第二特征直接接触,或第一和第二特征通过中间媒介间接接触。而且,第一特征在第二特征“之上”、“上方”和“上面”可是第一特征在第二特征正上方或斜上方,或仅仅表示第一特征水平高度高于第二特征。第一特征在第二特征“之下”、“下方”和“下面”可以是第一特征在第二特征正下方或斜下方,或仅仅表示第一特征水平高度小于第二特征。
在本说明书的描述中,参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中,对上述术语的示意性表述不必须针对的是相同的实施例或示例。而且,描述的具体特征、结构、材料或者特点可以在任一个或多个实施例或示例中以合适的方式结合。此外,在不相互矛盾的情况下,本领域的技术人员可以将本说明书中描述的不同实施例或示例以及不同实施例或示例的特征进行结合和组合。
尽管上面已经示出和描述了本发明的实施例,可以理解的是,上述实施例是示例性的,不能理解为对本发明的限制,本领域的普通技术人员在本发明的范围内可以对上述实施例进行变化、修改、替换和变型。

Claims (10)

1.一种基于频域调制的Fourier鬼成像方法,其特征在于,包括以下步骤:
通过频域调制得到二值条纹,以根据所述二值条纹生成灰度条纹;
通过所述灰度条纹采集图像频谱的幅值和相位;
根据所述图像频谱的幅值和相位恢复出频谱;以及
通过Fourier反变换得到成像图像。
2.根据权利要求1所述的基于频域调制的Fourier鬼成像方法,其特征在于,滤除所述二值条纹中高次频率点,以得到对应于所述二值条纹的正弦条纹。
3.根据权利要求2所述的基于频域调制的Fourier鬼成像方法,其特征在于,所述对应于所述二值条纹的正弦条纹的获取步骤包括:
获取二值图像;
通过在频域截取多个像素为半径的圆得到所述正弦条纹。
4.根据权利要求3所述的基于频域调制的Fourier鬼成像方法,其特征在于,根据垂直方向的灰度值判断所述正弦条纹是否满足预设条件。
5.根据权利要求2所述的基于频域调制的Fourier鬼成像方法,其特征在于,通过空间光调制器DMD频域调制所述二值条纹以得到所述正弦条纹。
6.一种基于频域调制的Fourier鬼成像系统,其特征在于,包括:
调制模块,用于通过频域调制得到二值条纹,以根据所述二值条纹生成灰度条纹;
采集模块,用于通过所述灰度条纹采集图像频谱的幅值和相位;
恢复模块,用于根据所述图像频谱的幅值和相位恢复出频谱;以及
成像模块,用于通过Fourier反变换得到成像图像。
7.根据权利要求6所述的基于频域调制的Fourier鬼成像系统,其特征在于,还包括:
滤除模块,用于滤除所述二值条纹中高次频率点,以得到对应于所述二值条纹的正弦条纹。
8.根据权利要求7所述的基于频域调制的Fourier鬼成像系统,其特征在于,所述滤除模块还用于获取二值图像,并且通过在频域截取多个像素为半径的圆得到所述正弦条纹。
9.根据权利要求8所述的基于频域调制的Fourier鬼成像系统,其特征在于,所述滤除模块还用于根据垂直方向的灰度值判断所述正弦条纹是否满足预设条件。
10.根据权利要求7所述的基于频域调制的Fourier鬼成像系统,其特征在于,通过空间光调制器DMD频域调制所述二值条纹以得到所述正弦条纹。
CN201611070198.6A 2016-11-28 2016-11-28 基于频域调制的快速Fourier鬼成像方法及系统 Active CN106530223B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201611070198.6A CN106530223B (zh) 2016-11-28 2016-11-28 基于频域调制的快速Fourier鬼成像方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201611070198.6A CN106530223B (zh) 2016-11-28 2016-11-28 基于频域调制的快速Fourier鬼成像方法及系统

Publications (2)

Publication Number Publication Date
CN106530223A CN106530223A (zh) 2017-03-22
CN106530223B true CN106530223B (zh) 2020-01-10

Family

ID=58355107

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201611070198.6A Active CN106530223B (zh) 2016-11-28 2016-11-28 基于频域调制的快速Fourier鬼成像方法及系统

Country Status (1)

Country Link
CN (1) CN106530223B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107241550B (zh) * 2017-06-26 2019-10-29 清华大学 基于时分复用的快速单像素频域成像方法
CN109239466B (zh) * 2018-06-15 2020-12-04 北京理工大学 一种基于经典鬼成像的网络测试系统
CN110132175B (zh) * 2019-05-30 2021-04-20 北京理工大学 基于幅度调制的单像素相位成像方法和装置
CN110809102B (zh) * 2019-10-11 2020-10-30 北京理工大学 基于二值调制的成像加速方法及装置
CN111474145A (zh) * 2020-03-19 2020-07-31 清华大学 单像素荧光和相位成像系统及方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101419668A (zh) * 2008-12-11 2009-04-29 潘玉利 基于二维傅立叶变换的水泥路面图像刻槽去除的方法
WO2011087083A1 (ja) * 2010-01-14 2011-07-21 国立大学法人富山大学 データ処理方法、データ処理装置、及びデータ処理プログラム
CN105023027A (zh) * 2014-04-18 2015-11-04 大连恒锐科技股份有限公司 基于多次反馈机制的现场鞋底痕迹花纹图像检索方法
CN105023247A (zh) * 2015-06-24 2015-11-04 南京理工大学 抑制单臂鬼成像离焦模糊的二次成像方法
CN105738073A (zh) * 2016-02-03 2016-07-06 中国科学院国家空间科学中心 一种在空间频率域进行像素响应函数测量的方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9454809B2 (en) * 2012-09-25 2016-09-27 The Board Of Trustees Of The University Of Illinois Phase derivative microscopy module having specified amplitude mask

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101419668A (zh) * 2008-12-11 2009-04-29 潘玉利 基于二维傅立叶变换的水泥路面图像刻槽去除的方法
WO2011087083A1 (ja) * 2010-01-14 2011-07-21 国立大学法人富山大学 データ処理方法、データ処理装置、及びデータ処理プログラム
CN105023027A (zh) * 2014-04-18 2015-11-04 大连恒锐科技股份有限公司 基于多次反馈机制的现场鞋底痕迹花纹图像检索方法
CN105023247A (zh) * 2015-06-24 2015-11-04 南京理工大学 抑制单臂鬼成像离焦模糊的二次成像方法
CN105738073A (zh) * 2016-02-03 2016-07-06 中国科学院国家空间科学中心 一种在空间频率域进行像素响应函数测量的方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
4-f光学架构频域纯相位调制压缩成像;张成 等;《北京邮电大学学报》;20140630;第37卷(第3期);第13-17页 *
基于多散斑图和三维计算的应用型鬼成像方案研究;仲亚军;《中国优秀硕士学位论文全文数据库 信息科技辑》;20160515;第I138-1225页 *
基于混合散斑图的压缩鬼成像方法研究;周成 等;《光学学报》;20160930;第36卷(第9期);第99-105页 *
多尺度压缩感知鬼成像系统;张娜 等;《量子电子学报》;20150731;第32卷(第4期);第385-390页 *

Also Published As

Publication number Publication date
CN106530223A (zh) 2017-03-22

Similar Documents

Publication Publication Date Title
CN106530223B (zh) 基于频域调制的快速Fourier鬼成像方法及系统
CN108370406B (zh) 摄像装置
US8773550B2 (en) Range measurement using multiple coded apertures
CN110057319B (zh) 一种大范围反射率变化物体的高速三维测量方法
US9471728B2 (en) Metrics for designing a plenoptic imaging system
KR20110124701A (ko) 감쇠 패턴을 포함하는 마스크를 이용한 광 필드 영상 처리 장치 및 방법
US10206580B2 (en) Full-field OCT system using wavelength-tunable laser and three-dimensional image correction method
CN109708842B (zh) 一种基于单像素成像的相机镜头点扩散函数测量方法
JP2011523057A (ja) 構造化された照射および均一な照射の両方を用いて光学的に切片化された画像を生成するためのシステムおよび方法
JP2007256171A (ja) ミリ波画像処理装置及びミリ波画像処理方法
CN110264540B (zh) 一种并行单像素成像方法
CN108873565B (zh) 摄像装置
JP6716375B2 (ja) 回転速度計測装置
Mercier et al. Combining visibilities from the giant meterwave radio telescope and the Nancay radio heliograph-High dynamic range snapshot images of the solar corona at 327 MHz
JP2003141529A (ja) 画像再構成法及び画像再構成装置
US10887504B2 (en) Distance measurement device and distance measurement method
US20220042923A1 (en) Structured illumination optical inspection platform for transparent materials
CN115546285B (zh) 基于点扩散函数解算的大景深条纹投影三维测量方法
CN113030962A (zh) 机载太赫兹合成孔径雷达及成像方法
JP2538435B2 (ja) 縞位相分布解析方法および縞位相分布解析装置
CN108445458B (zh) 一种合成孔径雷达轨道误差消除方法和装置
JP2001174235A (ja) フーリエ変換を用いた縞解析方法
JP6807281B2 (ja) 撮像装置、撮像システム、及び撮像方法
CN104463863A (zh) 基于时间外差投影的运动干涉场的标定方法与系统
JP6748504B2 (ja) Mtf測定用チャート

Legal Events

Date Code Title Description
C06 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