一种基于太赫兹衍射图分解的太赫兹波数字全息成像方法技术

本发明专利技术公开了一种基于太赫兹衍射图分解的太赫兹波数字全息成像方法，包括太赫兹全息图记录，全息图预传播，预传播衍射面复振幅图像分解和多层面复振幅图像重建合成过程。其中：(1)太赫兹全息图记录过程：在太赫兹数字全息成像系统中采集处理得到大孔径归一化太赫兹全息图。(2)全息图预传播过程：自由传播特定距离获得预传播衍射面的复振幅图像。(3)预传播衍射面复振幅图像分解过程：计算预传播衍射面复振幅图像前景面和后景面阈值，制作二值化掩膜图像，获得预传播衍射面的子复振幅图像。(4)多层面复振幅图像重建合成过程：重建各物平面的子复振幅图像，合成得到太赫兹波数字全息多层面成像结果。

A THz Wave Digital Holographic Imaging Method Based on Terahertz Diffraction Pattern Decomposition

【技术实现步骤摘要】
一种基于太赫兹衍射图分解的太赫兹波数字全息成像方法
本专利技术涉及太赫兹波数字全息多层面成像方法
，具体为一种基于太赫兹衍射图分解的太赫兹波数字全息成像方法。
技术介绍
太赫兹波段因处于红外到微波的过渡波段而兼具两者的特性，使之在成像领域应用前景广阔，而数字全息技术因其成像结构简单和灵活的数据处理方式，对特定目标样品的多个聚焦层面成像有着独特优势。在传统成像过程中，存在样品分布在多个聚焦距离而无法获得所需多聚焦层面探测结果的问题，太赫兹数字全息技术可以充分发挥两者的优势和特点，成为一种在太赫兹波段的实现多聚焦层面成像技术。本专利技术旨在提出一种基于衍射图分解的太赫兹波数字全息分层成像方法。这是一项在太赫兹波段的新型成像技术，其依靠计算机对太赫兹全息图的数字处理过程，通过单次成像记录样品的太赫兹波全息图即可获得样品包含不同深度方向聚焦面的振幅和相位信息。相比传统太赫兹波成像手段，该专利技术无需机械调节聚焦面位置以及多次曝光测量，即可实现样品的多层面复振幅探测成像。结合太赫兹波本身的选择穿透性和无损安全性，该专利技术是一种能够穿透非金属非极性掩盖物，单次记录曝光即可实现不同深度聚焦面信息的成像技术，在无损探测和材料表征领域会有广泛应用。
技术实现思路
本专利技术的目的在于提供一种基于太赫兹衍射图分解的太赫兹波数字全息成像方法，以解决上述
技术介绍
中提出的现有技术不足，获得目标样品不同深度聚焦面的复振幅信息，实现样品的多层面复振幅成像。为实现上述目的，本专利技术提供如下技术方案：一种基于太赫兹衍射图分解的太赫兹波数字全息成像方法，包括以下步骤：太赫兹全息图记录过程，全息图预传播过程，预传播衍射面复振幅图像分解过程，多层面复振幅图像重建合成过程：(1)太赫兹全息图记录过程。将深度方向上包含多层面信息的样品放置于太赫兹数字全息成像系统中，太赫兹波经过样品调制，携带样品的多层面信息，与不经过样品调制的太赫兹参考波干涉形成太赫兹全息图。为提高系统成像分辨率，搭载太赫兹探测器的二维平移台在垂直于太赫兹波照明方向的记录平面内按照给定扫描顺序进行二维移动，采集不同位置的太赫兹子全息图。移除样品后，搭载太赫兹探测器的二维平移台按照太赫兹子全息图记录过程中相同的扫描顺序，记录太赫兹子背景图像。将相同位置采集的对应太赫兹子全息图和太赫兹子背景图相除，得到子归一化全息图，按照扫描顺序进行子归一化全息图图像拼接，得到大孔径归一化太赫兹全息图。(2)全息图预传播过程。将步骤(1)中的大孔径归一化太赫兹全息图在计算机中通过自由空间传播算法，沿物体平面方向预传播至衍射平面P，获得预传播衍射面的振幅和相位分布，即预传播衍射面复振幅图像U(m，n)。该衍射面的传播距离为dp。定义样品不同分层面距离太赫兹探测器记录平面的最大距离和最小距离分别为dmax和dmin，则dp的取值定义为(3)预传播衍射面复振幅图像分解过程。计算提取预传播衍射面复振幅图像U(m，n)的振幅和相位图像，采用最大类间方差法进行前景面和后景面的最大方差阈值T的判定计算，基于得到的阈值T进行预传播衍射面复振幅图像的二值化运算，得到二值化掩膜图像K(m，n)，可表示为对二值化掩膜图像K(m，n)完成边缘均值模糊化处理，得到预传播衍射面的掩膜图像K′(m，n)，通过非零值区域识别，对掩膜图像K′(m，n)中的x个孤立非零值区域图像进行逐一提取，并将其补零值至原掩膜图像K′(m，n)的像素尺寸大小，得到x个只包含单个孤立非零值图像区域的子掩膜图像K′i(i＝1，2，...，x)。将得到的x个子掩膜图像K′i分别与预传播衍射面的复振幅图像U相乘，获得x个预传播衍射面的子复振幅图像Ui(i＝1，2，...，x)。Ui＝U×K'i,i＝1,2,...,x(4)多层面复振幅图像重建合成过程。通过自动聚焦算法计算预传播衍射面的各个子复振幅图像Ui的聚焦传播距离，并通过自由空间传播算法进行传播，重建得到各个物平面复振幅图像Oi(i＝1，2，...，x)。将各个物平面复振幅图像Oi分别与各自对应的子掩膜图像K′i相乘得到图像O′i，计算各个物平面复振幅图像Oi的平均值并与做取反运算的对应子掩膜图像(1-K′i)相乘得到图像O″i，将图像O′i和图像O″i相加合成目标样品的多层面复振幅成像结果O，该步骤可表示为优选的，所述雪崩二极管太赫兹波源的一侧输出有太赫兹波，且雪崩二极管太赫兹波源朝向太赫兹波的一侧设置有样品，所述太赫兹波热释电探测器安装于太赫兹波的末端，且太赫兹波热释电探测器的下方固定有探测器二维平移台。优选的，所述雪崩二极管太赫兹波源、样品和太赫兹波热释电探测器的结构为太赫兹同轴成像系统。优选的，所述预传播衍射面复振幅图像分解过程，针对振幅型样品，提取预传播衍射面的振幅分布图像，针对相位型样品，提取预传播衍射面的相位分布图像。本专利技术提供了一种基于太赫兹衍射图分解的太赫兹波数字全息成像方法，与现有技术相比，本专利技术的有益效果如下：本专利技术提供了一种基于衍射图分解的太赫兹波数字全息分层成像方法，通过该专利技术的方法和系统，对太赫兹全息图进行数字处理，可非接触、无损获得目标样品不同深度聚焦面的复振幅信息，实现样品在太赫兹波段的多层面复振幅成像，实用性广泛。附图说明图1为本专利技术一种基于太赫兹衍射图分解的太赫兹波数字全息成像方法的流程图；图2为本专利技术一种基于太赫兹衍射图分解的太赫兹波数字全息分层成像系统的示意图。图中：1、雪崩二极管太赫兹波源；1a、太赫兹波；2、样品；3、探测器二维平移台；4、太赫兹波热释电探测器。具体实施方式下面将结合本专利技术实施例中的数据，对本专利技术实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本专利技术一部分实施例，而不是全部的实施例。请参阅图1-2，本专利技术提供一种技术方案：一种基于太赫兹衍射图分解的太赫兹波数字全息成像方法，包括以下步骤：(1)太赫兹全息图记录过程。将深度方向上包含多层面信息的样品放置于太赫兹数字全息成像系统中，太赫兹波经过样品调制，携带样品的多层面信息，与不经过样品调制的太赫兹参考波干涉形成太赫兹全息图。为提高系统成像分辨率，搭载太赫兹探测器的二维平移台在垂直于太赫兹波照明方向的记录平面内按照给定扫描顺序进行二维移动，采集不同位置的太赫兹子全息图。移除样品后，搭载太赫兹探测器的二维平移台按照太赫兹子全息图记录过程中相同的扫描顺序，记录太赫兹子背景图像。将相同位置采集的对应太赫兹子全息图和太赫兹子背景图相除，得到子归一化全息图，按照扫描顺序进行子归一化全息图图像拼接，得到大孔径归一化太赫兹全息图。(2)全息图预传播过程。将步骤(1)中的大孔径归一化太赫兹全息图在计算机中通过自由空间传播算法，沿物体平面方向预传播至衍射平面P，获得预传播衍射面的振幅和相位分布，即预传播衍射面复振幅图像U(m，n)。该衍射面的传播距离为dp。定义样品不同分层面距离太赫兹探测器记录平面的最大距离和最小距离分别为dmax和dmin，则dp的取值定义为(3)预传播衍射面复振幅图像分解过程。针对振幅型样品，提取预传播衍射面的振幅分布图像，针对相位型样品，提取预传播衍射面的相位分布图像。计算提取预传播衍射面复振幅图像U(m，n)的振幅和相位图像，采用最大类间方差法进行前景面和后景面的最大方差阈本文档来自技高网...

