专利名称：一种用于下视阵列天线合成孔径雷达的三维聚焦成像方法
技术领域：
本发明涉及合成孔 径雷达的信息获取与处理技术领域，尤其是一种用于下视阵列 天线合成孔径雷达的三维聚焦成像方法，能够获取载体正下方成像区域的三维雷达图像。
背景技术：
下视阵列天线合成孔径雷达是一种新的合成孔径雷达三维成像技术，如图1所 示。其中，Oxyz为下视阵列天线合成孔径雷达及其成像区域所在的空间直角坐标系，PiUi, Yi, Zi)为成像区域第i个点目标及其直角坐标。Qxr θ为下视阵列天线合成孔径雷达及其 成像区域所在的柱面坐标系(方位向、斜距向和辐角向)，如图2所示，PiO^ri, Θ,)为成 像区域第i个点目标及其柱面坐标。Lx、Wy分别为下视阵列天线合成孔径雷达方位向合成 孔径长度和地距向成像幅宽。通过在地距向(即方位向-高程向平面的法线方向)上分布 线性阵列天线，使得在地距向上形成了异于方位向合成孔径的第二个真实孔径地距向真 实孔径，L为线阵长度，也即地距向真实孔径长度，从而形成了位于水平面(即方位向-地 距向平面)上的一个二维孔径，通过发射带宽信号，实现了对成像区域在方位向、地距向和 高程向上的三维分辨成像。1999年以来，国际上在下视阵列天线合成孔径雷达三维成像理论和方法方面开展 了一些研究工作。1999 年，德国宇航院(DLR)的 Christoph H. Gierull (C. Gierull. On a Concept for an Airborne Downward-Looking Imaging Radar. ΑΕ , No. 6, pp295_304, 1999.)分析了目前机载单通道合成孔径雷达的成像观测方式主要有正侧视和斜视，其共 同点是在同一时刻只能对载机飞行路线一侧的场景目标进行成像，这是由于单通道合成孔 径雷达存在着由飞行两侧距离对称而引起的左/右模糊问题，无法区分目标来自飞行航迹 的左侧还是右侧，因而单通道合成孔径雷达不具备对载机正前方和正下方区域进行高分 辨率成像的能力，如图3所示。为解决机载合成孔径雷达在载机正下方存在的成像盲区， Gierull将线性阵列天线引入到合成孔径雷达成像中，提出了下视成像雷达的概念，并给出 了下视成像系统的几何模型，将合成孔径和真实孔径综合在一起，方位向采用合成孔径原 理，地距向上引入线性阵列天线构成一个真实孔径，组成了一个二维孔径。由于Gierull提 出的下视成像雷达系统主要用于说明下视阵列天线合成孔径雷达成像的基本原理，因而采 用单频脉冲信号而不是带宽信号，实现了二维成像。2006 年，德国 FGAN-FHR 的 J. Klare 等人(J. Klare, A. Brenner, J. Ender. A new Airborne Radar for 3D Imaging-Image Formation using the ARTINO Principle-. EUSAR2006, Dresden, Germany, 2006 ；M. Wei β , J.Ender, 0. Peters, et al. An airborne Radar for Three Dimensional Imaging and Observation-technical realisation and status of ARTIN0. EUSAR2006，Dresden，Germany，2006.)提出了构建采用天底观 测的机载三维成像雷达(Airborne Radar for Three-dimensional Imaging and Nadir Observation, ARTIN0)系统。ARTINO原理基于沿机翼分布的一组稀疏天线阵，将发射天线 阵分成两部分，分别置于在机翼两端，将接收天线阵置于在机翼中间位置。通过距离向脉冲压缩技术、方位向合成孔径原理及地距向波束形成运算形成最终的三维图像。ARTINO系统 可直接生成真实的三维分辨单元而不是干涉图像，这使得通过旋转三维图像可以从不同的 角度来观察场景存在着可行性。通过在平行的飞行路线获得雷达数据，可以将在地距向上 的几幅图像连在一起，保证了可以得到各个方向上的大面积三维图像。FGAN-FHR将ARTINO 系统安装在小型UAV上，给出技术实现方案。2007年，J. Klare等人(Jens Klare，Delphine Cerutti-Maori, Andreas Brenner, et al. Image quality analysis of the vibrating sparse MIMO antenna array of the airborne 3D imaging radar ARTINO. IGARSS2007, Barcelona, Spain, 2007.)考虑了 ARTINO系统的载机机翼振动对三维聚焦成像的影响，并 给出 了补偿方案。2008 年，J. Klare 等人(Jens Klare,Andreas Brenner, Joachim Ender. Impact of Platform Attitude Disturbances on the 3D Imaging Quality of the UAV ARTIN0. EUSAR2008, Friedrichshafen, Germany, 2008.)进一步考虑了载机姿态的变化(包 括偏航、俯仰和横滚)对三维聚焦成像的影响，并考虑了补偿方案。上述已开展的研究中，主要存在三个问题首先，ARTINO系统所采用的成像 算法 中没有考虑距离徙动校正问题，由于下视阵列天线合成孔径雷达中存在着由方位向合成孔 径和地距向真实孔径组成的一个二维孔径，使得斜距向快时间t与方位向快时间u以及地 距向快时间υ均存在耦合，使得基于匹配滤波的成像算法计算更为复杂；其次，成像算法 中地距向上的压缩采用波束形成来实现，当天线阵元数目较多时，运算效率较低；最后，地 距向采样间隔，也即天线阵元间距d要求不大于λ。/2，使得线阵中包含的天线阵元数目大 为增加，不利于实际应用。
发明内容
本发明的目的在于提供一种用于下视阵列天线合成孔径雷达的三维聚焦成像方 法，该方法考虑了二维距离徙动校正，提高地距向压缩效率同时放宽地距向采样间隔的要 求，解决了现有的三维聚焦成像方法没有进行距离徙动校正的问题和地距向压缩运算效率 较低，下视阵列天线合成孔径雷达地距向采样间隔要求过小等不利于实际应用的问题，为了实现所述目的，本发明的技术解决方案是一种用于下视阵列天线合成孔径雷达的三维聚焦成像方法，包括步骤如下步骤Sl 下视阵列天线合成孔径雷达系统采集到的机下成像区域的原始回波信 号为S(l(t，u，u)，斜距向压缩单元将原始回波信号沿斜距向进行傅立叶变换，生成第1信号 为S1 (ft，u, υ),采用斜距向匹配滤波器H1 (ft)对第1信号进行斜距向聚焦成像处理，生成 第2信号为S2(ft，u，υ )，然后对第2信号进行斜距向傅里叶逆变换，生成第3信号为S3 (t， u, u)，其中，斜距向快时间域为t，斜距向频率域为ft，方位向时间域为u和地距向时间域 为u ；步骤S2 方位向距离徙动校正单元将生成第3信号沿方位向进行傅里叶变换，生 成第4信号为S4 (t，fu，u )，其中，方位向频率域为fu，然后对第4信号沿斜距向进行傅里叶 变换，生成第5信号为S5(ft，fu，u)，采用距离徙动校正因子H2 (ft)对第5信号进行方位向 距离徙动校正，生成第6信号为S6(ft，fu，u)，然后对第6信号进行斜距向傅里叶逆变换， 生成第7信号为S7(t，fu，υ)；步骤S3 方位向压缩单元采用方位向匹配滤波器H3(ft)对第7信号进行方位向聚焦成像处理，生成第8信号为S8 (t,fu, υ),然后对第8信号进行方位向傅里叶逆变换，生成 第9信号为S9(t，U，υ)；步骤S4 地距向距离徙动校正单元将第9信号沿斜距向进行傅里叶变换，生成第 10信号Sltl (ft，u, U)，对第10信号进行地距向距离徙动校正，生成第11信号为S11 (ft，u, u)，对第11信号进行斜距向傅里叶逆变换，生成第12信号为S12(t，u，υ)；步骤S5 地距向压缩单元首先采用解斜参考函数H6(U)对第12信号进行解斜处 理，生成第13信号为S13 (t，u, u)，然后对第13信号进行地距向傅里叶变换，生成第14信 号为S14(u，r，θ )，该信号即为获得的成像区域在柱面坐标系中包含幅度和相位信息的下视 阵列天线合成孔径雷达三维聚焦图像，其中，成像区域中第i个目标在柱面坐标系中方位 向坐标值为Xi，斜距向坐标值为^，辐角向坐标值为θ ρ所述的用于下视阵列天线合成孔径雷达的三维聚焦成像方法，其在成像几何中， 下视阵列天线合成孔径雷达系统中阵列天线阵元间隔即系统沿地距向的采样间隔d小于 等于A。/[2sin(i3/2)]，其中λ。为雷达发射信号的载 波波长，β为发射天线的地距向波束宽度。所述的用于下视阵列天线合成孔径雷达的三维聚焦成像方法，其所述采样间隔d 小于等于队，其中De为发射天线在地距向上的尺寸。所述的用于下视阵列天线合成孔径雷达的三维聚焦成像方法，其所述采样间隔d 小于等于\^/4五2+</(2%)，其中H为系统平台飞行高度，Wy为系统在地距向上的成像幅
觅ο所述的用于下视阵列天线合成孔径雷达的三维聚焦成像方法，其所述下视阵列天 线合成孔径雷达系统中，阵列天线采用单发多收、收发分置的线阵结构，其中单个发射天线 位于线阵的中央，接收天线阵元沿线阵等间隔分布。所述的用于下视阵列天线合成孔径雷达的三维聚焦成像方法，其在成像几何中， 地距向的真实孔径中心位于下视阵列天线合成孔径雷达成像区域的地距向中心。所述的用于下视阵列天线合成孔径雷达的三维聚焦成像方法，其在成像几何中， 成像区域的地距向成像幅宽大于下视阵列天线合成孔径雷达系统地距向线阵真实孔径长度。所述的用于下视阵列天线合成孔径雷达的三维聚焦成像方法，其所述地距向距离 徙动校正单元处理步骤是步骤S41 将第9信号沿斜距向进行傅里叶变换，生成第10信号Sltl (ft，U，υ)；步骤S42 采用距离徙动校正因子扎(&)对第10信号进行地距向第一次距离徙动 校正，生成第101信号为Sltll (ft，u，U)，将第10信号中与斜距位置有关的距离徙动项移除；步骤S43 采用距离徙动校正因子H5(ft)对第101信号进行地距向第二次距离徙 动校正，生成第11信号为S11 (ft，u, U)，将第101信号中与辐角位置有关的距离徙动项移 除；步骤S44 对第11信号进行斜距向傅里叶逆变换，生成第12信号为S12 (t，u，υ )。本发明的有益效果本发明为一种用于下视阵列天线合成孔径雷达的三维聚焦 成像方法，该方法将采集到的下视阵列天线合成孔径雷达原始回波数据变换到斜距向频率 域，进行斜距向匹配滤波；然后将斜距向匹配滤波后的回波信号变换到方位向频率域，选择斜距参考平面进行方位向距离徙动校正，将校正后的信号进行方位向匹配滤波以及方位向 傅里叶逆变换；将斜距向和方位向匹配滤波后的信号进行地距向距离徙动校正，然后进行 斜距向傅里叶逆变换，得到每个通道所成的二维雷达图像；将所有通道所得到的二维图像 进行地距向解斜和地距向傅里叶变换，重建出成像区域在柱面坐标系中的三维雷达图像； 通过坐标变换，得到成像区域在空间直角坐标系中的三维雷达图像。本发明考虑了已有的三维成像方法中未考虑的二维距离徙动校正问题，增强了本 方法适用的普遍性。本发明将已有的三维成像方法中的地距向波束形成运算变换为一次复乘和傅里 叶变换即获取三维分辨图像，缩短成像所需的时间。本发明的适用条件为下视阵列天线合成孔径雷达沿地距向采样间隔，即阵列天线 阵元间距d小于等于X。/[2sin(i3/2)]，从公式中可以看出本发明的地距向采样间隔要大 于现有技术的地距向采样间隔λ。/2，从而减少了下视阵列天线合成孔径雷达中阵列天线 所需要的天线阵元数目，有利于实际应用。
图1是本发明采用的下视阵列天线合成孔径雷达三维成像几何示意图； 图2是本发明采用的下视阵列天线合成孔径雷达成像几何中地距向_高程向剖面 图；图3是机载合成孔径雷达成像区域示意图；图4是本发明采用的下视阵列天线合成孔径雷达三维聚焦成像处理流程图；图5是本发明采用的地距向距离徙动校正单元处理流程图；图6是本发明实施例的仿真输入7个点目标的空间位置分布图；图7是本发明实施例的包含7个点目标的下视阵列天线合成孔径雷达三维成像聚 焦图像；图8是本发明实施例的图7在辐角向θ = 0°剖面上的5个点目标分布图；图9是本发明实施例的图7在方位向U= IOm剖面上的5个点目标分布图；图10是本发明实施例的图7在斜距向r = 490m剖面上的5个点目标分布图。
具体实施例方式下面结合附图详细说明本发明技术方案中所涉及的各个细节问题。应指出的是， 所描述的实施例仅旨在便于对本发明的理解，而对其不起任何限定作用。如图1所示，Oxyz为下视阵列天线合成孔径雷达及其成像区域所在的空间直角坐 标系，成像区域位于下视阵列天线合成孔径雷达飞行航线的正下方区域，并且关于飞行航 迹，即χ轴方向对称，PiO^yi, Zi)为成像区域第i个点目标及其直角坐标。Qxre为下视 阵列天线合成孔径雷达及其成像区域所在的柱面坐标系，垂直于方位向，即χ轴方向的地 距向-高程向平面如图2所示，PiO^ri, Θ》为成像区域第i个点目标及其柱面坐标。点 SfePiO^ri, Θ,)的空间直角坐标与柱面坐标之间的转换关系为 式中，H为下视阵列天线合成孔径雷达载机飞行高度。Lx为下视阵列天线合成孔 径雷达方位向合成孔径长度，L为下视阵列天线合成孔径雷达地距向真实孔径长度，即阵列 天线长度，Wy为下视阵列天线合成孔径雷达地距向成像幅宽，其中地距向成像幅宽Wy大于 下视阵列天线合成孔径雷达系统地距向真实孔径长度L。在下视阵列天线合成孔径雷达系统中，阵列天线采用单发多收、收发分置的线阵 结构，其中单个发射天线位于线阵的中央，也即地距向的真实孔径中心Q处，位于下视阵列 天线合成孔径雷达成像区域的地距向中心0附近。接收天线阵元沿线阵等间隔分布，\为 下视阵列天线合成孔径雷达中第η个接收阵列天线阵元所对应的地距向坐标，其中η = 1， 2，…，N，d为接收阵列天线阵元间距。该下视阵列天线合成孔径雷达成像几何中，通过下视 阵列天线合成孔径雷达的方位向(χ轴)运动形成了平行于航迹向(方位向)的方位向合 成孔径；通过沿地距向(y轴)上分布线性阵列天线形成地距向真实孔径，进而形成的二维 合成孔径平面位于水平面(即方位向-地距向平面)。发射天线发射的线性调频信号P (t) 为 式中，j表示单位虚数，f。表示雷达发射信号的载频，Y表示雷达发射线性调频 信号的调频率，T为脉冲持续时间，rect(.)为矩形窗函数。则线性调频信号的带宽B = Y · T0对于观测场景内后向散射系数为Si的目标PiO^yi, Zi)，其接收到的理想回波信 号为 式中，代^，力= < ㈠+ <卜⑷为观测场景内目标PiOci, Zi)与下视阵列天 线合成孔径雷达之间的往返距离，其中 为目标PiO^yi, Zi)与发射天线之间的距离。 为目标PiO^yi, Zi)与接收天线之间的距离。根据泰勒公式，Ri(u，υ)可近似为 如图4本发明提出的用于下视阵列天线合成孔径雷达的三维聚焦成像方法的流 程图所示，对于下视阵列天线合成孔径雷达采集的数据，其中，输入为下视阵列天线合成 孔径雷达原始回波，经过三维成像聚焦处理后的可获得下视阵列天线合成孔径雷达三维图 像，具体实施步骤如下步骤Sl 下视阵列天线合成孔径雷达系统采集到的机下成像区域的原始回波信 号为s^t，U, U)，斜距向压缩单元将原始回波信号沿斜距向进行傅立叶变换，生成第1信 号为Sjft，!!，u)，表达式为 生成第2信号为S2 (ft，u，u)，表达式为: 然后对第2信号进行斜距向傅里叶逆变换，生成第3信号为S3 (t，u, υ )，表达式 为 其中，斜距向快时间域为t，斜距向频率域为ft，方位向时间域为U和地距向时间 域为u ；步骤S2 方位向距离徙动校正单元将生成第3信号沿方位向进行傅里叶变换，生 成第4信号为S4(t，fu，u)，表达式为:
f 1 f2)^4 ( , / , = <5, · rect. exp jtt^- . exp (-^2τγ/Λ) · (11)
/^y [jrY
\ 1 X XJ\1 X /
Rl TT V2 - 2vr sin θ.) 4π 、sine j B t----- ‘ exp -j----- ■ exp —j — r
CΛ v.Λ
乂)Vc / \ Cj其中，方位向频率域为fu，方位向调频率Yx = 2/U。ri)，方位频率域中的距离徙 动项R' i(fu，u)表达式为:
然后对第4信号沿斜距向进行傅里叶变换，生成第5信号为S5 (ft，fu，υ )，表达式 为 采用距离徙动校正因子H2(ft)对第5信号进行方位向距离徙动校正，H2 (ft)表达 式为 生成第6信号为S6 (ft，fu，υ )，表达式为 然后对第6信号进行斜距向傅里叶逆变换，生成第7信号为S7(t，fu，u)，表达式 为 sine\B t - 2r H---ν sin Θ. C \ · exp -j--1-L
' 2r' /Λr.步骤S3 方位向压缩单元采用方位向匹配滤波器H3(ft)对第7信号进行方位向聚
焦成像处理，H3 (ft)表达式为
f ) f)H11 (f) = rect · exp -jtt^-(17)
3 v i!7 L7
\ X X /\,生成第8信号为S8 (t，fu，υ )，表达式为
Λ 471、S8 {t, fu， 二 Si . rect -γ- . exp(-;2τγ/Λ)· exp‘ (18)
,1X X JVc ν
ν2) i τγ ν2 — 2vr sin θsine \B t - 2r. + ---ν sin θ. C · exp -j--、-L
1 2r.1 /Λτ
L \ι//JJVc1)然后对第8信号进行方位向傅里叶逆变换，生成第9信号为S9 (t，u, u)，表达式 为
式中，Bx为方位向多普勒带宽，表达式为Bx= YxLx0步骤S4 地距向距离徙动校正单元步骤如图5所示，具体实施过程为步骤S41 将第9信号沿斜距向进行傅里叶变换，生成第10信号Sltl (ft，u，υ )，表 达式为 步骤S42 采用距离徙动校正因子扎(&)对第10信号进行地距向第一次距离徙动
校正，H4(ft)表达式为 生成第101信号为Sltll (ft，u，u)，表达式为: 将第10信号中与斜距位置有关的距离徙动项移除；步骤S43 采用距离徙动校正因子H5(ft)对第101信号进行地距向第二次距离徙 动校正，H5(ft)表达式为 生成第11信号为S11 (ft，u，u)，表达式为: 将第101信号中与辐角位置有关的距离徙动项移除；步骤S44 对第11信号进行斜距向傅里叶逆变换，生成第12信号为S12 (t，u，υ), 表达式为
步骤S5 地距向压缩单元首先采用解斜参考函数H6(U)对第12信号进行解斜处
S, Η6(υ)表达式为 生成第13信号为S13(t，u，U)，表达式为 然后对第13信号进行地距向傅里叶变换，生成第14信号为S14 (u，r, θ )，表达式 为 式中，B, = 2B/C。该信号即为获得的成像区域在柱面坐标系(方位向、斜距向和 辐角向)中包含幅度和相位信息的下视阵列天线合成孔径雷达三维聚焦图像，其中，成像 区域中第i个目标在柱面坐标系中方位向坐标值为Xi，斜距向坐标值为ri;辐角向坐标值为
Qi;下视阵列天线合成孔径雷达中阵列天线的阵元间距d由系统在地距向上空间采 样定理决定，阵元之间相距不能太近，否则阵元天线之间的互藕效应会严重影响系统的方 向性能，但也不能相距太远，否则系统在辐角向上会出现角度模糊问题。由式(28)可以看出，第14 信号 S14(x，r，θ)中的 sin [ ji Nd (sin θ-sin θ》/ λ J / sin [Jid(Sine-Sinei)AJ 是关于 sin θ 的周期函数，当满足 π d(sin θ main-sin θ》/ 入。-^^。^^--^!^》/^。= Ill JI时，除在辐角θ main出现主瓣峰值夕卜，还会在eamb出 现栅瓣，式中J G 0表示栅瓣的级数。因而可计算出目标辐角与模糊角的关系为 下视阵列天线合成孔径雷达辐角向成像范围由发射天线地距向波束宽度β决 定，为[-β/2，β/2]，如果θ amb落入到该范围内，则最终的成像结果会在eamb位置处出现栅瓣，此时系统辐角向存在角度模糊问题。为避免辐角模糊问题，式(29)应满足 由于下视阵列天线合成孔径雷达俯辐角向成像范围[_β/2，β/2]在[_π/2， η/2]之内，该范围内辐角θ的正弦函数sine是单调变化的，因此，如果第一级栅瓣(即
1 = 1)角度位置落在辐角向成像范围[-β/2，β/2]之外，则级数更高的栅瓣(即|1| > 1)角度位置也会落在辐角向成像范围之外。因此，当地距向空间采样间隔d满足下面条件 时，系统在辐角向上不会出现角度模糊问题。 由于发射天线地距向波束宽度β e
，满足I 2 sinf，并且β与发射
Δι/L天线地距向尺寸De之间的关系为β = A。/De，故式(31)中d需要满足的约束条 件可另写为 由图1所示，由如下几何关系:R ^ ^JH2 + W2 ，Sinf =孟，代入式(31)可得
地距向空间采用间隔d与跨航向成像幅宽Wy以及载体飞行高度H之间的关系为 以位于三维空间中的点目标为例，雷达载波波长为λ。= 8mm，线性调频信号带宽 B = 375MHz，载体飞行高度H = 500m,阵列天线实孔径长度L = 6m,发射天线地距向波束宽 度β = 6°，由式(31)可以计算出阵列天线阵元间距d = 7. 64cm，则阵列天线中由1个发 射天线和28个接收天线组成。成像区域共有7个点目标，其空间位置如图6所示，其后向 散射系数均为δ =1。则通过本发明的三维聚焦成像方法得到的7个点目标的成像结果 如图7-10所示，其中图7是包含7个点目标的下视阵列天线合成孔径雷达三维成像归一化 后-3dB轮廓显示结果，方位向和斜距向的坐标单位均为m，辐角向的坐标单位为。，坐标轴 上的数字表示相对下视阵列天线合成孔径雷达成像几何中的柱面坐标位置，通过与图6对 比可以看到，通过上述处理步骤后，能够精确地重建出7个点目标的空间位置和幅度信息。 图8是图7在辐角向θ =0°剖面(方位向-斜距向平面)上的5个点目标分布图，由于仿 真中各点目标的幅度一致，因此动态显示图像也一致，同时，各点目标的位置也跟实际目标 的位置保持一致。图9是图7在斜距向U= IOm剖面(斜距向-辐角向平面)上的5个点 目标分布图，由于仿真中各点目标的幅度一致，因此动态显示图像也一致，同时，各点目标的位置也跟实际目标的位置保持一致。图10是图7在斜距向r = 490m剖面(方位向-辐 角向平面)上的5个点目标分布图，由于仿真中各点目标的幅度一致，因此动态显示图像也 一致，同时，各点目标的位置也跟实际目标的位置保持一致。
本发明考虑了已有的三维成像方法中未考虑的二维距离徙动校正问题，增强了本 方法适用的普遍性。本发明将已有的三维成像方法中的地距向波束形成运算变换为一次复乘和傅里 叶变换即获取三维分辨图像，缩短成像所需的时间。本发明使得下视阵列天线合成孔径雷达沿地距向采样间隔，即阵列天线阵元间距 d小于等于A。/[2sin(i3/2)]也能适用，从公式中可以看出本发明的地距向采样间隔要大 于现有技术的地距向采样间隔λ。/2，从而减少了下视阵列天线合成孔径雷达中阵列天线 所需要的天线阵元数目，有利于实际应用。本发明上述的方法，已经在计算机上应用MATLAB软件得到验证，图7_图10是在 MATLAB软件环境下使用本发明所述方法获得的结果。该方法在计算机或专用设备上主要通 过5个程序模块来实现，如图4所示，所述斜距向压缩单元、方位向距离徙动校正单元、方位 向压缩单元、地距向距离徙动校正单元和地距向压缩单元，分别完成步骤Si、步骤S2、步骤 S3、步骤S4和步骤S5的功能。斜距向压缩单元，输入为下视阵列天线合成孔径雷达原始回波数据，输出为处于 斜距向时间域、方位向时间域和地距向时间域的第3信号；方位向距离徙动校正单元，输入为第3信号，输出为处于斜距向时间域、方位向频 率域和地距向时间域的第7信号；方位向压缩单元，输入为第7信号，输出为处于斜距向时间域、方位向时间域和地 距向时间域的第9信号；地距向距离徙动校正单元，输入为第9信号，输出为处于斜距向时间域、方位向时 间域和地距向时间域的第12信号；地距向压缩单元，输入为第12信号，输出为柱面坐标系中聚焦在方位向、斜距向 和辐角向的包含幅度和相位信息的下视阵列天线合成孔径雷达三维图像。以上所述，仅为本发明中的具体实施方式
，但本发明的保护范围并不局限于此，任 何熟悉该技术的人在本发明所揭露的技术范围内，可理解想到的变换或替换，都应涵盖在 本发明的包含范围之内，因此，本发明的保护范围应该以权利要求书的保护范围为准。
权利要求
一种用于下视阵列天线合成孔径雷达的三维聚焦成像方法，其特征在于，包括步骤如下步骤S1下视阵列天线合成孔径雷达系统采集到的成像区域原始回波信号为S0(t，u，υ)，斜距向压缩单元将原始回波信号沿斜距向进行傅立叶变换，生成第1信号为S1(ft，u，υ)，采用斜距向匹配滤波器H1(ft)对第1信号进行斜距向聚焦成像处理，生成第2信号为S2(ft，u，υ)，然后对第2信号进行斜距向傅里叶逆变换，生成第3信号为S3(t，u，υ)，其中，斜距向快时间域为t，斜距向频率域为ft，方位向时间域为u和地距向时间域为υ；步骤S2方位向距离徙动校正单元将生成的第3信号沿方位向进行傅里叶变换，生成第4信号为S4(t，fu，υ)，其中，方位向频率域为fu，然后对第4信号沿斜距向进行傅里叶变换，生成第5信号为S5(ft，fu，υ)，采用距离徙动校正因子H2(ft)对第5信号进行方位向距离徙动校正，生成第6信号为S6(ft，fu，υ)，然后对第6信号进行斜距向傅里叶逆变换，生成第7信号为S7(t，fu，υ)；步骤S3方位向压缩单元采用方位向匹配滤波器H3(ft)对第7信号进行方位向聚焦成像处理，生成第8信号为S8(t，fu，υ)，然后对第8信号进行方位向傅里叶逆变换，生成第9信号为S9(t，u，υ)；步骤S4地距向距离徙动校正单元将第9信号沿斜距向进行傅里叶变换，生成第10信号S10(ft，u，υ)，对第10信号进行地距向距离徙动校正，生成第11信号为S11(ft，u，υ)，对第11信号进行斜距向傅里叶逆变换，生成第12信号为S12(t，u，υ)；步骤S5地距向压缩单元首先采用解斜参考函数H6(υ)对第12信号进行解斜处理，生成第13信号为S13(t，u，υ)，然后对第13信号进行地距向傅里叶变换，生成第14信号为S14(u，r，θ)，该信号即为获得的成像区域在柱面坐标系中包含幅度和相位信息的下视阵列天线合成孔径雷达三维聚焦图像，其中，成像区域中第i个目标在柱面坐标系中方位向坐标值为xi，斜距向坐标值为ri，辐角向坐标值为θi。
2.根据权利要求1中所述的用于下视阵列天线合成孔径雷达的三维聚焦成像方法，其 特征在于，在成像几何中，下视阵列天线合成孔径雷达系统中阵列天线阵元间隔即系统沿 地距向的采样间隔d小于等于X。/[2sin (0/2)]，其中入。为雷达发射信号的载波波长， 3为发射天线的地距向波束宽度。
3.根据权利要求2中所述的用于下视阵列天线合成孔径雷达的三维聚焦成像方法，其 特征在于，所述采样间隔d小于等于De，其中De为发射天线在地距向上的尺寸。
4.根据权利要求2中所述的用于下视阵列天线合成孔径雷达的三维聚焦成像方法，其特征在于，所述采样间隔d小于等于√4II2+W2y/(2Wy)，其中H为系统平台飞行高度，ffy为系统在地距向上的成像幅宽。
5.根据权利要求1中所述的用于下视阵列天线合成孔径雷达的三维聚焦成像方法，其 特征在于，所述下视阵列天线合成孔径雷达系统中，阵列天线采用单发多收、收发分置的线 阵结构，其中单个发射天线位于线阵的中央，接收天线阵元沿线阵等间隔分布。
6.根据权利要求1中所述的用于下视阵列天线合成孔径雷达的三维聚焦成像方法，其 特征在于，在成像几何中，地距向的真实孔径中心位于下视阵列天线合成孔径雷达成像区 域的地距向中心。
7.根据权利要求1中所述的用于下视阵列天线合成孔径雷达的三维聚焦成像方法，其 特征在于，在成像几何中，成像区域的地距向成像幅宽大于下视阵列天线合成孔径雷达系 统地距向线阵真实孔径长度。
8.根据权利要求1中所述的用于下视阵列天线合成孔径雷达的三维聚焦成像方法，其 特征在于，所述步骤S4中，地距向距离徙动校正单元的处理步骤是步骤S41 将第9信号沿斜距向进行傅里叶变换，生成第10信号Sltl(ft，u，u)； 步骤S42 采用距离徙动校正因子扎(&)对第10信号进行地距向第一次距离徙动校 正，生成第101信号为S皿(ft，u, u),将第10信号中与斜距位置有关的距离徙动项移除； 步骤S43 采用距离徙动校正因子氏(&)对第101信号进行地距向第二次距离徙动校 正，生成第11信号为Sn (ft，u, u),将第101信号中与辐角位置有关的距离徙动项移除； 步骤S44 对第11信号进行斜距向傅里叶逆变换，生成第12信号为S12 (t，u，u )。
全文摘要
一种用于下视阵列天线合成孔径雷达的三维聚焦成像方法，涉及雷达技术，包括A)将采集到的下视阵列天线合成孔径雷达原始回波数据进行斜距向压缩；B)将A)步所得信号进行方位向距离徙动校正；C)将B)步所得信号进行方位向压缩E)将斜距向和方位向压缩后的信号进行地距向距离徙动校正F)将E)步所得信号进行地距向解斜和地距向傅里叶变换，重建出成像区域在柱面坐标系中的三维雷达图像。本发明方法使用较少的天线阵元来重建出成像区域的三维雷达图像，考虑了已有的三维成像方法中未考虑的距离徙动校正问题，并将已有的三维成像方法中的地距向波束形成运算变为一次复乘和傅里叶变换，即获取三维分辨图像，缩短了成像所需时间。
文档编号G01S13/90GK101866001SQ20091008199
公开日2010年10月20日 申请日期2009年4月15日 优先权日2009年4月15日
发明者吴一戎, 杜磊, 洪文, 王彦平 申请人:中国科学院电子学研究所