CN114252878B - 一种基于逆合成孔径雷达对运动目标进行成像及横向定标的方法 - Google Patents
一种基于逆合成孔径雷达对运动目标进行成像及横向定标的方法 Download PDFInfo
- Publication number
- CN114252878B CN114252878B CN202111558921.6A CN202111558921A CN114252878B CN 114252878 B CN114252878 B CN 114252878B CN 202111558921 A CN202111558921 A CN 202111558921A CN 114252878 B CN114252878 B CN 114252878B
- Authority
- CN
- China
- Prior art keywords
- target
- image
- imaging
- time
- distance
- 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
Links
- 238000003384 imaging method Methods 0.000 title claims abstract description 42
- 238000000034 method Methods 0.000 title claims abstract description 32
- 238000002592 echocardiography Methods 0.000 claims abstract description 10
- 238000004458 analytical method Methods 0.000 claims abstract description 9
- 238000012545 processing Methods 0.000 claims abstract description 9
- 238000010586 diagram Methods 0.000 claims abstract description 6
- 230000009466 transformation Effects 0.000 claims abstract description 6
- 239000011159 matrix material Substances 0.000 claims description 12
- 238000005070 sampling Methods 0.000 claims description 5
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 230000003321 amplification Effects 0.000 claims description 2
- 238000003199 nucleic acid amplification method Methods 0.000 claims description 2
- 238000013178 mathematical model Methods 0.000 claims 1
- 238000004422 calculation algorithm Methods 0.000 abstract description 29
- 230000001133 acceleration Effects 0.000 description 10
- 238000004088 simulation Methods 0.000 description 10
- 208000018989 iridogoniodysgenesis Diseases 0.000 description 6
- 230000005540 biological transmission Effects 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 238000009825 accumulation Methods 0.000 description 1
- 230000001427 coherent effect Effects 0.000 description 1
- 230000006835 compression Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000013519 translation Methods 0.000 description 1
Abstract
一种基于逆合成孔径雷达对运动目标进行成像及横向定标的方法,属于雷达技术领域。该方法包括:通过对几何观测模型图分析,构建雷达与目标之间数学表达式,推算出定标尺度因子;采用时频分析的方法对采集回的雷达信号分别在距离向和方位向上进行处理,得到二维目标成像;取两组相邻的等长回波,对两组回波的同一时刻进行成像,然后对两幅图像进行图像旋转相关处理,检索出最佳角速度;将估计的目标旋转参数代入尺度因子公式中,最终完成对目标的横向定标。本发明将瞬时距离—瞬时多普勒算法与图像旋转相关变换相结合,降低了估计旋转参数的误差,更好地对目标进行横向定标。
Description
技术领域
本发明属于雷达技术领域,具体涉及一种基于逆合成孔径雷达对运动目标进行成像及横向定标的方法。
背景技术
合成孔径雷达和逆合成孔径雷达都可以在任何条件下工作,它们具备优越的遥感,高分辨率成像等能力。与合成孔径雷达相比,逆合成孔径雷达在目标的运动未知的情况下扮演重要角色。传统的逆合成孔径雷达成像采用的是距离—多普勒算法,它将回波数据在距离向上和方位向上分别进行傅里叶变换,可以得到一幅聚焦良好的二维图像。
为了能更好的提取成像目标的尺寸信息,我们需要对成像后的图像进行定标处理。要想对图像正确进行定标,首先需要获知的是该图像的距离和方位向分辨率。由于已知雷达发射信号的参数信息,所以可以直接计算出图像的距离分辨率,而方位向分辨率取决于相参积累时间内目标相对雷达视线的有效转角,所以对图像的定标处理就变成了对其方位向分辨率的求解问题,也就是横向定标问题。但是逆合成孔径成像往往针对的是非合作机动目标,所以它的旋转参数是未知的,因此,要想对图像定标,就需要通过一些方法得到目标的旋转参数。目前对逆合成孔径雷达成像的定标方法有很多种。第一种方法是基于图像质量来估计旋转速度。这类方法虽然稳定性较高,但是需要反复迭代搜索最佳有效转速估计值,计算量较大,效率较低。第二种方法是通过提取选定的突出散射体的相位历史来获得旋转参数,这种方法对信号的旁瓣很敏感。还有一种方法是基于图像旋转相关理论估计旋转速度,这种方法采用快速傅里叶变换完成图像平移和旋转,计算效率高。
发明内容
本发明提供了一种基于逆合成孔径雷达对运动目标进行成像及横向定标的方法。为了避免二维图像模糊以及目标高速性的问题,该方案采用瞬时距离—瞬时多普勒算法来代替传统的距离多普勒算法进行成像,在距离向和方位向上采用时频分析来代替传统的傅里叶变换,并且对于匀速旋转的目标,采用图像旋转相关方法估计旋转参数。同时,对于非匀速旋转的物体,采用最小二乘法拟合的方式求得初始角速度和加速度。该方案对散射体的位置没有要求,而且可以明显减小估计的旋转参数的误差。
基于逆合成孔径雷达对运动目标进行成像及横向定标的方法,如图1所示,主要包括如下步骤:
步骤1:通过对几何观测模型图分析,构建雷达与目标之间数学表达式,推算出定标尺度因子;
步骤2:采用时频分析的方法对采集回的雷达信号在距离方向上进行压缩,经过时间采样之后,得到一维距离像;经过运动补偿,同样采用时频分析的方法对雷达信号在方位向上进行压缩,经过时间采样之后,得到二维目标成像;
步骤3:当目标是匀速运动时,取两组相邻的等长回波,对两组回波的同一时刻进行成像,然后对两幅图像进行图像旋转相关处理,检索出最佳角速度;当目标是匀加速运动时,取两组相邻的等长回波,分别在同一时刻对其进行成像。一共取20个时刻,将每一时刻对应的角速度都描绘在坐标轴上,对其进行最小二乘拟合,拟合出的斜线的斜率即为目标旋转加速度,斜线与纵轴的截距即为目标的角初速度。
步骤4:将估计的目标角速度或角加速度代入到中,得到横向分辨率的尺度因子,也就是图像域中的一个单位长度与目标实际尺寸之间的线性转换关系,最终完成对目标的横向定标。
本发明所述的算法的特点:
(1)本发明应用的算法可以使二维目标成像结果更清晰,并采用平均梯度、图像熵这两个指标验证了该算法的有效性。
(2)本发明所提出方案将瞬时距离—瞬时多普勒算法与图像旋转相关算法相结合,降低了估计旋转参数的误差,更好地对目标进行横向定标。
附图说明
图1:算法流程框图;
图2:几何观测模型图;
图3:瞬时距离—瞬时多普勒算法流程图;
图4:目标散射点模型图;
图5:目标成像结果图:(a)采用距离多普勒算法(b)采用瞬时距离—瞬时多普勒算法;
图6:图像旋转相关算法流程图;
图7:目标在1.98s时二维成像结果图:(a)第一组回波信号(b)第二组回波信号;
图8:相关系数曲线图;
图9:横向定标之后的目标图像;
图10:目标在1.98s时二维成像结果图:(a)第一组回波信号(b)第二组回波信号;
图11:相关系数曲线图;
图12:最小二乘拟合结果图;
图13:横向定标之后的目标图像。
具体实施方式
实施例1:
本发明提供了一种基于逆合成孔径雷达成像及横向定标的方法,算法流程图如图1所示,具体包括以下步骤:
步骤1:在理想情况下,目标相当于放置在转台上,围绕转台中心转动,几何模型如图2所示。假设目标围绕原点O以速度ω0做匀速转动,Q点代表目标上其中的任意一个散射点,rq代表原点O与散射点Q的距离,θ0代表OQ与x轴的夹角,ro表示原点O与雷达的距离,由于ro远远大于成像目标的尺寸,所以根据图2的几何模型我们可以得到散射点Q与雷达的距离为:
其中,xq=rq cosθo、yq=rqsinθo分别为Q点的横纵坐标,tm是慢时间。对应的多普勒频率可以表示为:
其中λ代表雷达的波长,从而对应的横向分辨单元为:
上式中,T表示成像时间,其对应的横向分辨率为:
ISAR图像的距离向与r(tm)对应,对应的距离分辨单元数目为
因此可以求得距离向的分辨率为
综上,我们用ηr、ηa分别代表距离分辨率和横向分辨率的定标尺度因子,所谓尺度因子,也就是图像域中的一个单位长度与目标实际尺寸之间的线性转换关系。
联立上式(3)(5)可以得到
其中
其中Yo表示ro对应的距离单元数目,矩阵S为伸缩矩阵,矩阵R为旋转矩阵。
由此可得,两个不同时刻tm1和tm2成的两幅ISAR图像有如下关系:
整理上式可得:
其中,H是旋转变换矩阵,它代表tm1时刻的图像经过放大、旋转、缩小可以变换成tm2时刻的图像,θd为tm1与tm2时刻之间目标旋转的角度,其中H(θd)可以展开为:
我们将tm1时刻的ISAR图像记为f1(Y1,X1),将tm2时刻的ISAR图像记为f2(Y2,X2),图像f1(Y1,X1)经过旋转矩阵H(θd)变换后的图像记为F1(Y,X),由此可得图像F1(Y,X)与f2(Y2,X2)的相关系数为:
其中
F1(Y,X)=H(ω)f1(Y1,X1) (14)
由此,可以在一定范围内对ω进行搜索,得到使相关系数S(ω)最大的转速值ω作为目标的转速估计值即:
最后,把估计的目标参数代入到公式(4)、(6)中,便可求得尺度因子。
步骤2:优化后的算法流程图如图3所示,在处理雷达接收信号时,采用时频分析在距离维上对信号进行距离压缩,取各次脉冲时频分布同一采样点时刻的切片,把这些切片排列在一起。在经过运动补偿之后,在信号的每个距离单元上对信号进行时频分析,这样就会得到时间—多普勒二维数据,再联合距离维,形成时间—距离—多普勒三维图像仿真,沿着每一瞬时的距离多普勒切片就是该时刻对应的目标的二维向。
步骤3:当目标处于匀速运动时,选取两组相邻的等长回波,对两组回波的同一时刻进行成像,然后对两幅图像进行图像旋转相关处理,检索出最佳角速度;当目标处于匀加速运动时,如图6所示,选取两组等长的雷达回波信号,分别在同一时刻对其进行成像。比如,对于第一组回波信号在t1时刻的成像结果记为t2时刻的成像结果记为以此类推,在tn时刻的成像结果记为对于第二组回波信号我们采取同样的方式,将第二组回波信号在t1时刻的成像结果记为在tn时刻的成像结果记为我们依次把t1,t2......tn时刻的每组图像,即进行图像旋转相关处理,可以得到t1,t2......tn时刻对应的角速度ω1,ω2......ωn。我们把每一时刻对应的角速度都描绘在坐标轴上,对其进行最小二乘拟合之后,可以得到一条斜线,这条斜线的斜率就是估计的旋转加速度,斜线与y轴的截距就是估计的角初速度。有一点需要注意,为了确保估计的准确性,我们通常不选取最初时刻和最末时刻的图像。
步骤4:将估计的旋转参数代入到公式(4)、(6)中,计算得到横向定标的尺度因子,在MATLAB仿真平台上完成对目标的横向定标。
本发明的效果可通过以下仿真实验进一步说明:
(1)仿真条件
本实施例的各实施步骤均在MATLAB2016仿真平台上进行。
(2)仿真内容
本节给出了仿真处理的结果,证明了所提出算法的有效性。
仿真1:为了验证成像算法的有效性,我们绘制了一个目标模型,如图4所示。我们分别采用RD算法和IRID算法对目标进行成像,成像结果如图5所示。根据成像结果,我们发现相对于RD算法,IRID算法成像效果明显更好,杂波更少。这是因为在各个瞬时散射点的多普勒值是固定的,因此不会造成图像的模糊。这种方法反映了ISAR图像随时间的变化规律,极大降低了对多普勒补偿的要求,同时也能解决方位角依赖的发散性散射点造成的模糊。根据图像熵和平均梯度的定义,我们可以发现当图像越清晰,它的图像熵越小,平均梯度越大。针对图5中采用RD算法和IRID算法得到的ISAR图像,我们对它们进行图像熵和平均梯度的计算,结果如表1所示。评价结果验证了IRID算法优于RD算法这一结论。
表1:成像质量比较数据
仿真2:针对匀速旋转的目标进行仿真。将图4中的散射点模型作为目标模型,所有散射点的散射强度相同,假设目标在转盘上以0.073rad/s的角速度进行匀速旋转。下面是一系列参数设置:雷达带宽500MHz,载波频率1GHz,脉冲重复时间0.02s,发射脉冲宽度50us,目标与雷达距离1km。我们选取两组相邻的回波信号进行ISAR成像,每组回波信号有300个脉冲。为了确保实验精度,不选取最初时刻和最末时刻的图像。两组回波信号在1.98s的ISAR图像如图7所示。这里我们采用的成像方法是IRID算法,由于目标在匀速转动,所以我们可以看到两幅图像之间有一个明显的转角。为了估计目标的旋转参数,这里采用图像旋转相关算法。根据图8可以看到,曲线的峰值对应的横坐标为0.06763,也就是说,采用该算法估计的旋转角速度为0.06763rad/s。我们将这个结果带入到公式(4)和公式(6)中,得到定标因子,最终完成对目标的定标,定标结果如图9所示。
仿真3:针对匀加速旋转的目标进行仿真。同样采用图4的目标模型。假设目标以0.0349rad/s的角初速度和0.00698rad/s2的角加速度在转盘上进行匀加速旋转。下面是一系列参数设置:雷达带宽500MHz,载波频率1GHz,脉冲重复时间0.02s,脉冲宽度500us,发射脉冲宽度50us,目标与雷达距离1km。同样选取两组等长的相邻回波数据,每组回波数据包含300个脉冲。在1.98s时,采用IRID算法得到ISAR图像如图10所示,经过图像旋转相关处理之后,如图11所示,可以估计此刻的角速度为0.06392rad/s。该方案共选取了二十个不同的时刻,通过图像旋转相关得到不同时刻对应的角速度,并将每一时刻对应的角速度描绘在坐标轴上,如图12所示,其中横轴代表时间,纵轴代表角速度,并采用最小二乘拟合法对这条灰色的折线进行拟合,得到一条黑色的直线。这条线在y轴上的截距就是旋转的初始角速度,斜率就是旋转加速度,即,该运动目标的初始角速度为0.0304rad/s,旋转加速度为0.0066rad/s2。与原始设定值的误差分别为12.89%和5.44%。图13是对目标定标的结果,图像的畸变得到了有效校正,并且横向定标精度高。
Claims (5)
1.一种基于逆合成孔径雷达对运动目标进行成像及横向定标的方法,其特征在于:
步骤1:通过对几何观测模型图分析,构建雷达与目标之间数学表达式,推算出定标尺度因子;
步骤2:采用时频分析的方法对采集回的雷达信号在距离方向上进行压缩,经过时间采样之后,得到一维距离像;经过运动补偿,同样采用时频分析的方法对雷达信号在方位向上进行压缩,经过时间采样之后,得到二维目标成像;
步骤3:取两组相邻的等长回波,对两组回波的同一时刻进行成像,然后对两幅图像进行图像旋转相关处理,检索出最佳角速度;
步骤4:将估计的目标旋转参数代入尺度因子公式中,最终完成对目标的横向定标。
2.根据权利要求1所述的一种基于逆合成孔径雷达对运动目标进行成像及横向定标的方法,其特征在于,根据雷达与目标之间距离构建数学模型,计算出定标尺度因子,具体过程为:
在理想情况下,目标相当于放置在转台上,围绕转台中心转动;假设目标围绕原点O以速度ω0做匀速转动,Q点代表目标上其中的任意一个散射点,rq代表原点O与散射点Q的距离,θ0代表OQ与x轴的夹角,ro表示原点O与雷达的距离,由于ro远远大于成像目标的尺寸,所以可以得到散射点Q与雷达的距离为:
其中,xq=rqcosθo、yq=rqsinθo分别为Q点的横纵坐标,tm是慢时间;对应的多普勒频率可以表示为:
其中λ代表雷达的波长,从而对应的横向分辨单元为:
上式中,T表示成像时间,其对应的横向分辨率为:
图像的距离向与r(tm)对应,对应的距离分辨单元数目为
因此可以求得距离向的分辨率为
综上,我们用ηr、ηa分别代表距离分辨率和横向分辨率的定标尺度因子,所谓尺度因子,也就是图像域中的一个单位长度与目标实际尺寸之间的线性转换关系。
3.根据权利要求1所述的一种基于逆合成孔径雷达对运动目标进行成像及横向定标的方法,其特征在于,在距离向和方位向上分别采取时频分析的方法来代替傅里叶变换,以此得到目标二维图像。
4.根据权利要求3所述的一种基于逆合成孔径雷达对运动目标进行成像及横向定标的方法,其特征在于,将成像结果进行图像旋转相关变换,检索出最佳角速度,具体过程为:
联立上式(3)(5)可以得到
其中
其中Yo表示ro对应的距离单元数目,矩阵S为伸缩矩阵,矩阵R为旋转矩阵;由此可得,两个不同时刻tm1和tm2成的两幅图像有如下关系:
整理上式可得:
其中,H是旋转变换矩阵,代表tm1时刻的图像经过放大、旋转、缩小可以变换成tm2时刻的图像,θd为tm1与tm2时刻之间目标旋转的角度,其中H(θd)可以展开为:
我们将tm1时刻的图像记为f1(Y1,X1),将tm2时刻的图像记为f2(Y2,X2),图像f1(Y1,X1)经过旋转矩阵H(θd)变换后的图像记为F1(Y,X),由此可得图像F1(Y,X)与f2(Y2,X2)的相关系数为:
其中
F1(Y,X)=H(ω)f1(Y1,X1) (14)
由此,可以在一定范围内对ω进行搜索,得到使相关系数S(ω)最大的转速值ω作为目标的转速估计值即:
。
5.根据权利要求4所述的一种基于逆合成孔径雷达对运动目标进行成像及横向定标的方法,其特征在于,将估计的目标旋转参数代入到公式(4)、(6)中,便可求得尺度因子,最终完成横向定标。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111558921.6A CN114252878B (zh) | 2021-12-20 | 一种基于逆合成孔径雷达对运动目标进行成像及横向定标的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111558921.6A CN114252878B (zh) | 2021-12-20 | 一种基于逆合成孔径雷达对运动目标进行成像及横向定标的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114252878A CN114252878A (zh) | 2022-03-29 |
CN114252878B true CN114252878B (zh) | 2024-07-16 |
Family
ID=
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101498788A (zh) * | 2008-02-01 | 2009-08-05 | 清华大学 | 一种逆合成孔径雷达的目标转角估计和横向定标方法 |
CN108318879A (zh) * | 2017-12-29 | 2018-07-24 | 哈尔滨工业大学 | 基于iaa谱估计技术的isar图像横向定标方法 |
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101498788A (zh) * | 2008-02-01 | 2009-08-05 | 清华大学 | 一种逆合成孔径雷达的目标转角估计和横向定标方法 |
CN108318879A (zh) * | 2017-12-29 | 2018-07-24 | 哈尔滨工业大学 | 基于iaa谱估计技术的isar图像横向定标方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111142105B (zh) | 复杂运动目标isar成像方法 | |
CN109856635B (zh) | 一种csar地面动目标重聚焦成像方法 | |
CN110146889B (zh) | 基于最优回波子区域选择的大转角isar成像方法 | |
CN109738894B (zh) | 一种用于大视场合成孔径雷达的高斜视多角度成像方法 | |
CN110568434B (zh) | 一种多通道匀加速sar动目标二维速度估计方法 | |
CN109407070B (zh) | 一种高轨平台地面运动目标检测方法 | |
CN101620272B (zh) | 一种逆合成孔径雷达的目标转速估计方法 | |
CN103728619B (zh) | 基于变重频技术的机载大斜视条带sar成像方法 | |
Peng et al. | Inverse synthetic aperture radar rotation velocity estimation based on phase slope difference of two prominent scatterers | |
CN109031299B (zh) | 低信噪比条件下基于相位差分的isar平动补偿方法 | |
Zhou et al. | Analysis of wide-angle radar imaging | |
CN111722227B (zh) | 基于近似观测矩阵的聚束sar压缩感知成像方法 | |
CN114966693B (zh) | 基于深度学习的机载舰船目标isar精细化成像方法 | |
CN111681272B (zh) | 一种基于奇异性功率谱的sar图像处理方法 | |
CN111060909B (zh) | 一种机载雷达斜前视超分辨成像方法 | |
CN114325704B (zh) | 一种基于波数谱拼接的合成孔径雷达快速时域成像方法 | |
CN109917383B (zh) | 基于回波降采样积累的低信噪比isar成像方法 | |
CN109188436B (zh) | 适用于任意平台轨迹的高效双基sar回波生成方法 | |
Kusk et al. | SAR focusing of P-band ice sounding data using back-projection | |
CN112415512B (zh) | 基于进退法和黄金分割法的sar运动目标聚焦方法 | |
CN111624602A (zh) | 一种旋转FrFT域构建SAR方位向分辨率成像方法 | |
CN114252878B (zh) | 一种基于逆合成孔径雷达对运动目标进行成像及横向定标的方法 | |
CN117169882A (zh) | 船载雷达海浪信息反演方法 | |
Xingyu et al. | Approach for ISAR imaging of near-field targets based on coordinate conversion and image interpolation | |
CN112305541A (zh) | 一种基于抽样序列长度约束条件下的sar成像方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant |