一种基于LAB颜色空间映射的幅度相位时频图表达方法技术

本发明专利技术公开了一种基于LAB颜色空间映射的幅度相位时频图表达方法，包括步骤如下：对时域信号进行短时傅里叶变换后取幅度值和相位值，从而获得初始的幅度矩阵和相位矩阵；对初始的幅度矩阵进行分段线性压扩处理，得到第一幅度矩阵；对初始的相位矩阵做差分正弦余弦变换，得到第一相位矩阵、第二相位矩阵，获得一组由第一幅度矩阵、第一相位矩阵、第二相位矩阵构成的三通道数据；将获得的三通道数据映射到LAB域，获得LAB颜色空间下的三通道时频域矩阵L、A、B；将三通道时频域矩阵L、A、B转换到RGB域生成可视化的时频图。本发明专利技术生成的时频图能够在具备幅度谱和相位谱信息的同时，抑制噪声对时频图的影响，强化时频图中有用信号部分的色彩表达，提升时频图可视化效果。提升时频图可视化效果。提升时频图可视化效果。

【技术实现步骤摘要】
一种基于LAB颜色空间映射的幅度相位时频图表达方法


[0001]本专利技术涉及无线电信号处理
，更具体的，涉及一种基于LAB颜色空间映射的幅度相位时频图表达方法。

技术介绍

[0002]时频图是通过对时域信号作时频变换(如短时傅立叶变换、小波变换等)所获得的时频域二维矩阵进行平面可视化绘制得到的二维图像。时频图描述了信号频谱成分的幅值和相位随时间变化的关系，其颜色及纹理的变化规律是对信号的一种极为有效的表达。在(通信、雷达、导航等)无线电信号处理领域，采用检测与估计理论结合人工智能等技术对无线电信号时频图进行建模分析，可以用于解决无线电信号的检测、识别、降噪、分离、参数估计等大量信号处理问题。最典型的时频图是只携带幅度(或能量)信息的幅度(或能量)谱，它的平面化二维图像是单层(单通道)的灰度图。另一种常见的时频图则是仅携带相位信息的相位谱单通道灰度图。此外，也有将幅度谱、能量谱和相位谱进行一同或部分组合得到RGB表达模式下的(双通道)伪彩色图或(三通道)彩色图。
[0003]目前主流的时频图生成方法是对时域信号进行时频变换后得到的复数矩阵进行求幅度、能量和相位从而得到幅度矩阵、能量矩阵和相位矩阵，通过直接将这三个矩阵中的元素值映射为灰度值来得到单幅度谱、能量谱或相位谱的时频灰度图。或者是将这三个矩阵分别独立地量化后映射到RGB三通道，以获得RGB模式的时频彩色图。
[0004]现有的生成单幅度谱(能量谱)或相位谱的时频灰度图仅考虑了信号时频变换结果的幅度(能量)信息或相位信息。而有的方法虽然通过映射幅度(能量)和相位信息到RGB三通道以获得时频彩色图，但幅度(能量)及相位信息的量化和映射通常是独立分开的，并且是生硬地将幅度及相位与RGB的三个通道进行对应。RGB三通道间的对等关系导致时频图的呈现缺乏物理意义上的可解释性，此外也没有考虑到对时频图的噪声抑制及可视化呈现效果的优化。

技术实现思路

