CN104076361B - 无人机载战场监视雷达超分辨广域成像方法 - Google Patents
无人机载战场监视雷达超分辨广域成像方法 Download PDFInfo
- Publication number
- CN104076361B CN104076361B CN201410317682.9A CN201410317682A CN104076361B CN 104076361 B CN104076361 B CN 104076361B CN 201410317682 A CN201410317682 A CN 201410317682A CN 104076361 B CN104076361 B CN 104076361B
- Authority
- CN
- China
- Prior art keywords
- data
- battlefield
- radar
- super
- aperture
- 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.)
- Expired - Fee Related
Links
- 238000003384 imaging method Methods 0.000 title claims abstract description 19
- 238000013213 extrapolation Methods 0.000 claims abstract description 36
- 238000012545 processing Methods 0.000 claims abstract description 24
- 238000000034 method Methods 0.000 claims abstract description 14
- 238000012544 monitoring process Methods 0.000 claims description 55
- 239000011159 matrix material Substances 0.000 claims description 14
- 238000012937 correction Methods 0.000 claims description 8
- 230000035485 pulse pressure Effects 0.000 claims description 8
- 230000006835 compression Effects 0.000 claims description 6
- 238000007906 compression Methods 0.000 claims description 6
- 238000005286 illumination Methods 0.000 claims description 2
- 230000017105 transposition Effects 0.000 claims description 2
- 241001272567 Hominoidea Species 0.000 claims 2
- 238000004088 simulation Methods 0.000 description 9
- 238000005516 engineering process Methods 0.000 description 6
- 238000007796 conventional method Methods 0.000 description 5
- 238000012986 modification Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 238000001514 detection method Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 102000016550 Complement Factor H Human genes 0.000 description 1
- 108010053085 Complement Factor H Proteins 0.000 description 1
- 238000009825 accumulation Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000001427 coherent effect Effects 0.000 description 1
- 230000021615 conjugation Effects 0.000 description 1
- 238000002592 echocardiography Methods 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
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
- G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
-
- 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
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/40—Means for monitoring or calibrating
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
本发明属于雷达超分辨广域成像技术领域,特别涉及无人机载战场监视雷达超分辨广域成像方法。该无人机载战场监视雷达超分辨广域成像方法包括以下步骤:利用战场监视雷达的每个照射区域的地面监视图像组合成广域地面监视图像;生成战场监视雷达的每个照射区域的地面监视图像包括以下步骤:战场监视雷达对对应的原始回波数据进行预处理,得出对应的中心补偿后数据;将对应的中心补偿后数据在方位向进行孔径外推,得出对应的孔径外推后数据;将对应的孔径外推后数据在方位向进行幅相估计法超分辨处理,得出对应的方位超分辨处理数据;根据对应的方位超分辨处理数据,得出对应的子图像。
Description
技术领域
本发明属于雷达超分辨广域成像技术领域,特别涉及无人机载战场监视雷达超分辨广域成像方法,可用于实时性的地面场景的快速广域成像与监视。
背景技术
合成孔径雷达(SAR)通过对地面场景进行长时间观测获得较高的方位分辨率。但由于这种体制不具备快速的数据重访能力,因此无法对远距离、大范围内的地面场景信息进行监视。与之相比,广域运动目标检测技术(Wide-area GMTI,Wide-area Ground MovingTarget Indication),采用较窄的天线波束和方位扫描的工作方式,可以实现对地面大场景的快速成像。由于其较高的采访率,可以对广域场景(200km×200km)反复进行探测,成为战场监视的最重要的工作方式。例如美国的JSTARS(Joint Surveillance Target AttackRadar System,联合星),德国的PAMIR(Phased Array Multifunctional Imaging Radar,多功能相控阵成像雷达),英国的ASTOR(Airborne Stand-off Radar机载防区外雷达)(即欧洲五国联合研制的SOSTAR(Stand-off Surveillance and Target Acquisition RadarDemonstrator防区外监视和目标获取雷达)***,均具备广域监视能力。JSTARS的广域监视能力在海湾战争中达到了很好的验证。
利用Wide-area GMTI模式进行广域监视时,通常采用多普勒波束锐化技术(DBS,Doppler Beam Sharpening)获得大范围的地面场景图像,其测绘带宽度时传统的SAR模式难以达到的。多普勒波束锐化技术是将一个天线的真实波束***成若干个窄的子波束,由于各子波束中心处目标相对雷达的径向速度不同,从而造成各子波束回波间的多普勒频差,此时在频域设置一组窄带滤波器组,且该滤波器的中心和带宽分别对应各个子波束的中心和带宽,这样就可以实现多普勒分割,有效改善方位分辨力。然而,由于广域GMTI的要求对数据的高重访率及其对场景的广域监视能力,因此雷达天线在进行扫面时,对某一方位角度的照射时间时非常短暂的,因此相干积累时间非常有限,这也大大限制的了多普勒波束锐化技术的成像分辨率。
发明内容
本发明的目的在于提出无人机载战场监视雷达超分辨广域成像方法。本发明能够将雷达方位分辨率提高一倍,在工程上实现对战场的实时监视能力。
为实现上述技术目的,本发明采用如下技术方案予以实现。
无人机载战场监视雷达超分辨广域成像方法包括以下步骤:针对战场监视雷达的每个照射区域,生成战场监视雷达的对应照射区域的地面监视图像,利用战场监视雷达的每个照射区域的地面监视图像组合成广域地面监视图像;生成战场监视雷达的每个照射区域的地面监视图像包括以下步骤:
步骤1,当战场监视雷达扫描到对应照射区域时,战场监视雷达接收对应的原始回波数据;战场监视雷达对对应的原始回波数据依次进行距离向脉冲压缩、距离走动校正、中心补偿,得出对应的中心补偿后数据;
步骤2,将对应的中心补偿后数据在方位向进行2倍的孔径外推,得出对应的孔径外推后数据;
步骤3,将对应的孔径外推后数据在方位向进行幅相估计法超分辨处理,得出对应的方位超分辨处理数据;
步骤4,针对对应的方位超分辨处理数据,在距离-多普勒域提取对应多普勒带宽内的数据;然后根据提取的数据,得出对应的子图像,对应子图像中每个像素点的幅度值为:提取的数据中对应位置的数据的取值。
本发明的特点和进一步改进在于:
所述步骤1的具体子步骤为:
(1.1)当战场监视雷达扫描到对应照射区域时,战场监视雷达接收对应的原始回波数据,然后对对应的原始回波数据进行距离向脉冲压缩处理,得出对应的距离脉压后数据;
(1.2)对对应的距离脉压后数据进行距离走动校正,得出对应的距离走动校正后数据;
(1.3)对对应的原始回波数据进行多普勒中心估计,得出对应的多普勒中心频率fd;
(1.4)根据对应的多普勒中心频率fd,对对应的距离走动校正后数据进行中心补偿,得出对应的中心补偿后数据。
在步骤1中,对应的中心补偿后数据表示为X,
X=[x(m,0),x(m,1),...,x(m,n),...,x(m,Na-1)]
其中,x(m,n)表示第m个距离单元第n个脉冲经过中心补偿后的脉压结果,m代表快时间,n代表慢时间,m为整数且m为0至M-1,M为战场监视雷达的距离单元数目,n为整数且n取0至Na-1,Na为战场监视雷达的方位向脉冲数目,Na为偶数;
所述步骤2的具体子步骤为:
(2.1)利用Burg算法求取对应的中心补偿后数据的孔径外推系数a(k);k为整数且k取0至Na-1;
(2.2)对对应的中心补偿后数据在前向进行孔径外推,得出对应的前向孔径外推数据,所述对应的前向孔径外推数据X1为:
其中,m代表快时间,n1代表慢时间,m为自然数且m取1至M,n1为自然数且n1取P+Na/2-1至P,P=Na/2,k'取1至P;
对对应的中心补偿后数据在后向进行孔径外推,得出对应的后向孔径外推数据,所述对应的后向孔径外推数据X2为:
其中,m代表快时间,n2代表慢时间,m为自然数且m取1至M,M为距离单元数目,n2为自然数且n2取P+Na/2至P+Na-1,P=Na/2,k'取1至P;上标H表示取共轭转置;
(2.3)根据对应的中心补偿后数据X、对应的前向孔径外推数据X1和对应的后向孔径外推数据X2,得出对应的孔径外推后数据X'
令则对应的孔径外推后数据X'为:
X'=[x(m,0),x(m,1),...,x(m,2Na-1)]。
所述步骤3的具体子步骤为:
(3.1)根据设定的有限冲击响应滤波器的阶数,对对应的孔径外推后数据X'进行重新排列,得出2Na-L+1组重新排列后数据,L为设定的FIR滤波器的阶数,L=Na;在2Na-L组重新排列后数据中,第l组数据Sm,l为:
Sm,l=[x(m,l),x(m,l+1),...,x(m,l+L-1)]T
其中,l取0至2Na-L,上标T表示矩阵或向量的转置;
(3.2)根据每组重新排列后数据,得出对应的数据序列的协方差矩阵
其中,N=2Na,上标H表示取共轭转置;
(3.3)估计对应的噪声和干扰的L维协方差矩阵
其中,N=2Na,上标H表示取共轭转置,说明w为APES超分辨处理后每组重新排列后数据的角频率;
(3.4)计算出对应的匹配滤波器系数
其中,上标-1表示矩阵的逆,β(w)=[1 exp(jw)…exp(j(L-1)w)]T,上标T表示矩阵或向量的转置,β(w)表示对应的归一化导向矢量;
(3.5)根据以下公式得出不同频率处的数据幅度估计值:
其中,w为APES超分辨处理后每组重新排列后数据的角频率,表示角频率w处的数据幅度估计值,则对应的方位超分辨处理数据为:子步骤(3.5)得出的不同频率处的数据幅度估计值。
利用战场监视雷达的每个照射区域的地面监视图像组合成广域地面监视图像的过程为:在得出战场监视雷达的每个照射区域的子图像之后,在直角坐标系下,按顺序完成各个子图像的拼接,得出广域地面监视图像。
本发明的有益效果为:本发明由于充分利用雷达孔径外推后的回波数据,在不引起数据失真的前提下将天线主波束内的回波数据量增加一倍,相当于将多普勒波束锐化技术的锐化比和方位分辨率都提高了一倍。同时,在进行方位向的多普勒分析时,采用APES(Amplitude and Phase Estimation幅相估计法)超分辨算法,这又大大提高了方位谱的分辨能力,因此理论上,该发明可以将传统的多普勒波束锐化技术横向分辨率提高一倍。实验结果表明,本发明能够应用在机械扫描雷达上,且能够得到较为理想的广域地面监视图像。
附图说明
图1为本发明的无人机载战场监视雷达超分辨广域成像方法的流程图;
图2为仿真实验中利用传统方法得出的广域地面监视图像的扇形图;
图3为仿真实验中利用本发明得出的广域地面监视图像的扇形图;
图4为图2的局部放大结果示意图;
图5为图3的局部放大结果示意图。
具体实施方式
下面结合附图对本发明作进一步说明:
参照图1,为本发明的无人机载战场监视雷达超分辨广域成像方法的流程图。该无人机载战场监视雷达超分辨广域成像方法包括以下步骤:当无人机飞行时,无人机上的战场监视雷达的照射区域发生连续性的变化;针对战场监视雷达的每个照射区域,生成战场监视雷达的对应照射区域的地面监视图像,利用战场监视雷达的每个照射区域的地面监视图像组合成广域地面监视图像;生成战场监视雷达的每个照射区域的地面监视图像包括以下步骤:
步骤1,当战场监视雷达(机械扫描雷达)扫描到对应照射区域时,战场监视雷达接收对应的原始回波数据;战场监视雷达对对应的原始回波数据依次进行距离向脉冲压缩、距离走动校正、中心补偿,得出对应的中心补偿后数据。
其具体子步骤为:
(1.1)当战场监视雷达(机械扫描雷达)扫描到对应照射区域时,战场监视雷达接收对应的原始回波数据,然后对对应的原始回波数据进行距离向脉冲压缩处理(匹配滤波处理),得出对应的距离脉压后数据,此时实现目标的距离向分离。
(1.2)对对应的距离脉压后数据进行距离走动校正,得出对应的距离走动校正后数据。具体地说,为减少雷达平台(无人机)运动所引起的目标跨距离单元走动现象,需要进行距离走动校正。对对应的距离脉压后数据在距离频域内乘以距离走动校正因子H(f,t),得出对应的距离走动校正后数据,H(f,t)为:
其中,R(t)≈R0-v·t·cosθ0,f为经过傅里叶变换后的距离频率,t为慢时间,R0为战场监视雷达与点目标之间的对应初始斜距,θ0为战场监视雷达的波束视线与载机运动方向的夹角,v为无人机(载机平台)的速度,c为光速。
(1.3)对对应的原始回波数据进行多普勒中心估计,得出对应的多普勒中心频率fd。具体地说,采用相关法对对应的原始回波数据进行多普勒中心估计。
(1.4)根据对应的多普勒中心频率fd,对对应的距离走动校正后数据进行中心补偿,得出对应的中心补偿后数据。具体地说,将对应的距离走动校正后数据乘以对应的多普勒中心补偿函数,得出对应的中心补偿后数据,对应的多普勒中心补偿函数为exp(-j2πfdk/fr),其中,k为整数且k取0至Na-1,Na为方位向脉冲数目,Na为偶数,fr为战场监视雷达的脉冲重复频率。
对应的中心补偿后数据X为:X=[x(m,0),x(m,1),...,x(m,n),...,x(m,Na-1)],x(m,n)表示第m个距离单元第n个脉冲经过中心补偿后的脉压结果,m代表快时间,n代表慢时间,m为自然数且m为0至M-1,M为距离单元数目,n为自然数且n取0至Na-1。
步骤2,将对应的中心补偿后数据在方位向进行2倍的孔径外推,得出对应的孔径外推后数据,在对应的孔径外推后数据中,对应的方位向脉冲数目变为2Na。
其具体子步骤为:
(2.1)利用Burg算法求取对应的中心补偿后数据的孔径外推系数a(k),在Burg算法中,以前向预测误差和后向预测误差的平均功率最小为准则求取孔径外推系数。
(2.2)对对应的中心补偿后数据在前向进行孔径外推,得出对应的前向孔径外推数据,所述对应的前向孔径外推数据X1为:
其中,m代表快时间,n1代表慢时间,m为自然数且m取1至M,n1为自然数且n1取P+Na/2-1至P,P=Na/2,k'取1至P;x(m,n)表示对应的中心补偿后数据。
对对应的中心补偿后数据在后向进行孔径外推,得出对应的后向孔径外推数据,所述对应的后向孔径外推数据X2为:
其中,m代表快时间,n2代表慢时间,m为自然数且m取1至M,M为距离单元数目,n2为自然数且n2取P+Na/2至P+Na-1,P=Na/2,k'取1至P;x(m,n)表示对应的中心补偿后数据,上标H表示取共轭转置。
(2.3)将对应的中心补偿后数据X、对应的前向孔径外推数据X1和对应的后向孔径外推数据X2进行合并,得出对应的孔径外推后数据X'。
令则对应的孔径外推后数据X'为:
X'=[x(m,0),x(m,1),...,x(m,2Na-1)]
步骤3,将对应的孔径外推后数据在方位向进行APES(Amplitude and PhaseEstimation幅相估计法)超分辨处理,得出对应的方位超分辨处理数据。
其具体子步骤为:
(3.1)根据设定的FIR(Finite Impulse Response,有限冲击响应滤波器)滤波器的阶数,对对应的孔径外推后数据X'进行重新排列,得出2Na-L+1组重新排列后数据,L为设定的FIR滤波器的阶数,本发明实施例中L=Na。在2Na-L组重新排列后数据中,第l组数据Sm,l为:
Sm,l=[x(m,l),x(m,l+1),...,x(m,l+L-1)]T
其中,l取0至2Na-L。
(3.2)根据每组重新排列后数据,得出对应的数据序列的协方差矩阵
其中,N=2Na,上标H表示取共轭转置。
(3.3)估计对应的噪声和干扰的L维协方差矩阵
其中,N=2Na,上标H表示取共轭转置,说明w为APES超分辨处理后每组重新排列后数据的角频率。
(3.4)计算出对应的匹配滤波器系数
其中,上标-1表示矩阵的逆,β(w)=[1 exp(jw)…exp(j(L-1)w)]T,上标T表示矩阵或向量的转置,β(w)表示对应的归一化导向矢量。
(3.5)实际应用中,在频率轴上利用谱峰搜索求出不同频率处的数据幅度估计值,即根据以下公式得出不同频率处的数据幅度估计值:
其中,w为APES超分辨处理后每组重新排列后数据的角频率,表示角频率w处的数据幅度估计值。
步骤4,针对对应的方位超分辨处理数据,在距离-多普勒域提取对应多普勒带宽内的数据;然后根据提取的数据,得出对应的子图像,对应子图像中每个像素点的幅度值为:提取的数据中对应位置的数据的取值。
在得出战场监视雷达的每个照射区域的子图像之后,在直角坐标系下(以第一幅子图为基准),按顺序完成各个子图像的拼接,得出广域地面监视图像。
本发明的效果可通过以下仿真实验进一步说明:
1)仿真实验环境和内容
仿真实验环境:MATLAB R2009b,Intel(R)Pentium(R)2CPU 2.7GHz,Window 7旗舰版。实验内容:将机载机械扫描雷达录取的回波数据,在仿真环境下,应用本发明和传统方法进行成像,得出对应的广域地面监视图像。
参照图2,为仿真实验中利用传统方法得出的广域地面监视图像的扇形图,参照图3,为仿真实验中利用本发明得出的广域地面监视图像的扇形图。然后在图2和图3中,选定共同区域进行放大,参照图4,为图2的局部放大结果示意图,参照图5,为图3的局部放大结果示意图。
2)仿真实验结果
从图2和图3可以看出,相较于传统方法,本发明的成像结果中,河流,村落,农田及周围地貌等细节信息更加丰富,说明本发明可以应用于机械扫描雷达。
从图4和图5可以看出,相较于传统方法,强散射点目标不仅得到了很好的聚焦,农田的边界信息也更加清晰,说明本发明可明显提高机载雷达的横向分辨率。
图2至图5的熵可以用下表进行说明:
图像的熵在一定程度上可以反映出反应图像的聚焦程度的好坏,熵越小,说明图像的聚焦效果越好。从表中可以看出,无论扇形图还是局部放大图,本发明都要比传统方法低1.5dB左右。
显然,本领域的技术人员可以对本发明进行各种改动和变型而不脱离本发明的精神和范围。这样,倘若本发明的这些修改和变型属于本发明权利要求及其等同技术的范围之内,则本发明也意图包含这些改动和变型在内。
Claims (3)
1.无人机载战场监视雷达超分辨广域成像方法,其特征在于,包括以下步骤:针对战场监视雷达的每个照射区域,生成战场监视雷达的对应照射区域的地面监视图像,利用战场监视雷达的每个照射区域的地面监视图像组合成广域地面监视图像;生成战场监视雷达的每个照射区域的地面监视图像包括以下步骤:
步骤1,当战场监视雷达扫描到对应照射区域时,战场监视雷达接收对应的原始回波数据;战场监视雷达对对应的原始回波数据依次进行距离向脉冲压缩、距离走动校正、中心补偿,得出对应的中心补偿后数据;
在步骤1中,对应的中心补偿后数据表示为X,
X=[x(m,0),x(m,1),...,x(m,n),...,x(m,Na-1)]
其中,x(m,n)表示第m个距离单元第n个脉冲经过中心补偿后的脉压结果,m代表快时间,n代表慢时间,m为自然数且m为0至M-1,M为战场监视雷达的距离单元数目,n为自然数且n取0至Na-1,Na为战场监视雷达的方位向脉冲数目,Na为偶数;
步骤2,将对应的中心补偿后数据在方位向进行2倍的孔径外推,得出对应的孔径外推后数据;
所述步骤2的具体子步骤为:
(2.1)利用Burg算法求取对应的中心补偿后数据的孔径外推系数a(k);k为整数且k取0至Na-1;
(2.2)对对应的中心补偿后数据在前向进行孔径外推,得出对应的前向孔径外推数据,所述对应的前向孔径外推数据X1为:
其中,m代表快时间,n1代表慢时间,m为自然数且m取1至M,n1为自然数且n1取P+Na/2-1至P,P=Na/2,k'取1至P;
对对应的中心补偿后数据在后向进行孔径外推,得出对应的后向孔径外推数据,所述对应的后向孔径外推数据X2为:
其中,m代表快时间,n2代表慢时间,m为自然数且m取1至M,M为距离单元数目,n2为自然数且n2取P+Na/2至P+Na-1,P=Na/2,k'取1至P;上标H表示取共轭转置;
(2.3)根据对应的中心补偿后数据X、对应的前向孔径外推数据X1和对应的后向孔径外推数据X2,得出对应的孔径外推后数据X'
令则对应的孔径外推后数据X'为:
X'=[x(m,0),x(m,1),...,x(m,2Na-1)];
步骤3,将对应的孔径外推后数据在方位向进行幅相估计法超分辨处理,得出对应的方位超分辨处理数据;
所述步骤3的具体子步骤为:
(3.1)根据设定的有限冲击响应滤波器的阶数,对对应的孔径外推后数据X'进行重新排列,得出2Na-L+1组重新排列后数据,L为设定的FIR滤波器的阶数,L=Na;在2Na-L组重新排列后数据中,第l组数据Sm,l为:
Sm,l=[x(m,l),x(m,l+1),...,x(m,l+L-1)]T
其中,l取0至2Na-L,上标T表示矩阵或向量的转置;
(3.2)根据每组重新排列后数据,得出对应的数据序列的协方差矩阵
其中,N=2Na,上标H表示取共轭转置;
(3.3)估计对应的噪声和干扰的L维协方差矩阵
其中,N=2Na,上标H表示取共轭转置,说明w为APES超分辨处理后每组重新排列后数据的角频率;
(3.4)计算出对应的匹配滤波器系数
其中,上标-1表示矩阵的逆,β(w)=[1 exp(jw) … exp(j(L-1)w)]T,上标T表示矩阵或向量的转置,β(w)表示对应的归一化导向矢量;
(3.5)根据以下公式得出不同频率处的数据幅度估计值:
其中,w为APES超分辨处理后每组重新排列后数据的角频率,表示角频率w处的数据幅度估计值,则对应的方位超分辨处理数据为:子步骤(3.5)得出的不同频率处的数据幅度估计值;
步骤4,针对对应的方位超分辨处理数据,在距离-多普勒域提取对应多普勒带宽内的数据;然后根据提取的数据,得出对应的子图像,对应子图像中每个像素点的幅度值为:提取的数据中对应位置的数据的取值。
2.如权利要求1所述的无人机载战场监视雷达超分辨广域成像方法,其特征在于,所述步骤1的具体子步骤为:
(1.1)当战场监视雷达扫描到对应照射区域时,战场监视雷达接收对应的原始回波数据,然后对对应的原始回波数据进行距离向脉冲压缩处理,得出对应的距离脉压后数据;
(1.2)对对应的距离脉压后数据进行距离走动校正,得出对应的距离走动校正后数据;
(1.3)对对应的原始回波数据进行多普勒中心估计,得出对应的多普勒中心频率fd;
(1.4)根据对应的多普勒中心频率fd,对对应的距离走动校正后数据进行中心补偿,得出对应的中心补偿后数据。
3.如权利要求1所述的无人机载战场监视雷达超分辨广域成像方法,其特征在于,利用战场监视雷达的每个照射区域的地面监视图像组合成广域地面监视图像的过程为:在得出战场监视雷达的每个照射区域的子图像之后,在直角坐标系下,按顺序完成各个子图像的拼接,得出广域地面监视图像。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410317682.9A CN104076361B (zh) | 2014-07-04 | 2014-07-04 | 无人机载战场监视雷达超分辨广域成像方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410317682.9A CN104076361B (zh) | 2014-07-04 | 2014-07-04 | 无人机载战场监视雷达超分辨广域成像方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104076361A CN104076361A (zh) | 2014-10-01 |
CN104076361B true CN104076361B (zh) | 2017-03-08 |
Family
ID=51597747
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410317682.9A Expired - Fee Related CN104076361B (zh) | 2014-07-04 | 2014-07-04 | 无人机载战场监视雷达超分辨广域成像方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104076361B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108614250B (zh) * | 2018-05-18 | 2022-02-22 | 西安电子科技大学 | 机载战场监视雷达的广域dbs图像拼接暗条纹校正方法 |
CN117237199B (zh) * | 2023-11-15 | 2024-01-26 | 中国科学院长春光学精密机械与物理研究所 | 基于无人机航拍生成模拟gmti雷达图像的方法 |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5448241A (en) * | 1994-05-26 | 1995-09-05 | Hughes Aircraft Company | Terrain height radar |
GB2430722B (en) * | 2001-09-26 | 2007-08-22 | Mbda Uk Ltd | A guidance system |
US7646327B1 (en) * | 2008-07-30 | 2010-01-12 | The United States Of America As Represented By The Secretary Of The Navy | Synthetic aperture radar motion estimation method |
CN102176018B (zh) * | 2011-03-15 | 2013-01-23 | 西安电子科技大学 | 机械扫描雷达的多普勒波束锐化快速成像方法 |
CN102778680B (zh) * | 2012-06-06 | 2014-03-19 | 西安电子科技大学 | 基于参数化的匀加速运动刚体群目标成像方法 |
CN103885058B (zh) * | 2014-02-26 | 2016-03-02 | 西安电子科技大学 | 一种利用稀疏统计特性的机载雷达前视超分辨成像方法 |
-
2014
- 2014-07-04 CN CN201410317682.9A patent/CN104076361B/zh not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CN104076361A (zh) | 2014-10-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103091674B (zh) | 基于hrrp序列的空间目标高分辨成像方法 | |
CN108051809B (zh) | 基于Radon变换的运动目标成像方法、装置及电子设备 | |
CN104931966B (zh) | 一种基于dcs算法的星载视频sar成像处理方法 | |
Park et al. | Cross-range scaling algorithm for ISAR images using 2-D Fourier transform and polar mapping | |
CN104698459B (zh) | 一种应用于缺失数据的条带sar压缩感知成像方法 | |
CN108427115B (zh) | 合成孔径雷达对运动目标参数的快速估计方法 | |
CN111142105A (zh) | 复杂运动目标isar成像方法 | |
CN110568434B (zh) | 一种多通道匀加速sar动目标二维速度估计方法 | |
CN114545411B (zh) | 一种基于工程实现的极坐标格式多模高分辨sar成像方法 | |
CN110082764B (zh) | 基于稳健正则化层析方法的sar图像成像方法 | |
CN109597075B (zh) | 一种基于稀疏阵列的成像方法及成像装置 | |
CN111880180A (zh) | 一种高分辨运动舰船sar成像的自聚焦方法 | |
CN109031299B (zh) | 低信噪比条件下基于相位差分的isar平动补偿方法 | |
CN102901966B (zh) | 基于解斜和楔石处理的动目标瞬时距离多普勒成像方法 | |
Yu et al. | Ground moving target motion parameter estimation using Radon modified Lv's distribution | |
CN109143236B (zh) | 适用于复杂飞行轨迹的双基聚束sar大场景成像方法 | |
CN108845318B (zh) | 基于Relax算法的星载高分宽幅成像方法 | |
Lu et al. | Resolution enhancement for forwarding looking multi-channel SAR imagery with exploiting space–time sparsity | |
CN104076361B (zh) | 无人机载战场监视雷达超分辨广域成像方法 | |
CN103278819A (zh) | 基于滑动接收窗的机载高分辨率斜视聚束sar成像方法 | |
CN106526544B (zh) | 基于高超声速平台的mimosar杂波抑制方法 | |
CN106772373B (zh) | 针对任意地面运动目标的sar成像方法 | |
CN105929381B (zh) | 一种机载sal振动估计方法 | |
Li et al. | Doppler keystone transform for SAR imaging of moving targets | |
CN113238229B (zh) | 一种geo星机双基sar无模糊成像方法 |
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: 20170308 Termination date: 20170704 |
|
CF01 | Termination of patent right due to non-payment of annual fee |