专利名称：一种非均匀分布的多基线合成孔径雷达三维成像方法
技术领域：
本发明涉及雷达三维成像技术领域，是一种针对非均勻分布的多基线合成孔径雷 达三维成像的处理方法。
背景技术：
多基线合成孔径雷达(SAR)通过雷达平台在高度方向上的多次重复飞行实现， 采集得到了不同视角下对目标沿高度方向的多次采样，构成了高度向合成孔径，具有高度 方向上的分辨能力；结合方位向上的孔径合成和距离向上的脉冲压缩，实现了对观测目 标的三维成像。2000年，德国学者A.Reigber等人，在论文《First Demonstration of Airborne SAR TomographyUsing Multibaseline L-Band Data》巾 冲f 了SAR H会H 成像的解析建模研究，并提出了使用谱估计方法实现多基线SAR三维成像的处理方法。在 后续石if 究中，F. Lomdardini、S. Guillaso 等人分别在论文《Adaptive spectral estimation for multibaseline SAR tomography with airborne L-band data》、《Polarimetric SAR Tomography》中，基于谱估计实现多基线SAR三维成像的原理，引入了 Capon、MUSIC等现 代信号处理方法进行高精度的高度向成像，获取了高分辨率的三维成像结果。这些处理方 法都是在基线均勻分布的情况下进行的，它是非均勻的多基线SAR数据实现三维成像的基 础。多基线SAR三维成像首先对每条轨迹采集的数据进行二维成像，获取目标在不同观测 视角下得到的二维聚焦数据向量，然后对高度向合成孔径的采样数据进行空间谱估计，得 到目标沿高度方向的图像，结合每次观测得到的目标距离-方位二维图像从而实现了对目 标的三维成像。一般的多基线合成孔径雷达三维成像的观测几何模型如图1所示。x为载机飞行 的方位向，y为地距方向，z为垂直高度方向。对目标场景在不同高度位置上共进行了 N次 观测，每次观测在正侧视条带式的观测模式下进行，总基线跨度为L，其与水平方向的夹角 为识。定义最低位置处的观测为参考位置，H是它的平台高度，参考位置处的波束中心视角为 9，其中心波束方向为参考斜距方向r，定义s为正交于载机飞行方向x和参考斜距方向r 的仰角高度方向，则建立了三维成像坐标系x-r-s，坐标原点位于参考平台位置处。第n次飞行观测所获取的目标二维聚焦信号可以表示为
f 2r (x' r' ’)) 「；r1 「万“r (x\r',s') = a(x\r\s') Qxp -j2m "v ‘， .sine -{x-xn) sine -(r-r )
IA J Ldx 」 LA 」以第一次观测获取的图像为参考图像，各幅图像进行配准，得到第n次观测的目 标二维聚焦数据，表示为fn(x',r',s') = a(x',r\s')-expf-j4^r"^X 'r 'S ) -sine —(x-x') sine
Ix ) L^x 」 LA .其中，a(x' , r'，s')为目标点的复散射系数，假定目标各向同性，其幅度和相 位在高度位置N次观测中保持不变；X为波长，rn(x' ,r' , s')为目标点与第n次观测 之间的距离，^^和、分别是二维聚焦成像后方位向和斜距向上的分辨率；x' , r' , s'分别为目标在方位向、距离向和高度向的坐标值。根据菲涅耳近似，rn(x' ,r' ,s')表示为 其中，ln为第n次观测与参考位置之间的距离。对不同高度位置采集得到的二维聚焦信号进行去除线性调频项的解斜处理，相位
因子为 经过解斜处理，第n次观测得到的聚焦数据表示为 此时，沿高度向采集信号的空间频率为 它只与目标散射点在仰角高度方向的位置有关，而与高度向数据采集位置无关。 因此通过傅里叶变换求取高度向合成孔径数据的空间谱，就可得到与目标高度位置有关的 函数表达式，实现对目标散射点沿高度方向的重建 当高度向的各次平行观测等间隔分布时，对目标沿高度向是均勻采样的，因此对 采集得到的高度向合成孔径数据直接使用快速傅里叶变换(FFT)的方法就可以求取空间 谱，实现对目标沿高度方向的重建。在实际的多基线SAR数据获取中，由于载机多次飞行无法保持相等的间隔，因 此采集的高度向合成孔径数据是非均勻的。针对非均勻的多基线合成孔径雷达数据，为 实现对目标的三维成像,2008年F. Lombardini等人在论文《3-D SAR Tomography :The Multibaseline Sector InterpolationApproalch》中使用部分插值的方法，重建未采集到 的数据，补充为均勻采样的数据，再利用已有的三维成像方法进行处理。2007年，F.Meglio 等人在论文《Three dimensional SAR image focusing from non-uniform samples》中也 利用非均勻采样定理将非均勻的采样数据重建为均勻的数据，然后使用基于FFT的方法重 建三维图像。这些方法以插值方法为基础，而根据高度向合成孔径数据的特点构造合适的 插值因子是非常困难的，难以得到精确的高度向成像结果；利用这些方法所得到的三维成 像仿真实验结果并不理想，对实际非均勻多基线SAR数据的处理还没有开展。在现代信号处理领域，关于非均勻采样信号的空间谱估计，已经由许多研究者在研究，提出了许多非均勻采样信号空间谱重建的方法。2005年，Yanwei Wang、Jian Li、 P. Stoica((Spectral Analysis of Signal, The Missing Data Case)) , J^Hj
了使用缺失数据幅度相位估计(Missing-data Amplitude and Phase Estimation,MAPES) 方法对非均勻采样的数据进行空间谱估计，得到了较好的谱估计结果。
发明内容
本发明的目的在于提出一种非均勻分布的多基线合成孔径雷达三维成像方法，以 降低由于非均勻采样所带来的高程模糊，获得高分辨率的三维成像结果。为实现上述目的，本发明的技术解决方案是一种非均勻分布的多基线合成孔径雷达三维成像方法，是通过雷达平台在不同高 度位置的多次重复飞行采集多基线合成孔径雷达的数据，各次飞行轨迹之间是非等间隔分 布的，获取的多基线合成孔径雷达数据是非均勻采样的；多基线合成孔径雷达各次飞行的轨迹保持平行，各次飞行是正侧视、斜视、条带式 或聚束式观测模式成像，能够单独生成观测场景的二维图像；各次轨迹的排列方向，是沿垂 直高度方向排列，或沿水平方向排列，或沿与水平方向有一夹角的方向排列；在基线非均勻分布的情况下，使用缺失数据的幅度和相位估计方法进行对非均勻 采样的高度向数据求取空间谱，得目标沿高度方向的分辨率图像，结合二维图像，实现对目 标沿高度方向的三维成像。所述的非均勻分布的多基线合成孔径雷达三维成像方法，其包括如下步骤步骤A 对每次飞行采集的目标场景原始回波数据进行二维成像，生成观测场景 距离-方位平面的图像序列，得到高度向每次观测下目标的二维聚焦复数据；步骤B 对二维图像序列进行配准，以最低轨道位置处的图像为参考图像，其余图 像经过位置移动与其对准；步骤C:对每幅图像进行相位补偿的解斜处理，去除高度向采样信号中的线性调 频项，使得频率值与目标高度位置线性相关；步骤D:提取目标由各次飞行观测得到的并经过解析处理后的二维复数据向量， 它是目标高度向合成孔径的非等间隔采样数据；步骤E 对每个距离_方位分辨单元内的目标高度向非均勻采样数据，使用缺失 数据幅度相位估计方法对非均勻的数据向量求取空间谱，得到目标沿高度方向的分辨率图 像；步骤F 将步骤E中空间谱估计得到的目标高度向图像，结合步骤A生成的观测场 景距离_方位平面的图像，得到目标在斜距_方位_仰角高度三维空间内的三维分辨率图 像；步骤G 斜地高程坐标转换，将目标三维成像结果的斜距_方位_仰角高度坐标， 转换为地距_方位_垂直高度方向坐标。所述的非均勻分布的多基线合成孔径雷达三维成像方法，其所述步骤E中得到目 标沿高度方向的分辨率图像，包括如下步骤步骤E1 由已观测的非均勻分布的多基线SAR数据，得到目标沿高度向观测数据 所在的非均勻采样的位置；根据这些非均勻采样的位置，以最小采样距离为单位间隔，重新划分采样位置，得到采样间隔相等的数据采集位置，从而重建出均勻采样时多基线SAR数 据位置分布；步骤E2 求取高度向合成孔径数据空间谱和噪声协方差矩阵的初始值；即根据重 新划分出的均勻采样位置，将采样位置处未观测到的数据置零，得到目标高度向合成孔径 数据向量，通过快速傅里叶变换求取初始值；步骤E3 根据计算得到的空间谱和噪声协方差矩阵值，求取以目标高度向已知数 据为条件的未知数据的数学期望和协方差矩阵；步骤E4:利用已求的数学期望和协方差矩阵，更新高度向合成孔径数据的空间谱 和噪声协方差矩阵；步骤E5 判断估计的空间谱是否收敛，当前空间谱与前一空间谱的功率之差小于 预设门限时，结果收敛，执行下一步；反之，如果功率差大于门限值，结果不收敛，返回步骤 E3 ；步骤E6 将估计得到的空间谱作为目标沿高度向的分辨率图像。所述的非均勻分布的多基线合成孔径雷达三维成像方法，其所述步骤E1，是将高 度向合成孔径数据分割成L个相互重叠的数据段、hM ...h^f ,1 = 1，. . . L，其中， h为目标高度向合成孔径数据第1个采样数据，K表示每个数据段中共有K个元素；同时得 到数据分段后各数据段中的已观测数据gl和未观测数据P ！；所述的非均勻分布的多基线合成孔径雷达三维成像方法，其所述步骤E2，是根据 采样数据的位置求取观测数据和缺失数据与均勻采样数据之间的酉变换矩阵Si、S2，它们之 间的关系表示为g,=片(/)6，约=对(/)6，其中，1为分割后的第1个数据段。所述的非均勻分布的多基线合成孔径雷达三维成像方法，其所述步骤E2中的高 度向合成孔径数据的空间谱初始值通过对整个高度向合成孔径数据进行快速傅里叶变换 计算得到，表示为；噪声协方差矩阵初始值的计算通过下式得到幺…一礼⑷+，其中，L为分割得到的
数据段的个数，《为频率值，a(co)为导向矢量，表示为a (co) = ..…(K_1)U]。所述的非均勻分布的多基线合成孔径雷达三维成像方法，其所述步骤E3中，求取 各个数据段的数学期望h和协方差矩阵&的公式为b^Sl^a^)q,_} (co)eJa"+st2 (/)a_, ns[ (/)[《(/)a_, (co)s{ (/)]" (gl - s{ (i)a(co)r )^/=52r(/)a_1(^)<s2r(/)
_2] -Sr2 (m (co)s[ (/)[《(/)么(Oj1 (m (0;其中，‘(咚)为第i_l次迭代之后得到的高度向合成孔径数据空间谱。么—i为第 i"l次迭代之后得到噪声协方差矩阵。所述的非均勻分布的多基线合成孔径雷达三维成像方法，其所述 步骤E4中，更新高度向合成孔径数据的空间谱以及噪声协方差矩阵公式
为高度向合成孔径数据的空间谱为q{o>)= H) [-,) i ；(，噪声协方
差矩阵为=々(⑴)+ [《(《) (仍)其中，
l i=ol i=0l /=oZl = SjDgi+SjDbp所述的非均勻分布的多基线合成孔径雷达三维成像方法，其所述步骤E5中，执行 下一步的收敛条件由所求得的空间谱表示为
k-]2 k-\2
一ZIU吟)1jt=0_k^Q_< 1
— - —lv
k=0其中，々,(咚)为第i次迭代之后的空间谱，I,(吟)为第i_l次迭代之后的空间谱， K为每个数据段中元素的个数。所述的非均勻分布的多基线合成孔径雷达三维成像方法，其所述步骤E5中，预设 门限设为10_3，空间谱迭代收敛的公式表示为其中，Pi为第i迭代计算的空间谱功率，Ph为第i_l次迭代计算的空间谱功率。本发明方法的有益效果是针对非均勻的高度向合成孔径采样数据，使用缺失数 据幅度和相位估计(Missing-data Amplitude and PhaseEstimation，MAPES)方法，不进行 插值运算，直接使用采集得到的非均勻数据，以观测数据的数学期望最大化为标准，通过迭 代运算的方法估计高度向合成孔径数据的空间谱，对非均勻采样的高度向合成孔径数据进 行空间谱的精确重建，从而获取目标的高分辨率三维成像。
图1是非均勻的多基线SAR三维成像观测几何示意图；图2是本发明方法的非均勻多基线SAR三维成像处理流程图；图3是本发明方法的非均勻数据高度向成像处理流程图；图4是本发明方法的非均勻数据高度向成像具体实施步骤图。
具体实施例方式本发明的一种非均勻分布的多基线合成孔径雷达三维成像方法，具体步骤是步骤A 对每次飞行采集的目标场景原始回波数据进行二维成像，生成观测场景 距离-方位平面的图像序列，得到高度向每次观测下目标的二维聚焦复数据；步骤B 对二维图像序列进行配准，以最低轨道位置处的图像为参考图像，其余图 像经过位置移动与其对准；步骤C:对每幅图像进行相位补偿的解斜处理，去除高度向采样信号中的线性调
9频项，使得频率值与目标高度位置线性相关；步骤D:提取目标由各次飞行观测得到的并经过解析处理后的二维复数据向量， 它是目标高度向合成孔径的非等间隔采样数据；步骤E 对每个距离_方位分辨单元内的目标高度向非均勻采样数据，使用缺失数 据幅度禾口相位估计(Missing-data Amplitude and PhaseEstimation, MAPES)方法对非均 勻的数据向量求取空间谱，得到目标沿高度方向的分辨率图像；步骤F 结合步骤A得到的目标距离_方位二维图像，与空间谱估计得到的目标高 度向图像，得到目标在斜距_方位_仰角高度三维空间内的三维分辨率图像。步骤G 斜地高程坐标转换，将目标三维成像结果的斜距_方位_仰角高度坐标， 转换为地距_方位_垂直高度方向坐标。在本发明的方法中，多基线合成孔径雷达数据通过载机多次重复飞行获得，构成 对目标观测的高度向合成孔径；同时由于多次飞行轨迹的非间隔分布，对高度向合成孔径 的采样是非均勻的，高度向非均勻采样的位置由载机飞行系统记录。每次飞行的观测模式相同，并生成相同分辨率的二维图像，观测模式可以是条 带式或聚束式、正侧视或斜视；对场景的二维成像由传统的SAR 二维成像算法实现，如距 离-多普勒算法(R-D)、线性调变算法(CS)、距离徙动算法(RMA)等。在本发明的方法中，序列图像的配准需要达到亚像素级的精配准，可以采用干涉 合成孔径雷达技术中的配准方法，与参考图像分别进行两两配准，最后实现对图像序列的 全部配准。在本发明的方法中，第n幅图像进行解斜处理时的相位补偿因子为 W =eXp|y'，/ 2cOS2p-^)j淇中，人为波长，r’为目标参考斜距值，1 为观测位置相对
于参考平台的距离，9为观测视角，炉为基线排列方向与水平方向夹角。在本发明的方法中，对目标的高度向成像使用谱估计的方法，并使用缺失数据幅 度估计方法对非均勻采样的高度向合成孔径数据进行空间谱重建，实现目标的高程成像。本发明中使用非均勻的高度向合成孔径采样数据实现对目标高程成像的具体步 骤为步骤E1 根据已观测的多基线SAR数据高度向采样位置，重新划分出高度向均勻 采样时的位置分布；步骤E2 将重新划分的采样位置处未观测到的数据置零，得到目标高度向合成孔 径数据向量，通过快速傅里叶变换的方法求取高度向合成孔径数据空间谱和噪声协方差矩 阵的初始值；步骤E3 根据计算得到的空间谱和噪声协方差矩阵值，求取以已知数据为条件的 未知数据的数学期望和协方差矩阵；步骤E4:利用已求的数学期望和协方差矩阵，更新高度向合成孔径数据的空间谱 和噪声协方差矩阵；步骤E5 判断估计的空间谱是否收敛，当前空间谱功率与前一空间谱功率差小于 预设门限时，结果收敛，执行下一步；反之，返回步骤E3 ；步骤E6 将估计得到的空间谱作为目标沿高度向的分辨率图像。
在本发明的方法中，高度向合成孔径数据空间谱和噪声协方差矩阵初始估计值表 示为q0( ) = FFT[hm],Q0(co) = E [e ( ) eH( )]，其中，FFT表示快速傅里叶变换，E表示 求取数学期望；gn，n = 1，. . . .，N是观测得到的高度向合成孔径非均勻采样数据，由此重建 得到的均勻采样数据为hm，m = 1，. . . M，即将未观测到的数据置零，同时缺失的数据向量表 示为P i，i = 1，. . .，M-N。高度向合成孔径采样数据可以表示为hm = q ((0) eJ"m+effl ( )，m = 1，...，M，co e
为导向矢量，Si、S2分别为观测数据和缺失数
据与均勻采样数据之间的酉变换矩阵，它们之间的关系表示为。在本发明的方法中，空间谱迭代计算的收敛条件表示为
2 K-l2
- ZM崎)j
、10_3，其中， 为第i次迭代之后的空间谱，为第i_l
次迭代之后的空间谱。下面结合附图，详细说明本发明的一种非均勻分布的多基线合成孔径雷达三维成 像方法，所涉及的各个细节问题。应指出的是，所描述的实施例仅旨在便于对本发明的理 解，而对其不起任何限定作用。利用采集的多基线SAR回波数据，通过对目标的二维聚焦成像、解斜处理，得到对 目标观测的高度向合成孔径数据。由于高度方向上非等间隔分布的观测位置，得到了目标 高度向合成孔径非均勻的采样数据，对其经过解斜处理，得到以上分析中的非均勻采样数 据gn，n = 1，....，N。假定高度向合成孔径完整的均勻采样数据为hm，m= 1，...M，观测数 据gn，n = 1，....，N是从hm中非均勻采样得到，其中缺失的数据向量为^，i = 1，...， M-N。根据谱估计方法实现多基线SAR三维成像的原理，高度向合成孔径数据的空间谱即为 目标沿高度方向的分辨率图像。对于非均勻采样的数据向量，我们使用缺失数据的幅度相 位估计(Missing-data Amplitude and Phase Estimation，MAPES)方法求取原始数据向量 的空间谱。高度向合成孔径采样数据表示为hm = q(o)eJ"m+em(w), m = 1，.，M，w G
其中，q( )表示频率中的幅度复数值，即为所要估计的空间谱，em( )表示噪声。将数据向量h分割成L个长度为KX1的相互重叠的子向量，表示为h,=[h, hM …矣1 = 1, L，L = N-K+lo第1个子向量数据可以表示为+司( )，其中，a(co)= [le,..^—1)"]为导向矢量，司= ..、“—)]为零均值的高斯随机向量， Q{co) 二 E[es (co)^ ( )]为未知的噪声协方差矩阵。非均勻数据的谱估计是利用已观测的采样信号估计任意频率位置《处的复数值 冷( )。缺失数据幅度相位估计(MAPES)方法通过最大化已知数据的条件数学期望，来求解 空间谱的值。对每个子向量的条件数学期望表示为+1”(⑷ |咖)，00))其中，p为求取概率密度函数。令a = {q( )，Q( )}，利用已估计出的新数据向 量进行迭代运算来求解该数学期望的最大化问题，则空间谱求解的最大化问题表示为a' = arg max E^lnp(g,ju\a)\g, a^}计算迭代运算的初始条件，即将未观测到的数据置零，使用快速傅里叶变换(FFT) 的方法对已观测的数据求取空间谱，并利用空间谱初始值求取噪声协方差矩阵初始值q,(co) = FFT[hm}4⑷二+￡p-⑷一]—^(似)《⑷eJalT
L i=Q利用空间谱和噪声矩阵求取观测数据的数学期望b和协方差矩阵K，高度向合 成孔径的每个子向量中缺失数据谱密度函数在已观测数据为条件下满足高斯分布，即 ^\g a'-x ~CN{b Kt)，其中CN表示高斯分布，。它的数学期望和协方差表示为b^E^g^}，A = cov{约|g,，}使用计算的b和K来更新空间谱和噪声矩阵a = {q(co),Q(co)}的数值。
制 a" {co)S-\co)Z{co)_)、"( —)e( ) = + (oi)a{(o)-Z{(o)\[q} (a)a(o})-Z(w)]H其中j ⑷=钱(O + jl^f-Z ⑷ ⑷，=，
Zl = Si (1) gl+S2 (1) bp Si、S2分别为观测数据和缺失数据与均勻采样数据之间的u变换矩阵，
^Ag^slii)^
表不为 cr.A,。将得到的空间谱和噪声矩阵再去求取b和K，然后重复进行迭代运算，当空间谱的 功率与前一空间谱功率差小于特定的门限时，即
结果收敛，得到最终收敛的空间谱估计结果W )。本发明所提出的非均勻多基线SAR数据三维成像方法中高程成像是本发明的核 心内容，其具体实施步骤如图4所示，主要有步骤S1 将高度向合成孔径数据分割成L个相互重叠的数据段， ht=[h, hM...hUK_,f ,1 = 1，...L，并得到分段后已观测数据gl和未观测数据同时根 据观测位置重建均勻采样的位置。步骤S2 根据采样数据的位置求取观测数据和缺失数据与均勻采样数据之间的U 变换失巨阵 &、&，&= Si (l)^, Ml= St2 (I) h,。步骤S3 高度向合成孔径数据的空间谱初始值，和噪声协方差矩阵初始值q。(《 ) = FFT [gn]，& ⑷=+ L B 1 M历)] —孔(历)a (历)^ T 导向
矢量 a(co) = [1^..』劍”。步骤S4 由空间谱和噪声协方差矩阵的值，求取各个数据段的数学期望和协方差 矩阵bl^ST2{l)a(co)ql-X{co)ej011+《(/)eM (co)s( (/)[X (z)e'-1 —K (i)T (/)aHr (^)-Si (/)纩(co)Sf (/)[《(/)eM HSl (/)『'Si (/)纩(co)St2 (/)步骤S5 求取高度向合成孔径数据的空间谱以及噪声协方差矩阵
“、aH (co)S-Uco)Z(o))q(co)= H\ ； ~ ,v. ； ； ！ ( )=々(《) +[象[么其中，
^(^)41^(oks"(o41-z(^)(^)‘半，
L i=o厶 /=oL /=oZl = sjDgi+sjDiv步骤S6 判断收敛条件，当前空间谱功率与前一空间谱功率差小于预设门限时， 结果收敛，执行下一步，即
K-\2 K-12^-—~巴-< 10"3
人丨2
Zl^Ml
k=0反之，返回步骤S4，重复上述步骤；
13
步骤S7 将迭代收敛后的空间谱结果，作为目标高度向的分辨率图像。获取目标高度向的分辨率图像之后，结合每条轨道生成的目标二维图像，得到目 标在斜距-方位-仰角高度方向的三维成像结果P(x，r，s)；通过斜距地距以及高程转换，
“1 _
将目标成像结果转换到地距_方位-垂直高度三维坐标系中，即= 力，其
sin 6
中，9为参考视角，得到易于直观显示的目标三维成像结果P(x，y，z)。本发明的方法，已经在计算机上应用MATLAB软件得到验证，并且对微波暗室采集 非均勻多基线SAR数据进行了三维成像处理，方法的有效性得到了验证。以上所述，仅为本发明中的具体实施方式
，但本发明的保护范围并不局限于此，任 何熟悉该技术的人在本发明所揭露的技术范围内，可理解想到的变换或替换，都应涵盖在 本发明的包含范围之内，因此，本发明的保护范围应以权利要求书的保护范围为准。
权利要求
一种非均匀分布的多基线合成孔径雷达三维成像方法，是通过雷达平台在不同高度位置的多次重复飞行采集多基线合成孔径雷达的数据，各次飞行轨迹之间是非等间隔分布的，获取的多基线合成孔径雷达数据是非均匀采样的；其特征在于多基线合成孔径雷达各次飞行的轨迹保持平行，各次飞行是正侧视、斜视、条带式或聚束式观测模式成像，能够单独生成观测场景的二维图像；各次轨迹的排列方向，是沿垂直高度方向排列，或沿水平方向排列，或沿与水平方向有一夹角的方向排列；在基线非均匀分布的情况下，使用缺失数据的幅度和相位估计方法进行对非均匀采样的高度向数据求取空间谱，得目标沿高度方向的分辨率图像，结合二维图像，实现对目标沿高度方向的三维成像。
2.如权利要求1所述的非均勻分布的多基线合成孔径雷达三维成像方法，其特征在 于包括如下步骤步骤A 对每次飞行采集的目标场景原始回波数据进行二维成像，生成观测场景距 离-方位平面的图像序列，得到高度向每次观测下目标的二维聚焦复数据；步骤B 对二维图像序列进行配准，以最低轨道位置处的图像为参考图像，其余图像经 过位置移动与其对准；步骤C:对每幅图像进行相位补偿的解斜处理，去除高度向采样信号中的线性调频项， 使得频率值与目标高度位置线性相关；步骤D 提取目标由各次飞行观测得到的并经过解析处理后的二维复数据向量，它是 目标高度向合成孔径的非等间隔采样数据；步骤E 对每个距离-方位分辨单元内的目标高度向非均勻采样数据，使用缺失数据幅 度相位估计方法对非均勻的数据向量求取空间谱，得到目标沿高度方向的分辨率图像；步骤F 将步骤E中空间谱估计得到的目标高度向图像，结合步骤A生成的观测场景距 离_方位平面的图像，得到目标在斜距_方位_仰角高度三维空间内的三维分辨率图像；步骤G 斜地高程坐标转换，将目标三维成像结果的斜距_方位_仰角高度坐标，转换 为地距_方位_垂直高度方向坐标。
3.如权利要求2所述的非均勻分布的多基线合成孔径雷达三维成像方法，其特征在 于所述步骤E中得到目标沿高度方向的分辨率图像，包括如下步骤步骤El 由已观测的非均勻分布的多基线SAR数据，得到目标沿高度向观测数据所在 的非均勻采样的位置；根据这些非均勻采样的位置，以最小采样距离为单位间隔，重新划分 采样位置，得到采样间隔相等的数据采集位置，从而重建出均勻采样时多基线SAR数据位 置分布；步骤E2 求取高度向合成孔径数据空间谱和噪声协方差矩阵的初始值；即根据重新划 分出的均勻采样位置，将采样位置处未观测到的数据置零，得到目标高度向合成孔径数据 向量，通过快速傅里叶变换求取初始值；步骤E3 根据计算得到的空间谱和噪声协方差矩阵值，求取以目标高度向已知数据为 条件的未知数据的数学期望和协方差矩阵；步骤E4 利用已求的数学期望和协方差矩阵，更新高度向合成孔径数据的空间谱和噪 声协方差矩阵；步骤E5 判断估计的空间谱是否收敛，当前空间谱与前一空间谱的功率之差小于预设门限时，结果收敛，执行下一步；反之，如果功率差大于预设门限值，结果不收敛，返回步骤 E3 ；步骤E6 将估计得到的空间谱作为目标沿高度向的分辨率图像。
4.如权利要求3所述的非均勻分布的多基线合成孔径雷达三维成像方法，其特 征在于所述步骤E1，是将高度向合成孔径数据分割成L个相互重叠的数据段，即 Til^hl hM ...hUKJ ,I = 1, ... L，其中，Ill为目标高度向合成孔径数据第1个采样数据，K 表示每个数据段中共有K个元素；同时得到数据分段后各数据段中的已观测数据gl和未观 测数据μ ρ
5.如权利要求3所述的非均勻分布的多基线合成孔径雷达三维成像方法，其特征在 于所述步骤Ε2，是根据采样数据的位置求取观测数据和缺失数据与均勻采样数据之间的 酉变换矩阵S” S2，它们之间的关系表示为g,=片(/μ,, M1=St2(I)H1，其中，1为分割后的第 1个数据段。
6.如权利要求3所述的非均勻分布的多基线合成孔径雷达三维成像方法，其特征在 于所述步骤Ε2中的高度向合成孔径数据的空间谱初始值通过对整个高度向合成孔径数 据进行快速傅里叶变换计算得到，表示为纟。(《)；噪声协方差矩阵初始值的计算通过下式得 到a…)=+ΣΡ - ^—Μ ), &(iyMw)W Γ，其中，L为分割得到的数据段的个数，ω为频率值，a(co)为导向矢量，表示为a(co) = ..…(κ_1)ω]。
7.如权利要求3所述的非均勻分布的多基线合成孔径雷达三维成像方法，其特征在 于所述步骤Ε3中，求取各个数据段的数学期望Id1和协方差矩阵K1的公式为δ,=《(/)+)‘ (f+ Sl (/)0_(/)[片(l)Q,_, (ω)S^ (Ζ)]" (gl - (1)α(ω)ξ^ )K1 =St2 (1)^)^(1)—Sl (/)β-, Hs! (/)[《(I)Q1, (ω)Ξ[ (/)]" S[ (/)α_, —K (0 ；其中，U饵)为第i_l次迭代之后得到的高度向合成孔径数据空间谱。I1为第H次 迭代之后得到噪声协方差矩阵。
8.如权利要求3所述的非均勻分布的多基线合成孔径雷达三维成像方法，其特征在 于所述步骤E4中，更新高度向合成孔径数据的空间谱以及噪声协方差矩阵的计算公式 为Λ/ 、 αΗ (ω)S'1 (ω)Z(ω)高度向合成孔径数据的空间谱为g —二 Λ [^1)；；, >噪声协方差矩阵为 ΚΗ外…+|>( ) ( )-2( )][奴《>(6))-ζ(…]"； 其中，= 补Κ Ο + Σν,-ΖΗΖ ，=，L /=0L /=0L· ；=0Z1 = S1 (DgfS2 (I)Iv
9.如权利要求3所述的非均勻分布的多基线合成孔径雷达三维成像方法，其特征在 于所述步骤E5中，执行下一步的收敛条件由所求得的空间谱表示为 其中，么(饵)为第i次迭代之后的空间谱Jm (饵)为第i-Ι次迭代之后的空间谱，K为 每个数据段中元素的个数。
10.如权利要求3所述的非均勻分布的多基线合成孔径雷达三维成像方法，其特征在 于所述步骤E5中，预设门限设为10_3，空间谱迭代收敛的公式表示为|Pi-Pi-1|/Pi-1≤10-3其中，Pi为第i迭代计算的空间谱功率，Pi^1为第i-ι次迭代计算的空间谱功率。
全文摘要
本发明公开了一种非均匀分布的多基线合成孔径雷达三维成像方法，涉及三维成像技术，首先将各次飞行观测得到的原始回波数据通过二维聚焦得到单视复图像，然后将单视复图像序列配准，获取观测目标在不同视角下的非均匀采样数据；针对这些非均匀采样的数据，经过解斜处理进行相位调制；然后使用基于缺失数据幅度相位估计的方法，通过最大化观测数据数学期望的迭代运算，估计出原始均匀采样数据的空间谱，从而实现对目标高度方向的成像；结合每条轨道二维成像中得到的目标二维图像，完成对目标的三维成像。本发明基于幅度相位估计方法进行高度向成像，降低了由于多基线非均匀分布所造成的高程模糊，获取清晰的高分辨率的目标三维成像结果。
文档编号G01S13/90GK101893710SQ20091008416
公开日2010年11月24日 申请日期2009年5月20日 优先权日2009年5月20日
发明者吴一戎, 洪文, 王彦平, 王斌, 谭维贤 申请人:中国科学院电子学研究所