【技术保护点】
1.一种基于太赫兹衍射图分解的太赫兹波数字全息成像方法，包括以下步骤：太赫兹全息图记录过程，全息图预传播过程，预传播衍射面复振幅图像分解过程，多层面复振幅图像重建合成过程：(1)太赫兹全息图记录过程。将深度方向上包含多层面信息的样品放置于太赫兹数字全息成像系统中，太赫兹波经过样品调制，携带样品的多层面信息，与不经过样品调制的太赫兹参考波干涉形成太赫兹全息图。为提高系统成像分辨率，搭载太赫兹探测器的二维平移台在垂直于太赫兹波照明方向的记录平面内按照给定扫描顺序进行二维移动，采集不同位置的太赫兹子全息图。移除样品后，搭载太赫兹探测器的二维平移台按照太赫兹子全息图记录过程中相同的扫描顺序，记录太赫兹子背景图像。将相同位置采集的对应太赫兹子全息图和太赫兹子背景图相除，得到子归一化全息图，按照扫描顺序进行子归一化全息图图像拼接，得到大孔径归一化太赫兹全息图。(2)全息图预传播过程。将步骤(1)中的大孔径归一化太赫兹全息图在计算机中通过自由空间传播算法，沿物体平面方向预传播至衍射平面P，获得预传播衍射面的振幅和相位分布，即预传播衍射面复振幅图像U(m，n)。该衍射面的传播距离为dp。定义样品不同分层面距离太赫兹探测器记录平面的最大距离和最小距离分别为dmax和dmin，则dp的取值定义为...

