3-D SAR波数域快速成像的方法技术

本发明专利技术提供了一种3-D?SAR波数域快速成像的方法。该方法只对感兴趣区域而不是整个三维观测区域场景进行重建，从而避免了数据量的激增，能够进行快速成像处理。

【技术实现步骤摘要】
3-D SAR波数域快速成像的方法
本专利技术涉及雷达成像和信号处理
，尤其是机载下视阵列3-DSAR和多基线3-D SAR工作模式中跨航向孔径长度远小于跨航向测绘带宽情况下的一种3-D SAR波数域快速成像的方法。
技术介绍
机载阵列下视3-D SAR、多基线3_D SAR能够获得观测区域场景三维散射信息，能够克服常规侧视SAR中的阴影、叠掩等问题，因此具有广泛的应用前景。目前3-D SAR成像处理主要是先在波传播方向和航迹向使用常规侧视SAR成像处理方法中的RD、CS、ω-k等方法进行处理，然后在跨航向使用SPECAN、波束形成、压缩感知等方法进行处理。这类成像处理方法对距离历程进行泰勒展开近似处理，仅保留泰勒展开中的一次项和二次项，将三维成像转换成波传播和航迹向二维处理与跨航向一维处理。由于忽略了距离历程中的三次及三次以上的高次项，因此对三维观测区域场景不能够精确重建，尤其是雷达到观测区域距离较短及观测区域测绘带较大的时候，距离历程的近似导致的重建误差更大。三维波数域方法能够完全补偿场景中的距离徙动，重建精度高，但是该类方法要求回波采集时综合孔径长度不小于成像区域以防止FFT出现卷绕，在机载下视阵列3-DSAR中，跨航向阵列天线长度一般为几米，而跨航向幅宽为几百米到几千米，因此需要对回波数据进行大量补零，补零会带来内存需求和计算量的激增。若直接使用三维波数域方法对整个观测区域场景进行重建，无论是数据量还是计算量都非常大，从而限制了该方法在机载下视阵列3-D SAR中的应用。
技术实现思路
(一 )要解决的技术问题为解决上述的一个或多个问题，本专利技术提供了一种3-D SAR波数域快速成像的方法。( 二 )技术方案根据本专利技术的一个方面，提供了一种3-D SAR波数域快速成像的方法。该方法包括:步骤A，从3-D SAR回波数据中提取波传播-航迹向二维数据和波传播-跨航向二维数据；步骤B，对波传播-航迹向二维数据和波传播-跨航向二维数据分别采用二维成像处理方法进行重建，得到波传播-航迹向二维图像和波传播-跨航向二维图像，其中σ为雷达后向散射系数；步骤C，由波传播-航迹向二维图像和波传播-跨航向二维图像反映出的观测区域场景目标散射特性，搜寻得到整个三维观测区域场景中感兴趣区域；步骤D，对3-DSAR回波数据进行匹配滤波，并沿波传播向设置距离门限，使感兴趣区域落在波传播方向距离门限内，感兴趣区域之外的区域落在波传播方向距离门限之外，仅保留距离门限内的数据；步骤Ε，对航迹向空域、跨航向空域、波传播向空域信号沿航迹向、跨航向和波传播向进行三维FFT，得到航迹向、跨航向和波传播向三维波数域信号；步骤F，将波传播方向、航迹向和跨航向三维波数域信号进行场景中心距离徙动校正和残余距离徙动校正；以及步骤G，对完成场景中心距离徙动校正和残余距离徙动校正的波传播方向、航迹向和跨航向三维波数域信号沿波传播方向、航迹向和跨航向进行三维IFFT，得到精确重建的感兴趣区域三维图像。(三)有益效果从上述技术方案可以看出，本专利技术3-D SAR波数域快速成像的方法只对感兴趣区域而不是整个三维观测区域场景进行重建，从而避免了数据量的激增，能够进行快速成像处理。【附图说明】图1为机载下视阵列3-D SAR成像场景的立体示意图；图2A至图2C分别为图1所示成像场景沿波传播-轨迹向、轨迹向-跨航向、波传播-跨航向的剖面示意图；图3为机载下视阵列3-D SAR波束覆盖区与感兴趣区域的示意图；图4为本专利技术实施例3-D SAR波数域快速成像方法的流程图；图5A和图5B为地基雷达实验场景目标放置位置沿不同视角的视图。图6A至图6F为地基雷达实验数据采用本专利技术方法的处理结果，其中:图6A为Y方向-Z方向二维波数域重建和ROI选取；图6B为X方向-Z方向二维波数域重建和ROI选取；图6C为ROI三维波数域重建结果；图6D为ROI三维重建结果在X-Y平面最大值投影；图6E为ROI三维重建结果在X-Z平面最大值投影；[0021 ]图6F为ROI三维重建结果在Y-Z平面最大值投影。【具体实施方式】为使本专利技术的目的、技术方案和优点更加清楚明白，以下结合具体实施例，并参照附图，对本专利技术进一步详细说明。需要说明的是，在附图或说明书描述中，相似或相同的部分都使用相同的图号。附图中未绘示或描述的实现方式，为所属
中普通技术人员所知的形式。本专利技术3D SAR波数域快速成像方法适用于机载下视阵列3-D SAR、多基线3_D SAR等成像模式中跨航向孔径长度远小于跨航向幅宽(Le/[2H sin(0e/2)]≤0.01，Le为跨航向孔径长度，一般为几米；2H sin(0e/2)为跨航向幅宽，一般为千米量级，其中H为载机飞行高度，为天线在跨航向波束宽度)的情况下，由于观测区域场景在整个三维空间中呈现稀疏分布，直接对观测区域场景中感兴趣区域(Region of Interest,简称R0I)进行重建，既能降低数据量又能减少运算量。图1为机载下视阵列3-D SAR成像场景的立体示意图。请参照图1，航迹向平行于载机平台运动方向，将其定义为X轴；跨航向平行于阵列向，将其定义为Y轴；波传播向自上而下垂直于XY平面，将其定义为Z轴。X轴、Y轴和Z轴相交于坐标原点O。P为观测区域场景目标。图2A至图2C分别为图1所示成像场景沿波传播_轨迹向、轨迹向_跨航向、波传播-跨航向的剖面示意图。从图2A至图2C可看出，整个三维观测区域场景在不同平面上的投影呈现不同的目标分布特征。目标在不同平面上投影后的分布特征，可在进行ROI确定时作为一个重要的依据和参考。图3为机载下视阵列3-D SAR波束覆盖区与感兴趣区域的示意图。请参照图3中左图，与雷达相作用的观测区域场景只是数字表面模型部分和电磁波穿透部分，从整个三维空间分布来看，观测区域场景是稀疏分布的。而感兴趣区域只占整个三维观测区域场景中很小的一部分，如图3右图所示。在本专利技术的一个示例性实施例中，提供了一种3-D SAR波数域快速成像的方法。图4为本专利技术实施例3-D SAR波数域快速成像方法的流程图。如图4所示，本实施例包括如下步骤:步骤A:从3-D SAR回波数据j中提取波传播-航迹向二维数据和波传播-跨航向二维数据；针对3-D SAR获取的三维回波数据，其中\代表航迹向采用位置，7?代表跨航向采样位置，；代表波传播方向快时间采样，从三维回波数据乂中沿跨航向中心提取波传播-航迹向(即X-Z方向)二维回波数据，其中，ync代表跨航向采样中心，从三维回波数据中沿航迹向中心提取波传播-跨航向(即Y-Z方向)二维回波数据，其中，Xm。代表航迹向采样中心。步骤B:对波传播-航迹向二维数据和波传播-跨航向二维数据P)分别采用二维成像处理方法进行重建，重建后得到波传播-航迹向二维图像O (χ,ζ)和波传播-跨航向二维图像σ (y，z)，其中σ为雷达后向散射系数；由于航迹向合成孔径由载机运动获得，一般比较长，因此对于波传播-航迹向二维数据的重建一般选取二维波数域方法、极坐标方法和后向投影方法。由于跨航向阵列长度较短，由于需要沿跨航向大量补零以防止FFT卷绕，不太适合选取二维波数域方法进行重建，因此对于波传播-跨航向二维图像一般选取极本文档来自技高网...

