CN102721965A - 用于合成孔径雷达成像处理的全向旁瓣抑制方法 - Google Patents

用于合成孔径雷达成像处理的全向旁瓣抑制方法 Download PDF

Info

Publication number
CN102721965A
CN102721965A CN2012101898639A CN201210189863A CN102721965A CN 102721965 A CN102721965 A CN 102721965A CN 2012101898639 A CN2012101898639 A CN 2012101898639A CN 201210189863 A CN201210189863 A CN 201210189863A CN 102721965 A CN102721965 A CN 102721965A
Authority
CN
China
Prior art keywords
orientation
distance
grid
imaging region
expression
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
Application number
CN2012101898639A
Other languages
English (en)
Other versions
CN102721965B (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.)
Beihang University
Original Assignee
Beihang 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 Beihang University filed Critical Beihang University
Priority to CN 201210189863 priority Critical patent/CN102721965B/zh
Publication of CN102721965A publication Critical patent/CN102721965A/zh
Application granted granted Critical
Publication of CN102721965B publication Critical patent/CN102721965B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明公开了一种用于合成孔径雷达成像处理的全向旁瓣抑制方法，该方法首先将合成孔径雷达获取的二维回波信号排列为一维数据，然后确定成像平面区域位置，构造成像平面对应的测量矩阵，最后采用线性最小二乘方法，结合测量矩阵对一维回波数据进行处理，从而实现了合成孔径雷达成像结果的全向旁瓣抑制，显著降低了能量泄漏。从而使得处理后的星载回波数据的成像结果的能量泄露最小，进而达到抑制旁瓣的目的。本发明对目标散射强度没有特殊要求，可以对任意场景目标进行成像处理，适用范围广泛。

Description

