CN106772365A - 一种基于贝叶斯压缩感知的多径利用穿墙雷达成像方法 - Google Patents
一种基于贝叶斯压缩感知的多径利用穿墙雷达成像方法 Download PDFInfo
- Publication number
- CN106772365A CN106772365A CN201611052793.7A CN201611052793A CN106772365A CN 106772365 A CN106772365 A CN 106772365A CN 201611052793 A CN201611052793 A CN 201611052793A CN 106772365 A CN106772365 A CN 106772365A
- Authority
- CN
- China
- Prior art keywords
- wall
- compressed sensing
- bayes
- radar
- multipath
- 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
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/887—Radar or analogous systems specially adapted for specific applications for detection of concealed objects, e.g. contraband or weapons
- G01S13/888—Radar or analogous systems specially adapted for specific applications for detection of concealed objects, e.g. contraband or weapons through wall detection
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了一种基于贝叶斯压缩感知的多径利用穿墙雷达成像方法。本发明采用单站单收发天线频率步进雷达***在平行于墙面多个位置发射脉冲串并获取回波信号,通过随机测量矩阵对回波信号进行压缩采样获得相对少量的测量信号,然后联合依据墙体多径反射特性设计好的字典构建压缩感知模型,使用贝叶斯压缩感知方法将墙后场景信息从测量信号中重构出来,最后利用重构出的信息完成穿墙雷达成像。本发明降低了成本,缩减了记录时间,同时贝叶斯压缩感知技术的使用使得计算速度更快、所需观测数据量更少,结合多径利用方法得到的图像精度更高且有效地抑制了虚影的出现。
Description
技术领域
本发明属于穿墙雷达成像技术领域,具体涉及一种通过运用贝叶斯压缩感知(Bayesian Compressive Sensing,BCS)原理并利用墙体多径反射现象处理雷达接收信号从而达到穿墙成像目的的方法。
背景技术
近年来,穿墙雷达成像技术在军、民用领域都收到了越来越广泛的关注。在建筑拆除、灾后救援和城市反恐作战等诸多场景中,都存在着对建筑物内或障碍物后的隐蔽目标探测的难题。传统的穿墙雷达***往往需要采用多个观测天线作为物理大孔径或者通过移动单个观测天线至多个观测位置合成一个大孔径来获得较大的距离分辨率。这样,***成本高且记录数据时间长。同时,巨大的原始数据采样量会极大地增加数据传输和存储负担。
针对上述问题,引入近些年兴起的压缩感知(Compressive Sensing,CS)技术是有效的解决办法。根据压缩感知原理,如果一个信号是稀疏的,在远小于奈奎斯特(Nyquist)采样率的条件下,用随机采样获取信号的离散样本,然后通过非线性重建算法就可以重构原信号。一般情况下,墙后场景信息在成像领域通常可被视为绝对稀疏的,因此CS技术能够将墙后场景信息从少量的测量样本中恢复出来。
尽管基于CS的成像方法能够有效地减少数据采样时间、数据量以及节省信号带宽,然而穿墙雷达成像技术的应用场景往往存在低信噪比的特点,并且在墙体多径反射现象的影响下,传统CS技术的成像质量不高,会出现虚假目标(或称“虚影”)。
发明内容
本发明的目的在于提供一种基于贝叶斯压缩感知的多径利用穿墙雷达成像方法,通过少量的测量数据得到精确、高分辨率、低噪声和虚影少的墙后场景图像。
实现本发明目的的技术解决方案为:一种基于贝叶斯压缩感知的多径利用穿墙雷达成像方法,采用单站单收发天线频率步进雷达***发射脉冲串并获取回波信号,然后通过基于贝叶斯压缩感知的多径利用方法处理回波信号并完成成像,具体实施步骤如下:
步骤1、将穿墙雷达置于初始位置即距墙体外墙面垂直距离dr处,整个方法共需要M个穿墙雷达的观测位置,该M个观测位置距墙体外墙面垂直距离均为dr;
步骤2、穿墙雷达的发射天线发射一串频率步进信号,该串频率步进信号包含N个频率步进脉冲,接收天线分别接收各频率下的回波信号zn,m,n∈{1,2,...,N},m∈{1,2,...,M};接收天线在第m个穿墙雷达观测位置获取的总回波信号表示为Zm=[z1,m,z2,m,...,zN,m]T;
步骤3、通过高斯随机测量矩阵Φm对回波信号Zm做压缩采样,得到测量向量Ym进行存储;
步骤4、沿平行于墙壁方向移动穿墙雷达至下一个观测位置,重复步骤2、3,继续移动穿墙雷达,直至获得M组测量信号为止;
步骤5、对成像区域进行网格化分,基于接收信号的相关域信息以及各网格像素与各观测位置间的传播时延构建M个字典Ψm;
步骤6、基于M组测量向量、M个字典和M个测量矩阵建立压缩感知模型;
步骤7、根据压缩感知模型,通过贝叶斯压缩感知算法将墙后图像信息矢量s重构,从而完成墙后场景成像。
本发明与现有技术相比,其显著优点为:(1)所需观测天线数量少、数据存储空间小、传输时间短,从而降低了成本,提高了效率;(2)BCS技术的使用使得该方法计算速度更快、所需观测数据量更少;(3)结合多径利用算法,得到的图像精度更高且有效地抑制了虚影的出现。
附图说明
图1为本发明基于贝叶斯压缩感知的多径利用穿墙雷达成像方法示意图。
图2为本发明基于贝叶斯压缩感知的多径利用穿墙雷达成像方法流程图。
图3为本发明基于贝叶斯压缩感知的多径利用穿墙雷达成像方法构成字典时成像区域网格划分示意图。
图4为本发明基于贝叶斯压缩感知的多径利用穿墙雷达成像方法成像多径几何关系模型图,其中,(a)是目标回波直接反射模型图,(b)是目标回波墙体内反射模型图,(c)是目标回波侧边内墙反射模型图。
图5为本发明基于贝叶斯压缩感知的多径利用穿墙雷达成像方法的实例效果图。
具体实施方式
下面结合附图和具体实施例进一步说明本发明基于贝叶斯压缩感知的多径利用穿墙雷达成像方法。
贝叶斯压缩感知算法(Bayesian Compressive Sensing,BCS)是用贝叶斯方法来处理压缩感知的信号重构问题,从统计的角度分析恢复问题,给原始信号赋予一个先验分布,由观测信号、稀疏基字典和观测矩阵恢复出原始稀疏信号的一个后验分布。与传统CS方法只能提供成像场景散射系数的点估计不同,该方法能够提供散射系数完整的后验密度估计,更加适合杂波环境;由于考虑了压缩采样过程中遇到的加性噪声,因而具有更好的鲁棒性;基于BCS方法得到的图像比基于传统CS方法得到的图像相比更加稀疏,且精度更高,计算速度更快,所需观测数量更少。因此,BCS方法相较于传统CS更加适用于穿墙雷达成像实际应用环境。
本发明是基于贝叶斯压缩感知的多径利用穿墙雷达成像方法,采用单站单收发天线频率步进雷达***发射脉冲串并获取回波信号,然后通过基于贝叶斯压缩感知的多径利用方法处理回波信号并完成成像。本发明应用场景构造如图1所示,按图2所示,具体实施步骤如下:
步骤1、按图1所示,将穿墙雷达置于初始位置即距墙体外墙面垂直距离dr处,整个方法共需要M个穿墙雷达的观测位置,该M个观测位置距墙体外墙面垂直距离均为dr;
步骤2、穿墙雷达的发射天线发射一串频率步进信号,该串频率步进信号包含N个频率步进脉冲,接收天线分别接收各频率下的回波信号zn,m,n∈{1,2,...,N},m∈{1,2,...,M};接收天线在第m个穿墙雷达观测位置获取的总回波信号表示为Zm=[z1,m,z2,m,...,zN,m]T。
所述频率步进信号的脉冲周期为Tr,载频起始频率为f1,频率步长为Δf,频率步进数为N,第n个脉冲信号频率为fn,信号带宽fB=fN-f1,雷达分辨率即为c/2fB,c为光在空气中的传播速度。这些参数可根据场景的实际需求进行设置。
所述接收天线接受的回波信号zn,m包含目标反射墙面多径反射及噪声vn,m:
其中,目标反射表示为:
其中,P表示场景中存在的目标总数,表示第p个目标的反射系数,τm,p表示第m个观测位置到第p个目标的往返传播时延,fn为频率步进信号中第n个脉冲信号频率;
墙面多径反射表示为:
其中,R表示墙反射路径总数,表示墙反射系数,βr表示与第r条路径相关的衰变系数,表示第r条反射路径的往返传播时延。
步骤3、通过高斯随机测量矩阵Φm对回波信号Zm做压缩采样,得到测量向量Ym进行存储,表达式为:
Ym=ΦmZm (4)
其中,Φm为L×N维的高斯随机矩阵,且L<<N,Zm为N×1的一维列向量,因此Ym是L×1的一维列向量。由于L<<N,Ym相较于Zm可以大大节省存储空间并降低传输所需时间。
步骤4、按图1所示,沿平行于墙壁方向移动穿墙雷达至下一个观测位置,重复步骤2、3,继续移动穿墙雷达,直至获得M组测量信号为止;
步骤5、对成像区域进行网格化分,基于接收信号的相关域信息以及各网格像素与各观测位置间的传播时延构建M个字典Ψm,具体为:
如图3所示,将墙后感兴趣的成像区域离散化为Q个网格的像素,网格距离单元的大小由成像***决定,通常被设为雷达分辨率。
考虑到回波信号是目标直接反射和各路径多径反射的叠加,这样,所有图像信息归入一个RQ×1的一维列向量s。该向量中的元素sq为一个加权指示函数,sq代表了第q点的反射系数大小,定义为:
R表示墙反射路径总数,则s为RQ×1的一维列向量;
Zm视为每个像素产生的回波信息叠加,利用每个像素点的特征设计字典Ψm,这样Zm和s之间的关系如下式表示:
Zm=Ψms (6)
其中,其中,字典Ψm由针对目标直接反射的子字典Ψm,t及针对各路径反射的子字典组成,即
Ψm,t=[ψm,t(n,q)]是一个N×Q矩阵,其第(n,q)个元素表达式为:
其中,τm,q表示第m个位置天线到第q个像素的理论往返传播时延,fn为频率步进信号中第n个脉冲信号频率;
是一个N×Q矩阵,表示第q个像素沿第r条路径反射的子字典,第(n,q)个元素表达式为:
其中,表示第m个观测位置到第q个像素沿第r条路径反射的理论往返传播时延;
τm,q、可根据图4所示穿墙雷达成像多径几何关系模型计算得出。对于图4(a)、(b)所示的模型,τm,q可计算为:
其中,lair1是电磁波从发射天线到前墙面的传播距离,lair2是电磁波从后墙面第q个像素的传播距离,lwall是电磁波在墙体内部的传播距离,c和v分别是电磁波在空气中和在墙体内的传播速度。lair1、lair2、lwall与入射角χmq、折射角以及第m个收发天线坐标(xm,-dr)、第q个像素坐标(xq,yq)之间的关系可写成如下表达式:
其中,式(10)中的dw为墙壁厚度,式(11)中的W为电磁波在墙体内往返的次数。墙体介电常数ε已知,根据Snell定律,有如下表达式:
联立式(9)-(13)即可求得τm,q的值。
计算图4(c)所示侧边内墙反射为例传播时延时,可将第q个像素对应侧边内墙做轴对称等效至其虚影,即坐标(2xw-xq,yq)处,然后用等效后的坐标替代原坐标(xq,yq)代入上述方法求解方程组即可得到
实际应用中,多径反射子字典的数量可以根据场景需求进行设置。一般情况下,只需考虑前墙体内一次反射和侧边内墙反射,其他路径反射往往由于衰减极大而可以忽略。
步骤6、基于M组测量向量、M个字典和M个测量矩阵建立压缩感知模型,具体为:
按照式(14)方式用M个高斯随机测量矩阵Φm构成总测量矩阵Φ,用M个字典Ψm构成总字典Ψ,用M个测量向量Ym构成总测量向量Y,这样根据压缩感知模型Y、Φ、Ψ及s的关系可由式(15)表示:
Y=ΦΨs=Θs (15)
其中,Θ为测量矩阵Φ和字典矩阵Ψ的积,称为传感矩阵;通过贝叶斯压缩感知方法求解Y=Θs获得墙后图像信息s。
步骤7、根据压缩感知模型,通过贝叶斯压缩感知算法将墙后图像信息矢量s重构,从而完成墙后场景成像。
从贝叶斯方法角度考虑压缩感知信号重构问题,图像信息矢量s具有稀疏先验特性,其稀疏度用K表示。通过计算观测信号Y在传感矩阵Θ下得到的重构图像信息矢量权重的后验概率。
Y为总测量向量、Θ为测量矩阵Φ和字典矩阵Ψ的积、s为图像信息矢量、ne为均值为0、方差为的高斯白噪声,贝叶斯压缩感知模型由式(15)修正重写为:
Y=Θs+ne (16)
Y的高斯似然模型为:
其中,K为图像信息矢量s的稀疏度;
至此,压缩感知问题转化为带有s的稀疏约束线性回归问题。Θ已知,则剩下工作是估计系数向量s和噪声方差
利用相对向量机(Relative Vector Machine,RVM)假设先验并估计参数。RVM采用层次先验模型促进稀疏性并满足共轭先验,对于s中的每个元素,假设均值为0的高斯先验:
其中,αi表示高斯密度函数的精确度(方差的倒数),αi先验为
对超参数α进行边缘分布计算得到s的先验为:
其中,符合Student-t分布,调整置信区间(a,b)使该分布在si=0获得峰值,所以此先验促进s的稀疏。
按照上述方法,为选择Γ(α0|c,d)先验。在假设式(19)先验的情况下,对于给定的Y和Θ,根据贝叶斯公式,s的后验是多变量的高斯分布,s的均值和方差分别为:
μ=α0∑ΘTY (20)
Σ=(α0ΘTΘ+Α)-1 (21)
其中,Α=diag(α1,α2,…αRQ)。参数α=(α1,α2,…αRQ)T和参数α0也要进行估计,通过第II类极大似然即Type-II ML(Type-II Maximum-Likelihood)逼近过程,计算α和α0得到:
其中,μi是依据式(20)计算的第i次的s均值,Σii是根据式(21)计算的第i次的s方差。
对s和α,α0进交替行迭代计算,直到收敛为止。
照此方法计算出矢量s后,按照s对应每个网格以Q为周期进行R次赋值,就能得到墙后场景的成像结果。
实施例1
下面举一实例说明,该实例通过Matlab软件仿真。例子中,仿真参数设置如下:墙壁厚度dw=0.15m,介电常数ε=6。雷达天线距离前墙面0.2m,以0.1m的间隔扫描从0m到5m间的感兴趣成像区域,即共有51个天线位置。发射信号为2~3GHz,1GHz带宽的步进频率信号(201个脉冲信号,步长5MHz)。1GHz带宽提供了0.15m的分辨率。将5m×5m的感兴趣成像区域划分为25×25的网格,网格数为625。考虑目标直接反射、墙体内一次反射和侧边内墙反射三条路径,分别就每条路径计算网格参数构建字典。测量数L设为10。三个目标分别位于(1.7m,0.9m),(2.5m,4.5m)和(4.1m,2.7m)。
仿真结果如图5所示。由图5可见,三个目标都被准确重构,并且除目标周围有少量可见噪点外,无其他虚影产生。由此可以证明本发明一种基于贝叶斯压缩感知的多径利用穿墙雷达成像方法的可行性和准确性。
本发明有效地解决了传统穿墙雷达SAR成像方法所需观测天线数量多、数据存储空间大、传输时间长的问题,从而降低成本,缩减记录时间。同时,相比于传统压缩感知方法,贝叶斯压缩感知技术的使用使得本发明中的方法计算速度更快、所需观测数据量更少,结合多径利用方法得到的图像精度更高且有效地抑制了虚影的出现。
Claims (7)
1.一种基于贝叶斯压缩感知的多径利用穿墙雷达成像方法,其特征在于,采用单站单收发天线频率步进雷达***发射脉冲串并获取回波信号,然后通过基于贝叶斯压缩感知的多径利用方法处理回波信号并完成成像,具体实施步骤如下:
步骤1、将穿墙雷达置于初始位置即距墙体外墙面垂直距离dr处,整个方法共需要M个穿墙雷达的观测位置,该M个观测位置距墙体外墙面垂直距离均为dr;
步骤2、穿墙雷达的发射天线发射一串频率步进信号,该串频率步进信号包含N个频率步进脉冲,接收天线分别接收各频率下的回波信号zn,m,n∈{1,2,...,N},m∈{1,2,...,M};接收天线在第m个穿墙雷达观测位置获取的总回波信号表示为Zm=[z1,m,z2,m,...,zN,m]T;
步骤3、通过高斯随机测量矩阵Φm对回波信号Zm做压缩采样,得到测量向量Ym进行存储;
步骤4、沿平行于墙壁方向移动穿墙雷达至下一个观测位置,重复步骤2、3,继续移动穿墙雷达,直至获得M组测量信号为止;
步骤5、对成像区域进行网格化分,基于接收信号的相关域信息以及各网格像素与各观测位置间的传播时延构建M个字典Ψm;
步骤6、基于M组测量向量、M个字典和M个测量矩阵建立压缩感知模型;
步骤7、根据压缩感知模型,通过贝叶斯压缩感知算法将墙后图像信息矢量s重构,从而完成墙后场景成像。
2.根据权利要求1所述的基于贝叶斯压缩感知的多径利用穿墙雷达成像方法,其特征在于,步骤2中所述频率步进信号的脉冲周期为Tr,载频起始频率为f1,频率步长为Δf,频率步进数为N,第n个脉冲信号频率为fn,信号带宽fB=fN-f1,雷达分辨率即为c/2fB,c为光在空气中的传播速度。
3.根据权利要求1所述的基于贝叶斯压缩感知的多径利用穿墙雷达成像方法,其特征在于,步骤2中所述接收天线接受的回波信号zn,m包含目标反射墙面多径反射及噪声vn,m:
其中,目标反射表示为:
其中,P表示场景中存在的目标总数,表示第p个目标的反射系数,τm,p表示第m个观测位置到第p个目标的往返传播时延,fn为频率步进信号中第n个脉冲信号频率;
墙面多径反射表示为:
其中,R表示墙反射路径总数,表示墙反射系数,βr表示与第r条路径相关的衰变系数,表示第r条反射路径的往返传播时延。
4.根据权利要求1所述的基于贝叶斯压缩感知的多径利用穿墙雷达成像方法,其特征在于,步骤3中所述通过高斯随机测量矩阵Φm对回波信号Zm做压缩采样,得到测量向量Ym进行存储,表达式为:
Ym=ΦmZm (4)
其中,Φm为L×N维的高斯随机矩阵,且L<<N,Zm为N×1的一维列向量,因此Ym是L×1的一维列向量。
5.根据权利要求1所述的基于贝叶斯压缩感知的多径利用穿墙雷达成像方法,其特征在于,步骤5所述对成像区域进行网格化分,基于接收信号的相关域信息以及各网格像素与各观测位置间的传播时延构建M个字典Ψm,具体为:
将墙后成像区域离散化为Q个网格的像素,网格距离单元的大小设为雷达分辨率,这样,所有图像像素信息归入一个一维列向量s,该向量中的元素sq为一个加权指示函数,sq代表了第q点的反射系数大小,定义为:
R表示墙反射路径总数,则s为RQ×1的一维列向量;
Zm视为每个像素产生的回波信息叠加,利用每个像素点的特征设计字典Ψm,这样Zm和s之间的关系如下式表示:
Zm=Ψms (6)
其中,字典Ψm由针对目标直接反射的子字典Ψmt及针对各路径反射的子字典组成,即
Ψm,t=[ψm,t(n,q)]是一个N×Q矩阵,其第(n,q)个元素表达式为:
其中,τm,q表示第m个位置天线到第q个像素的理论往返传播时延,fn为频率步进信号中第n个脉冲信号频率;
是一个N×Q矩阵,表示第q个像素沿第r条路径反射的子字典,第(n,q)个元素表达式为:
其中,表示第m个观测位置到第q个像素沿第r条路径反射的理论往返传播时延;
τm,q、根据墙壁特性及穿墙雷达成像多径几何关系计算得出,墙壁的介电常数、厚度分别为ε,dw。
6.根据权利要求1所述的基于贝叶斯压缩感知的多径利用穿墙雷达成像方法,其特征在于,步骤6所述基于M组测量向量、M个字典和M个测量矩阵建立压缩感知模型,具体为:
按照式(9)方式用M个高斯随机测量矩阵Φm构成总测量矩阵Φ,用M个字典Ψm构成总字典Ψ,用M个测量向量Ym构成总测量向量Y,这样压缩感知模型Y、Φ、Ψ及s的关系由式(10)表示:
Y=ΦΨs=Θs (10)
其中,Θ为测量矩阵Φ和字典矩阵Ψ的积,称为传感矩阵;通过贝叶斯压缩感知方法求解Y=Θs获得墙后图像信息s。
7.根据权利要求1所述的基于贝叶斯压缩感知的多径利用穿墙雷达成像方法,其特征在于,步骤7所述根据压缩感知模型,通过贝叶斯压缩感知算法将墙后图像信息矢量s重构,从而完成墙后场景成像,具体为:
Y为总测量向量、Θ为测量矩阵Φ和字典矩阵Ψ的积、s为图像信息矢量、ne为均值为0、方差为的高斯白噪声,贝叶斯压缩感知模型如下式所示:
Y=Θs+ne (11)
Y的高斯似然模型为:
其中,K为图像信息矢量s的稀疏度;
利用相对向量机假设先验并估计参数,相对向量机采用层次先验模型促进稀疏性并满足共轭先验,对于s中的每个元素si,假设均值为0的高斯先验:
其中,αi表示高斯密度函数的精确度,αi先验为
对超参数α进行边缘分布计算得到s的先验为:
其中,符合Student-t分布,调整置信区间(a,b)使该分布在si=0获得峰值,所以此先验促进s的稀疏;
按照上述方法,为选择Γ(α0|c,d)先验;在假设式(14)先验的情况下,对于给定的Y和Θ,根据贝叶斯公式,s的后验是多变量的高斯分布,s均值和方差分别为:
μ=α0ΣΘTY (15)
Σ=(α0ΘTΘ+Α)-1 (16)
其中,Α=diag(α1,α2,…αRQ),参数α=(α1,α2,…αRQ)T和参数α0也要进行估计,通过第II类极大似然即Type-II ML逼近过程,计算α和α0得到:
其中,μi是依据式(15)计算的第i次的s均值,Σii是根据式(16)计算的第i次的s方差;
对s和α、α0进交替行迭代计算,直到收敛为止;
计算出矢量s后,按照s对应每个网格以Q为周期进行R次赋值,就能得到墙后场景的成像结果。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611052793.7A CN106772365B (zh) | 2016-11-25 | 2016-11-25 | 一种基于贝叶斯压缩感知的多径利用穿墙雷达成像方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611052793.7A CN106772365B (zh) | 2016-11-25 | 2016-11-25 | 一种基于贝叶斯压缩感知的多径利用穿墙雷达成像方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106772365A true CN106772365A (zh) | 2017-05-31 |
CN106772365B CN106772365B (zh) | 2019-07-12 |
Family
ID=58912549
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201611052793.7A Expired - Fee Related CN106772365B (zh) | 2016-11-25 | 2016-11-25 | 一种基于贝叶斯压缩感知的多径利用穿墙雷达成像方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106772365B (zh) |
Cited By (24)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107479043A (zh) * | 2017-08-15 | 2017-12-15 | 电子科技大学 | 基于多成像字典的合成孔径穿墙雷达多径假目标剔除方法 |
CN107861125A (zh) * | 2017-11-01 | 2018-03-30 | 河海大学 | 基于广义帕累托分布的贝叶斯压缩感知成像方法 |
CN107894591A (zh) * | 2017-09-30 | 2018-04-10 | 沈阳航空航天大学 | 基于压缩感知的穿墙雷达衍射层析成像方法 |
CN107942314A (zh) * | 2017-11-22 | 2018-04-20 | 中南大学 | 基于lasso特征提取的多普勒穿墙雷达定位方法 |
CN108872980A (zh) * | 2018-06-19 | 2018-11-23 | 电子科技大学 | 一种基于窄带***的自适应穿墙成像方法 |
CN108896990A (zh) * | 2018-05-10 | 2018-11-27 | 桂林电子科技大学 | 一种利用耦合模式字典学习的建筑物墙体成像方法和装置 |
CN109085572A (zh) * | 2018-09-05 | 2018-12-25 | 西安电子科技大学昆山创新研究院 | 隧道内基于多径利用毫米波雷达的运动目标跟踪方法 |
CN109283530A (zh) * | 2018-09-14 | 2019-01-29 | 浙江大学 | 一种利用压缩感知提高微波成像线性度的方法 |
CN109471097A (zh) * | 2018-11-14 | 2019-03-15 | 湖南华诺星空电子技术有限公司 | 一种穿墙雷达信号优化处理方法及装置 |
CN109696671A (zh) * | 2018-12-29 | 2019-04-30 | 南京汇君半导体科技有限公司 | 基于组结构与相关性学习的高分辨穿墙雷达成像方法 |
CN109799499A (zh) * | 2019-01-28 | 2019-05-24 | 沈阳航空航天大学 | 一种穿墙雷达墙体参数估计方法 |
CN109917361A (zh) * | 2019-04-02 | 2019-06-21 | 电子科技大学 | 一种基于双基雷达的三维未知场景成像方法 |
CN110275158A (zh) * | 2018-03-15 | 2019-09-24 | 南京理工大学 | 基于贝叶斯压缩感知的宽带雷达回波信号参数估计方法 |
CN110568445A (zh) * | 2019-08-30 | 2019-12-13 | 浙江大学 | 一种轻量化卷积神经网络的激光雷达与视觉融合感知方法 |
CN111045006A (zh) * | 2019-12-23 | 2020-04-21 | 电子科技大学 | 一种基于多成像字典融合的拐角隐蔽目标成像方法 |
CN111624632A (zh) * | 2020-06-30 | 2020-09-04 | 青岛杰瑞自动化有限公司 | 一种卫星导航信号捕获方法及装置 |
CN111665500A (zh) * | 2020-06-12 | 2020-09-15 | 沈阳航空航天大学 | 一种基于单比特压缩感知的脉冲穿墙雷达成像方法 |
CN112444809A (zh) * | 2020-10-28 | 2021-03-05 | 浙江大学 | 一种基于智能反射面的分布式计算成像方法 |
CN112698326A (zh) * | 2020-11-30 | 2021-04-23 | 浙江华消科技有限公司 | 雷达穿墙探测对象的方法、雷达探测设备及电子装置 |
CN113109807A (zh) * | 2021-03-09 | 2021-07-13 | 桂林电子科技大学 | 基于压缩感知的频率分集阵列雷达地下目标三维成像方法 |
CN113965881A (zh) * | 2021-09-28 | 2022-01-21 | 浙江大学 | 遮挡效应下的毫米波一体化通信与感知方法 |
CN115291185A (zh) * | 2022-10-09 | 2022-11-04 | 南京理工大学 | 一种雷达目标的参数检测方法、装置及电子设备 |
CN117630924A (zh) * | 2023-11-30 | 2024-03-01 | 中国科学院空天信息创新研究院 | 一种基于先验信息的穿墙雷达三维成像方法和装置 |
US11979200B2 (en) | 2020-07-23 | 2024-05-07 | Camero-Tech Ltd. | System and a method for extracting low-level signals from hi-level noisy signals |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20050220310A1 (en) * | 2004-03-30 | 2005-10-06 | Mcgrath William R | Technique and device for through-the-wall audio surveillance |
CN103197289A (zh) * | 2013-04-03 | 2013-07-10 | 电子科技大学 | 一种适用于穿墙雷达成像的多径幻象抑制方法 |
-
2016
- 2016-11-25 CN CN201611052793.7A patent/CN106772365B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20050220310A1 (en) * | 2004-03-30 | 2005-10-06 | Mcgrath William R | Technique and device for through-the-wall audio surveillance |
CN103197289A (zh) * | 2013-04-03 | 2013-07-10 | 电子科技大学 | 一种适用于穿墙雷达成像的多径幻象抑制方法 |
Non-Patent Citations (2)
Title |
---|
V.H. TANG等: "Enhanced through-the-wall radar imaging using Bayesian compressive sensing", 《PROC. SPIE 8717, COMPRESSIVE SENSING ‖, 87170I》 * |
WU QISONG等: "Robust Multipath Exploitation Radar Imaging in Urban Sensing Based on Bayesian Compressive Sensing", 《IEEE, ASILOMAR》 * |
Cited By (37)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107479043A (zh) * | 2017-08-15 | 2017-12-15 | 电子科技大学 | 基于多成像字典的合成孔径穿墙雷达多径假目标剔除方法 |
CN107479043B (zh) * | 2017-08-15 | 2020-10-30 | 电子科技大学 | 基于多成像字典的合成孔径穿墙雷达多径假目标剔除方法 |
CN107894591A (zh) * | 2017-09-30 | 2018-04-10 | 沈阳航空航天大学 | 基于压缩感知的穿墙雷达衍射层析成像方法 |
CN107861125A (zh) * | 2017-11-01 | 2018-03-30 | 河海大学 | 基于广义帕累托分布的贝叶斯压缩感知成像方法 |
CN107861125B (zh) * | 2017-11-01 | 2021-09-03 | 河海大学 | 基于广义帕累托分布的贝叶斯压缩感知成像方法 |
CN107942314A (zh) * | 2017-11-22 | 2018-04-20 | 中南大学 | 基于lasso特征提取的多普勒穿墙雷达定位方法 |
CN107942314B (zh) * | 2017-11-22 | 2021-06-04 | 中南大学 | 基于lasso特征提取的多普勒穿墙雷达定位方法 |
CN110275158A (zh) * | 2018-03-15 | 2019-09-24 | 南京理工大学 | 基于贝叶斯压缩感知的宽带雷达回波信号参数估计方法 |
CN108896990A (zh) * | 2018-05-10 | 2018-11-27 | 桂林电子科技大学 | 一种利用耦合模式字典学习的建筑物墙体成像方法和装置 |
CN108896990B (zh) * | 2018-05-10 | 2022-06-03 | 桂林电子科技大学 | 一种利用耦合模式字典学习的建筑物墙体成像方法和装置 |
CN108872980B (zh) * | 2018-06-19 | 2022-06-03 | 电子科技大学 | 一种基于窄带***的自适应穿墙成像方法 |
CN108872980A (zh) * | 2018-06-19 | 2018-11-23 | 电子科技大学 | 一种基于窄带***的自适应穿墙成像方法 |
CN109085572A (zh) * | 2018-09-05 | 2018-12-25 | 西安电子科技大学昆山创新研究院 | 隧道内基于多径利用毫米波雷达的运动目标跟踪方法 |
CN109283530A (zh) * | 2018-09-14 | 2019-01-29 | 浙江大学 | 一种利用压缩感知提高微波成像线性度的方法 |
CN109283530B (zh) * | 2018-09-14 | 2020-08-14 | 浙江大学 | 一种利用压缩感知提高微波成像线性度的方法 |
CN109471097B (zh) * | 2018-11-14 | 2022-04-22 | 湖南华诺星空电子技术有限公司 | 一种穿墙雷达信号优化处理方法及装置 |
CN109471097A (zh) * | 2018-11-14 | 2019-03-15 | 湖南华诺星空电子技术有限公司 | 一种穿墙雷达信号优化处理方法及装置 |
CN109696671A (zh) * | 2018-12-29 | 2019-04-30 | 南京汇君半导体科技有限公司 | 基于组结构与相关性学习的高分辨穿墙雷达成像方法 |
CN109799499A (zh) * | 2019-01-28 | 2019-05-24 | 沈阳航空航天大学 | 一种穿墙雷达墙体参数估计方法 |
CN109799499B (zh) * | 2019-01-28 | 2023-04-28 | 沈阳航空航天大学 | 一种穿墙雷达墙体参数估计方法 |
CN109917361A (zh) * | 2019-04-02 | 2019-06-21 | 电子科技大学 | 一种基于双基雷达的三维未知场景成像方法 |
CN109917361B (zh) * | 2019-04-02 | 2023-04-25 | 电子科技大学 | 一种基于双基雷达的三维未知场景成像方法 |
CN110568445A (zh) * | 2019-08-30 | 2019-12-13 | 浙江大学 | 一种轻量化卷积神经网络的激光雷达与视觉融合感知方法 |
CN111045006A (zh) * | 2019-12-23 | 2020-04-21 | 电子科技大学 | 一种基于多成像字典融合的拐角隐蔽目标成像方法 |
CN111045006B (zh) * | 2019-12-23 | 2022-10-14 | 电子科技大学 | 一种基于多成像字典融合的拐角隐蔽目标成像方法 |
CN111665500A (zh) * | 2020-06-12 | 2020-09-15 | 沈阳航空航天大学 | 一种基于单比特压缩感知的脉冲穿墙雷达成像方法 |
CN111624632B (zh) * | 2020-06-30 | 2023-04-11 | 青岛杰瑞自动化有限公司 | 一种卫星导航信号捕获方法及装置 |
CN111624632A (zh) * | 2020-06-30 | 2020-09-04 | 青岛杰瑞自动化有限公司 | 一种卫星导航信号捕获方法及装置 |
US11979200B2 (en) | 2020-07-23 | 2024-05-07 | Camero-Tech Ltd. | System and a method for extracting low-level signals from hi-level noisy signals |
CN112444809A (zh) * | 2020-10-28 | 2021-03-05 | 浙江大学 | 一种基于智能反射面的分布式计算成像方法 |
CN112444809B (zh) * | 2020-10-28 | 2023-04-25 | 浙江大学 | 一种基于智能反射面的分布式计算成像方法 |
CN112698326A (zh) * | 2020-11-30 | 2021-04-23 | 浙江华消科技有限公司 | 雷达穿墙探测对象的方法、雷达探测设备及电子装置 |
CN113109807A (zh) * | 2021-03-09 | 2021-07-13 | 桂林电子科技大学 | 基于压缩感知的频率分集阵列雷达地下目标三维成像方法 |
CN113965881A (zh) * | 2021-09-28 | 2022-01-21 | 浙江大学 | 遮挡效应下的毫米波一体化通信与感知方法 |
CN115291185A (zh) * | 2022-10-09 | 2022-11-04 | 南京理工大学 | 一种雷达目标的参数检测方法、装置及电子设备 |
CN117630924A (zh) * | 2023-11-30 | 2024-03-01 | 中国科学院空天信息创新研究院 | 一种基于先验信息的穿墙雷达三维成像方法和装置 |
CN117630924B (zh) * | 2023-11-30 | 2024-06-07 | 中国科学院空天信息创新研究院 | 一种基于先验信息的穿墙雷达三维成像方法和装置 |
Also Published As
Publication number | Publication date |
---|---|
CN106772365B (zh) | 2019-07-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106772365A (zh) | 一种基于贝叶斯压缩感知的多径利用穿墙雷达成像方法 | |
CN109116311B (zh) | 基于知识辅助稀疏迭代协方差估计的杂波抑制方法 | |
US8193967B2 (en) | Method and system for forming very low noise imagery using pixel classification | |
CN105911544B (zh) | 一种基于压缩感知技术的扩频穿墙雷达成像方法 | |
US7796829B2 (en) | Method and system for forming an image with enhanced contrast and/or reduced noise | |
CN107340518B (zh) | 一种用于信号缺失下的isar雷达成像方法 | |
CN105353355B (zh) | 一种基于稀疏重构和投影成像的多基地雷达多目标定位方法 | |
CN108845316A (zh) | 一种基于压缩感知理论的雷达稀疏探测方法 | |
CN105467388A (zh) | 低信噪比下压缩感知雷达高分辨成像设备及其成像方法 | |
CN108226927A (zh) | 基于加权迭代最小稀疏贝叶斯重构算法的sar成像方法 | |
CN104076360B (zh) | 基于压缩感知的二维sar稀疏目标成像方法 | |
CN109696672B (zh) | 一种基于空间结构关联性的高分辨率穿墙雷达成像方法 | |
CN112198506A (zh) | 一种超宽带穿墙雷达学习成像的方法、装置、***和可读存储介质 | |
CN111796272A (zh) | 穿墙雷达人体图像序列的姿态实时识别方法及计算机设备 | |
CN108896990A (zh) | 一种利用耦合模式字典学习的建筑物墙体成像方法和装置 | |
CN109298417B (zh) | 一种基于雷达信号处理的建筑物内部结构探测方法及装置 | |
Ren et al. | 3D Imaging Algorithm for Down‐Looking MIMO Array SAR Based on Bayesian Compressive Sensing | |
CN109856636A (zh) | 曲线合成孔径雷达自适应三维成像方法 | |
CN106291500B (zh) | 基于fpga的片上***及其目标检测方法 | |
Dai et al. | Scattering simulation and reconstruction of a 3-D complex target using downward-looking step-frequency radar | |
CN108983192B (zh) | 基于gps辐射源的雷达运动目标参数估计方法 | |
CN112835006A (zh) | 一种基于帧间积累的跟踪雷达海上小目标检测方法及*** | |
Wei et al. | Sparse reconstruction for linear array SAR 3-D imaging based on Bayesian estimation | |
CN116256701B (zh) | 一种基于深度学习的探地雷达互扰波抑制方法及*** | |
Khamlichi et al. | Using CFAR algorithm to further improve a combined through-wall imaging method |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190712 Termination date: 20201125 |
|
CF01 | Termination of patent right due to non-payment of annual fee |