【技术保护点】
一种3?D?SAR波数域快速成像的方法，其特征在于，包括：步骤A，从3?D?SAR回波数据中提取波传播?航迹向二维数据和波传播?跨航向二维数据；步骤B，对波传播?航迹向二维数据和波传播?跨航向二维数据分别采用二维成像处理方法进行重建，得到波传播?航迹向二维图像和波传播?跨航向二维图像，其中σ为雷达后向散射系数；步骤C，由波传播?航迹向二维图像和波传播?跨航向二维图像反映出的观测区域场景目标散射特性，搜寻得到整个三维观测区域场景中感兴趣区域；步骤D，对3?D?SAR回波数据进行匹配滤波，并沿波传播向设置距离门限，使感兴趣区域落在波传播方向距离门限内，感兴趣区域之外的区域落在波传播方向距离门限之外，仅保留距离门限内的数据；步骤E，对航迹向空域、跨航向空域、波传播向空域信号沿航迹向、跨航向和波传播向进行三维FFT，得到航迹向、跨航向和波传播向三维波数域信号；步骤F，将波传播方向、航迹向和跨航向三维波数域信号进行场景中心距离徙动校正和残余距离徙动校正；以及步骤G，对完成场景中心距离徙动校正和残余距离徙动校正的波传播方向、航迹向和跨航向三维波数域信号沿波传播方向、航迹向和跨航向进行三维IFFT，得到精确重建的感兴趣区域三维图像。...

