弹载大斜视小孔径多通道SAR的成像方法与流程

弹载大斜视小孔径多通道sar的成像方法
技术领域
1.本发明属于雷达信号处理技术领域，具体涉及一种弹载大斜视小孔径多通道sar的成像方法。


背景技术：

2.sar(synthetic aperture radar，合成孔径雷达)是一种实用的成像雷达，能全天候、全天时地工作，现已广泛应用于地面场景观测、海洋遥感等诸多方面。随着合成孔径雷达的发展，sar也被应用于许多新的平台，如导弹和战斗机等，由于轨道多样、提前获取图像和实时成像等特点，大斜视sar也被广泛应用在高速机动平台上。
3.在实际分辨率要求不高的情况下，可以利用一小段孔径内的回波信号进行成像以提高实时性，但小孔径sar应用于宽幅场景时，方位向与距离向的模糊将会严重影响成像效果。apc(azimuth phase coding，方位相位编码)可以有效地解决距离模糊的问题，但由于滤波器只能将多普勒谱移动到脉冲重复频率(prf，pulse repetition frequency)之内，所以压缩效果很大程度上取决于重采样率，这也会影响成像处理的实时性。
4.与常规小孔径sar成像相比，弹载大斜视小孔径多通道sar则适用于机动平台的高机动性、宽幅成像的情况。多通道sar成像算法的研究和应用主要集中在星载和直线轨迹机载方面。这种算法对于多通道sar数据的距离模糊，首先通过将各通道信号都补偿到相位中心，之后结合数字波束形成的方法，建立导向矩阵，通过重构多普勒谱完成多通道信号联合解模糊，然后再通过匹配滤波完成聚焦成像。
5.但是，结合机动平台大斜视小孔径的情况下，多通道sar中的通道指向向量是空变的且与平台速度向量的夹角也是空变的，这将导致各通道信号谱的非线性和多普勒谱混叠的问题。


技术实现要素：