用于合成孔径雷达成像处理的全向旁瓣抑制方法
技术领域
本发明涉及一种合成孔径雷达的信号处理方法，更特别地说，是指一种用于合成孔径雷达成像处理的全向旁瓣抑制方法。
背景技术
合成孔径雷达(Synthetic Aperture Radar，SAR)系统是一种有源的遥感设备，它主动向目标发射电磁波，利用目标散射特性作用下接收到回波信号进行成像。星载合成孔径雷达是一个非常复杂的系统，系统参数众多且相互关联、相互制约。在距离向，雷达发射线性调频信号，经过对回波信号进行脉冲压缩，从而得到较高的距离向分辨率；在方位向，通过卫星的雷达平台在前进过程波束与目标相对位置的变化，形成等效线性阵列天线，从而提高了方位分辨率。
影响星载SAR的图像质量的关键因素之一就是方位模糊性问题。这与星载SAR的脉冲工作方式有关，其脉冲重复频率(Pulse Repetitive Frequency，PRF)的选择与方位模糊紧密相关。由于多普勒效应，SAR所收到的回波沿方位向可以等效为一个线性调频(linear frequency modulation，LFM)信号，并且这个LFM信号是离散的，抽样频率为PRF。方位向LFM信号是被方位向的天线方向图加权的，但是天线方向图都存在着旁瓣且处于较高频段上。
1999年10月哈尔滨工业大学出版社出版、刘永坦编著的《雷达成像技术》指出，合成孔径雷达(SAR，Synthetic Aperture Radar)安装在运动平台上，按照一定的重复频率发射、接收脉冲，形成回波信号。SAR系统的结构框图如图1所示，SAR系统包括有星上雷达系统、卫星平台及数据下传系统和地面系统三部分，对合成孔径雷达成像处理是在地面系统中完成的。地面系统通过地面接收站接收卫星平台及数据下传系统下发的回波信号，该回波信号经SAR信号处理器进行成像处理，获得SAR图像；所述SAR图像存储于备档操作系统中。
2001年第49卷《IEEE Transaction on Antennas and Propagation》XuXiaojian、Narayanan R.M.发表的《《Range sidelobe suppression technique forcoherent ultra wide-band random noise radar imaging》指出，SAR信号处理器对回波信号分别沿距离向和方位向进行匹配滤波，获取成像结果。由于成像过程中二维频域的支撑域有限，使得合成孔径雷达的冲激响应在距离向和方位向上均为辛格函数，其包络具有逐渐递减的旁瓣。这种现象的存在使得在成像结果中后向散射强度较弱的目标被淹没在相邻强散射目标的旁瓣中，降低了合成孔径雷达系统对多个相邻目标的分辨能力。
1989年11月科学出版社出版、张澄波编著的《综合孔径雷达原理、系统分析与应用》中介绍了基于幅度加权的旁瓣抑制方法。该方法能够在一定程度上压低旁瓣电平，但会导致成像分辨率变差和信噪比降低。2009年第7期《信号处理》王建、周智敏、宋千等发表的《SAR图像二维旁瓣自适应抑制技术》提出了空变切趾滤波技术，采用非线性方法得到较低的旁瓣电平。该方法对强点目标的旁瓣抑制效果显著，但是不适用于弱目标，并且运算量较大。
最小二乘法(least square method)是一种数学优化技术，它通过最小化误差的平方寻找数据的最佳函数匹配。2006年9月科学出版社出版、魏木生编著的《广义最小二乘问题理论和计算》中给出了线性最小二乘法概念，当给定的矩阵
Figure BDA00001745888000021
非奇异，则对任何向量
Figure BDA00001745888000022
和线性方程组Ax＝y总有唯一解x＝A-1y，其中
Figure BDA00001745888000024
表示所有m×l复元素矩阵的实体，
Figure BDA00001745888000025
和
Figure BDA00001745888000026
分别表示所有m、n个复元素向量的实体。
发明内容
本发明的目的是提出一种适用于合成孔径雷达成像处理的全向旁瓣抑制方法，该方法首先将合成孔径雷达获取的二维回波信号排列为一维数据，然后确定成像平面区域位置，构造成像平面对应的测量矩阵，最后采用线性最小二乘方法，结合测量矩阵对一维回波数据进行处理，从而实现了合成孔径雷达成像结果的全向旁瓣抑制，显著降低了能量泄漏。
本发明是基于现有的SAR信号处理器对接收到的回波信号进行全向旁瓣抑制，具体的全向旁瓣抑制方法包括有下列处理步骤：
第一步：回波信号的二维转一维
第二步：成像区域划分
第三步：构建观测向量
第四步：观测向量构建
第五步：构建成像区域对应的观测矩阵
第六步：获取成像区域散射系数
第七步：向量矩阵转换
在本发明中，经SAR信号处理器采用全向旁瓣抑制处理后的图像，应用的图像评价参数包括有方位向能量泄漏比
Figure BDA00001745888000031
距离向能量泄漏比
Figure BDA00001745888000032
和二维能量泄漏比
Figure BDA00001745888000033
在合成孔径雷达旁瓣抑制处理中，应用本发明的方法具有如下优点：
①将回波信号由二维排列为一维，构建成像区域对应的观测矩阵，采用线性最小二乘法进行处理，消弱了合成孔径雷达沿方位向、距离向以及交叉向的旁瓣。
②与基于幅度加权的旁瓣抑制方法相比，本发明在抑制旁瓣的同时，不会引起合成孔径雷达空间分辨率的损失。
③与空变切趾滤波技术相比，本发明对目标的散射强度没有特殊要求，适用于各种类型的成像区域。
附图说明
图1是传统SAR系统的结构框图。
图1A是本发明全向旁瓣抑制方法的处理流程图。
图2是BP成像方法对点目标回波信号处理后得到幅值归一化三维图像。
图3是本发明方法对点目标回波信号处理后得到幅值归一化三维图像。
图4是BP成像方法对点目标回波信号处理后得到幅值归一化的方位向剖面图。
图5是BP成像方法对点目标回波信号处理后得到幅值归一化的距离向剖面图。
图6是本发明方法对点目标回波信号处理后得到幅值归一化的方位向剖面图。
图7是本发明方法对点目标回波信号处理后得到幅值归一化的距离向剖面图。
具体实施方式
下面将结合附图和仿真实例对本发明做进一步的详细说明。
参见图1A所示，本发明是基于现有的SAR信号处理器对接收到的回波信号进行全向旁瓣抑制，具体的全向旁瓣抑制方法包括有下列处理步骤：
第一步：回波信号的二维转一维
由于SAR信号处理器接到的回波数据是以二维矩阵进行排列的，为了实现对回波数据形成的成像区域进行处理，本发明采用对二维回波数据E二维回波进行一维回波数据E一维回波的重新排列。
在本发明中，合成孔径雷达的二维回波数据E二维回波采用矩阵形式表示为：
Figure BDA00001745888000041
在本发明中，合成孔径雷达的一维回波数据E一维回波采用矩阵形式表示为：
N方位表示沿方位向的采样点数；N距离表示沿距离向的采样点数；
a1-1表示沿方位向上采集到的第1个方位时刻的距离向的第1个采样点；
a1-2表示沿方位向上采集到的第1个方位时刻的距离向的第2个采样点；
Figure BDA00001745888000043
表示沿方位向上采集到的第1个方位时刻的距离向的第N距离个采样点；
a2-1表示沿方位向上采集到的第2个方位时刻的距离向的第1个采样点；
a2-2表示沿方位向上采集到的第2个方位时刻的距离向的第2个采样点；
Figure BDA00001745888000051
表示沿方位向上采集到的第2个方位时刻的距离向的第N距离个采样点；
表示沿方位向上采集到的第N方位个方位时刻的距离向的第1个采样点；
Figure BDA00001745888000053
表示沿方位向上采集到的第N方位个方位时刻的距离向的第2个采样点；
Figure BDA00001745888000054
表示沿方位向上采集到的第N方位个方位时刻的距离向的第N距离个采样点。
第二步：成像区域划分
在本发明中，SAR信号处理器对接到的回波数据采用网格剖分进行成像区域划分，得到带网格成像区域；即一个成像区域可以划分为D方位×D距离个带网格成像区域。
在带网格成像区域中，方位向的网格数目中任意一个网格点记为Ia，简称为方位向网格点Ia；
在带网格成像区域中，距离向的网格数目
Figure BDA00001745888000056
中任意一个网格点记为Ir，简称为距离向网格点Ir；
因此，方位向网格点Ia与距离向网格点Ir的网格中心坐标记为
Figure BDA00001745888000057
D方位表示方位向的网格数目，且
Figure BDA00001745888000058
网格数目向上取整；
D距离表示距离向的网格数目，且
Figure BDA00001745888000059
网格数目向上取整；
λ表示雷达系统波长；
Rmin表示下传天线到成像区域的最短距离；
La表示星上雷达系统中方位向天线尺寸；
V表示卫星平台飞行时的等效速度；
PRT表示脉冲重复时间间隔；
τp表示发射脉冲的时间宽度；
Fs表示距离向脉冲采样率；
c表示光速。
在本发明中，对成像区域进行网格划分，划分后的成像区域(带网格成像区域)中的所有点能够满足星载合成孔径雷达全孔径照射。
第三步：构建观测向量
在卫星平台飞行的任意时刻tk下对第二步骤中获得的网格中心坐标W网格进行斜距计算，该斜距
Figure BDA00001745888000061
在tk时刻下，第二步骤中获得的网格中心坐标W网格对应的观测向量为
Figure BDA00001745888000062
ek-1表示在tk时刻下，网格中心坐标W网格的距离向的第1个采样点观测值；
ek-2表示在tk时刻下，网格中心坐标W网格的距离向的第2个采样点观测值；
ek-n表示在tk时刻下，网格中心坐标W网格的距离向的第n个采样点观测值；
Figure BDA00001745888000063
表示在tk时刻下，网格中心坐标W网格的距离向的第N距离个采样点观测值。
在本发明中，距离向的第n个采样点观测值
Figure BDA00001745888000064
Nbegin≤n≤Nend；
n表示网格中心坐标W网格的距离向的采样点；
Nbegin表示距离向起点，且
Figure BDA00001745888000065
Nend表示距离向终点，且
Figure BDA00001745888000066
i表示虚数单位；
b表示发射信号调频率；
Fs表示距离向脉冲采样率；
τp表示发射脉冲的时间宽度；
Figure BDA00001745888000067
表示网格中心坐标W网格的斜距；
c表示光速；
Rmin表示下传天线到成像区域的最短距离。
第四步：观测向量构建
因为回波数据和成像区域散射点之间的相关性，为了得到带网格成像区域中每个网格点的散射系数，需要构建在全部运动时间内的方位向第Ia个、距离向第Ir个网格中心对应的观测向量
Figure BDA00001745888000071
所述网格中心对应的观测向量中
Figure BDA00001745888000073
表示
Figure BDA00001745888000074
的转置；
Figure BDA00001745888000075
表示的转置；表示
Figure BDA00001745888000078
的转置；表示
Figure BDA000017458880000710
的转置。
Figure BDA000017458880000711
表示在t1时刻成像区域的方位向第Ia个网格、距离向第Ir个网格中心点对应的观测向量；
Figure BDA000017458880000712
表示在t2时刻成像区域的方位向第Ia个网格、距离向第Ir个网格中心点对应的观测向量；
Figure BDA000017458880000713
表示在tk时刻成像区域的方位向第Ia个网格、距离向第Ir个网格中心点对应的观测向量；
Figure BDA000017458880000714
表示在时刻成像区域的方位向第Ia个网格、距离向第Ir个网格中心点对应的观测向量。
第五步：构建成像区域对应的观测矩阵
因为回波数据和成像区域散射点之间的相关性，为了得到带网格成像区域中所有网格点的散射系数，需要在本发明中构建带网格成像区域对应的观测矩阵φ。
所述观测矩阵
Figure BDA00001745888000081
中
Figure BDA00001745888000082
表示在全部运动时间内成像区域的方位向第1个网格、距离向第1个网格中心点的观测向量；
表示在全部运动时间内成像区域的方位向第1个网格、距离向第2个网格中心点的观测向量；
Figure BDA00001745888000084
表示在全部运动时间内成像区域的方位向第1个网格、距离向第D距离个网格中心点的观测向量；
Figure BDA00001745888000085
表示在全部运动时间内成像区域的方位向第2个网格、距离向第1个网格中心点的观测向量；
Figure BDA00001745888000086
表示在全部运动时间内成像区域的方位向第2个网格、距离向第2个网格中心点的观测向量；
Figure BDA00001745888000087
表示在全部运动时间内成像区域的方位向第2个网格、距离向第D距离个网格中心点的观测向量；
表示在全部运动时间内成像区域的方位向第D方位个网格、距离向第1个网格中心点的观测向量；
Figure BDA00001745888000089
表示在全部运动时间内成像区域的方位向第D方位个网格、距离向第2个网格中心点的观测向量；
Figure BDA000017458880000810
表示在全部运动时间内成像区域的方位向第D方位个网格、距离向第D距离个网格中心点的观测向量。
第六步：获取成像区域散射系数
依据第一步的一维回波数据E一维回波和第五步的观测矩阵φ按照线性最小二乘法求解E一维散射＝(φH×φ)-1×φH×E一维回波，从而得到成像区域散射系数，即一维散射系数E一维散射。其中，φH表示φ的共轭转置。
所述一维散射系数
Figure BDA00001745888000091
S1-1表示成像区域的方位向第1个网格、距离向第1个网格中心点的散射系数；
S1-2表示成像区域的方位向第1个网格、距离向第2个网格中心点的散射系数；
Figure BDA00001745888000092
表示成像区域的方位向第1个网格、距离向第D距离个网格中心点的散射系数；
S2-1表示成像区域的方位向第2个网格、距离向第1个网格中心点的散射系数；
S2-2表示成像区域的方位向第2个网格、距离向第2个网格中心点的散射系数；
表示成像区域的方位向第2个网格、距离向第D距离个网格中心点的散射系数；
Figure BDA00001745888000094
表示成像区域的方位向第D方位个网格、距离向第1个网格中心点的散射系数；
Figure BDA00001745888000095
表示成像区域的方位向第D方位个网格、距离向第2个网格中心点的散射系数；
Figure BDA00001745888000096
表示成像区域的方位向第D方位个网格、距离向第D距离个网格中心点的散射系数。
在本发明中，利用最小二乘法进行成像处理，该方法的优点是能够准确估计图像中目标的散射系数。
第七步：向量矩阵转换
将第六步的一维散射系数E一维散射由一维向量转换为二维散射矩阵E二维散射，并对矩阵EF求绝对值，得到成像区域的成像结果Iabs。
所述二维散射矩阵
Figure BDA00001745888000101
所述成像区域的成像结果
Figure BDA00001745888000102
在某型机载SAR中，将本发明的扩展非线性变标方法下载至信号处理系统中，然后对合成孔径雷达前斜视子孔径成像进行成像处理。
采用本发明的全向旁瓣抑制方法进行仿真处理，其参数设置如表1：
表1成像参数设置
  雷达波长λ   0.03m
  方位向天线长度La   1m
  脉冲宽度τp   1μs
  发射信号调频率b   150MHz/μs
  采样率Fs   160MHz
  脉冲重复时间间隔PRT   0.0033s
  最短距离Rmin   7km
  飞行速度V   160m/s
  方位向采样点N方位   512
  距离向采样点N距离   256
 成像区域方位向网格点D方位   164
 成像区域距离向网格点D距离   96