【技术特征摘要】
1.一种基于太赫兹衍射图分解的太赫兹波数字全息成像方法，包括以下步骤：太赫兹全息图记录过程，全息图预传播过程，预传播衍射面复振幅图像分解过程，多层面复振幅图像重建合成过程：(1)太赫兹全息图记录过程。将深度方向上包含多层面信息的样品放置于太赫兹数字全息成像系统中，太赫兹波经过样品调制，携带样品的多层面信息，与不经过样品调制的太赫兹参考波干涉形成太赫兹全息图。为提高系统成像分辨率，搭载太赫兹探测器的二维平移台在垂直于太赫兹波照明方向的记录平面内按照给定扫描顺序进行二维移动，采集不同位置的太赫兹子全息图。移除样品后，搭载太赫兹探测器的二维平移台按照太赫兹子全息图记录过程中相同的扫描顺序，记录太赫兹子背景图像。将相同位置采集的对应太赫兹子全息图和太赫兹子背景图相除，得到子归一化全息图，按照扫描顺序进行子归一化全息图图像拼接，得到大孔径归一化太赫兹全息图。(2)全息图预传播过程。将步骤(1)中的大孔径归一化太赫兹全息图在计算机中通过自由空间传播算法，沿物体平面方向预传播至衍射平面P，获得预传播衍射面的振幅和相位分布，即预传播衍射面复振幅图像U(m，n)。该衍射面的传播距离为dp。定义样品不同分层面距离太赫兹探测器记录平面的最大距离和最小距离分别为dmax和dmin，则dp的取值定义为(3)预传播衍射面复振幅图像分解过程。计算提取预传播衍射面复振幅图像U(m，n)的振幅和相位图像，采用最大类间方差法进行前景面和后景面的最大方差阈值T的判定计算，基于得到的阈值T进行预传播衍射面复振幅图像的二值化运算，得到二值化掩膜图像K(m，n)，可表示为对二值化掩膜图像K(m，n)完成边缘均值模糊化处理，得到预传播衍射面的掩膜图像K′(m，n)，通过非零值区域识别，对掩膜图像K′(m，n)中的x个孤立非零值区域图像进行逐一...

【专利技术属性】
技术研发人员：黄昊翀，杨雍青，仇佩瑶，郝思博，马媛媛，郑志远，
申请(专利权)人：中国地质大学北京，
类型：发明
国别省市：北京,11