CN108614250B - 机载战场监视雷达的广域dbs图像拼接暗条纹校正方法 - Google Patents

机载战场监视雷达的广域dbs图像拼接暗条纹校正方法 Download PDF

Info

Publication number
CN108614250B
CN108614250B CN201810481692.4A CN201810481692A CN108614250B CN 108614250 B CN108614250 B CN 108614250B CN 201810481692 A CN201810481692 A CN 201810481692A CN 108614250 B CN108614250 B CN 108614250B
Authority
CN
China
Prior art keywords
doppler
airborne
range
data
battlefield
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.)
Active
Application number
CN201810481692.4A
Other languages
English (en)
Other versions
CN108614250A (zh
Inventor
王彤
刘嘉昕
刘程
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Xidian University
Original Assignee
Xidian University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Xidian University filed Critical Xidian University
Priority to CN201810481692.4A priority Critical patent/CN108614250B/zh
Publication of CN108614250A publication Critical patent/CN108614250A/zh
Application granted granted Critical
Publication of CN108614250B publication Critical patent/CN108614250B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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
    • G01S7/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/40Means for monitoring or calibrating
    • G01S7/4004Means for monitoring or calibrating of parts of a radar system
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Systems 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/88Radar or analogous systems specially adapted for specific applications
    • G01S13/89Radar or analogous systems specially adapted for specific applications for mapping or imaging

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Physics & Mathematics (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • General Physics & Mathematics (AREA)
  • Electromagnetism (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明公开了一种机载战场监视雷达的广域DBS图像拼接暗条纹校正方法,思路为:确定机载战场监视雷达,在机载战场监视雷达成像周期中共设置I个波位,机载战场监视雷达在I个波位发射I×T个脉冲信号后接收回波信号,并根据接收的回波信号,得到机载战场监视雷达I个波位的距离‑多普勒数据后据进行融合,得到机载战场监视雷达成像周期的距离‑多普勒融合数据;计算机载战场监视雷达成像周期的距离‑多普勒融合数据均值;然后得到增益校正处理后机载战场监视雷达成像周期I个波位的距离‑多普勒融合数据后得到暗条纹校正后的广域DBS图像,所述暗条纹校正后的广域DBS图像为机载战场监视雷达的广域DBS图像拼接暗条纹校正结果。

Description

机载战场监视雷达的广域DBS图像拼接暗条纹校正方法
技术领域
本发明属于雷达技术领域,特别涉及一种机载战场监视雷达的广域DBS图像拼接暗条纹校正方法,即一种机载战场监视雷达的广域多普勒波束锐化(DBS)图像拼接暗条纹校正方法,适用于消除机载战场监视雷达的广域DBS图像拼接暗条纹。
背景技术
机载战场监视雷达工作在广域地面动目标指示(GMTI)模式下时,机载战场监视雷达天线采用周期性的方位向扫描或方位-俯仰二维扫描的方式工作,从而可对远距离宽广范围的地面场景成像并实时监视地面运动目标;多普勒波束锐化(DBS)成像技术是机载战场监视雷达工作在广域GMTI模式下时经常配合使用的一项关键成像技术,DBS成像技术的概念最早于1951年被美国Goodyear公司提出,他们将这种利用多普勒滤波来提高方位向分辨力的成像方法命名为“多普勒锐化”;虽然它的成像分辨率不及聚焦成像的SAR成像技术,但是DBS成像技术处理方式较为简单,对处理器的运算负荷要求较低,可以在短时间内完成对广域地面场景的成像,为后续对地面运动目标的检测、跟踪与打击提供实时的地形参考,因此其在战场监视、战场侦察、地形匹配制导等军事领域均有十分重要的应用价值。
单波位DBS图像仅可显示主波束范围所对应的窄扇形地面区域,虽可观察到地面的部分场景,但是由于范围过窄,因此无法正确识别该场景所对应的地形信息;为实现广域地面区域的大场景测绘,需要将多幅单波位DBS图像准确拼接起来,图像拼接工作就是在给定的地面坐标系中将每幅单波位DBS图像放置于正确的位置上;然而,机载战场监视雷达在实际工作中,由于风速、气流等因素的影响,使得载机无法完全按照规定航路维持理想的匀速运动,这将导致雷达的距离-多普勒数据无法与地面上的实际坐标点配准,从而使得拼接后的图像出现错位、重影等引起图像质量下降的问题。2012年,西安电子科技大学胡瑞贤提出了一种基于惯导数据运动补偿的DBS图像拼接算法,该算法利用惯导数据对载机的运动状态实时补偿,可有效解决由载机非理想运动引起的雷达距离-多普勒数据与实际地面坐标失配的问题。
在实际工程中运用上述基于惯导数据运动补偿的DBS图像拼接算法时,一些拼接后的广域图像上会出现宽窄不等且分布不规则的暗条纹,这些暗条纹不仅影响了图像质量,严重情况下甚至也干扰了地形信息的正确识别。
发明内容
针对上述现有技术存在的不足,本发明的目的在于提出一种机载战场监视雷达的广域DBS图像拼接暗条纹校正方法,该种一种机载战场监视雷达的广域DBS图像拼接暗条纹校正方法通过分析相邻波位回波数据的距离-多普勒关系,依此构造出反映扫描周期回波功率与多普勒变化关系的曲线,根据该曲线对每个波位的距离-多普勒数据进行适度增益补偿,从而实现对拼接暗条纹区域回波功率的校正。
为达到上述技术目的,本发明采用如下技术方案予以实现。
一种机载战场监视雷达的广域DBS图像拼接暗条纹校正方法,包括以下步骤:
步骤1,确定机载战场监视雷达,在机载战场监视雷达成像周期中共设置I个波位,机载战场监视雷达在每一个波位都发射数目为T的脉冲信号;其中,I和T分别为大于1的正整数;
步骤2,机载战场监视雷达在I个波位发射I×T个脉冲信号后接收回波信号,并根据接收的回波信号,得到机载战场监视雷达I个波位的距离-多普勒数据;
步骤3,对机载战场监视雷达I个波位的距离-多普勒数据进行融合,得到机载战场监视雷达成像周期的距离-多普勒融合数据Pall
步骤4,计算机载战场监视雷达成像周期的距离-多普勒融合数据Pall中每一列距离-多普勒融合数据的均值,进而得到机载战场监视雷达成像周期的距离-多普勒融合数据均值;
步骤5,根据机载战场监视雷达成像周期的距离-多普勒融合数据均值,得到增益校正处理后机载战场监视雷达成像周期I个波位的距离-多普勒融合数据P′all
步骤6,根据增益校正处理后机载战场监视雷达成像周期中I个波位的距离-多普勒融合数据P′all,得到暗条纹校正后的广域DBS图像,所述暗条纹校正后的广域DBS图像为机载战场监视雷达的广域DBS图像拼接暗条纹校正结果。
本发明的有益效果:针对拼接后的广域图像上出现分布不规则且宽窄不等的暗条纹情况,本发明创新性地提出了一种广域DBS图像拼接暗条纹校正方法,该种方法能够有效地消除拼接后的广域图像上出现的暗条纹,使成像质量获得了明显提升。
附图说明
下面结合附图和具体实施方式对本发明作进一步详细说明。
图1为本发明的一种机载战场监视雷达的广域DBS图像拼接暗条纹校正方法流程图;
图2a为未使用本发明方法直接对实测数据成像拼接后得到的广域DBS图;
图2b为采用本发明方法进行暗条纹消除后获得的广域DBS图。
具体实施方式
参照图1,为本发明的一种机载战场监视雷达的广域DBS图像拼接暗条纹校正方法流程图;其中所述机载战场监视雷达的广域DBS图像拼接暗条纹校正方法,包括以下步骤:
步骤1,确定机载战场监视雷达,机载战场监视雷达发射波长为λ的电磁波;该机载战场监视雷达置于其载机上,记为战场监视雷达载机;设置战场监视雷达载机始终在平行于地面且高度为H的平面内以速度v飞行,以战场监视雷达载机在地面上的投影为原点o建立xoyz坐标系,其中x轴与战场监视雷达载机飞行速度同方向,y轴为战场监视雷达载机飞行速度方向的垂直方向,z轴为经过战场监视雷达载机机身中心且垂直于地面的方向,且xoy平面为地平面;选取地平面上任意一点,记为点目标P,点目标P的坐标为(xp,yp),xp表示点目标P在地平面中的x轴坐标,yp表示点目标P在地平面中的y轴坐标;点目标P的方位角为θ,点目标的俯仰角为
Figure BDA0001666027990000031
战场监视雷达载机飞行速度方向与点目标P的空间锥角为ψ。
在机载战场监视雷达成像周期中共设置I个波位,波位指机载战场监视雷达波束在空间的驻留位置;战场监视雷达载机在飞行过程中,机载战场监视雷达在一个波位持续发射数目为T的脉冲信号,然后切换至下一个波位。
机载战场监视雷达在I个波位发射I×T个脉冲信号后接收回波信号,将机载战场监视雷达在I个波位发射I×T个脉冲信号后接收的回波信号,记为接收机载战场监视雷达的接收信号,所述接收机载战场监视雷达的接收信号是I个R×T维矩阵,其中第i个R×T维矩阵为Di,其表达式为:
Figure BDA0001666027990000032
其中,i=1,2,…,I,R表示机载战场监视雷达的接收信号包括的距离单元总个数,T表示机载战场监视雷达在每个波位发射脉冲信号的总个数,Di(ri,ti)表示第i个波位、第ri个距离单元、第ti个脉冲信号的接收信号,ri=1,2,…,R,ti=1,2,…,T,I和T分别为大于1的正整数。
步骤2,对第i个R×T维矩阵Di依次做距离压缩、距离徙动校正、方位向快速傅里叶变换FFT处理,将第i个R×T维矩阵Di变换为距离-多普勒数据,进而得到第i个波位的距离-多普勒数据Pi,第i个波位的距离-多普勒数据Pi为R×F维矩阵,其具体表达式如下:
Figure BDA0001666027990000041
其中,Pi(ri,fi)表示第i个波位、第ri个距离单元、第fi个多普勒单元的距离-多普勒数据,ri=1,2,…,R,R表示机载战场监视雷达的接收信号包括的距离单元总个数,fi=1,2,…,F,F表示第i个波位的距离-多普勒数据Pi包括的多普勒单元总个数,且F与机载战场监视雷达在每个波位发射脉冲信号的总个数T取值相等。
步骤3,令i的值分别取1至I,重复执行步骤2,进而分别得到第1个波位的距离-多普勒数据P1至第I个波位的距离-多普勒数据PI,记为机载战场监视雷达I个波位的距离-多普勒数据。
步骤4,对机载战场监视雷达I个波位的距离-多普勒数据进行融合,其具体子步骤为:
4.1取出第i'个波位的距离-多普勒数据Pi',所述第i'个波位的距离-多普勒数据Pi'包括R×F个距离-多普勒单元,其中第i'个距离-多普勒单元为(ri',fi'),ri'=1,2,…,R,R表示机载战场监视雷达的接收信号包括的距离单元总个数,fi'=1,2,…,F,F表示第i'个波位的距离-多普勒数据Pi'包括的多普勒单元总个数,且F与机载战场监视雷达在每个波位发射脉冲信号的总个数T取值相等;i'=1,2,…,I-1,I表示机载战场监视雷达的成像周期中设置的波位总个数,i'的初始值为1。
4.2计算第i'个距离-多普勒单元(ri',fi')在第i'+1个波位的距离-多普勒数据Pi'+1中所处的距离-多普勒单元(ri'+1,fi'+1),其具体计算方法为:
Figure BDA0001666027990000051
Figure BDA0001666027990000052
其中,v表示战场监视雷达载机的飞行速度,T表示机载战场监视雷达在每个波位发射脉冲信号的总个数,λ表示机载战场监视雷达发射的电磁波波长。
然后,将所述距离-多普勒单元(ri'+1,fi'+1)处的距离-多普勒数据,记为第i'+1个波位距离-多普勒单元(ri'+1,fi'+1)的距离-多普勒数据Pi'+1(ri'+1,fi'+1)。
4.3比较第i'个波位、第ri'个距离单元、第fi'个多普勒单元的距离-多普勒数据Pi'(ri',fi')与第i'+1个波位距离-多普勒单元(ri'+1,fi'+1)的距离-多普勒数据Pi'+1(ri'+1,fi'+1)的大小,选取其中最大值作为第i'个距离-多普勒单元和第i'+1个距离多普勒单元的融合数据。
4.4令ri'分别取1至R,令fi'分别取1至F,重复执行子步骤4.2-4.3遍历第i'个波位的距离-多普勒数据Pi'中的R×F个距离-多普勒单元,完成第i'个波位的距离-多普勒数据Pi'与第i'+1个波位的距离-多普勒数据Pi'+1的融合,得到第i'个波位的距离-多普勒数据Pi'与第i'+1个波位的距离-多普勒数据Pi'+1的融合数据,然后将ri'的初始化为1,将fi'的值初始化1。
4.5令i'的值分别取1至I-1,重复执行子步骤4.2-4.4,直到得到第1个波位的距离-多普勒数据P1与第2个波位的距离-多普勒数据P2的融合数据至第I-1个波位的距离-多普勒数据PI-1与第I个波位的距离-多普勒数据PI的融合数据,记为机载战场监视雷达成像周期的距离-多普勒融合数据Pall,机载战场监视雷达成像周期的距离-多普勒融合数据Pall是Fall列向量,其具体表达式为:
Figure BDA0001666027990000061
其中,Pall_f表示机载战场监视雷达成像周期的距离-多普勒融合数据Pall的第f列距离-多普勒融合数据,Pall_f=[Pall(1,f),…,Pall(r,f),…,Pall(R,f)]T',上标T'表示转置操作,Pall(r,f)表示机载战场监视雷达成像周期的距离-多普勒融合数据Pall中第r个距离单元、第f个多普勒单元处的距离-多普勒融合数据,r=1,2,…,R,R表示机载战场监视雷达的接收信号包括的距离单元总个数,f=1,2,…,Fall,Fall表示机载战场监视雷达成像周期的距离-多普勒融合数据Pall包括的多普勒单元总个数。
步骤5,计算机载战场监视雷达成像周期的距离-多普勒融合数据Pall中每一列距离-多普勒融合数据的均值,进而得到机载战场监视雷达成像周期的距离-多普勒融合数据均值。
步骤5的具体子步骤如下:
5.1计算机载战场监视雷达成像周期的距离-多普勒融合数据Pall的第f列距离-多普勒融合数据Pall_f的均值pavg_f,其计算表达式为:
Figure BDA0001666027990000062
其中,f的初始值为1。
5.2令f的值分别取1至Fall,重复执行子步骤5.1,进而分别得到机载战场监视雷达成像周期的距离-多普勒融合数据Pall的第1列距离-多普勒融合数据Pall_1的均值pavg_1至机载战场监视雷达成像周期的距离-多普勒融合数据Pall的第Fall列距离-多普勒融合数据
Figure BDA0001666027990000063
的均值
Figure BDA0001666027990000064
记为机载战场监视雷达成像周期的距离-多普勒融合数据均值,然后将f的值初始化为1。
步骤6,根据机载战场监视雷达成像周期的距离-多普勒融合数据均值,得到增益校正处理后机载战场监视雷达成像周期I个波位的距离-多普勒融合数据P′all
步骤6的子步骤为:
6.1对机载战场监视雷达成像周期的距离-多普勒融合数据Pall的第f列距离-多普勒融合数据Pall_f的均值pavg_f进行增益校正处理,得到增益校正处理后的第f列距离-多普勒融合数据P′all_f,其计算公式为:
P′all_f=Pall_f⊙(Pall_f-pavg_f)
其中,⊙表示矩阵点乘。
6.2令f的值分别取1至Fall,重复执行子步骤6.1,进而分别得到增益校正处理后的第1列距离-多普勒融合数据P′all_1至增益校正处理后的第Fall列距离-多普勒融合数据
Figure BDA0001666027990000071
记为增益校正处理后机载战场监视雷达成像周期I个波位的距离-多普勒融合数据P′all,其表达式为:
Figure BDA0001666027990000072
步骤7,利用DBS成像技术和基于惯导数据运动补偿的图像拼接技术对增益校正处理后机载战场监视雷达成像周期中I个波位的距离-多普勒融合数据P′all依次进行成像与拼接处理,以达到暗条纹校正目的,进而得到暗条纹校正后的广域DBS图像,所述暗条纹校正后的广域DBS图像为机载战场监视雷达的广域DBS图像拼接暗条纹校正结果。
本发明的效果可由以下仿真实验作进一步说明:
(一)仿真条件:
对一组机载战场监视雷达在广域GMTI模式下获取的实测回波数据进行仿真处理,该组机载战场监视雷达的工作频率都为X波段,信号带宽约为25MHz,采用双通道接收回波数据,发射脉冲数目为1024,距离单元数目为9216,机载战场监视雷达成像周期中共有61个波位,每个机载战场监视雷达天线的扫描角范围为[-60°,60°],俯仰角的变化范围为[2°,4°],每个机载战场监视雷达载机的飞行高度约为8400m。
(二)仿真结果分析:
图2a为未使用本发明方法直接对实测数据成像拼接后得到的广域DBS图像,图2b为采用本发明方法进行暗条纹消除后获得的广域DBS图。
对比图2a与图2b可以看出,图2a中出现的不规则分布暗条纹在图2b的成像结果中被有效地改善,使用本发明方法校正后的广域DBS图像质量相比于原图像获得了明显提升。
显然,本领域的技术人员可以对本发明进行各种改动和变型而不脱离本发明的精神和范围;这样,倘若本发明的这些修改和变型属于本发明权利要求及其等同技术的范围之内,则本发明也意图包含这些改动和变型在内。

Claims (4)

1.一种机载战场监视雷达的广域DBS图像拼接暗条纹校正方法,其特征在于,包括以下步骤:
步骤1,确定机载战场监视雷达,在机载战场监视雷达成像周期中共设置I个波位,机载战场监视雷达在每一个波位都发射数目为T的脉冲信号;其中,I和T分别为大于1的正整数;
步骤2,机载战场监视雷达在I个波位发射I×T个脉冲信号后接收回波信号,并根据接收的回波信号,得到机载战场监视雷达I个波位的距离-多普勒数据;
在所述步骤2中,所述机载战场监视雷达在I个波位发射I×T个脉冲信号后接收回波信号,还包括:
将机载战场监视雷达在I个波位发射I×T个脉冲信号后接收的回波信号,记为接收机载战场监视雷达的接收信号,所述接收机载战场监视雷达的接收信号是I个R×T维矩阵,其中第i个R×T维矩阵为Di,其表达式为:
Figure FDA0003299735060000011
其中,i=1,2,...,I,R表示机载战场监视雷达的接收信号包括的距离单元总个数,T表示机载战场监视雷达在每个波位发射脉冲信号的总个数,Di(ri,ti)表示第i个波位、第ri个距离单元、第ti个脉冲信号的接收信号,ri=1,2,...,R,ti=1,2,...,T,I和T分别为大于1的正整数;
在所述步骤2中,所述机载战场监视雷达I个波位的距离-多普勒数据,其得到过程为:
2.1对第i个R×T维矩阵Di依次做距离压缩、距离徙动校正、方位向快速傅里叶变换FFT处理,将第i个R×T维矩阵Di变换为距离-多普勒数据,进而得到第i个波位的距离-多普勒数据Pi,第i个波位的距离-多普勒数据Pi为R×F维矩阵,其具体表达式如下:
Figure FDA0003299735060000012
其中,Pi(ri,fi)表示第i个波位、第ri个距离单元、第fi个多普勒单元的距离-多普勒数据,ri=1,2,...,R,R表示机载战场监视雷达的接收信号包括的距离单元总个数,fi=1,2,...,F,F表示第i个波位的距离-多普勒数据Pi包括的多普勒单元总个数,且F与机载战场监视雷达在每个波位发射脉冲信号的总个数T取值相等;
2.2令i的值分别取1至I,重复执行2.1,进而分别得到第1个波位的距离-多普勒数据P1至第I个波位的距离-多普勒数据PI,记为机载战场监视雷达I个波位的距离-多普勒数据;
步骤3,对机载战场监视雷达I个波位的距离-多普勒数据进行融合,得到机载战场监视雷达成像周期的距离-多普勒融合数据Pall
在所述步骤3中,所述机载战场监视雷达成像周期的距离-多普勒融合数据Pall,其得到过程为:
3.1取出第i′个波位的距离-多普勒数据Pi′,所述第i′个波位的距离-多普勒数据Pi′包括R×F个距离-多普勒单元,其中第i′个距离-多普勒单元为(ri′,fi′),ri′=1,2,...,R,R表示机载战场监视雷达的接收信号包括的距离单元总个数,fi′=1,2,...,F,F表示第i′个波位的距离-多普勒数据Pi′包括的多普勒单元总个数,且F与机载战场监视雷达在每个波位发射脉冲信号的总个数T取值相等;i′=1,2,…,I-1,I表示机载战场监视雷达的成像周期中设置的波位总个数,i′的初始值为1;
3.2计算第i′个距离-多普勒单元(ri′,fi′)在第i′+1个波位的距离-多普勒数据Pi′+1中所处的距离-多普勒单元(ri′+1,fi′+1),其具体计算方法为:
Figure FDA0003299735060000021
Figure FDA0003299735060000022
其中,v表示战场监视雷达载机的飞行速度,T表示机载战场监视雷达在每个波位发射脉冲信号的总个数,λ表示机载战场监视雷达发射的电磁波波长;
然后,将所述距离-多普勒单元(ri′+1,fi′+1)处的距离-多普勒数据,记为第i′+1个波位距离-多普勒单元(ri′+1,fi′+1)的距离-多普勒数据Pi′+1(ri′+1,fi′+1);
3.3比较第i′个波位、第ri′个距离单元、第fi′个多普勒单元的距离-多普勒数据Pi′(ri′,fi′)与第i′+1个波位距离-多普勒单元(ri′+1,fi′+1)的距离-多普勒数据Pi′+1(ri′+1,fi′+1)的大小,选取其中最大值作为第i′个距离-多普勒单元和第i′+1个距离多普勒单元的融合数据;
3.4令ri′分别取1至R,令fi′分别取1至F,重复执行子步骤3.2-3.3遍历第i′个波位的距离-多普勒数据Pi′中的R×F个距离-多普勒单元,完成第i′个波位的距离-多普勒数据Pi′与第i′+1个波位的距离-多普勒数据Pi′+1的融合,得到第i′个波位的距离-多普勒数据Pi′与第i′+1个波位的距离-多普勒数据Pi′+1的融合数据,然后将ri′的初始化为1,将fi′的值初始化1;
3.5令i′的值分别取1至I-1,重复执行子步骤3.2-3.4,直到得到第1个波位的距离-多普勒数据P1与第2个波位的距离-多普勒数据P2的融合数据至第I-1个波位的距离-多普勒数据PI-1与第I个波位的距离-多普勒数据PI的融合数据,记为机载战场监视雷达成像周期的距离-多普勒融合数据Pall,机载战场监视雷达成像周期的距离-多普勒融合数据Pall是Fall列向量,其具体表达式为:
Figure FDA0003299735060000031
其中,Pall_f表示机载战场监视雷达成像周期的距离-多普勒融合数据Pall的第f列距离-多普勒融合数据,Pall_f=[Pall(1,f),...,Pall(r,f),...,Pall(R,f)]T′,上标T′表示转置操作,Pall(r,f)表示机载战场监视雷达成像周期的距离-多普勒融合数据Pall中第r个距离单元、第f个多普勒单元处的距离-多普勒融合数据,r=1,2,...,R,R表示机载战场监视雷达的接收信号包括的距离单元总个数,f=1,2,...,Fall,Fall表示机载战场监视雷达成像周期的距离-多普勒融合数据Pall包括的多普勒单元总个数;
步骤4,计算机载战场监视雷达成像周期的距离-多普勒融合数据Pall中每一列距离-多普勒融合数据的均值,进而得到机载战场监视雷达成像周期的距离-多普勒融合数据均值;
步骤5,根据机载战场监视雷达成像周期的距离-多普勒融合数据均值,得到增益校正处理后机载战场监视雷达成像周期I个波位的距离-多普勒融合数据P′all
步骤6,根据增益校正处理后机载战场监视雷达成像周期中I个波位的距离-多普勒融合数据P′all,得到暗条纹校正后的广域DBS图像,所述暗条纹校正后的广域DBS图像为机载战场监视雷达的广域DBS图像拼接暗条纹校正结果。
2.如权利要求1所述的一种机载战场监视雷达的广域DBS图像拼接暗条纹校正方法,其特征在于,在步骤4中,所述机载战场监视雷达成像周期的距离-多普勒融合数据均值,其得到过程为:
4.1计算机载战场监视雷达成像周期的距离-多普勒融合数据Pall的第f列距离-多普勒融合数据Pall_f的均值pavg_f,其计算表达式为:
Figure FDA0003299735060000041
其中,f的初始值为1;
4.2令f的值分别取1至Fall,重复执行子步骤4.1,进而分别得到机载战场监视雷达成像周期的距离-多普勒融合数据Pall的第1列距离-多普勒融合数据Pall_1的均值pavg_1至机载战场监视雷达成像周期的距离-多普勒融合数据Pall的第Fall列距离-多普勒融合数据
Figure FDA0003299735060000042
的均值pavg_Fall,记为机载战场监视雷达成像周期的距离-多普勒融合数据均值,然后将f的值初始化为1。
3.如权利要求2所述的一种机载战场监视雷达的广域DBS图像拼接暗条纹校正方法,其特征在于,在步骤5中,所述增益校正处理后机载战场监视雷达成像周期I个波位的距离-多普勒融合数据P′all,其得到过程为:
5.1对机载战场监视雷达成像周期的距离-多普勒融合数据Pall的第f列距离-多普勒融合数据Pall_f的均值pavg_f进行增益校正处理,得到增益校正处理后的第f列距离-多普勒融合数据P′all_f,其计算公式为:
P′all_f=Pall_f⊙(Pall_f-pavg_f)
其中,⊙表示矩阵点乘;
5.2令f的值分别取1至Fall,重复执行子步骤5.1,进而分别得到增益校正处理后的第1列距离-多普勒融合数据P′all_1至增益校正处理后的第Fall列距离-多普勒融合数据P′all_Fall,记为增益校正处理后机载战场监视雷达成像周期I个波位的距离-多普勒融合数据P′all,其表达式为:
Figure FDA0003299735060000051
4.如权利要求1所述的一种机载战场监视雷达的广域DBS图像拼接暗条纹校正方法,其特征在于,在步骤6中,所述暗条纹校正后的广域DBS图像,具体是利用DBS成像技术和基于惯导数据运动补偿的图像拼接技术对增益校正处理后机载战场监视雷达成像周期中I个波位的距离-多普勒融合数据P′all依次进行成像与拼接处理,以达到暗条纹校正目的,进而得到的暗条纹校正后的广域DBS图像。
CN201810481692.4A 2018-05-18 2018-05-18 机载战场监视雷达的广域dbs图像拼接暗条纹校正方法 Active CN108614250B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810481692.4A CN108614250B (zh) 2018-05-18 2018-05-18 机载战场监视雷达的广域dbs图像拼接暗条纹校正方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810481692.4A CN108614250B (zh) 2018-05-18 2018-05-18 机载战场监视雷达的广域dbs图像拼接暗条纹校正方法

Publications (2)

Publication Number Publication Date
CN108614250A CN108614250A (zh) 2018-10-02
CN108614250B true CN108614250B (zh) 2022-02-22

Family

ID=63663544

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810481692.4A Active CN108614250B (zh) 2018-05-18 2018-05-18 机载战场监视雷达的广域dbs图像拼接暗条纹校正方法

Country Status (1)

Country Link
CN (1) CN108614250B (zh)

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5237329A (en) * 1991-07-08 1993-08-17 Deutsche Forschungsanstalt Fur Luft- Und Raumfahrt E.V. Method of correcting range migration in image generation in synthetic aperture radar
US6836285B1 (en) * 1999-09-03 2004-12-28 Arete Associates Lidar with streak-tube imaging,including hazard detection in marine applications; related optics
US7535412B1 (en) * 2007-01-25 2009-05-19 The United States Of America As Represented By The Secretary Of The Navy Single pulse imaging radar system and method
CN103941257A (zh) * 2014-04-11 2014-07-23 哈尔滨工程大学 一种基于波数能量谱的导航雷达图像反演海面风向的方法
CN104076361A (zh) * 2014-07-04 2014-10-01 西安电子科技大学 无人机载战场监视雷达超分辨广域成像方法
CN106707278A (zh) * 2017-01-18 2017-05-24 北京无线电测量研究所 一种基于稀疏表示的多普勒波束锐化成像方法及装置

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2002071096A2 (en) * 2001-03-02 2002-09-12 Massachusetts Institute Of Technology High-definition imaging apparatus and method
WO2011083300A1 (en) * 2010-01-11 2011-07-14 Bae Systems Plc Pulse radar range profile motion compensation
US10073177B2 (en) * 2014-11-14 2018-09-11 Massachusetts Institute Of Technology Methods and apparatus for phased array imaging

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5237329A (en) * 1991-07-08 1993-08-17 Deutsche Forschungsanstalt Fur Luft- Und Raumfahrt E.V. Method of correcting range migration in image generation in synthetic aperture radar
US6836285B1 (en) * 1999-09-03 2004-12-28 Arete Associates Lidar with streak-tube imaging,including hazard detection in marine applications; related optics
US7535412B1 (en) * 2007-01-25 2009-05-19 The United States Of America As Represented By The Secretary Of The Navy Single pulse imaging radar system and method
CN103941257A (zh) * 2014-04-11 2014-07-23 哈尔滨工程大学 一种基于波数能量谱的导航雷达图像反演海面风向的方法
CN104076361A (zh) * 2014-07-04 2014-10-01 西安电子科技大学 无人机载战场监视雷达超分辨广域成像方法
CN106707278A (zh) * 2017-01-18 2017-05-24 北京无线电测量研究所 一种基于稀疏表示的多普勒波束锐化成像方法及装置

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
《Synthetic Aperture Ultrasound Fourier Beamformation Using Virtual Sources》;Elahe Moghimirad et.al;《IEEE TRANSACTIONS ON ULTRASONICS,FERROELECTRICS,AND FREQUENCY CONTROL》;20161231;第63卷(第12期);第2018-2030页 *
《分布式小卫星SAR实现全孔径分辨率的信号处理》;李真芳 等;《电子学报》;20031231;第31卷(第12期);第1800-1803页 *

Also Published As

Publication number Publication date
CN108614250A (zh) 2018-10-02

Similar Documents

Publication Publication Date Title
CN108459321B (zh) 基于距离-方位圆模型的大斜视高分辨率sar成像方法
CN111443339A (zh) 一种双基sar空变校正成像方法、装置、设备及存储介质
CN107340434B (zh) 基于内定标数据的星载合成孔径雷达天线方向图获取方法
CN105137408B (zh) 一种最优天线方向图选取的雷达角超分辨方法
CN107390181B (zh) 一种基于多波束扫描的雷达高分辨成像方法
CN113567982B (zh) 一种基于混合范数的方位向周期采样数据稀疏sar成像方法及装置
CN105223572B (zh) 一种基于pfa算法的正前视双基sar成像处理方法
CN111352107B (zh) 基于多通道数字和差的单脉冲跟踪与成像方法
CN111880180A (zh) 一种高分辨运动舰船sar成像的自聚焦方法
CN111025256A (zh) 一种机载雷达的微弱生命体征信号的检测方法及系统
CN109655819B (zh) 一种基于实孔径多普勒波束锐化的杂波抑制三维成像方法
Bishop et al. Video-SAR using higher order Taylor terms for differential range
CN106019280B (zh) 基于距离多普勒校正的fmcw sar成像方法及装置
CN108363046B (zh) 基于双接收机椭圆扫描的无源雷达测向方法及系统
CN111208514B (zh) 一种曲线运动轨迹SAR的切比雪夫斜距模型和Chirp Scaling成像方法
CN108107427A (zh) 基于超分辨技术的机载/弹载阵列雷达前视成像方法
CN110261852B (zh) 一种频率分集逆合成孔径雷达成像系统及方法
CN108562898B (zh) 一种前侧视sar的距离和方位两维空变自聚焦方法
CN112255596B (zh) 一种基于空间频率插值的星载sar欺骗干扰信号生成方法
CN112649806B (zh) 一种mimo雷达近场三维成像方法
CN108614250B (zh) 机载战场监视雷达的广域dbs图像拼接暗条纹校正方法
CN117192553A (zh) 一种多航过sar三维成像自聚焦方法
Al-Ibadi et al. DEM extraction of the basal topography of the Canadian archipelago ICE caps via 2D automated layer-tracker
CN113640794B (zh) 一种mimo-sar三维成像自聚焦方法
CN113419240B (zh) 基于双通道sar的动目标检测方法,双通道sar及存储介质

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