图2、图4、图5是采用传统的BP成像方法(见朱国富2001年在信息处理的第十七卷第五期发表的《超宽带LFM信号的BP成像算法》)对点目标回波信号处理后得到的成像结果。
图3、图6、图7是采用本发明的全向旁瓣抑制方法对点目标回波信号处理后得到的成像结果。
将图3、图6、图7分别于图2、图4、图5进行相比较，可以看出本发明比传统BP成像算法的方位向、距离向和二维能量泄漏比值小很多。评估指标结果列于表2。可见，本发明方法实现了全向旁瓣抑制。
表2点目标评估指标结果
Figure BDA00001745888000111
在本发明中，关键的SAR图像参数有距离和方位分辨率、峰值旁瓣比、积分旁瓣比等。由于本发明抑制旁瓣效果非常好，几乎没有任何旁瓣(如图3所示)，因而常规的图像评价参数不再适合本发明方法处理后的图像，特此定义四个新评价参数。本发明中方位向能量泄漏比
Figure BDA00001745888000112
距离向能量泄漏比和二维能量泄漏比其中，Pa表示方位向成像结果中除峰值所在位置外所有像素点的能量之和，Pap表示方位向成像结果中峰值的能量；Pr表示距离向成像结果中除峰值所在位置外所有像素点的能量之和，Prp表示距离向成像结果中峰值的能量；Pe表示二维成像结果中除峰值所在位置外所有像素点的能量之和，Pep表示二维成像结果中峰值的能量。