[0005]本专利技术为了解决以上现有技术在时频图生成过程中，对信号时频变换的幅度信息和相位信息利用不充分、不合理等问题，提供了一种基于LAB颜色空间映射的幅度相位时频图表达方法，其使生成的时频图能够在具备幅度谱和相位谱信息的同时，抑制噪声对时频图的影响，强化时频图中有用信号部分的色彩表达，提升时频图可视化效果。
[0006]为实现上述本专利技术目的，采用的技术方案如下：
[0007]一种基于LAB颜色空间映射的幅度相位时频图表达方法，所述的方法包括步骤如下：
[0008]S1：对时域信号进行短时傅里叶变换后取幅度值和相位值，获得初始的幅度矩阵和相位矩阵；
[0009]S2：对初始的幅度矩阵进行分段线性压扩处理，得到第一幅度矩阵；对初始的相位
矩阵做差分正弦余弦变换，得到第一相位矩阵、第二相位矩阵，从而获得一组由第一幅度矩阵、第一相位矩阵、第二相位矩阵构成的三通道数据；
[0010]S3：将获得的三通道数据映射到LAB域，获得LAB颜色空间下的三通道时频域矩阵L、A、B；
[0011]S4：将三通道时频域矩阵L、A、B转换到RGB域生成可视化的时频图。
[0012]优选地，S1，具体地，
[0013]将无线电的时域信号表示为其中L是信号的样本长度；时域信号的短时傅立叶变换表示为复数时频矩阵其中T和F对应于时间帧和频率点的数量；
[0014]则时域信号的时频矩阵所对应初始的幅度矩阵和相位矩阵分别表示为：
[0015][0016]P＝∠S∈[
‑
π,π]T
×
F
[0017]其中，∠为求弧度角运算、M表示初始的幅度矩阵、P表示初始的相位矩阵。
[0018]进一步地，对初始的幅度矩阵进行分段线性压扩处理，得到第一幅度矩阵，具体如下：
[0019]对初始的幅度矩阵M的矩阵元素进行统计，选取第一阈值M
max
使得M
max
不小于98％初始的幅度矩阵M的矩阵元素值；
[0020]对于初始的幅度矩阵M的矩阵元素进行限幅的归一化，
[0021][0022]其中，m
t,f
为初始的幅度矩阵M的矩阵元素，t∈{1,
…
,T},f∈{1,
…
,F}；m
′
t,f
为归一化后幅度矩阵M
normal
∈[0,1]T
×
F
的矩阵元素；
[0023]根据噪声估计的结果确定划分时频变换幅度值中噪声成分和有用信号成分的第二阈值M
th
；
[0024]将第二阈值以下判定为噪声，对噪声做线性压缩处理；第二阈值以上判定为有用信号，对有用信号做线性扩展处理；由此得到分段线性压扩处理后的第一幅度矩阵。
[0025]再进一步地，所述的分段线性压扩处理的函数公式表达式如下：
[0026][0027]其中，k∈(0,1)为线性压缩部分的斜率，m
″
t,f
表示第一幅度矩阵M
compand
中第t行f列的矩阵元素；M
compand
∈[0,1]T
×
F
表示归一化后幅度矩阵M
normal
进行分段线性压扩后的第一幅度矩阵。
[0028]再进一步地，对初始的相位矩阵做差分正弦余弦变换，具体如下：
[0029]采用在频率维度对初始的相位矩阵P进行一阶前向差分变换来增大目标信号与噪声在相位谱图上的区分度，获得差分相位矩阵ΔP∈[
‑
2π,2π]T
×
F
，即
[0030]Δp
t,f
＝p
t,f+1
‑
p
t,f
[0031]其中，Δp
t,f
为差分相位矩阵ΔP的矩阵元素，p
t,f
为初始相位矩阵P的矩阵元素，t∈{1,
…
,T},f∈{1,
…
,F}；
[0032]对差分相位矩阵进行三角变换以消除
‑
2π、+2π的弧度角跳变；所述的三角变换后的相位矩阵表示为：
[0033]P
sin
＝sin(ΔP)∈[
‑
1,1]T
×
F
[0034]P
cos
＝cos(ΔP)∈[
‑
1,1]T
×
F
[0035]其中，cos(
·
)表示三角余弦变换，sin(
·
)表示三角正弦变换，P
sin
表示通过三角正弦变换后得到的第一相位矩阵、P
cos
表示通过三角余弦变换后得到的第二相位矩阵。
[0036]再进一步地，将获得的三通道数据映射到LAB域，具体如下：
[0037]结合L、A、B三通道，以及第本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于LAB颜色空间映射的幅度相位时频图表达方法，其特征在于：所述的方法包括步骤如下：S1：对时域信号进行短时傅里叶变换后取幅度值和相位值，从而获得初始的幅度矩阵和相位矩阵；S2：对初始的幅度矩阵进行分段线性压扩处理，得到第一幅度矩阵；对初始的相位矩阵做差分正弦余弦变换，得到第一相位矩阵、第二相位矩阵，从而获得一组由第一幅度矩阵、第一相位矩阵、第二相位矩阵构成的三通道数据；S3：将获得的三通道数据映射到LAB域，获得LAB颜色空间下的三通道时频域矩阵L、A、B；S4：将三通道时频域矩阵L、A、B转换到RGB域生成可视化的时频图。2.根据权利要求2所述的基于LAB颜色空间映射的幅度相位时频图表达方法，其特征在于：S1，具体地，将无线电的时域信号表示为其中L是信号的样本长度；时域信号的短时傅立叶变换表示为复数时频矩阵其中T和F对应于时间帧和频率点的数量；则时域信号的时频矩阵所对应初始的幅度矩阵和相位矩阵分别表示为：P＝∠S∈[
‑
π,π]
T
×
F
其中，∠为求弧度角运算、M表示初始的幅度矩阵、P表示初始的相位矩阵。3.根据权利要求2所述的基于LAB颜色空间映射的幅度相位时频图表达方法，其特征在于：对初始的幅度矩阵进行分段线性压扩处理，得到第一幅度矩阵，具体如下：对初始的幅度矩阵M的矩阵元素进行统计，选取第一阈值M
max
使得M
max
不小于98％初始的幅度矩阵M的矩阵元素值；对于初始的幅度矩阵M的矩阵元素进行限幅的归一化，其中，m
t,f
为初始的幅度矩阵M的矩阵元素，t∈{1,
…
,T},f∈{1,
…
,F}；m
′
t,f
为归一化后幅度矩阵M
nomal
∈[0，1]
T
×
F
的矩阵元素；根据噪声估计的结果确定划分时频变换幅度值中噪声成分和有用信号成分的第二阈值M
th
；将第二阈值以下判定为噪声，对噪声做线性压缩处理；第二阈值以上判定为有用信号，对有用信号做线性扩展处理；由此得到分段线性压扩处理后的第一幅度矩阵。4.根据权利要求3所述的基于LAB颜色空间映射的幅度相位时频图表达方法，其特征在于：所述的分段线性压扩处理的函数公式表达式如下：
其中，k∈(0，1)为线性压缩部分的斜率，m
″
t，f
表示第一幅度矩阵M
compand
中第t行f列的矩阵元素；M
compand
∈[0，1]
T
×
F
表示归一化后幅度矩阵M
nomal
进行分段线性压扩的第一幅度矩阵。5.根据权利要求4所述的基于LAB颜色空间映射的幅度相位时频图表达方法，其特征在于：对初始的相位矩阵做差分正弦余弦变换，得到第一相位矩阵、第二相位矩阵，具体如下：采用在频率维度对初始的相位矩阵P进行一阶前向差分变换来增大目标信号与噪声在相位谱图上的区分度，获得差分相位矩阵ΔP∈[
‑
2π，2π]
T
×
F
，即Δp
t，f
＝p
t，f+1
‑
p
t，f
其中，Δp
t，f
为差分相位矩阵ΔP的矩阵元素，p
t，f
为初始相位矩阵P的矩阵元素，t∈{1，...，T}，f∈{1，...，F}；对差分相位矩阵进行三角变换以消除
‑
2π、+2π的弧度角跳变；所述的三角变换后的相位矩阵表示为：P
sin
＝sin(ΔP)∈[
‑
1，1]
T
×
F
P
cos...

【专利技术属性】
技术研发人员：李扬清，张志龙，宋沅隆，洪江涛，陈章烜，陈啸宇，高秋雅，唐辉敏，
申请(专利权)人：北京邮电大学，
类型：发明
国别省市：