CN110849841A - 基于暗场成像的反射式太赫兹成像方法及装置 - Google Patents
基于暗场成像的反射式太赫兹成像方法及装置 Download PDFInfo
- Publication number
- CN110849841A CN110849841A CN201911152543.4A CN201911152543A CN110849841A CN 110849841 A CN110849841 A CN 110849841A CN 201911152543 A CN201911152543 A CN 201911152543A CN 110849841 A CN110849841 A CN 110849841A
- Authority
- CN
- China
- Prior art keywords
- lens
- imaging
- terahertz
- target
- reflected
- 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
Links
- 238000003384 imaging method Methods 0.000 title claims abstract description 123
- 238000001228 spectrum Methods 0.000 claims abstract description 8
- 230000003595 spectral effect Effects 0.000 claims description 14
- 230000000694 effects Effects 0.000 claims description 4
- 230000005540 biological transmission Effects 0.000 claims description 3
- 230000000903 blocking effect Effects 0.000 claims description 3
- 230000005855 radiation Effects 0.000 claims description 3
- 238000005516 engineering process Methods 0.000 abstract description 3
- 238000004088 simulation Methods 0.000 description 4
- 238000012634 optical imaging Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 238000002310 reflectometry Methods 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 239000004615 ingredient Substances 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000000149 penetrating effect Effects 0.000 description 1
- 230000001902 propagating effect Effects 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 238000002834 transmittance Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/25—Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
- G01N21/31—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
- G01N21/35—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light
- G01N21/3581—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light using far infrared light; using Terahertz radiation
- G01N21/3586—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light using far infrared light; using Terahertz radiation by Terahertz time domain spectroscopy [THz-TDS]
Landscapes
- Physics & Mathematics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Health & Medical Sciences (AREA)
- Toxicology (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
一种基于暗场成像的反射式太赫兹成像方法及装置,该方法包括准直的太赫兹波束经分束镜反射后得到反射波束;反射波束辐射到成像目标并经成像目标反射后得到第一波束;第一波束传播距离f后经过第一透镜调制后得到第二波束;第二波束传播距离f后经过频谱面的波束挡板后得到第三波束;第三波束传播距离f后经过第二透镜调制后得到第四波束;第四波束经第二透镜调制后传播距离f被太赫兹相机接收成像。本发明将暗场成像技术引入到太赫兹成像中,结合暗场成像可实现高对比度成像的特点,从而实现太赫兹高对比度成像。
Description
技术领域
本发明涉及太赫兹准光成像技术领域,尤其涉及一种基于暗场成像的反射式太赫兹成像方法及装置。
背景技术
太赫兹波对生物组织具有一定的穿透能力,且单光子能量低,不会引起生物组织电离,在生物医学成像上非常安全,适用于生物医学成像。但是,太赫兹生物医学图像往往对比度很低,对于生物组织来说,由于其含水量较高使得太赫兹反射率很低,太赫兹对生物组织的透射率远小于反射率,因此,亟需研究反射式太赫兹高对比度成像方法。暗场成像技术被广泛应用于光学成像中,相对于普通成像方法,暗场成像通过提取成像目标的高频信息从而实现成像对比度增强。将暗场成像技术引入到反射式太赫兹成像领域中,可以实现太赫兹高对比度成像。因此,研究反射式太赫兹暗场成像方法具有重要的研究价值。
发明内容
有鉴于此,本发明的主要目的之一在于提出一种基于暗场成像的反射式太赫兹成像方法及装置,以期至少部分地解决上述技术问题中的至少之一。
为了实现上述目的,作为本发明的一个方面,提供了一种基于暗场成像的反射式太赫兹成像方法,包括:
准直的太赫兹波束经分束镜反射后得到反射波束;
反射波束辐射到成像目标并经成像目标反射后得到第一波束;
第一波束传播距离f后经过第一透镜调制后得到第二波束;
第二波束传播距离f后经过频谱面的波束挡板后得到第三波束;
第三波束传播距离f后经过第二透镜调制后得到第四波束;
第四波束经第二透镜调制后传播距离f被太赫兹相机接收成像。
作为本发明的另一个方面,还提供了一种成像装置,用于执行如上所述的反射式太赫兹成像方法,包括:
分束镜,用于将太赫兹波束反射到成像目标上并透射经目标反射的太赫兹波;
第一透镜,用于调制经过成像目标反射后的第一波束;
波束挡板,用于遮挡经过第一透镜的第二波束中心;
第二透镜,用于调制经过波束挡的第三波束;以及
太赫兹相机,用于接收经过第二透镜后的第四波束并成像。
基于上述技术方案可知,本发明的基于暗场成像的反射式太赫兹成像方法及装置相对于现有技术至少具有以下优势之一:
1、基于暗场成像的太赫兹高对比度成像方法将暗场成像技术引入到太赫兹成像中,结合暗场成像可实现高对比度成像的特点,从而实现太赫兹高对比度成像;
2、采用反射式成像相比透射式成像,可以获得更高的信噪比;
3、利用太赫兹相机进行接收,可以实现对成像目标的实时成像。
附图说明
图1为本发明实施例的基于暗场成像的反射式太赫兹高对比度成像方法的步骤示意图;
图2为本发明实施例的基于暗场成像的反射式太赫兹高对比度成像装置的光路结构示意图;
图3为本发明实施例的基于暗场成像的反射式太赫兹高对比度成像方法的仿真成像结果。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明白,以下结合具体实施例,并参照附图,对本发明作进一步的详细说明。
目前,传统的太赫兹生物医学成像方法大多成像对比度低,成像质量差。本发明旨在将光学成像中的暗场成像引入到反射式太赫兹成像中,以提高成像对比度。相对于传统成像方法,暗场成像在频谱面加入波束挡板,用于滤除目标的低频信息,从而提取出目标的高频信息。由于低频信息代表背景,高频信息代表目标的细节,因此,暗场成像可以抑制背景,提取出目标细节信息,从而实现高对比度成像。成像系统中利用太赫兹相机作为接收,可以实现对目标实时高对比度成像。由于太赫兹波长相对于光波要长的多,容易发生较强的衍射效应,因此实验系统元件尺寸和间距应远大于波长。
本发明公开了一种基于暗场成像的反射式太赫兹成像方法,包括:
准直的太赫兹波束经分束镜反射后得到反射波束;
反射波束辐射到成像目标并经成像目标反射后得到第一波束;
第一波束传播距离f后经过第一透镜调制后得到第二波束;
第二波束传播距离f后经过频谱面的波束挡板后得到第三波束;
第三波束传播距离f后经过第二透镜调制后得到第四波束;
第四波束经第二透镜调制后传播距离f被太赫兹相机接收成像。
其中,所述第一波束传播距离f过程中经所述分束镜透射后到达第一透镜;
其中,所述第二波束中心被波束挡板遮挡。
其中,所述第一透镜和第二透镜均为凸透镜,第一透镜和第二透镜的焦距均为f。
其中,经所述成像目标反射后的第一波束起始位置的场分布E1(x1,y1)为:
E1(x1,y1)=E0(x1,y1)t(x1,y1);
其中,t(x1,y1)为目标的反射函数,E0(x1,y1)为辐射到成像目标的场分布,(x1,y1)为目标平面的坐标。
其中,所述到达第一透镜的场分布满足:
F{E2(x2,y2)}=F{E1(x1,y1)}H(fx,fy);
其中,第二波束到达所述频谱面的场分布E3(x3,y3)为:
其中,所述第二波束到达频谱面的场分布展开为:
其中,第一项为没有目标时频谱面上的场分布,即低频成分;第二项是由衍射效应引起的,代表目标的细节信息,即高频成分;
由于波束挡板遮挡了波束的中心,即低频成分,保留了高频成分;波束未被遮挡的部分场分布为:
其中,所述第三波束到达第二透镜的场分布E4(x4,y4)为:
F{E4(x4,y4)}=F{E′3(x3,y3)}H(fx,fy);
其中,fx=x4/λf,fy=y4/λf,(x4,y4)为第二透镜所在平面的坐标。
其中,所述第四波束到达太赫兹相机表面的场分布E5(x5,y5)为:
其中,(x5,y5)为太赫兹相机所在平面的坐标;
其中,所述太赫兹相机接收的为光强,其分布I为:
本发明该公开了一种成像装置,用于执行如上所述的反射式太赫兹成像方法,包括:
分束镜,用于将太赫兹波束反射到成像目标上并透射经目标反射的太赫兹波;
第一透镜,用于调制经过成像目标反射后的第一波束;
波束挡板,用于遮挡经过第一透镜的第二波束中心;
第二透镜,用于调制经过波束挡的第三波束;以及
太赫兹相机,用于接收经过第二透镜后的第四波束并成像;
其中,所述波束挡板设置在第二透射波束中心位置;
其中,所述成像目标、第一透镜、螺旋相位板、第二透镜、太赫兹相机等距且平行设置;
其中,所述分束镜设置在成像目标与第一透镜之间;
其中,所述太赫兹波束经过分束镜反射后的波束传播方向与成像目标表面垂直。
以下通过具体实施例结合附图对本发明的技术方案做进一步阐述说明。需要注意的是,下述的具体实施例仅是作为举例说明,本发明的保护范围并不限于此。
如图1和图2所示,本实施例的基于暗场成像的太赫兹高对比度成像方法,包括以下步骤:
S1、准直的太赫兹波束经分束镜反射后辐射到成像目标。
所述辐射到成像目标的场分布记为E0(x1,y1),(x1,y1)为目标平面的坐标。
S2、所述反射波束经成像目标反射后得到第一波束;
所述经成像目标反射后的场分布为:
E1(x1,y1)=E0(x1,y1)t(x1,y1);
其中,t(x1,y1)为目标的反射函数。
S3、所述第一波束传播距离f后到达第一透镜。
所述到达第一透镜的场分布满足:
F{E2(x2,y2)}=F{E1(x1,y1)}H(fx,fy)
S4、所述第一波束经第一透镜调制后传播距离f到达频谱面。
第一波束经过第一透镜后形成第二波束;到达频谱面的场分布为:
S5、所述第二波束中心被位于频谱面的波束挡板遮挡,第二波束其余部分未被遮挡。
到达频谱面的场分布可以改写为:
其中,第一项为没有目标时频谱面上的场分布,即低频成分;第二项是由衍射效应引起的,代表目标的细节信息,即高频成分。波束挡板遮挡了波束的中心,即低频成分,保留了高频成分。波束未被遮挡的部分场分布为:
S6、所述第二波束未被遮挡的部分传播距离f后到达第二透镜;第二波束经过波束挡板后形成第三波束;
第三波束到达第二透镜的场分布为:
F{E4(x4,y4)}=F{E′3(x3,y3)}H(fx,fy)
其中,fx=x4/λf,fy=y4/λf,(x4,y4)为第二透镜所在平面的坐标。
S7、所述第三波束波束经第二透镜调制后传播距离f被太赫兹相机接收。
到达太赫兹相机表面的场分布为:
其中,(x5,y5)为太赫兹相机所在平面的坐标。
所述太赫兹相机接收的为光强,其分布为:
图3中(a)图、(c)图、(e)图分别为仿真目标,目标图像中的字母“THz”的反射率均为1,背景的反射率分别为0.6、0.9、0.99,由对比度计算公式:(最大值-最小值)/(最大值+最小值)可知,目标图像的对比度分别25%、5.26%、0.5%,即对比度依此降低,当对比度为0.5%时,人眼无法识别出目标;(b)图、(d)图、(f)图、为对应的(a)图、(c)图、(e)图太赫兹螺旋相衬成像仿真结果,可以看出,仿真结果均能较好地提取出目标的高频信息,三个图像的对比度均为100%,从而实现对低对比度目标的高对比度成像。综上,本发明的基于螺旋相衬成像的太赫兹高对比度成像方法,能实现高对比度成像,采用太赫兹相机作为接收,可以实现实时成像。
本实施例中第一透镜和第二透镜均为凸透镜,第一透镜和第二透镜的焦距均为f。第一透镜的目的是对含有目标信息的太赫兹波束进行傅里叶变换到频谱面然后在频谱面滤波。暗场成像在频谱面加中心挡板以滤除低频分量,凸显高频分量。第二透镜的目的是对在频谱面滤波后的波束进行傅里叶变换,还原出目标滤波后的图像。
综上,本发明的基于暗场成像的太赫兹高对比度成像方法,能实现高对比度成像,采用太赫兹相机作为接收,可以实现实时成像。
除非有所知名为相反之意,本说明书及所附权利要求中的数值参数是近似值,能够根据通过本公开的内容所得的所需特性改变。具体而言,所有使用于说明书及权利要求中表示组成的含量、反应条件等等的数字,应理解为在所有情况中是受到“约”的用语所修饰。一般情况下,其表达的含义是指包含由特定数量在一些实施例中±10%的变化、在一些实施例中±5%的变化、在一些实施例中±1%的变化、在一些实施例中±0.5%的变化。
再者,“包含”不排除存在未列在权利要求中的元件或步骤。位于元件之前的“一”或“一个”不排除存在多个这样的元件。
说明书与权利要求中所使用的序数例如“第一”、“第二”、“第三”等的用词,以修饰相应的元件,其本身并不意味着该元件有任何的序数,也不代表某一元件与另一元件的顺序、或是制造方法上的顺序,该些序数的使用仅用来使具有某命名的一元件得以和另一具有相同命名的元件能做出清楚区分。
以上所述的具体实施例,对本发明的目的、技术方案和有益效果进行了进一步详细说明,应理解的是,以上所述仅为本发明的具体实施例而已,并不用于限制本发明,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (10)
1.一种基于暗场成像的反射式太赫兹成像方法,包括:
准直的太赫兹波束经分束镜反射后得到反射波束;
反射波束辐射到成像目标并经成像目标反射后得到第一波束;
第一波束传播距离f后经过第一透镜调制后得到第二波束;
第二波束传播距离f后经过频谱面的波束挡板后得到第三波束;
第三波束传播距离f后经过第二透镜调制后得到第四波束;
第四波束经第二透镜调制后传播距离f被太赫兹相机接收成像。
2.根据权利要求1所述的反射式太赫兹成像方法,其特征在于,
所述第一波束传播距离f过程中经所述分束镜透射后到达第一透镜;
所述第二波束中心被波束挡板遮挡。
3.根据权利要求1所述的反射式太赫兹成像方法,其特征在于,
所述第一透镜和第二透镜均为凸透镜,第一透镜和第二透镜的焦距均为f。
4.根据权利要求1所述的反射式太赫兹成像方法,其特征在于,
经所述成像目标反射后的第一波束起始位置的场分布E1(x1,y1)为:
E1(x1,y1)=E0(x1,y1)t(x1,y1);
其中,t(x1,y1)为目标的反射函数,E0(x1,y1)为辐射到成像目标的场分布,(x1,y1)为目标平面的坐标。
5.根据权利要求4所述的反射式太赫兹成像方法,其特征在于,
所述到达第一透镜的场分布满足:
F{E2(x2,y2)}=F{E1(x1,y1)}H(fx,fy);
其中,F{·}为傅里叶变换,fx=x2/λf,fy=y2/λf,j为虚数单位,λ为太赫兹波波长,(x2,y2)为第一透镜所在平面的坐标。
8.根据权利要求7所述的反射式太赫兹成像方法,其特征在于,
所述第三波束到达第二透镜的场分布E4(x4,y4)为:
F{E4(x4,y4)}=F{E′3(x3,y3)}H(fx,fy);
其中,fx=x4/λf,fy=y4/λf,(x4,y4)为第二透镜所在平面的坐标。
10.一种成像装置,用于执行如权利要求1-9任一项所述的反射式太赫兹成像方法,其特征在于,包括:
分束镜,用于将太赫兹波束反射到成像目标上并透射经目标反射的太赫兹波;
第一透镜,用于调制经过成像目标反射后的第一波束;
波束挡板,用于遮挡经过第一透镜的第二波束中心;
第二透镜,用于调制经过波束挡的第三波束;以及
太赫兹相机,用于接收经过第二透镜后的第四波束并成像;
其中,所述波束挡板设置在第二透射波束中心位置;
其中,所述成像目标、第一透镜、螺旋相位板、第二透镜、太赫兹相机等距且平行设置;
其中,所述分束镜设置在成像目标与第一透镜之间;
其中,所述太赫兹波束经过分束镜反射后的波束传播方向与成像目标表面垂直。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911152543.4A CN110849841A (zh) | 2019-11-20 | 2019-11-20 | 基于暗场成像的反射式太赫兹成像方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911152543.4A CN110849841A (zh) | 2019-11-20 | 2019-11-20 | 基于暗场成像的反射式太赫兹成像方法及装置 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN110849841A true CN110849841A (zh) | 2020-02-28 |
Family
ID=69603513
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911152543.4A Pending CN110849841A (zh) | 2019-11-20 | 2019-11-20 | 基于暗场成像的反射式太赫兹成像方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110849841A (zh) |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2005172774A (ja) * | 2003-12-05 | 2005-06-30 | Semiconductor Res Found | 反射光学特性によって物性を測定する装置および測定方法 |
CN103348235A (zh) * | 2011-02-10 | 2013-10-09 | 株式会社日立高新技术 | 异物检测装置和异物检测方法 |
CN105534481A (zh) * | 2016-01-21 | 2016-05-04 | 华中科技大学 | 一种频域光学相干层析装置与方法 |
CN105581773A (zh) * | 2014-10-21 | 2016-05-18 | 深圳先进技术研究院 | 一种太赫兹成像系统 |
CN110455834A (zh) * | 2019-07-23 | 2019-11-15 | 中国科学院上海光学精密机械研究所 | 基于光强传输方程的x射线单次曝光成像装置及方法 |
-
2019
- 2019-11-20 CN CN201911152543.4A patent/CN110849841A/zh active Pending
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2005172774A (ja) * | 2003-12-05 | 2005-06-30 | Semiconductor Res Found | 反射光学特性によって物性を測定する装置および測定方法 |
CN103348235A (zh) * | 2011-02-10 | 2013-10-09 | 株式会社日立高新技术 | 异物检测装置和异物检测方法 |
CN105581773A (zh) * | 2014-10-21 | 2016-05-18 | 深圳先进技术研究院 | 一种太赫兹成像系统 |
CN105534481A (zh) * | 2016-01-21 | 2016-05-04 | 华中科技大学 | 一种频域光学相干层析装置与方法 |
CN110455834A (zh) * | 2019-07-23 | 2019-11-15 | 中国科学院上海光学精密机械研究所 | 基于光强传输方程的x射线单次曝光成像装置及方法 |
Non-Patent Citations (2)
Title |
---|
王炯: "螺旋相位板参数对螺旋相衬成像系统性能影响研究", 《光学学报》 * |
王秉超: "《光学》", 31 December 1991 * |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Poon et al. | Fundamentals of optics | |
Badon et al. | Smart optical coherence tomography for ultra-deep imaging through highly scattering media | |
US11143555B2 (en) | Methods and devices for reducing spectral noise and spectrometry systems employing such devices | |
CN104783767B (zh) | 一种利用正交偏振光谱成像探测人体微循环的装置及方法 | |
US9377362B2 (en) | Systems and methods for high-contrast, near-real-time acquisition of terahertz images | |
Tchaka et al. | Chromaticity based smoke removal in endoscopic images | |
Yushkov et al. | Contrast enhancement in microscopy of human thyroid tumors by means of acousto-optic adaptive spatial filtering | |
WO2019175546A1 (en) | Optical device | |
Ritt et al. | Evaluation of protection measures against laser dazzling for imaging sensors | |
CN110793943A (zh) | 基于螺旋相衬成像的反射式太赫兹成像方法及装置 | |
CN110823832B (zh) | 基于暗场成像的太赫兹成像方法及装置 | |
CN110849841A (zh) | 基于暗场成像的反射式太赫兹成像方法及装置 | |
Huang et al. | Active imaging through dense fog by utilizing the joint polarization defogging and denoising optimization based on range-gated detection | |
JP6661791B2 (ja) | 防眩撮像カメラ及び方法 | |
Venkanna et al. | Edge imaging characteristics of aberrated coherent optical systems by edge masking of circular apertures | |
CN111061063B (zh) | 光瞳滤波远场超分辨成像系统及光瞳滤波器设计方法 | |
CN117288720A (zh) | 一种基于子区域均一化的非侵入式单帧宽谱散射成像系统及成像方法 | |
Holmes et al. | Image reconstructions with active illumination in strong-turbulence scenarios with single-frame blind deconvolution approaches | |
CN210534350U (zh) | 毫米波太赫兹成像设备 | |
Stepanov et al. | Impact of illumination spectrum and eye pigmentation on image quality from a fundus camera using transscleral illumination | |
Tahara | Adding dimensions with Lucy–Richardson–Rosen algorithm to incoherent imaging | |
Tajahuerce et al. | Phase-object fractional differentiation using Fourier plane filters | |
Liu et al. | Enhancement of low-quality reconstructed digital hologram images based on frequency extrapolation of large objects under the diffraction limit | |
Leng et al. | Speckle Reduction by Optical Pre-processing and Digital Filtering in Digital Holography. | |
Xu et al. | Image clarity of membrane diffractive imaging system using total-scale Retinex |
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 | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20200228 |
|
WD01 | Invention patent application deemed withdrawn after publication |