CN110823832A - 基于暗场成像的太赫兹成像方法及装置 - Google Patents
基于暗场成像的太赫兹成像方法及装置 Download PDFInfo
- Publication number
- CN110823832A CN110823832A CN201911152487.4A CN201911152487A CN110823832A CN 110823832 A CN110823832 A CN 110823832A CN 201911152487 A CN201911152487 A CN 201911152487A CN 110823832 A CN110823832 A CN 110823832A
- Authority
- CN
- China
- Prior art keywords
- terahertz
- lens
- transmission beam
- imaging
- imaging method
- 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.)
- Granted
Links
- 238000003384 imaging method Methods 0.000 title claims abstract description 94
- 230000005540 biological transmission Effects 0.000 claims abstract description 63
- 238000001228 spectrum Methods 0.000 claims abstract description 7
- 230000003595 spectral effect Effects 0.000 claims description 17
- 230000000694 effects Effects 0.000 claims description 4
- 230000000903 blocking effect Effects 0.000 claims description 3
- 238000005516 engineering process Methods 0.000 abstract description 3
- 238000004088 simulation Methods 0.000 description 5
- 238000012634 optical imaging Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 238000002834 transmittance Methods 0.000 description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 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
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);
其中,所述第二透射波束到达频谱面的场分布E3(x3,y3)为:
其中,所述第二透射波束到达频谱面的场分布展开为:
其中,第一项为没有目标时频谱面上的场分布,即低频成分;第二项是由衍射效应引起的,代表目标的细节信息,即高频成分;
由于波束挡板遮挡了第二透射波束的中心,即低频成分,保留了高频成分;所以第二透射波束未被遮挡的部分位于频谱面的场分布E′3(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)为第二透镜所在平面的坐标。
其中,所述第四透射波束到达太赫兹相机表面的场分布E4(x4,y4)为:
其中,(x5,y5)为太赫兹相机所在平面的坐标;
其中,所述太赫兹相机接收的为光强,其分布I为:
本发明还公开了一种成像装置,用于执行如上所述的太赫兹成像方法,包括:
第一透镜,用于调制经过成像目标透射后的第一透射波束;
波束挡板,用于遮挡经过第一透镜的第二透射波束中心;
第二透镜,用于调制经过波束挡板的第三透射波束;以及
太赫兹相机,用于接收经过第二透镜后的第四透射波束并成像;
其中,所述成像目标、第一透镜、波束挡板、第二透镜、太赫兹相机等距且平行设置;
其中,所述波束挡板设置在第二透射波束中心位置。
以下通过具体实施例结合附图对本发明的技术方案做进一步阐述说明。需要注意的是,下述的具体实施例仅是作为举例说明,本发明的保护范围并不限于此。
如图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);
S3、所述第二透射波束经第一透镜调制后传播距离f到达频谱面。
到达频谱面的场分布为:
其中,k为波数,为第一透镜的调制函数,(x3,y3)为频谱面的坐标。
S4、所述第二透射波束中心被位于频谱面的波束挡板遮挡,第二透射波束其余部分未被遮挡;
第二透射波束经过波束挡板后形成第三透射波束;
到达频谱面的场分布可以改写为:
其中,第一项为没有目标时频谱面上的场分布,即低频成分;第二项是由衍射效应引起的,代表目标的细节信息,即高频成分。波束挡板遮挡了波束的中心,即低频成分,保留了高频成分。第二透射波束未被遮挡的部分场分布为:
S5、所述第二透射波束未被遮挡的部分(即第三透射波束)传播距离f后到达第二透镜;
到达第二透镜的场分布为:
F{E4(x4,y4)}=F{E′3(x3,y3)}H(fx,fy);
其中,fx=x4/λf,fy=y4/λf,(x4,y4)为第二透镜所在平面的坐标。
S6、所述第三波束经第二透镜调制后传播距离f被太赫兹相机接收。第三透射波束经过第二透镜后得到第四透射波束。到达太赫兹相机表面的场分布为:
其中,(x5,y5)为太赫兹相机所在平面的坐标。
所述太赫兹相机接收的为光强,其分布为:
图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 (10)
1.一种基于暗场成像的太赫兹成像方法,包括:
准直的太赫兹波束经目标透射后得到第一透射波束;
第一透射波束传播距离f后经第一透镜调制后得到第二透射波束;
第二透射波束传播距离f后经过频谱面的波束挡板后得到第三透射波束;
第三透射波束传播距离f后经第二透镜调制后得到第四透射波束;
第四透射波束传播距离f后被太赫兹相机接收后成像。
2.根据权利要求1所述的太赫兹成像方法,其特征在于,
所述第二透射波束中心被波束挡板遮挡。
3.根据权利要求1所述的太赫兹成像方法,其特征在于,
所述第一透镜和第二透镜均为凸透镜,第一透镜和第二透镜的焦距均为f。
4.根据权利要求1所述的太赫兹成像方法,其特征在于,
所述第一透射波束起始位置的场分布E1(x1,y1)为:
E1(x1,y1)=E0(x1,y1)t(x1,y1);
其中,E0(x1,y1)为照射到目标的准直太赫兹场分布,t(x1,y1)为目标的透射函数,(x1,y1)为目标平面的坐标。
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 |
---|---|---|---|
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 true CN110823832A (zh) | 2020-02-21 |
CN110823832B 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) |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103743681A (zh) * | 2014-01-24 | 2014-04-23 | 中国工程物理研究院流体物理研究所 | 一种太赫兹光谱仪及太赫兹收发探头 |
CN110160984A (zh) * | 2019-01-08 | 2019-08-23 | 南开大学 | 一种基于超表面和铌酸锂混合结构的片上太赫兹传感增强器件 |
CN110455834A (zh) * | 2019-07-23 | 2019-11-15 | 中国科学院上海光学精密机械研究所 | 基于光强传输方程的x射线单次曝光成像装置及方法 |
-
2019
- 2019-11-20 CN CN201911152487.4A patent/CN110823832B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103743681A (zh) * | 2014-01-24 | 2014-04-23 | 中国工程物理研究院流体物理研究所 | 一种太赫兹光谱仪及太赫兹收发探头 |
CN110160984A (zh) * | 2019-01-08 | 2019-08-23 | 南开大学 | 一种基于超表面和铌酸锂混合结构的片上太赫兹传感增强器件 |
CN110455834A (zh) * | 2019-07-23 | 2019-11-15 | 中国科学院上海光学精密机械研究所 | 基于光强传输方程的x射线单次曝光成像装置及方法 |
Non-Patent Citations (1)
Title |
---|
王秉超等: "《光学》", 31 December 1991, 吉林大学出版社 * |
Also Published As
Publication number | Publication date |
---|---|
CN110823832B (zh) | 2021-03-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Li et al. | Underwater image enhancement by dehazing and color correction | |
US20230306558A1 (en) | Frequency domain-based method for removing periodic noise from reconstructed light field image | |
US20170294001A1 (en) | Method and system for reducing ringing artifacts in x-ray image | |
WO2014002744A1 (ja) | 画像処理装置、画像処理方法、及び画像処理プログラム | |
Lu et al. | Single underwater image descattering and color correction | |
CN110823832B (zh) | 基于暗场成像的太赫兹成像方法及装置 | |
CN101995731B (zh) | 一种基于液晶的直视型防眩光成像仪及防眩光成像方法 | |
CN109993703A (zh) | 多尺度图像降噪方法及装置 | |
CN111192208A (zh) | 一种基于边窗滤波器的牙齿cr图像增强方法及装置 | |
CN110811615A (zh) | 基于螺旋相衬成像的太赫兹成像方法及装置 | |
Sandeep | Remote sensing image dehazing using guided filter | |
CN110793943A (zh) | 基于螺旋相衬成像的反射式太赫兹成像方法及装置 | |
CN110849841A (zh) | 基于暗场成像的反射式太赫兹成像方法及装置 | |
Cao et al. | Dehazing method through polarimetric imaging and multi-scale analysis | |
Venkanna et al. | Edge imaging characteristics of aberrated coherent optical systems by edge masking of circular apertures | |
CN117288720A (zh) | 一种基于子区域均一化的非侵入式单帧宽谱散射成像系统及成像方法 | |
CN117391972A (zh) | 一种基于同态滤波的水下偏振成像方法 | |
Zhou et al. | Inverse sinusoidal phase mask to extend the depth of field of incoherent imaging systems | |
Luo et al. | Endoscopic video defogging using luminance blending | |
Rao et al. | Contrast enhancement for the infrared vein image of leg based on the optical angular spectrum theory | |
Meng et al. | Learning based polarization image fusion under an alternative paradigm | |
Liu et al. | Multi-Focus Image Fusion Based on Non-Subsampled Shearlet Transform and Consistency Optimization | |
Liu et al. | Enhancement of low-quality reconstructed digital hologram images based on frequency extrapolation of large objects under the diffraction limit | |
Song et al. | Research on image enhancement algorithms based on Matlab | |
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 |