Claims (2)

1.一种用于合成孔径雷达成像处理的全向旁瓣抑制方法，其特征在于：是基于现有的SAR信号处理器对接收到的回波信号进行全向旁瓣抑制，具体的全向旁瓣抑制方法包括有下列处理步骤：
第一步：回波信号的二维转一维
在SAR信号处理器中将合成孔径雷达的二维回波数据E二维回波重新排列为一维回波数据E一维回波；
所述的二维回波数据E二维回波采用矩阵形式表示为
Figure FDA00001745887900011
所述的一维回波数据E一维回波采用矩阵形式表示为
Figure FDA00001745887900012
N方位表示沿方位向的采样点数；N距离表示沿距离向的采样点数；
a1-1表示沿方位向上采集到的第1个方位时刻的距离向的第1个采样点；
a1-2表示沿方位向上采集到的第1个方位时刻的距离向的第2个采样点；
表示沿方位向上采集到的第1个方位时刻的距离向的第N距离个采样点；
a2-1表示沿方位向上采集到的第2个方位时刻的距离向的第1个采样点；
a2-2表示沿方位向上采集到的第2个方位时刻的距离向的第2个采样点；
Figure FDA00001745887900021
表示沿方位向上采集到的第2个方位时刻的距离向的第N距离个采样点；
Figure FDA00001745887900022
表示沿方位向上采集到的第N方位个方位时刻的距离向的第1个采样点；
表示沿方位向上采集到的第N方位个方位时刻的距离向的第2个采样点；
Figure FDA00001745887900024
表示沿方位向上采集到的第N方位个方位时刻的距离向的第N距离个采样点；
第二步：成像区域划分
在SAR信号处理器中对回波数据形成的图像区域采用网格剖分进行成像区域划分，得到带网格成像区域；即一个成像区域可以划分为D方位×D距离个带网格成像区域；
在带网格成像区域中，方位向的网格数目
Figure FDA00001745887900025
中任意一个网格点记为Ia，简称为方位向网格点Ia；
在带网格成像区域中，距离向的网格数目
Figure FDA00001745887900026
中任意一个网格点记为Ir，简称为距离向网格点Ir；
因此，方位向网格点Ia与距离向网格点Ir的网格中心坐标记为
D方位表示方位向的网格数目，且
Figure FDA00001745887900028
网格数目向上取整；
D距离表示距离向的网格数目，且
Figure FDA00001745887900029
网格数目向上取整；
λ表示雷达系统波长；
Rmin表示下传天线到成像区域的最短距离；
La表示星上雷达系统中方位向天线尺寸；
V表示卫星平台飞行时的等效速度；
PRT表示脉冲重复时间间隔；
τp表示发射脉冲的时间宽度；
Fs表示距离向脉冲采样率；
c表示光速；
第三步：构建观测向量
在卫星平台飞行的任意时刻tk下对第二步骤中获得的网格中心坐标W网格进行斜距计算，该斜距
Figure FDA00001745887900031
在tk时刻下，第二步骤中获得的网格中心坐标W网格对应的观测向量为
ek-1表示在tk时刻下，网格中心坐标W网格的距离向的第1个采样点观测值；
ek-2表示在tk时刻下，网格中心坐标W网格的距离向的第2个采样点观测值；
ek-n表示在tk时刻下，网格中心坐标W网格的距离向的第n个采样点观测值；
Figure FDA00001745887900033
表示在tk时刻下，网格中心坐标W网格的距离向的第N距离个采样点观测值；
所述距离向的第n个采样点观测值
Figure FDA00001745887900034
Nbegin≤n≤Nend；
n表示网格中心坐标W网格的距离向的采样点；
Nbegin表示距离向起点，且
Figure FDA00001745887900035
Nend表示距离向终点，且
Figure FDA00001745887900036
i表示虚数单位；
b表示发射信号调频率；
Fs表示距离向脉冲采样率；
τp表示发射脉冲的时间宽度；
Figure FDA00001745887900037
表示网格中心坐标W网格的斜距；
c表示光速；
Rmin表示下传天线到成像区域的最短距离；
第四步：观测向量构建
因为回波数据和成像区域散射点之间的相关性，为了得到带网格成像区域中每个网格点的散射系数，需要构建在全部运动时间内的方位向第Ia个、距离向第Ir个网格中心对应的观测向量
Figure FDA00001745887900038
所述网格中心对应的观测向量
Figure FDA00001745887900041
中表示
Figure FDA00001745887900043
的转置；
Figure FDA00001745887900044
表示
Figure FDA00001745887900045
的转置；
Figure FDA00001745887900046
表示
Figure FDA00001745887900047
的转置；
Figure FDA00001745887900048
表示的转置；
表示在t1时刻成像区域的方位向第Ia个网格、距离向第Ir个网格中心点对应的观测向量；
Figure FDA000017458879000411
表示在t2时刻成像区域的方位向第Ia个网格、距离向第Ir个网格中心点对应的观测向量；
Figure FDA000017458879000412
表示在tk时刻成像区域的方位向第Ia个网格、距离向第Ir个网格中心点对应的观测向量；
Figure FDA000017458879000413
表示在
Figure FDA000017458879000414
时刻成像区域的方位向第Ia个网格、距离向第Ir个网格中心点对应的观测向量；
第五步：构建成像区域对应的观测矩阵
因为回波数据和成像区域散射点之间的相关性，为了得到带网格成像区域中所有网格点的散射系数，需要构建带网格成像区域对应的观测矩阵φ；
所述观测矩阵
Figure FDA000017458879000415
中
Figure FDA000017458879000416
表示在全部运动时间内成像区域的方位向第1个网格、距离向第1个网格中心点的观测向量；
Figure FDA000017458879000417
表示在全部运动时间内成像区域的方位向第1个网格、距离向第2个网格中心点的观测向量；
Figure FDA000017458879000418
表示在全部运动时间内成像区域的方位向第1个网格、距离向第D距离个网格中心点的观测向量；
表示在全部运动时间内成像区域的方位向第2个网格、距离向第1个网格中心点的观测向量；
Figure FDA000017458879000420
表示在全部运动时间内成像区域的方位向第2个网格、距离向第2个网格中心点的观测向量；
Figure FDA00001745887900051
表示在全部运动时间内成像区域的方位向第2个网格、距离向第D距离个网格中心点的观测向量；
Figure FDA00001745887900052
表示在全部运动时间内成像区域的方位向第D方位个网格、距离向第1个网格中心点的观测向量；
Figure FDA00001745887900053
表示在全部运动时间内成像区域的方位向第D方位个网格、距离向第2个网格中心点的观测向量；
Figure FDA00001745887900054
表示在全部运动时间内成像区域的方位向第D方位个网格、距离向第D距离个网格中心点的观测向量；
第六步：获取成像区域散射系数
依据第一步的一维回波数据E一维回波和第五步的观测矩阵φ按照线性最小二乘法求解E一维散射＝(φH×φ)-1×φH×E一维回波，从而得到成像区域散射系数，即一维散射系数E一维散射；其中，φH表示φ的共轭转置；
所述一维散射系数
Figure FDA00001745887900055
S1-1表示成像区域的方位向第1个网格、距离向第1个网格中心点的散射系数；
S1-2表示成像区域的方位向第1个网格、距离向第2个网格中心点的散射系数；
Figure FDA00001745887900056
表示成像区域的方位向第1个网格、距离向第D距离个网格中心点的散射系数；
S2-1表示成像区域的方位向第2个网格、距离向第1个网格中心点的散射系数；
S2-2表示成像区域的方位向第2个网格、距离向第2个网格中心点的散射系数；
Figure FDA00001745887900061
表示成像区域的方位向第2个网格、距离向第D距离个网格中心点的散射系数；
Figure FDA00001745887900062
表示成像区域的方位向第D方位个网格、距离向第1个网格中心点的散射系数；
表示成像区域的方位向第D方位个网格、距离向第2个网格中心点的散射系数；
Figure FDA00001745887900064
表示成像区域的方位向第D方位个网格、距离向第D距离个网格中心点的散射系数；
第七步：向量矩阵转换
将第六步的一维散射系数E一维散射由一维向量转换为二维散射矩阵E二维散射，并对矩阵EF求绝对值，得到成像区域的成像结果Iabs；
所述二维散射矩阵
Figure FDA00001745887900065
所述成像区域的成像结果
Figure FDA00001745887900066
2.根据权利要求1所述的用于合成孔径雷达成像处理的全向旁瓣抑制方法，其特征在于：SAR信号处理器采用全向旁瓣抑制处理后的图像，应用的图像评价参数包括有方位向能量泄漏比
Figure FDA00001745887900067
距离向能量泄漏比
Figure FDA00001745887900068
和二维能量泄漏比其中，Pa表示方位向成像结果中除峰值所在位置外所有像素点的能量之和，Pap表示方位向成像结果中峰值的能量；Pr表示距离向成像结果中除峰值所在位置外所有像素点的能量之和，Prp表示距离向成像结果中峰值的能量；Pe表示二维成像结果中除峰值所在位置外所有像素点的能量之和，Pep表示二维成像结果中峰值的能量。
CN 201210189863 2012-06-08 2012-06-08 用于合成孔径雷达成像处理的全向旁瓣抑制方法 Expired - Fee Related CN102721965B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN 201210189863 CN102721965B (zh) 2012-06-08 2012-06-08 用于合成孔径雷达成像处理的全向旁瓣抑制方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN 201210189863 CN102721965B (zh) 2012-06-08 2012-06-08 用于合成孔径雷达成像处理的全向旁瓣抑制方法

