CN110823832B - 基于暗场成像的太赫兹成像方法及装置 - Google Patents

基于暗场成像的太赫兹成像方法及装置 Download PDF

Info

Publication number
CN110823832B
CN110823832B CN201911152487.4A CN201911152487A CN110823832B CN 110823832 B CN110823832 B CN 110823832B CN 201911152487 A CN201911152487 A CN 201911152487A CN 110823832 B CN110823832 B CN 110823832B
Authority
CN
China
Prior art keywords
lens
terahertz
imaging
transmission beam
transmitted
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
CN201911152487.4A
Other languages
English (en)
Other versions
CN110823832A (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.)
Institute of Electronics of CAS
Original Assignee
Institute of Electronics of CAS
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 Institute of Electronics of CAS filed Critical Institute of Electronics of CAS
Priority to CN201911152487.4A priority Critical patent/CN110823832B/zh
Publication of CN110823832A publication Critical patent/CN110823832A/zh
Application granted granted Critical
Publication of CN110823832B publication Critical patent/CN110823832B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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、利用太赫兹相机进行接收，可以实现对成像目标的实时成像。
附图说明
图1为本发明实施例的基于暗场成像的太赫兹高对比度成像方法的步骤示意图；
图2为本发明实施例的基于暗场成像的太赫兹高对比度成像装置的光路结构示意图；
图3为本发明实施例的基于暗场成像的太赫兹高对比度成像方法的仿真成像结果。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明白，以下结合具体实施例，并参照附图，对本发明作进一步的详细说明。
目前，传统的太赫兹生物医学成像方法大多成像对比度低，成像质量差。本发明旨在将光学成像中的暗场成像引入到太赫兹成像中，以提高成像对比度。相对于传统成像方法，暗场成像在频谱面加入波束挡板，用于滤除目标的低频信息，从而提取出目标的高频信息。由于低频信息代表背景，高频信息代表目标的细节，因此，暗场成像可以抑制背景，提取出目标细节信息，从而实现高对比度成像。成像系统中利用太赫兹相机作为接收，可以实现对目标实时高对比度成像。由于太赫兹波长相对于光波要长的多，容易发生较强的衍射效应，因此实验系统元件尺寸和间距应远大于波长。
本发明公开了一种基于暗场成像的太赫兹成像方法，包括：
准直的太赫兹波束经目标透射后得到第一透射波束；
第一透射波束传播距离f后经第一透镜调制后得到第二透射波束；
第二透射波束传播距离f后经过频谱面的波束挡板后得到第三透射波束；
第三透射波束传播距离f后经第二透镜调制后得到第四透射波束；
第四透射波束传播距离f后被太赫兹相机接收后成像。
其中，所述第二透射波束中心被波束挡板遮挡。
其中，所述第一透镜和第二透镜均为凸透镜，第一透镜和第二透镜的焦距均为f。
其中，所述第一透射波束起始位置的场分布E1(x1，y1)为：
E1(x1，y1)＝E0(x1，y1)t(x1，y1)；
其中，E0(x1，y1)为照射到目标的准直太赫兹场分布，t(x1，y1)为目标的透射函数，(x1，y1)为目标平面的坐标。
其中，所述第一透射波束到达第一透镜的场分布E2(x2，y2)满足：
F{E2(x2，y2)}＝F{E1(x1，y1)}H(fx，fy)；
其中，F{·}为傅里叶变换，
Figure BDA0002281150990000031
fx＝x2/λf，fy＝y2/λf，j为虚数单位，λ为太赫兹波波长，(x2，y2)为第一透镜所在平面的坐标。
其中，所述第二透射波束到达频谱面的场分布E3(x3，y3)为：
Figure BDA0002281150990000041
其中，k为波数，
Figure BDA0002281150990000042
为第一透镜的调制函数，(x3，y3)为频谱面的坐标。
其中，所述第二透射波束到达频谱面的场分布展开为：
Figure BDA0002281150990000043
其中，第一项为没有目标时频谱面上的场分布，即低频成分；第二项是由衍射效应引起的，代表目标的细节信息，即高频成分；
由于波束挡板遮挡了第二透射波束的中心，即低频成分，保留了高频成分；所以第二透射波束未被遮挡的部分位于频谱面的场分布E′3(x3，y3)为：
Figure BDA0002281150990000044
其中，所述第三透射波束到达第二透镜的场分布E4(x4，y4)为：
F{E4(x4，y4)}＝F{E′3(x3，y3)}H(fx，fy)；
其中，fx＝x4/λf，fy＝y4/λf，(x4，y4)为第二透镜所在平面的坐标。
其中，所述第四透射波束到达太赫兹相机表面的场分布E4(x4，y4)为：
Figure BDA0002281150990000051
其中，(x5，y5)为太赫兹相机所在平面的坐标；
其中，所述太赫兹相机接收的为光强，其分布I为：
Figure BDA0002281150990000052
本发明还公开了一种成像装置，用于执行如上所述的太赫兹成像方法，包括：
第一透镜，用于调制经过成像目标透射后的第一透射波束；
波束挡板，用于遮挡经过第一透镜的第二透射波束中心；
第二透镜，用于调制经过波束挡板的第三透射波束；以及
太赫兹相机，用于接收经过第二透镜后的第四透射波束并成像；
其中，所述成像目标、第一透镜、波束挡板、第二透镜、太赫兹相机等距且平行设置；
其中，所述波束挡板设置在第二透射波束中心位置。
以下通过具体实施例结合附图对本发明的技术方案做进一步阐述说明。需要注意的是，下述的具体实施例仅是作为举例说明，本发明的保护范围并不限于此。
如图1和图2所示，本实施例的基于暗场成像的太赫兹高对比度成像方法，包括以下步骤：
S1、准直的太赫兹波束经目标透射后得到第一透射波束。
所述第一透射波束起始位置的场分布为：
E1(x1，y1)＝E0(x1，y1)t(x1，y1)；
其中，E0(x1，y1)为照射到目标的准直太赫兹场分布，t(x1，y1)为目标的透射函数，(x1，y1)为目标平面的坐标。
S2、所述第一透射波束传播距离f后到达第一透镜。
第一透射波束经过第一透镜后形成第二透射波束；
所述到达第一透镜的场分布满足：
F{E2(x2，y2)}＝F{E1(x1，y1)}H(fx，fy)；
其中，F{·}为傅里叶变换，
Figure BDA0002281150990000061
fx＝x2/λf，fy＝y2/λf，j为虚数单位，λ为太赫兹波波长，(x2，y2)为第一透镜所在平面的坐标。
S3、所述第二透射波束经第一透镜调制后传播距离f到达频谱面。
到达频谱面的场分布为：
Figure BDA0002281150990000071
其中，k为波数，
Figure BDA0002281150990000072
为第一透镜的调制函数，(x3，y3)为频谱面的坐标。
S4、所述第二透射波束中心被位于频谱面的波束挡板遮挡，第二透射波束其余部分未被遮挡；
第二透射波束经过波束挡板后形成第三透射波束；
到达频谱面的场分布可以改写为：
Figure BDA0002281150990000073
其中，第一项为没有目标时频谱面上的场分布，即低频成分；第二项是由衍射效应引起的，代表目标的细节信息，即高频成分。波束挡板遮挡了波束的中心，即低频成分，保留了高频成分。第二透射波束未被遮挡的部分场分布为：
Figure BDA0002281150990000074
S5、所述第二透射波束未被遮挡的部分(即第三透射波束)传播距离f后到达第二透镜；
到达第二透镜的场分布为：
F{E4(x4，y4)}＝F{E′3(x3，y3)}H(fx，fy)；
其中，fx＝x4/λf，fy＝y4/λf，(x4，y4)为第二透镜所在平面的坐标。
S6、所述第三波束经第二透镜调制后传播距离f被太赫兹相机接收。第三透射波束经过第二透镜后得到第四透射波束。到达太赫兹相机表面的场分布为：
Figure BDA0002281150990000081
其中，(x5，y5)为太赫兹相机所在平面的坐标。
所述太赫兹相机接收的为光强，其分布为：
Figure BDA0002281150990000082
图3为仿真结果，图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 (3)

