(19)中华 人民共和国 国家知识产权局 (12)发明 专利 (10)授权公告 号 (45)授权公告日 (21)申请 号 202111638247.2 (22)申请日 2021.12.3 0 (65)同一申请的已公布的文献号 申请公布号 CN 114004833 A (43)申请公布日 2022.02.01 (73)专利权人 首都师范大学 地址 100089 北京市海淀区西三环北路10 5 号 (72)发明人 李志龙　左剑　张存林　赵源萌　 杨墨轩　 (74)专利代理 机构 北京超凡宏宇专利代理事务 所(特殊普通 合伙) 11463 代理人 彭星 (51)Int.Cl. G06T 7/00(2017.01)G06T 7/13(2017.01) G06T 5/00(2006.01) G06T 5/10(2006.01) G06T 5/20(2006.01) 审查员 董泽华 (54)发明名称 复合材料太赫兹成像分辨力增强方法、 装 置、 设备和介质 (57)摘要 本申请提供了一种复合材料太赫兹成像分 辨力增强方法、 装置、 设备和介质， 该方法包括： 采用太赫兹成像系统， 对被测复合材料进行成像 检测， 得到被测 复合材料的第一一维距离像； 对 第一一维距离像去噪， 获得去噪后的第二一维距 离像； 将第二一维距离像、 细节增强后的第二一 维距离像分别作为第一引导图、 第二引导图， 对 第二一维距离像 分别进行引导滤波， 分别得到细 节层、 基础层； 对细节与基础层分别进行第一增 益和第二增益， 将进行第一增益后的细节层和进 行第二增益后的基础层叠加为第三一维距离像； 对第三一维距离像进行幅值 成像， 得到成像后的 结果图像。 本申请实施例通过上述方法， 能够在 有效抑制背景噪声的同时， 增强图像的边缘等细 节信息。 权利要求书4页 说明书18页 附图3页 CN 114004833 B 2022.04.01 CN 114004833 B 1.一种复合材 料太赫兹成像分辨力增强方法， 其特 征在于， 所述方法包括： 采用太赫兹成像系统， 对被测复合材料进行成像检测， 以得到针对所述被测复合材料 的第一一维距离像； 对所述第一 一维距离像进行高频去噪， 获得去噪后的第二 一维距离像； 将所述第二一维距离像作为第 一引导图， 对所述第 二一维距离像进行引导滤波得到细 节层， 以及， 将细节增强后的所述第二一 维距离像作为第二引导图， 对所述第二一 维距离像 进行引导滤波得到基础 层； 对所述细节层与所述基础层分别进行第 一增益和第 二增益， 将进行所述第 一增益后的 细节层和进行 所述第二增益后的基础 层叠加为第三 一维距离像； 对所述第三 一维距离像进行幅值成像， 得到成像后的结果图像。 2.根据权利要求  1 所述的方法， 其特征在于， 所述采用太赫兹成像系统， 对被测复合 材料进行成像 检测， 以得到针对所述被测复合材 料的第一 一维距离像， 包括： 通过所述太赫兹成像系统对所述被测复合材料进行成像检测， 以获取所述被测复合材 料在不同界面的反射回波的至少一个幅值信息； 根据所述至少一个幅值信息得到所述被测复合材 料在太赫兹时域的检测信号； 对所述检测信号进行快速傅里叶变换， 得到所述第 一一维距离像； 其中， 所述快速傅里 叶变换的公式为： ； 式中， 表征快速傅里叶变换， 为所述检测信号， 为所述第一 一维距离像。 3.根据权利要求  1 所述的方法， 其特征在于， 对所述第一一维距离像进行高频去噪， 获得去噪后的第二 一维距离像， 包括： 对所述第一一维距离像进行离散小波变换， 以根据小波母函数将所述第 一一维距离像 分解为多个尺度下 的第一变换结果； 其中， 所述第一一维距离像与所述第一变换结果的第 一关系为： ； 式中， 为预设的离散化伸缩指数， 为预设的离散化平移系数， t为太赫兹时域的检 测信号中的时间变量， 为所述第一变换结果， 为所述离散化伸缩指数为 的尺度 因子， 为所述离散化伸缩指数为离散化伸缩指数为 的尺度因子，   为小波母 函数， 为小波基函数， 所述小波基函数是基于所述离散化伸缩指数、 所述离散 化平移系数和所述时间变量下的小波母函数， 为所述小波基函数的共轭； 为 所述第一 一维距离像； 根据尺度函数将所述第一一维距离像分解为多个尺度下的第二变换结果； 其中， 所述 第一一维距离像与所述第二变换 结果的第二关系为：权　利　要　求　书 1/4 页 2 CN 114004833 B 2； 式中， 为所述第二变换结果， 为尺度函数， 为所述尺度函数经 过 倍伸缩和 倍平移之后得到的第一尺度函数； 为所述第一尺度函数的共 轭； 根据多分辨率分析方程确定多个尺度 下所述小波基函数与所述尺度函数的第 三关系； 其中， 所述第三关系为： ； 式中， 为平移倍数， 为所述尺度函数 在 尺度下进行 倍 平移后的第二尺度函数； 为所述小波基函数对应的第一小波系数； 为j尺度 下的小波函数与j+1尺度下 的所述尺度函数建立关系时的对应平移变化量； 小波函数包括 小波基函数、 小 波母函数； 根据所述第一关系、 所述第 二关系、 所述第三关系， 得到多个尺度 下所述第 一变换结果 与目标第二变换 结果的第四关系； 其中， 所述第四关系为： ； 式中， 为所述第二尺度函数对应的目标第二变换 结果； 根据所述第四关系， 确定多个尺度下每 个所述第一变换 结果对应的第一小 波系数； 根据预设阈值对由所有所述第一小波系数构成的第一小波系数集合进行硬阈值函数 收缩去噪， 以将 每个尺度下的所述第一小波系数更改为对应的第二小波系数； 其中， 通过以 下公式进行硬阈值 函数收缩去噪： ； 式中， 为所述预设阈值； 为所述第二小波系数； 为所述第 一小波系数的模； 所述预设阈值是多个尺度下 所述第一小 波系数的标准偏差； 确定所述第二小 波系数中值 不为零的至少一个目标第二小 波系数； 根据所述第四关系， 确定所述至少一个目标第 二小波系数中每个所述目标第 二小波系 数对应的目标第一变换 结果；权　利　要　求　书 2/4 页 3 CN 114004833 B 3