CN102129067A - 一种前向散射雷达地面运动目标信号建模及成像方法 - Google Patents
一种前向散射雷达地面运动目标信号建模及成像方法 Download PDFInfo
- Publication number
- CN102129067A CN102129067A CN2010105879799A CN201010587979A CN102129067A CN 102129067 A CN102129067 A CN 102129067A CN 2010105879799 A CN2010105879799 A CN 2010105879799A CN 201010587979 A CN201010587979 A CN 201010587979A CN 102129067 A CN102129067 A CN 102129067A
- Authority
- CN
- China
- Prior art keywords
- prime
- target
- centerdot
- ground
- receiver
- 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
Images
Landscapes
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明涉及一种运动目标前向散射信号建模及成像的信号处理方法,属于雷达信号处理技术领域。本发明所提出的方法考虑地面反射、大衍射角以及目标运动方向较为随机等实际情况,建立了精确的地面运动目标前向散射信号模型,修正了SISAR成像算法,并对目标像进行了更为精确的估计。修正后的成像算法在大衍射角、目标斜穿基线和多径干涉情况下依然可以获得精确的目标像,这些目标像与目标真实轮廓形状非常接近,可有效地对不同形状的目标进行分类识别。
Description
技术领域
本发明涉及一种运动目标前向散射信号建模及成像的信号处理方法,尤其涉及前向散射雷达地面运动目标信号建模及对运动目标侧影轮廓进行成像,属于雷达信号处理技术领域。
背景技术
前向散射雷达(FSR)是一种特殊的双基地雷达,其工作区域很窄,主要集中于双基地角大于140度的前向散射区域内。前向散射雷达具有许多独特的优势:前向散射的RCS迅速增加,而且与是否涂有吸波涂层无关,具有较强的反隐身能力;雷达结构简单,功率低;接收信号起伏小,具有较长的相关时间等等。
2000年提出的阴影逆合成孔径雷达(SISAR)理论为前向散射雷达运动目标成像与识别提供了新的途径。这种新的成像方法是利用目标的阴影信号对运动目标成像,也称为SISAR成像技术。获得的目标像是目标阴影轮廓的中线和高度差,它们的形状与目标轮廓的真实外形非常相像,可以直观地对目标进行分类识别。由于目标阴影信号的强弱与目标是否涂有吸波涂层无关,因此SISAR成像技术在识别隐身目标方面独具优势。
然而目前的SISAR成像算法是以空中运动目标为研究对象,基于如下假设和近似建立的:1、假设目标尺寸远大于工作波长,系统工作于光学区,进而采用光学近似的方法;2、假设天线波束很窄,指向空中目标,进而忽略地面反射采用自由空间模型;3、假设收发基线很长,相参积累时间很短,目标始终处于小衍射角区域,忽略了目标多普勒频率高次分量的影响;4、假设目标近似垂直穿越收发基线,忽略了目标斜穿基线的影响。当上述假设和近似条件不成立时,目前的SISAR成像算法将不再精确,甚至无法使用。
本发明的应用背景是利用前向散射雷达对地面运动目标进行监视与识别。由于前向散射雷达在低空目标和隐身目标检测识别等方面独特的优势,使得地面前向散射雷达系统具有广泛的应用需求。典型的地面前向散射雷达系统具有体积小、重量轻、功耗小等优点,非常便于地面交通工具或无人机(UAV)的运送、投放和布阵,可组成微型无线传感器网络,实现地面交通监视或战场状态感知。该系统对于常规小目标(如敌方士兵)以及具有隐身能力的地面运动目标(如隐身坦克或者军用卡车等),均可进行有效地探测和识别。此外,对于编队行进中的部队、车队以及坦克群等的有效识别,可以很好的了解敌军的行军路线、动态、实力等重要信息,为战场态势的感知提供重要信息,具有非常有价值的军事应用潜力。由于在地面系统中,传统成像算法的假设条件不再满足,故传统的SISAR成像算法不再适用。
发明内容
本发明的目的是对前向散射雷达系统地面运动目标前向散射雷达信号建模并利用前向散射雷达系统对地面运动目标侧影轮廓成像。解决大衍射角情况下或目标斜穿基线时SISAR成像结果失真的问题以及地面反射信号引起多径干涉致使无法成像的问题,以获得精确的地面运动目标侧影轮廓像。
一种前向散射雷达地面运动目标信号建模及成像方法包括如下具体步骤:
步骤1、建立地面前向散射雷达系统的几何结构模型,如图1所示。
xoy平面是地平面,发射机置于坐标系(x,y,z)的原点O,接收机在y轴B(0,L,0)处。发射机和接收机之间的距离为基线长度L。假设目标以速度v平行于xoy平面运动,与基线夹角为αh,βh分别表示从发射机和接收机方向观测的目标方位角。发射机到目标穿越基线位置的距离记为dT,接收机到目标穿越基线位置的距离记为dR。坐标系(x′,y′,z′)与坐标系(x,y,z)各轴相互平行,原点P(xp,yp,zp)为目标的中心。T,R,T′,R′分别表示发射机、接收机、镜像发射机和镜像接收机天线相位中心;发射机和镜像发射机天线相位中心的高度是hT;接收机和镜像接收机的天线高度分别为hR。M1和M2是地面反射点;θ是反射波擦地角。假设P′为目标上任意一点。r1和r2分别表示发射机和接收机天线相位中心到该点的斜距,和分别表示镜像发射机和镜像接收机天线相位中心到该点的斜距。
步骤2、对接收信号建模,在地面前向散射雷达系统中,接收到的信号认为是图1所示四条路径上传播信号的矢量和,这四条路径分别为:发射机到目标到接收机(以下称为路径1);发射机到地面反射点M1到目标再到接收机(以下称为路径2);发射机到目标到地面反射点M2到接收机(以下称为路径3);发射机到M1到目标到M2到接收机(以下称为路径4)。四条路径的长度分别为:r1+r2, 接收信号表示如公式(1)所示:
其中:
其中,和分别表示反射系数等于1时,四条路径上的接收信号;Γ为复地面反射系数,与地面特性、反射波擦地角、波长和天线的极化方式有关;A是发射信号幅度,λ是连续波信号的波长,k=2π/λ是波数,S是阴影孔径,rc1,rc2是目标中心点分别到发射机和接收机之间的距离。当目标以速度v平行于xoy平面运动,与基线夹角为时,目标的运动方程可以表示为公式(6)和公式(7):
步骤3、对接收信号进行相位补偿和傅里叶变换得到目标复轮廓函数,如公式(8)所示:
其中,Ts为相参积累时间,运动补偿参数和γ可分别表示如公式(9)和(10):
步骤4、多径干涉下目标目标复轮廓函数可表示为:
其中:
是路径1上由发射天线和接收天线高度的一次项和二次项引入的前向散射信号相位。α1,α2,α3,α4分别如公式(19)至(20)所示,可近似认为是四条路径上目标中心点俯仰向观测角αv,βv之和。
步骤5、在多径条件下估计目标轮廓高度差和中线。当地面为理想地面,即反射系数Γ≈-1时:
由目标复轮廓函数的相位信息提取目标的中线高度:
由目标中线高度及目标复轮廓函数幅度可以获得目标轮廓的高度差信息。
当地面为真实地面时,假设地面反射系数(其中:ρ为反射系数的幅度,为反射系数的相位)为简化复杂性,假设发射机和接收机的天线具有相同的高度,目标垂直穿越基线中点。且经路径4到接收机相位中心的前向散射信号,由于经两次地面反射幅度衰减很大,忽略不计,仅考虑经一次反射到达到接收机相位中心的多径信号。此时,可表示如下:
其中:
由于在真实地面的情况下,目标复轮廓函数幅度受多径信号的影响相对较小,目标复轮廓函数幅度与目标高度差的相关性很高。因此,初步粗略的认为 然后再由中线高度的估计值及地面反射率根据公式(25)修正目标轮廓高度差的估计值。通过迭代的方法最终可以获得精确的目标轮廓高度差和中线高度估计值。
有益效果:
本发明提出的一种前向散射雷达地面运动目标信号建模及成像方法的优点在于,考虑地面反射、大衍射角以及目标运动方向较为随机等实际情况,建立了精确的地面运动目标前向散射信号模型,修正了SISAR成像算法,并对目标像进行了更为精确的估计。修正后的成像算法在大衍射角、目标斜穿基线和多径干涉情况下依然可以获得精确的目标像,这些目标像与目标真实轮廓形状非常接近,可有效地对不同形状的目标进行分类识别。
附图说明
图1为地面FSR系统几何结构;
图2为前向散射雷达地面运动目标信号建模及成像方法流程图。
Claims (1)
1.一种前向散射雷达地面运动目标信号建模及成像方法,其特征在于包括如下具体步骤:
步骤1、建立地面前向散射雷达系统的几何结构模型;
xoy平面是地平面,发射机置于坐标系(x,y,z)的原点O,接收机在y轴B(0,L,0)处,发射机和接收机之间的距离为基线长度L,假设目标以速度v平行于xoy平面运动,与基线夹角为αh,βh分别表示从发射机和接收机方向观测的目标方位角,发射机到目标穿越基线位置的距离记为dT,接收机到目标穿越基线位置的距离记为dR,坐标系(x′,y′,z′)与坐标系(x,y,z)各轴相互平行,原点P(xp,yp,zp)为目标的中心,T,R,T′,R′分别表示发射机、接收机、镜像发射机和镜像接收机天线相位中心;发射机和镜像发射机天线相位中心的高度是hT;接收机和镜像接收机的天线高度分别为hR,M1和M2是地面反射点;θ是反射波擦地角,假设P′为目标上任意一点,r1和r2分别表示发射机和接收机天线相位中心到该点的斜距,和分别表示镜像发射机和镜像接收机天线相位中心到该点的斜距;
步骤2、对接收信号建模,在地面前向散射雷达系统中,接收到的信号认为是图1所示四条路径上传播信号的矢量和,这四条路径分别为:发射机到目标到接收机(以下称为路径1);发射机到地面反射点M1到目标再到接收机(以下称为路径2);发射机到目标到地面反射点M2到接收机(以下称为路径3);发射机到M1到目标到M2到接收机(以下称为路径4),四条路径的长度分别为:r1+r2,接收信号表示如公式(1)所示:
其中:
其中,和分别表示反射系数等于1时,四条路径上的接收信号;Γ为复地面反射系数,与地面特性、反射波擦地角、波长和天线的极化方式有关;A是发射信号幅度,λ是连续波信号的波长,k=2π/λ是波数,S是阴影孔径,rc1,rc2是目标中心点分别到发射机和接收机之间的距离,当目标以速度v平行于xoy平面运动,与基线夹角为时,目标的运动方程可以表示为公式(6)和公式(7):
步骤3、对接收信号进行相位补偿和傅里叶变换得到目标复轮廓函数,如公式(8)所示:
其中,Ts为相参积累时间,运动补偿参数和γ可分别表示如公式(9)和(10):
步骤4、多径干涉下目标目标复轮廓函数可表示为:
其中:
步骤5、在多径条件下估计目标轮廓高度差和中线,当地面为理想地面,即反射系数Γ≈-1时:
由目标复轮廓函数的相位信息提取目标的中线高度:
由目标中线高度及目标复轮廓函数幅度可以获得目标轮廓的高度差信息,
当地面为真实地面时,假设地面反射系数(其中:ρ为反射系数的幅度,为反射系数的相位)为简化复杂性,假设发射机和接收机的天线具有相同的高度,目标垂直穿越基线中点,且经路径4到接收机相位中心的前向散射信号,由于经两次地面反射幅度衰减很大,忽略不计,仅考虑经一次反射到达到接收机相位中心的多径信号,此时,可表示如下:
其中:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201010587979 CN102129067B (zh) | 2011-03-15 | 2011-03-15 | 一种前向散射雷达地面运动目标信号建模及成像方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201010587979 CN102129067B (zh) | 2011-03-15 | 2011-03-15 | 一种前向散射雷达地面运动目标信号建模及成像方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102129067A true CN102129067A (zh) | 2011-07-20 |
CN102129067B CN102129067B (zh) | 2013-07-03 |
Family
ID=44267218
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 201010587979 Expired - Fee Related CN102129067B (zh) | 2011-03-15 | 2011-03-15 | 一种前向散射雷达地面运动目标信号建模及成像方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102129067B (zh) |
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102981153A (zh) * | 2012-11-13 | 2013-03-20 | 北京理工大学 | 一种阴影逆合成孔径雷达无线电全息信号提取方法 |
CN103197293A (zh) * | 2013-04-15 | 2013-07-10 | 中国科学院电子学研究所 | 一种机载干涉sar多路径误差的提取与补偿方法 |
CN103675814A (zh) * | 2013-09-25 | 2014-03-26 | 中国科学院电子学研究所 | 基于圆周sar确定建筑物地平面高度的方法 |
CN104049253A (zh) * | 2014-06-16 | 2014-09-17 | 北京理工大学 | 一种前向散射雷达阴影逆合成孔径成像方法 |
CN107110954A (zh) * | 2014-12-23 | 2017-08-29 | 巴鲁夫公司 | 用于测量与目标相距的距离的接近传感器和方法 |
CN108562916A (zh) * | 2018-01-26 | 2018-09-21 | 中国科学院上海天文台 | 适用于复杂寒区冻融地表的前向gps多路径模型的建立方法 |
CN109358325A (zh) * | 2018-12-06 | 2019-02-19 | 西安电子科技大学 | 起伏地形背景下雷达高度表的地形反演方法 |
US10209344B2 (en) | 2015-03-12 | 2019-02-19 | Src, Inc. | Methods and systems for mitigating multipath propagation |
CN109490872A (zh) * | 2017-09-12 | 2019-03-19 | 德尔福技术有限责任公司 | 使用幅度和多个接收元件的单散射体测试 |
CN110174668A (zh) * | 2019-04-26 | 2019-08-27 | 中国矿业大学 | 一种矿井无源动目标轮廓识别方法 |
CN111208499A (zh) * | 2020-01-16 | 2020-05-29 | 深圳以正科技有限公司 | 基于声前向散射雷达的运动目标探测方法和系统 |
CN111896924A (zh) * | 2020-08-17 | 2020-11-06 | 扬州船用电子仪器研究所(中国船舶重工集团公司第七二三研究所) | 一种雷达海面低空探测性能仿真方法及仿真模型 |
CN113393509A (zh) * | 2021-08-16 | 2021-09-14 | 浙江光珀智能科技有限公司 | 一种激光雷达多径补偿方法 |
CN114039679A (zh) * | 2022-01-10 | 2022-02-11 | 中国人民解放军海军工程大学 | 一种低频正交天线信号检测方法及系统 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6756935B1 (en) * | 2003-01-31 | 2004-06-29 | The Boeing Company | Full polarization ground moving target indicator radar automatic target detection algorithm |
CN101881825A (zh) * | 2009-05-06 | 2010-11-10 | 中国人民解放军海军航空工程学院 | 利用延时的frft模之差的动目标检测器 |
-
2011
- 2011-03-15 CN CN 201010587979 patent/CN102129067B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6756935B1 (en) * | 2003-01-31 | 2004-06-29 | The Boeing Company | Full polarization ground moving target indicator radar automatic target detection algorithm |
CN101881825A (zh) * | 2009-05-06 | 2010-11-10 | 中国人民解放军海军航空工程学院 | 利用延时的frft模之差的动目标检测器 |
Non-Patent Citations (1)
Title |
---|
李晓良等: "《多极化前向散射RCS分析及其对目标分类识别的影响》", 《电子与信息学报》, vol. 23, no. 9, 30 September 2010 (2010-09-30) * |
Cited By (23)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102981153A (zh) * | 2012-11-13 | 2013-03-20 | 北京理工大学 | 一种阴影逆合成孔径雷达无线电全息信号提取方法 |
CN102981153B (zh) * | 2012-11-13 | 2014-07-02 | 北京理工大学 | 一种阴影逆合成孔径雷达无线电全息信号提取方法 |
CN103197293A (zh) * | 2013-04-15 | 2013-07-10 | 中国科学院电子学研究所 | 一种机载干涉sar多路径误差的提取与补偿方法 |
CN103197293B (zh) * | 2013-04-15 | 2014-02-26 | 中国科学院电子学研究所 | 一种机载干涉sar多路径误差的提取与补偿方法 |
CN103675814A (zh) * | 2013-09-25 | 2014-03-26 | 中国科学院电子学研究所 | 基于圆周sar确定建筑物地平面高度的方法 |
CN103675814B (zh) * | 2013-09-25 | 2016-03-23 | 中国科学院电子学研究所 | 基于圆周sar确定建筑物地平面高度的方法 |
CN104049253A (zh) * | 2014-06-16 | 2014-09-17 | 北京理工大学 | 一种前向散射雷达阴影逆合成孔径成像方法 |
CN107110954A (zh) * | 2014-12-23 | 2017-08-29 | 巴鲁夫公司 | 用于测量与目标相距的距离的接近传感器和方法 |
CN107110954B (zh) * | 2014-12-23 | 2020-12-22 | 巴鲁夫公司 | 用于测量与目标相距的距离的接近传感器和方法 |
US10209344B2 (en) | 2015-03-12 | 2019-02-19 | Src, Inc. | Methods and systems for mitigating multipath propagation |
CN109490872A (zh) * | 2017-09-12 | 2019-03-19 | 德尔福技术有限责任公司 | 使用幅度和多个接收元件的单散射体测试 |
CN108562916A (zh) * | 2018-01-26 | 2018-09-21 | 中国科学院上海天文台 | 适用于复杂寒区冻融地表的前向gps多路径模型的建立方法 |
WO2019144597A1 (zh) * | 2018-01-26 | 2019-08-01 | 中国科学院上海天文台 | 适用于复杂寒区冻融地表的前向gps多路径模型的建立方法 |
CN109358325A (zh) * | 2018-12-06 | 2019-02-19 | 西安电子科技大学 | 起伏地形背景下雷达高度表的地形反演方法 |
CN109358325B (zh) * | 2018-12-06 | 2022-10-04 | 西安电子科技大学 | 起伏地形背景下雷达高度表的地形反演方法 |
CN110174668A (zh) * | 2019-04-26 | 2019-08-27 | 中国矿业大学 | 一种矿井无源动目标轮廓识别方法 |
CN110174668B (zh) * | 2019-04-26 | 2022-12-20 | 中国矿业大学 | 一种矿井无源动目标轮廓识别方法 |
CN111208499A (zh) * | 2020-01-16 | 2020-05-29 | 深圳以正科技有限公司 | 基于声前向散射雷达的运动目标探测方法和系统 |
CN111896924A (zh) * | 2020-08-17 | 2020-11-06 | 扬州船用电子仪器研究所(中国船舶重工集团公司第七二三研究所) | 一种雷达海面低空探测性能仿真方法及仿真模型 |
CN111896924B (zh) * | 2020-08-17 | 2023-09-26 | 扬州船用电子仪器研究所(中国船舶重工集团公司第七二三研究所) | 一种雷达海面低空探测性能仿真方法及仿真模型 |
CN113393509A (zh) * | 2021-08-16 | 2021-09-14 | 浙江光珀智能科技有限公司 | 一种激光雷达多径补偿方法 |
CN114039679A (zh) * | 2022-01-10 | 2022-02-11 | 中国人民解放军海军工程大学 | 一种低频正交天线信号检测方法及系统 |
CN114039679B (zh) * | 2022-01-10 | 2022-04-01 | 中国人民解放军海军工程大学 | 一种低频正交天线信号检测方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN102129067B (zh) | 2013-07-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102129067B (zh) | 一种前向散射雷达地面运动目标信号建模及成像方法 | |
CN102967858B (zh) | 雷达前视超分辨成像方法 | |
CA2613800C (en) | Subsurface imaging radar | |
Zhu et al. | Application of feature extraction methods for landmine detection using the Wichmann/Niitek ground-penetrating radar | |
CN102778681B (zh) | 固定发射站双基地前视合成孔径雷达成像方法 | |
Paglieroni et al. | Imaging modes for ground penetrating radar and their relation to detection performance | |
Shishanov et al. | Height-finding for automotive THz radars | |
Vriesman et al. | An experimental analysis of rain interference on detection and ranging sensors | |
Watson et al. | Non-line-of-sight radar | |
Prager et al. | Application of ultra-wideband synthesis in software defined radar for UAV-based landmine detection | |
CN105652257A (zh) | 一种利用双站单阵元高频地波雷达探测船只目标的方法 | |
CN105044721A (zh) | 机载正前视扫描雷达角超分辨方法 | |
Iizuka et al. | MilliSign: mmWave-Based Passive Signs for Guiding UAVs in Poor Visibility Conditions | |
CN113376628B (zh) | 一种基于空间等离子体云的集群无人机超视距探测方法 | |
Pei et al. | mmCTD: Concealed Threat Detection for Cruise Ships Via Mmwave Radar | |
Hu et al. | An accurate SISAR imaging method of ground moving target in forward scatter radar | |
Wang et al. | Sparse time-frequency representation based feature extraction method for landmine discrimination | |
Latthe et al. | A review on ISAR imaging techniques for low RCS targets | |
Liu et al. | Spaceborne-airborne bistatic radar clutter modeling and analysis | |
Wang et al. | An imaging algorithm for obstacle detection of tunnel horizontal transporter based on millimeter wave radar | |
Gessel et al. | Three-Dimensional Emitter Selection Optimization for Passive GMTI | |
Sheng et al. | Method to reconnoiter pulse amplitude train for phased array radar based on Hausdorff distance | |
Fan et al. | Shallowly Buried Trash Detection in Sandy Land Based on IR-UWB Radar | |
Santos et al. | FMCW Radar Point Cloud Multiperson Tracking Using a Kalman Filter–Based Approach | |
Kitamura et al. | Range Migration Algorithm Using Doppler-Division Multiplexed Multiple-Input Multiple-Output Imaging |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20130703 Termination date: 20160315 |
|
CF01 | Termination of patent right due to non-payment of annual fee |