【技术特征摘要】
1.一种3-D SAR波数域快速成像的方法，其特征在于，包括: 步骤A，从3-D SAR回波数据中提取波传播-航迹向二维数据和波传播-跨航向二维数据； 步骤B，对波传播-航迹向二维数据和波传播-跨航向二维数据分别采用二维成像处理方法进行重建，得到波传播-航迹向二维图像和波传播-跨航向二维图像，其中σ为雷达后向散射系数； 步骤C，由波传播-航迹向二维图像和波传播-跨航向二维图像反映出的观测区域场景目标散射特性，搜寻得到整个三维观测区域场景中感兴趣区域； 步骤D，对3-D SAR回波数据进行匹配滤波，并沿波传播向设置距离门限，使感兴趣区域落在波传播方向距离门限内，感兴趣区域之外的区域落在波传播方向距离门限之外，仅保留距离门限内的数据； 步骤E，对航迹向空域、跨航向空域、波传播向空域信号沿航迹向、跨航向和波传播向进行三维FFT，得到航迹向、跨航向和波传播向三维波数域信号； 步骤F，将波传播方向、航迹向和跨航向三维波数域信号进行场景中心距离徙动校正和残余距离徙动校正；以及 步骤G，对完成场景中心距离徙动校正和残余距离徙动校正的波传播方向、航迹向和跨航向三维波数域信号沿波传播方向、航迹向和跨航向进行三维IFFT，得到精确重建的感兴趣区域三维图像。2.根据权利要求1所述的方法，其特征在于，所述步骤D包括: 子步骤D1，对航迹向空域、跨航向空域、波传播时域3-D SAR回波数据沿波传播方向进行FFT，变换到航迹向空域、跨航向空域、波传播频域； 子步骤D2，在波传播方向乘以匹配滤波器，沿波传播方向进行IFFT，完成匹配滤波，得到航迹向空域、跨航向空域、波传播时域...

【专利技术属性】
技术研发人员：王彦平，彭学明，谭维贤，洪文，吴一戎，
申请(专利权)人：中国科学院电子学研究所，
类型：发明
国别省市：