Publications (2)

Publication Number Publication Date
CN102721965A true CN102721965A (zh) 2012-10-10
CN102721965B CN102721965B (zh) 2013-08-21

Family

ID=46947781

Family Applications (1)

Application Number Title Priority Date Filing Date
CN 201210189863 Expired - Fee Related CN102721965B (zh) 2012-06-08 2012-06-08 用于合成孔径雷达成像处理的全向旁瓣抑制方法

Country Status (1)

Country Link
CN (1) CN102721965B (zh)

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103149457A (zh) * 2013-03-01 2013-06-12 西安电子科技大学 空间电磁强度分布分析方法
CN103323820A (zh) * 2013-06-19 2013-09-25 北京理工大学 一种抑制频率步进sar一维距离图像栅瓣的方法
CN104635230A (zh) * 2015-03-05 2015-05-20 北京航空航天大学 一种用于mimo-sar近场测量成像方位向旁瓣抑制的方法
CN106772276A (zh) * 2016-12-21 2017-05-31 南京信息工程大学 一种地球同步轨道圆迹sar水平面二维旁瓣抑制方法
CN108318891A (zh) * 2017-11-28 2018-07-24 西安电子科技大学 一种基于改进sva和cs的sal数据旁瓣的压低方法
CN109471097A (zh) * 2018-11-14 2019-03-15 湖南华诺星空电子技术有限公司 一种穿墙雷达信号优化处理方法及装置
CN110333489A (zh) * 2019-07-24 2019-10-15 北京航空航天大学 采用cnn与rsva组合对sar回波数据旁瓣抑制的处理方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US3950747A (en) * 1975-02-28 1976-04-13 International Telephone & Telegraph Corporation Optical processing system for synthetic aperture radar
JP2001091651A (ja) * 1999-09-22 2001-04-06 Toshiba Corp 合成開口レーダ
JP2001281332A (ja) * 2000-03-30 2001-10-10 Mitsubishi Electric Corp 合成開口レーダ装置および目標画像再生方法
CN101839982A (zh) * 2010-05-18 2010-09-22 中国人民解放军国防科学技术大学 一种合成孔径雷达图像旁瓣抑制方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US3950747A (en) * 1975-02-28 1976-04-13 International Telephone & Telegraph Corporation Optical processing system for synthetic aperture radar
JP2001091651A (ja) * 1999-09-22 2001-04-06 Toshiba Corp 合成開口レーダ
JP2001281332A (ja) * 2000-03-30 2001-10-10 Mitsubishi Electric Corp 合成開口レーダ装置および目標画像再生方法
CN101839982A (zh) * 2010-05-18 2010-09-22 中国人民解放军国防科学技术大学 一种合成孔径雷达图像旁瓣抑制方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
张先义 等: "切趾滤波旁瓣抑制技术在合成孔径雷达中的应用", 《电子与信息学报》 *
田旭文 等: "合成孔径雷达图像的自适应旁瓣抑制算法分析和实验研究", 《遥感技术与应用》 *