6.为了解决现有技术中存在的上述问题，本发明提供了一种弹载大斜视小孔径多通道sar的成像方法。本发明要解决的技术问题通过以下技术方案实现：
7.本发明提供了一种弹载大斜视小孔径多通道sar的成像方法，包括：
8.s1：建立多通道大斜视同心圆斜距模型；
9.s2：获取目标回波信号；
10.s3：对所述目标回波信号进行距离脉冲压缩处理和距离走动校正处理，得到第一校正信号；
11.s4：对所述第一校正信号进行场景中心非空变多通道相位差异项的补偿和空变多通道相位差异项的补偿，得到补偿信号；
12.s5：对所述补偿信号进行多通道信号矩阵重构和基于改进导向矩阵的多普勒谱重构，得到重构信号；
13.s6：对所述重构信号进行距离徙动校正处理和stolt插值差值处理，得到第二校正
信号；
14.s7：对所述第二校正信号进行方位时间谱压缩处理和频域非线性变标处理，得到二维时域信号；
15.s8：对所述二维时域信号进行方位向specan处理，得到聚焦信号；
16.s9：对所述聚焦信号进行几何形变插值校正处理，得到sar雷达成像图。
17.在本发明的一个实施例中，所述多通道大斜视同心圆斜距模型的瞬时斜距包含阵列中心斜距的四阶泰勒级数和多通道差异的二阶泰勒级数，在所述s1中：
18.设雷达平台以速度和恒定加速度沿着曲线移动了t
a
的合成孔径时间，雷达平台上各通道定义为ch
m
，m＝1,2,
…
,m，当雷达位于方位时间t
a
＝0时，中心通道ch1星下点为坐标原点o，中心通道波束指向为定义x轴与的夹角为偏航角θ
yaw
，q点是与p点相差方位角的点，其中
19.雷达平台在曲线上的c点时，定义其方位时间为t
a
，瞬时斜距的四阶麦克劳林级数为：
[0020][0021]
其中，表示雷达中心通道在t
a
时刻的瞬时斜距，t
a
表示方位时间，i表示展开阶数，r表示参考斜距，表示泰勒展开后各次项的系数；
[0022]
参考通道与第m个通道的瞬时斜距差表示为：
[0023][0024]
其中，d
m
表示第m个通道与参考通道的距离，r
d
(d
m
,t
a
,r)表示瞬时斜距差，σ0(d
m
,r)表示泰勒展开后的零次项系数，σ1(d
m
,r)表示泰勒展开后的一次项系数，σ2(d
m
,r)表示泰勒展开后的二次项系数；
[0025]
第m个通道的瞬时斜距表示为：
[0026][0027]
在本发明的一个实施例中，所述s2包括：
[0028]
获取目标回波信号的第一距离时域和方位时域信号a(t
r
,t
a
)，对所述第一距离时域和方位时域信号a(t
r
,t
a
)进行距离向的傅里叶变换，得到第一距离频域和方位时域信号a
′
(f
r
,t
a
)，其中，
[0029]
所述第一距离频域和方位时域信号a
′
(f
r
,t
a
)表示为：
[0030][0031]
其中，t
r
表示距离时间，t
a
表示方位时间，f
r
表示距离频率，w
r
(
·
)表示距离向频域
包络，w
a
(
·
)表示方位向时域包络，c表示光速，f
c
表示载波频率，γ表示调频率。
[0032]
在本发明的一个实施例中，所述s3包括：
[0033]
将所述第一距离频域和方位时域信号a
′
(f
r
,t
a
)分别与距离脉冲压缩函数h
rc
(f
r
,t
a
)和距离走动校正函数h
rwc
(f
r
,t
a
)进行点乘，得到第二距离频域和方位时域信号b(f
r
,t
a
)作为所述第一校正信号，其中，
[0034]
所述距离脉冲压缩函数h
rc
(f
r
,t
a
)和所述距离走动校正函数h
rwc
(f
r
,t
a
)分别表示为：
[0035][0036]
h
rwc
(f
r
,t
a
)＝exp(j4π(f
c
f
r
)k
10
(r
s
)t
a
/c)，
[0037]
其中，r
s
表示参考斜距，k
10
(r
s
)表示在r＝r
s
和时的值，j表示虚数单位。
[0038]
在本发明的一个实施例中，所述s4包括：
[0039]
s41：将所述第二距离频域和方位时域信号b(f
r
,t
a
)与场景中心非空变多通道相位差异函数h1(f
r
,t
a
)相乘，得到第三距离频域和方位时域信号c(f
r
,t
a
)，其中，
[0040]
所述场景中心非空变多通道相位差异函数h1(f
r
,t
a
)表示为：
[0041][0042]
其中，δr
m
(r
s
,t
a
)表示非空变的多通道参考斜距差异；
[0043]
s42：对所述第三距离频域和方位时域信号c(f
r
,t
a
)进行距离向逆傅里叶变换，得到第二距离时域和方位时域信号c
′
(t
r
,t
a
)，其中，
[0044]
所述第二距离时域和方位时域信号c
′
(t
r
,t
a
)表示为：
[0045][0046]
其中，r
′
m
＝r
m
‑
δr
m
(r
s
,t
a
)
‑
k1(r
s
)t
a
；
[0047]
s43：将所述第二距离时域和方位时域信号c
′
(t
r
,t
a
)与空变多通道相位差异函数h2(t
r
,t
a
)相乘，得到第三距离时域和方位时域信号c
″
(t
r
,t
a
)作为所述补偿信号，其中，
[0048]
所述空变多通道相位差异函数h2(t
r
,t
a
)表示为：
[0049][0050]
其中，δr
m
(r,t
a
)＝δr
m
(r,t
a
)
‑
δr
m
(r
s
,t
a
)表示目标非空变的多通道斜距差异，r＝r
‑
r
s
表示目标距离向与参考点的相对位置。
[0051]
在本发明的一个实施例中，所述s5包括：
[0052]
s51：对所述第三距离时域和方位时域信号c
″
(t
r
,t
a
)进行方位向傅里叶变换，得到第一距离时域和方位频域信号d(t
r
,f
a
)；
[0053]
s52：对所述第一距离时域和方位频域信号d(t
r
,f
a
)进行多通道信号矩阵重构处理，得到第二距离时域和方位频域信号d
′
(t
r
,f
b
)，其中，
[0054]
所述第二距离时域和方位频域信号d
′
(t
r
,f
b
)表示为：
[0055][0056]
其中，f
b
∈[
‑
prf/2,prf/2]表示重构后方位向基带频率，i∈[i
min
,i
max
]表示f
a
的模糊数，f
a
表示方位频率，h(r,d
m
,f
b
i
·
prf)表示通道转移函数，prf表示脉冲重复频率；
[0057]
s53：对所述第二距离时域和方位频域信号d
′
(t
r
,f
b
)进行基于改进导向矩阵的多普勒谱重构，得到第三距离时域和方位频域信号d
″
(t
r
,f
al
)作为所述重构信号，其中，
[0058]
所述第三距离时域和方位频域信号d
″
(t
r
,f
al
)表示为：
[0059]
d
″
(t
r
,f
al
)＝w
opt
(r,f
b
)
·
d
′
(t
r
,f
b
)，
[0060]
其中，f
al
表示经过多普勒谱重构后的方位频率，w
opt
(r,f
b
)表示基于数字波束形成的改进加权矩阵；
[0061]
所述基于数字波束形成的改进加权矩阵w
opt
(r,f
b
)表示为：
[0062][0063]
其中，r
‑1表示斜距矩阵的协方差的逆，表示导向矩阵，h表示共轭转置。
[0064]
在本发明的一个实施例中，所述s6包括：
[0065]
s61：对所述第三距离时域和方位频域信号d
″
(t
r
,f
al
)进行距离向傅里叶变换，得到第一距离频域和方位频域信号e(f
r
,f
al
)，其中，
[0066]
所述第一距离频域和方位频域信号e(f
r
,f
al
)表示为：
[0067][0068]
其中，f
al
表示经过多普勒谱重构后的方位频率，φ(r；f
r
,f
al
)表示信号剩余相位；
[0069]
s62：将所述第一距离频域和方位频域信号e(f
r
,f
al
)，与距离非空变的距离徙动校正函数h
rcmc
(f
r
,f
al
)相乘，得到第二距离频域和方位频域信号e
′
(f
r
,f
al
)，其中，
[0070]
所述距离非空变距离徙动校正函数h
rcmc
(f
r
,f
al
)表示为：
[0071]
h
rcmc
(f
r
,f
al
)＝exp(
‑
jφ0(r
s
；f
r
,f
al
))，
[0072]
其中，φ0(r
s
；f
r
,f
al
)表示非空变项；
[0073]
s63：对所述第二距离频域和方位频域信号e
′
(f
r
,f
al
)进行改进的stolt插值，得到第三距离频域和方位频域信号e
″
(f
r
,f
al
)，其中，
[0074]
所述改进的stolt插值操作表示为：
[0075][0076]
其中，f
r
′
表示stolt插值后新的距离频率，φ1(f
r
,f
al
)表示空变项；
[0077]
s64：对所述第三距离频域和方位频域信号e
″
(f
r
,f
al
)进行距离向逆傅里叶变换，得到第四距离时域和方位频域信号e
″′
(t
r
,f
al
)作为所述第二校正信号，其中，
[0078]
所述第四距离时域和方位频域信号e
″′
(t
r
,f
al
)表示为：
[0079][0080]
其中，ψ(r,f
al
)＝φ(r；f
r
＝0,f
al
)
‑
φ0(r
s
；f
r
＝0,f
al
)，表示与方位向聚焦相关的相位，λ表示载波波长。
[0081]
在本发明的一个实施例中，所述s7包括：
[0082]
s71：将所述第四距离时域和方位频域信号e
″′
(t
r
,f
al
)，与时间维谱压缩函数h
tsc
(f
al
)相乘，得到第五距离时域和方位频域信号f(t
r
,f
al
)，其中，
[0083]
所述时间维谱压缩函数h
tsc
(f
al
)表示为：
[0084]
h
tsc
(f
al
)＝exp(jφ0(r
s
；f
r
＝0,f
al
))；
[0085]
s72：将所述第五距离时域和方位频域信号f(t
r
,f
al
)与频域非线性变标函数h
ncs
(f
al
)相乘，得到第六距离时域和方位频域信号f
′
(t
r
,f
al
)，其中，
[0086]
所述频域非线性变标函数h
ncs
(f
al
)表示为：
[0087][0088]
其中，p表示三次项系数，q表示四次项的系数；
[0089]
s73：对所述第六距离时域和方位频域信号f
′
(t
r
,f
al
)进行方位向逆傅里叶变换，得到所述二维时域信号f
″
(t
r
,t
al
)，其中，t
al
表示重构后的方位时间。
[0090]
在本发明的一个实施例中，所述s8包括：
[0091]
s81：将所述二维时域信号f
″
(t
r
,t
al
)与方位向specan函数h
spe
(t
al
)相乘，得到第四距离时域和方位时域信号g(t
r
,t
al
)，其中，
[0092]
所述方位向specan函数h
spe
(t
al
)表示为：
[0093][0094]
其中，h
20
对应的相位，h
30
对应的相位，h
40
对应的相位；
[0095]
s82：对所述第四距离时域和方位时域信号g(t
r
,t
al
)进行方位向傅里叶变换，得到聚焦信号h(t
r
,f
al
)，其中，
[0096]
所述聚焦信号h(t
r
,f
al
)表示为：
[0097]
h(t
r
,f
al
)＝σ0sinc(b
r
(t
r
‑
2k0/c))sinc(t
a
(f
al
‑
h1/2π))，
[0098]
其中，σ0表示地面目标的后向散射系数，h1表示信号方位频域一次项的系数，b
r
表示距离相带宽，t
a
表示合成孔径时间，k0表示目标距离向聚焦位置。
[0099]
与现有技术相比，本发明的有益效果在于：
[0100]
1.本发明的弹载大斜视小孔径多通道sar的成像方法，以机动平台为背景，建立了多通道大斜视同心圆斜距模型，其对阵列中心采用四阶泰勒级数近似瞬时斜距，对多通道差异采用了二阶泰勒级数展开，减小了多通道sar成像过程中的参考斜距的误差，更好的实时适应载有多通道的高机动运动轨迹。
[0101]
2.本发明的弹载大斜视小孔径多通道sar的成像方法，通过距离频域补偿非空变多通道相位差异和距离时域补偿空变多通道相位差异的方法，解决了通道间相位差异带来的无法联合处理的问题。
[0102]
3.本发明的弹载大斜视小孔径多通道sar的成像方法，通过多通道信号矩阵重构
和基于改进导向矩阵的多普勒谱重构，解决了大斜视小孔径机动平台方位带宽不足的问题，完成多通道sar信号的重构，在机动平台大斜视sar成像中有良好的效果及应用价值。
[0103]
上述说明仅是本发明技术方案的概述，为了能够更清楚了解本发明的技术手段，而可依照说明书的内容予以实施，并且为了让本发明的上述和其他目的、特征和优点能够更明显易懂，以下特举较佳实施例，并配合附图，详细说明如下。
附图说明
[0104]
图1是本发明实施例提供的一种弹载大斜视小孔径多通道sar的成像方法的流程框图；
[0105]
图2是本发明实施例提供的一种弹载大斜视小孔径多通道sar的成像方法的流程示意图；
[0106]
图3是本发明实施例提供的一种大斜视多通道同心圆斜距模型的示意图；
[0107]
图4是本发明实施例提供的一种多通道相位差异校正效果的时频线分析示意图；
[0108]
图5是本发明实施例提供的一种方位时域谱压缩时频分析图；
[0109]
图6是本发明实施例提供的一种仿真的目标分布示意图；
[0110]
图7是本发明实施例提供的一种仿真目标距离徙动校正结果分析示意图；
[0111]
图8是本发明实施例提供的一种仿真目标时间域谱压缩过程分析示意图；
[0112]
图9是本发明实施例提供的一种仿真目标的等高线图。
具体实施方式
[0113]
为了进一步阐述本发明为达成预定发明目的所采取的技术手段及功效，以下结合附图及具体实施方式，对依据本发明提出的一种弹载大斜视小孔径多通道sar的成像方法进行详细说明。
[0114]
有关本发明的前述及其他技术内容、特点及功效，在以下配合附图的具体实施方式详细说明中即可清楚地呈现。通过具体实施方式的说明，可对本发明为达成预定目的所采取的技术手段及功效进行更加深入且具体地了解，然而所附附图仅是提供参考与说明之用，并非用来对本发明的技术方案加以限制。
[0115]
实施例一
[0116]
请集合参见图1和图2，图1是本发明实施例提供的一种弹载大斜视小孔径多通道sar的成像方法的流程框图；图2是本发明实施例提供的一种弹载大斜视小孔径多通道sar的成像方法的流程示意图。如图所示，本实施例的成像方法包括：
[0117]
s1：建立多通道大斜视同心圆斜距模型；
[0118]
在本实施例中，多通道大斜视同心圆斜距模型的瞬时斜距包含阵列中心斜距的四阶泰勒级数和多通道差异的二阶泰勒级数。
[0119]
请参见图3，图3是本发明实施例提供的一种大斜视多通道同心圆斜距模型的示意图。如图所示，本实施例的多通道大斜视同心圆斜距模型是一种机动轨迹的时频谱重构的大斜视成像几何模型，具体地，在s1中：
[0120]
设雷达平台在合成孔径时间t
a
内沿着曲线以速度矢量和恒定加速度矢量运动。在雷达平台上配备有m个通道的天线阵列，各通道定义为ch
m
，m＝1,2,
…
,m。天线上每
一个通道都可以收到信号，同时中心通道ch1发射线性调频信号，即ch1也作为信号发射通道(为了简化示意图，图2仅展示三个通道)。
[0121]
通道阵列的指向是时变的，当雷达位于方位时间t
a
＝0时，雷达平台位于位置b，中心通道ch1星下点为坐标原点o，中心通道波束指向为定义x轴与的夹角为雷达系统的偏航角θ
yaw
。p点是偏航线上任意一点，其到零时刻雷达平台位置的距离为r。q点为与p等距离环上的任意一点，故op＝oq，它们之间的方位角也就是，q点是与p点相差方位角的点，其中
[0122]
当雷达平台在方位时间t
a
时到达c点，信号发射通道ch1的瞬时斜距可以展开为四阶麦克劳林级数：
[0123][0124]
其中，表示雷达中心通道在t
a
时刻的瞬时斜距，t
a
表示方位时间，i表示展开阶数，r表示参考斜距，表示泰勒展开后各次项的系数，表示在t
a
＝0处的i阶导数；
[0125]
可见各次项系数是关于参考斜距r和方位角的二维空变参数，所以系数可以在时进行麦克劳林展开：
[0126][0127]
其中，i＝1,2,3,4，j＝0,1,
…4‑
i，k
ij
(r)表示在处的各阶导数。
[0128]
之后，第m个通道的瞬时斜距可以表示为：
[0129][0130]
其中，表示与通道间距同向的单位向量。
[0131]
由于平台轨迹和通道位置可以通过惯导系统获取，故信号发射通道ch1与第m个通道间的斜距差可以表示为方位时间t
a
的二阶多项式：
[0132][0133]
其中，d
m
表示第m个通道与参考通道的距离，r
d
(d
m
,t
a
,r)表示瞬时斜距差，σ0(d
m
,r)表示泰勒展开后的零次项系数，σ1(d
m
,r)表示泰勒展开后的一次项系数，σ2(d
m
,r)表示泰勒展开后的二次项系数；
[0134]
需要说明的是，在本实施例中，信号发射通道ch1作为参考通道。
[0135]
在本实施例中，系数σ
i
(d
m
,r),i＝0,1,2可以通过曲线拟合获得，故第m个通道的两步展开瞬时斜距为：
[0136][0137]
由上式可见，信号发射通道ch1与第m个通道间的斜距是距离耦合且方位空变的。
[0138]
s2：获取目标回波信号；
[0139]
在本实施例中，sar雷达在信号发射通道ch1处发射线性调频信号，并由各通道ch
m
接收目标回波信号，具体地，s2包括：
[0140]
获取目标回波信号的第一距离时域和方位时域信号a(t
r
,t
a
)，它可以表示为：
[0141][0142]
其中，w
r
(*)表示距离时域包络，w
a
(*)表示方位向时域包络，c表示光速，λ表示载波波长，j表示虚数单位，t
m
表示方位时间。
[0143]
对第一距离时域和方位时域信号a(t
r
,t
a
)进行距离向的傅里叶变换，得到第一距离频域和方位时域信号a
′
(f
r
,t
a
)，其中，
[0144]
第一距离频域和方位时域信号a
′
(f
r
,t
a
)表示为：
[0145][0146]
其中，t
r
表示距离时间，t
a
表示方位时间，f
r
表示距离频率，w
r
(
·
)表示距离向频域包络，w
a
(*)表示方位向时域包络，c表示光速，f
c
表示载波频率，γ表示调频率。
[0147]
需要说明的是，本实例中利用传统的重构手段，将斜距转化为：
[0148][0149]
其中，δr
m
(r,t
a
)表示多通道相位差异中需要补偿的距离变化且时变的部分。
[0150]
s3：对目标回波信号进行距离脉冲压缩处理和距离走动校正处理，得到第一校正信号；
[0151]
具体地，s3包括：
[0152]
将第一距离频域和方位时域信号a
′
(f
r
,t
a
)分别与距离脉冲压缩函数h
rc
(f
r
,t
a
)和距离走动校正函数h
rwc
(f
r
,t
a
)进行点乘，得到第二距离频域和方位时域信号b(f
r
,t
a
)作为第一校正信号，其中，
[0153]
距离脉冲压缩函数h
rc
(f
r
,t
a
)和距离走动校正函数h
rwc
(f
r
,t
a
)分别表示为：
[0154][0155]
其中，r
s
表示参考斜距，k
10
(r
s
)表示在r＝r
s
和时的值，j表示虚数单
位。
[0156]
需要说明的是，本实施例中，由于较大的斜视角，产生了很大的多普勒中心，同时多普勒谱宽也超过了方位波束宽度，使一个目标的包络将分布在很多个距离单元之内，因此，本实施例中提出距离走动校正函数来解决大斜视角带来的问题。
[0157]
s4：对第一校正信号进行场景中心非空变多通道相位差异项的补偿和空变多通道相位差异项的补偿，得到补偿信号；
[0158]
具体地，s4包括：
[0159]
s41：将第二距离频域和方位时域信号b(f
r
,t
a
)与场景中心非空变多通道相位差异函数h1(f
r
,t
a
)相乘，得到第三距离频域和方位时域信号c(f
r
,t
a
)，其中，
[0160]
场景中心非空变多通道相位差异函数h1(f
r
,t
a
)表示为：
[0161][0162]
其中，δr
m
(r
s
,t
a
)表示非空变的多通道参考斜距差异，即在r＝r
s
时的δr
m
(r,t
a
)；经过该步骤后，信号中的与距离无关的非空变多通道相位差异被消除。
[0163]
s42：对第三距离频域和方位时域信号c(f
r
,t
a
)进行距离向逆傅里叶变换，得到第二距离时域和方位时域信号c
′
(t
r
,t
a
)，其中，
[0164]
第二距离时域和方位时域信号c
′
(t
r
,t
a
)表示为：
[0165][0166]
其中，r
′
m
＝r
m
‑
δr
m
(r
s
,t
a
)
‑
k1(r
s
)t
a
；
[0167]
s43：将第二距离时域和方位时域信号c
′
(t
r
,t
a
)与空变多通道相位差异函数h2(t
r
,t
a
)相乘，得到第三距离时域和方位时域信号c
″
(t
r
,t
a
)作为补偿信号，其中，
[0168]
空变多通道相位差异函数h2(t
r
,t
a
)表示为：
[0169][0170]
其中，δr
m
(r,t
a
)＝δr
m
(r,t
a
)
‑
δr
m
(r
s
,t
a
)表示目标非空变的多通道斜距差异，r＝r
‑
r
s
表示目标距离向与参考点的相对位置。
[0171]
需要说明的是，由于搭载在机动平台的大斜视多通道小孔径sar具有很短的合成孔径时间，故单个目标的信号在距离走动校正之后会聚集在同一距离单元，故可以用空变多通道相位差异函数h2(t
r
,t
a
)来校正。通过步骤s4，距离变化且时变的多通道相位差异被补偿，故谱重构可以在单个距离单元内完成。
[0172]
s5：对补偿信号进行多通道信号矩阵重构和基于改进导向矩阵的多普勒谱重构，得到重构信号；
[0173]
具体地，s5包括：
[0174]
s51：对第三距离时域和方位时域信号c
″
(t
r
,t
a
)进行方位向傅里叶变换，得到第一距离时域和方位频域信号d(t
r
,f
a
)；
[0175]
具体地，在多通道相位差异补偿后，第m个通道的斜距可以表示为：
[0176][0177]
其中r
ref
(t
a
)表示在距离走动校正后的参考通道斜距，因此第m个通道的多普勒谱可以表示为：
[0178]
d
m
(t
r
,f
a
)＝s
ref
(r,f
a
)h(r,d
m
,f
a
)
ꢀꢀꢀ
(14)，
[0179]
多普勒谱中通道转移函数h(r,d
m
,f
a
)可以表示为：
[0180][0181]
可以看出导向矢量指数相位与方位频率f
a
具有线性关系，因此，在多通道相位差异补偿后，请参见图4，图4是本发明实施例提供的一种多通道相位差异校正效果的时频线分析示意图，其中，(a)图表示多通道相位差异校正前的时频图，(b)图表示多通道相位差异校正后的时频图，如图4中的(b)图所示，时频线曲线变成直线。
[0182]
s52：对第一距离时域和方位频域信号d(t
r
,f
a
)进行多通道信号矩阵重构处理，得到第二距离时域和方位频域信号d
′
(t
r
,f
b
)，其中，
[0183]
第二距离时域和方位频域信号d
′
(t
r
,f
b
)表示为：
[0184][0185]
其中，f
b
∈[
‑
prf/2,prf/2]表示重构后方位向基带频率，i∈[i
min
,i
max
]表示f
a
的模糊数，f
a
表示方位频率，h(r,d
m
,f
b
i
·
prf)表示通道转移函数，prf表示脉冲重复频率；
[0186]
那么，所有通道的多普勒谱相对于方位向基带频率f
b
可以表示为：
[0187][0188]
其中，s(r,f
b
)是m
×
1的信号向量，n(r,f
b
)是m
×
1的噪声向量。
[0189]
s53：对第二距离时域和方位频域信号d
′
(t
r
,f
b
)进行基于改进导向矩阵的多普勒谱重构，得到第三距离时域和方位频域信号d
″
(t
r
,f
al
)作为重构信号，其中，
[0190]
第三距离时域和方位频域信号d
″
(t
r
,f
al
)表示为：
[0191]
d
″
(t
r
,f
al
)＝w
opt
(r,f
b
)
·
d
′
(t
r
,f
b
)
ꢀꢀꢀ
(18)，
[0192]
其中，f
al
表示经过多普勒谱重构后的方位频率，w
opt
(r,f
b
)表示基于数字波束形成的改进加权矩阵；
[0193]
基于数字波束形成的改进加权矩阵w
opt
(r,f
b
)表示为：
[0194][0195]
其中，r
‑1表示斜距矩阵的协方差的逆，表示导向矩阵，h表示共轭转
置。
[0196]
在本实施例中，为了重构模糊的多普勒普，导向矩阵可以构造为：
[0197][0198]
其中a(r,f
b
)是m
×
1的向量，a
i
(r)是导向矢量表示为：
[0199][0200]
同时参考通道地模糊部分可以表示为：
[0201][0202]
将公式(16
‑
17)和(20
‑
22)组合，信号矩阵可以表示为：
[0203]
s(r,f
b
)＝a(r,f
b
)
·
s
ref
(r,f
b
) n(f
b
)
ꢀꢀꢀ
(23)，
[0204]
由于本实施例利用了多普勒谱重构的算法，通过将s
ref
(r,f
b
)的模糊部分进行补偿，获得了等价的单通道不模糊信号，与此同时方位向采样率也提高了n倍，即新的信号的prf
i
＝i
·
prf，其中i＝i
max
‑
i
min
1为多普勒模糊数。
[0205]
s6：对重构信号进行距离徙动校正处理和stolt插值差值处理，得到第二校正信号；
[0206]
具体地，s6包括：
[0207]
s61：对第三距离时域和方位频域信号d
″
(t
r
,f
al
)进行距离向傅里叶变换，得到第一距离频域和方位频域信号e(f
r
,f
al
)，其中，
[0208]
第一距离频域和方位频域信号e(f
r
,f
al
)表示为：
[0209][0210]
其中，f
al
表示经过多普勒谱重构后的方位频率，φ(r；f
r
,f
al
)表示信号剩余相位；
[0211]
另外，信号剩余相位φ(r；f
r
,f
al
)表示为：
[0212][0213]
其中，表示参考方位相位中心。
[0214]
在本实施例中，距离空变的距离徙动可以通过分块距离徙动校正来解决，所以在距离向处理中距离徙动的方位空变可以忽略，仅有线性距离空变需要考虑。因此，本实施例
斜距模型中的系数可以在r＝r
s
处泰勒展开为：
[0215][0216]
s62：将第一距离频域和方位频域信号e(f
r
,f
al
)，与距离非空变的距离徙动校正函数h
rcmc
(f
r
,f
al
)相乘，得到第二距离频域和方位频域信号e
′
(f
r
,f
al
)，其中，
[0217]
在本实施例中，为了解决距离徙动，将信号剩余相位φ(r；f
r
,f
al
)展开为两部分：
[0218]
φ(r；f
r
,f
al
)＝φ0(r
s
；f
r
,f
al
) φ1(f
r
,f
al
)r
ꢀꢀꢀ
(27)，
[0219]
其中，φ0(r
s
；f
r
,f
al
)表示非空变项，φ1(f
r
,f
al
)表示空变项。
[0220]
故在本实施例中，距离非空变距离徙动校正函数h
rcmc
(f
r
,f
al
)表示为：
[0221]
h
rcmc
(f
r
,f
al
)＝exp(
‑
jφ0(r
s
；f
r
,f
al
))
ꢀꢀꢀꢀ
(28)，
[0222]
s63：对第二距离频域和方位频域信号e
′
(f
r
,f
al
)进行改进的stolt插值，得到第三距离频域和方位频域信号e
″
(f
r
,f
al
)，其中，
[0223]
改进的stolt插值操作表示为：
[0224][0225]
其中，f
r
′
表示stolt插值后新的距离频率；
[0226]
s64：对第三距离频域和方位频域信号e
″
(f
r
,f
al
)进行距离向逆傅里叶变换，得到第四距离时域和方位频域信号e
″′
(t
r
,f
al
)作为第二校正信号，其中，
[0227]
第四距离时域和方位频域信号e
″′
(t
r
,f
al
)表示为：
[0228][0229]
其中，ψ(r,f
al
)表示与方位向聚焦相关的相位，λ表示载波波长，具体地，
[0230]
ψ(r,f
al
)＝φ(r；f
r
＝0,f
al
)
‑
φ0(r
s
；f
r
＝0,f
al
)
ꢀꢀꢀ
(31)。
[0231]
需要说明的是，由于目标共享位于相同同心圆上的聚焦在相同的距离单元上的参考倾斜，因此，可以根据每个距离单元更新方位向处理的参数。
[0232]
s7：对第二校正信号进行方位时间谱压缩处理和频域非线性变标处理，得到二维时域信号；
[0233]
具体地，s7包括：
[0234]
s71：将第四距离时域和方位频域信号e
″′
(t
r
,f
al
)，与时间维谱压缩函数h
tsc
(f
al
)相乘，得到第五距离时域和方位频域信号f(t
r
,f
al
)，
[0235]
具体地，对于小合成孔径数据而言，方位向应当利用时间维谱分析(specan)方式，将目标聚焦在方位频域。然而由于公式(28)中参考相位的叠加。请参见图5，图5是本发明实施例提供的一种方位时域谱压缩时频分析图，其中，(a)图表示不同目标的时频线旋转了一个角度α的结果图，(b)图表示时域谱压缩函数的时频图，(c)图表示时域谱压缩结果图，如图5中的(b)图所示，不同目标的时频线旋转了一个角度α，即产生了时间维谱扩张。如果没有补零操作，信号在方位时域将会产生混叠。因此提出时间维谱压缩函数来解决方位时域混叠的问题。
[0236]
其中，时间维谱压缩函数h
tsc
(f
al
)表示为：
[0237]
h
tsc
(f
al
)＝exp(jφ0(r
s
；f
r
＝0,f
al
))
ꢀꢀꢀ
(32)；
[0238]
在本实施例中，信号经过与时间维谱压缩函数相乘后，如图5中的(c)图所示，每个目标的时频谱都旋转了角度α，从而解决了方位时域谱混叠的问题。
[0239]
s72：将第五距离时域和方位频域信号f(t
r
,f
al
)与频域非线性变标函数h
ncs
(f
al
)相乘，得到第六距离时域和方位频域信号f
′
(t
r
,f
al
)。
[0240]
经过方位时域谱压缩后的信号相位可以表示为：
[0241][0242]
其中f
dr
＝2(k1‑
k
10
)/λ，表示方位空变的多普勒中心。为了校正多普勒中心的方位空变，本实施了提出了频域非线性变标函数h
ncs
(f
al
)，表示为：
[0243][0244]
其中，p表示三次项系数，q表示四次项的系数；
[0245]
经过方位向逆傅里叶变换，信号的相位在时域表示为：
[0246][0247]
其中，h1t
ai
表示相位的线性项，代表聚焦的位置，其他项高阶影响聚焦效果。为了获得频域非线性变标函数的具体系数，高阶项中的二次项系数又可以对方位角展开为：
[0248][0249]
当方位时间二次项中的空变项h
21
和h
22
为零时，可以得到三次项系数p和四次项的系数q，具体为：
[0250][0251]
s73：对第六距离时域和方位频域信号f
′
(t
r
,f
al
)进行方位向逆傅里叶变换，得到二维时域信号f
″
(t
r
,t
al
)，其中，t
al
表示重构后的方位时间。
[0252]
s8：对二维时域信号进行方位向specan处理，得到聚焦信号；
[0253]
具体地，s8包括：
[0254]
s81：将二维时域信号f
″
(t
r
,t
al
)与方位向specan函数h
spe
(t
al
)相乘，得到第四距离时域和方位时域信号g(t
r
,t
al
)，其中，
[0255]
方位向specan函数h
spe
(t
al
)表示为：
[0256][0257]
其中，h
20
对应的相位，h
30
对应的相位，h
40
对应的相位；
[0258]
s82：对第四距离时域和方位时域信号g(t
r
,t
al
)进行方位向傅里叶变换，得到聚焦信号h(t
r
,f
al
)，其中，
[0259]
聚焦信号h(t
r
,f
al
)表示为：
[0260]
h(t
r
,f
al
)＝σ0sinc(b
r
(t
r
‑
2k0/c))sinc(t
a
(f
al
‑
h1/2π))
ꢀꢀꢀ
(39)，
[0261]
其中，σ0表示地面目标的后向散射系数，h1表示信号方位频域一次项的系数，b
r
表示距离相带宽，t
a
表示合成孔径时间，k0表示目标距离向聚焦位置。
[0262]
s9：对聚焦信号进行几何形变插值校正处理，得到sar雷达成像图。
[0263]
本实施例的弹载大斜视小孔径多通道sar的成像方法，以机动平台为背景，建立了多通道大斜视同心圆斜距模型，其对阵列中心采用四阶泰勒级数近似瞬时斜距，对多通道差异采用了二阶泰勒级数展开，减小了多通道sar成像过程中的参考斜距的误差，更好的实时适应载有多通道的高机动运动轨迹。
[0264]
而且，通过距离频域补偿非空变多通道相位差异和距离时域补偿空变多通道相位差异的方法，解决了通道间相位差异带来的无法联合处理的问题。通过多通道信号矩阵重构和基于改进导向矩阵的多普勒谱重构，解决了大斜视小孔径机动平台方位带宽不足的问题，完成多通道sar信号的重构，在机动平台大斜视sar成像中有良好的效果及应用价值。
[0265]
实施例二
[0266]
本实施例对实施例一的弹载大斜视小孔径多通道sar的成像方法进行了仿真实验验证。
[0267]
1、仿真条件
[0268]
本实验仿真参数如表1所示，距离向和方位向宽度分别为3km和6km，点目标的分布如图6所示，图6是本发明实施例提供的一种仿真的目标分布示意图。
[0269]
表1 sar仿真参数
[0270][0271][0272]
2、仿真内容和结果分析
[0273]
仿真1：
[0274]
请参见图7，图7是本发明实施例提供的一种仿真目标距离徙动校正结果分析示意图。图7中左边的图表示常规方式的距离徙动校正的结果，右边的图为本实施例的距离徙动的距离徙动校正的结果，图7中(a)为目标a1的距离徙动校正结果，图7中(b)为目标b2的距离徙动校正结果，图7中(c)为目标c3的距离徙动校正结果。从图中可以看出，右边的结果图普遍比左边常规方式的结果图，在距离向聚焦效果更好，因此，证明了本实施例改进的omega
‑
k方式进行距离徙动校正的有效性。
[0275]
请参见图8，图8是本发明实施例提供的一种仿真目标时间域谱压缩过程分析示意图，图8中的(a)图为时间维谱压缩前的时频线，可以看出，目标在方位时域产生了混叠，本应分布在两边的目标点由于方位时间谱不足折叠到了中间。图8中的(b)图为时间维谱压缩后的时频线，可以看出，目标在方位时域的混叠被消除，各点的信号依次排布在是频谱内。
因此，证明了本实施例时间维谱压缩解决方位时域混叠的有效性。
[0276]
请参见图9，图9是本发明实施例提供的一种仿真目标的等高线图，图9展示了位于不同方位不同距离的九个点的二维剖面图。图9中(a)图从左到右依次表示位于同一距离环上的三个点a1，b1和c1的仿真结果剖面图；图9中(b)图从左到右依次表示位于场景中心距离环上的三个点a2，b2和c2的仿真结果剖面图；图9中(c)图从左到右依次表示位于远距离环上的三个点a3，b3和c3的仿真结果剖面图。
[0277]
从图中可以看出，本实施例所提出的方法获得的成像点目标质量很好，为了进一步评估所提出的聚焦算法的性能，本实施例将目标的测量参数列于表2中，峰值旁瓣比(pslr，peak sidelobe ratio)和积分旁瓣比(islr，integrated sidelobe ratio)的无窗口理想值分别为
‑
13.26db和
‑
10.8db，因此，本实施例所提出的方法的聚焦性能接近理想值。
[0278]
表2 挑选出的点目标聚焦性能评估
[0279][0280]
应当说明的是，在本文中，诸如第一和第二等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来，而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。而且，术语“包括”、“包含”或者任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的物品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素。在没有更多限制的情况下，由语句“包括一个
……”
限定的要素，并不排除在包括所述要素的物品或者设备中还存在另外的相同要素。
[0281]
以上内容是结合具体的优选实施方式对本发明所作的进一步详细说明，不能认定本发明的具体实施只局限于这些说明。对于本发明所属技术领域的普通技术人员来说，在不脱离本发明构思的前提下，还可以做出若干简单推演或替换，都应当视为属于本发明的保护范围。