1.一种基于暗场成像的太赫兹成像方法，包括：
准直的太赫兹波束经目标透射后得到第一透射波束；
第一透射波束传播距离f后经第一透镜调制后得到第二透射波束；
第二透射波束传播距离f后经过频谱面的波束挡板后得到第三透射波束；
第三透射波束传播距离f后经第二透镜调制后得到第四透射波束；
第四透射波束传播距离f后被太赫兹相机接收后成像；
其中，所述第一透镜和第二透镜均为凸透镜，第一透镜和第二透镜的焦距均为f；
其中，所述第一透射波束起始位置的场分布E1(x1，y1)为：
E1(x1，y1)＝E0(x1，y1)t(x1，y1)；
其中，E0(x1，y1)为照射到目标的准直太赫兹场分布，t(x1，y1)为目标的透射函数，(x1，y1)为目标平面的坐标；
所述第一透射波束到达第一透镜的场分布E2(x2，y2)满足：
F{E2(x2，y2)}＝F{E1(x1，y1)}H(fx，fy)；
其中，F{·}为傅里叶变换，
Figure FDA0002813454640000011
fx＝x2/( λf) ，fy＝y2/( λf) ，j为虚数单位，λ为太赫兹波波长，(x2，y2)为第一透镜所在平面的坐标；
所述第二透射波束到达频谱面的场分布E3(x3，y3)为：
Figure FDA0002813454640000021
其中，k为波数，
Figure FDA0002813454640000022
为第一透镜的调制函数，(x3，y3)为频谱面的坐标；
所述第二透射波束到达频谱面的场分布展开为：
Figure FDA0002813454640000023
其中，第一项为没有目标时频谱面上的场分布，即低频成分；第二项是由衍射效应引起的，代表目标的细节信息，即高频成分；
由于波束挡板遮挡了第二透射波束的中心，即低频成分，保留了高频成分，所以第二透射波束未被遮挡的部分位于频谱面的场分布E′3(x3，y3)为：
Figure FDA0002813454640000024
所述第三透射波束到达第二透镜的场分布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)为：
Figure FDA0002813454640000031
其中，(x5，y5)为太赫兹相机所在平面的坐标；
其中，所述太赫兹相机接收的为光强，其分布I为：
Figure FDA0002813454640000032
2.根据权利要求1所述的太赫兹成像方法，其特征在于，
所述第二透射波束中心被波束挡板遮挡。
3.一种成像装置，用于执行如权利要求1-2任一项所述的太赫兹成像方法，其特征在于，包括：
第一透镜，用于调制经过成像目标透射后的第一透射波束；
波束挡板，用于遮挡经过第一透镜的第二透射波束中心；
第二透镜，用于调制经过波束挡板的第三透射波束；以及
太赫兹相机，用于接收经过第二透镜后的第四透射波束并成像；
其中，所述成像目标、第一透镜、波束挡板、第二透镜、太赫兹相机等距且平行设置；
其中，所述波束挡板设置在第二透射波束中心位置。
CN201911152487.4A 2019-11-20 2019-11-20 基于暗场成像的太赫兹成像方法及装置 Active CN110823832B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911152487.4A CN110823832B (zh) 2019-11-20 2019-11-20 基于暗场成像的太赫兹成像方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911152487.4A CN110823832B (zh) 2019-11-20 2019-11-20 基于暗场成像的太赫兹成像方法及装置

Publications (2)

Publication Number Publication Date
CN110823832A CN110823832A (zh) 2020-02-21
CN110823832B true CN110823832B (zh) 2021-03-19

Family

ID=69558031

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911152487.4A Active CN110823832B (zh) 2019-11-20 2019-11-20 基于暗场成像的太赫兹成像方法及装置

Country Status (1)

Country Link
CN (1) CN110823832B (zh)

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103743681B (zh) * 2014-01-24 2017-01-25 中国工程物理研究院流体物理研究所 一种太赫兹光谱仪及太赫兹收发探头
CN110160984B (zh) * 2019-01-08 2021-12-24 南开大学 一种基于超表面和铌酸锂混合结构的片上太赫兹传感增强器件
CN110455834B (zh) * 2019-07-23 2021-11-05 中国科学院上海光学精密机械研究所 基于光强传输方程的x射线单次曝光成像装置及方法

Also Published As

Publication number Publication date
CN110823832A (zh) 2020-02-21

Similar Documents

Publication Publication Date Title
US10401294B2 (en) OAM microscope for edge enhancement of biomedical and condensed matter samples and objects
Li et al. Underwater image enhancement by dehazing and color correction
Poon et al. Fundamentals of optics
Tchaka et al. Chromaticity based smoke removal in endoscopic images
US10373296B2 (en) Method and system for reducing ringing artifacts in X-ray image
US20230306558A1 (en) Frequency domain-based method for removing periodic noise from reconstructed light field image
WO2014002744A1 (ja) 画像処理装置、画像処理方法、及び画像処理プログラム
CN110823832B (zh) 基于暗场成像的太赫兹成像方法及装置
CN101995731B (zh) 一种基于液晶的直视型防眩光成像仪及防眩光成像方法
CN109993703A (zh) 多尺度图像降噪方法及装置
CN110811615A (zh) 基于螺旋相衬成像的太赫兹成像方法及装置
CN110793943A (zh) 基于螺旋相衬成像的反射式太赫兹成像方法及装置
CN110849841A (zh) 基于暗场成像的反射式太赫兹成像方法及装置
Zhang et al. Gabor-based anisotropic diffusion for speckle noise reduction in medical ultrasonography
Venkanna et al. Edge imaging characteristics of aberrated coherent optical systems by edge masking of circular apertures
CN111192208B (zh) 一种基于边窗滤波器的牙齿cr图像增强方法及装置
Bedi et al. Ultrasound image despeckling and enhancement using modified multiscale anisotropic diffusion model in non-subsampled shearlet domain
CN115578304A (zh) 一种结合显著性区域检测的多波段图像融合方法和系统
Di Claudio et al. Predicting blur visual discomfort for natural scenes by the loss of positional information
KR20230102142A (ko) 1회 촬영으로 듀얼 에너지 영상 획득이 가능한 x-선 센싱 장치 및 영상처리 방법
Rao et al. Contrast enhancement for the infrared vein image of leg based on the optical angular spectrum theory
Song et al. Research on image enhancement algorithms based on Matlab
Liu et al. Enhancement of low-quality reconstructed digital hologram images based on frequency extrapolation of large objects under the diffraction limit
Yang et al. Designing Phase Masks for Under-Display Cameras
Wang et al. Nonlinear retinal image enhancement for vessel detection

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