Cited By (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103149457A (zh) * 2013-03-01 2013-06-12 西安电子科技大学 空间电磁强度分布分析方法
CN103323820A (zh) * 2013-06-19 2013-09-25 北京理工大学 一种抑制频率步进sar一维距离图像栅瓣的方法
CN103323820B (zh) * 2013-06-19 2015-07-15 北京理工大学 一种抑制频率步进sar一维距离图像栅瓣的方法
CN104635230A (zh) * 2015-03-05 2015-05-20 北京航空航天大学 一种用于mimo-sar近场测量成像方位向旁瓣抑制的方法
CN106772276A (zh) * 2016-12-21 2017-05-31 南京信息工程大学 一种地球同步轨道圆迹sar水平面二维旁瓣抑制方法
CN106772276B (zh) * 2016-12-21 2019-05-21 南京信息工程大学 一种地球同步轨道圆迹sar水平面二维旁瓣抑制方法
CN108318891A (zh) * 2017-11-28 2018-07-24 西安电子科技大学 一种基于改进sva和cs的sal数据旁瓣的压低方法
CN108318891B (zh) * 2017-11-28 2021-09-10 西安电子科技大学 一种基于改进sva和cs的sal数据旁瓣的压低方法
CN109471097A (zh) * 2018-11-14 2019-03-15 湖南华诺星空电子技术有限公司 一种穿墙雷达信号优化处理方法及装置
CN109471097B (zh) * 2018-11-14 2022-04-22 湖南华诺星空电子技术有限公司 一种穿墙雷达信号优化处理方法及装置
CN110333489A (zh) * 2019-07-24 2019-10-15 北京航空航天大学 采用cnn与rsva组合对sar回波数据旁瓣抑制的处理方法

Also Published As

Publication number Publication date
CN102721965B (zh) 2013-08-21

Similar Documents

Publication Publication Date Title
CN102721965B (zh) 用于合成孔径雷达成像处理的全向旁瓣抑制方法
CN106886021B (zh) 高分辨率星载sar成像质量提升方法
CN107976660B (zh) 弹载多通道雷达超低空目标分析与多径回波建模方法
CN105785327B (zh) 频率分集阵列合成孔径雷达高分辨宽测绘带成像方法
CN113253213B (zh) 高频地波雷达电离层回波空间分布特性测量方法及系统
Tarchi et al. SAR imaging with noise radar
CN101980048B (zh) 基于天线组阵技术的空间碎片地基雷达系统
CN103353595B (zh) 基于阵列内插压缩感知的米波雷达测高方法
RU2682661C1 (ru) Способ активной обзорной моноимпульсной радиолокации с инверсным синтезированием апертуры антенны
CN110967692B (zh) 一种成像方法
CN104181531B (zh) 一种基于相控阵雷达的三维关联成像方法
CN111649803B (zh) 基于垂直线性阵列的三维雷达物位计及其设计方法
CN107741586A (zh) 基于DBF‑TOPS加权的星载Ka InSAR信号处理方法
CN107942295B (zh) 一种前视阵列sar系统的稀疏天线
CN104635230A (zh) 一种用于mimo-sar近场测量成像方位向旁瓣抑制的方法
CN102736073B (zh) 一种通用模式下星载sar距离向模糊度的计算方法
CN107192992B (zh) 雷达校准的方法、校准器、雷达装置及系统
CN105137430A (zh) 一种前视阵列sar的回波稀疏获取及其三维成像方法
Zeng et al. Multi-static MIMO-SAR three dimensional deformation measurement system
CN102565772B (zh) 基于sar子孔径序列图像的海洋动态信息提取方法
CN107607945A (zh) 一种基于空间嵌入映射的扫描雷达前视成像方法
CN111190151A (zh) 扫描模式下多模态小卫星sar的系统参数设计及发射功率优化方法
CN102680956B (zh) 星载合成孔径雷达方位向第一模糊区回波信号能量抑制系统
CN104020465A (zh) 基于八单元小孔径圆阵天线的外辐射源雷达测角方法
CN102495408B (zh) 一种合成孔径雷达点阵目标图像数据的自动寻点方法

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: 20130821

CF01 Termination of patent right due to non-payment of annual fee