专利名称：斜视合成孔径雷达多普勒模糊数快速测定方法
技术领域：
本发明涉及合成孔径雷达(Synthetic Aperture Radar，SAR)技术领域，特别是涉 及斜视SAR的参数测定技术领域。
背景技术：
合成孔径雷达(SAR)具有穿透性、全天时和全天候工作能力的独特优点，目前已 被广泛应用在地球遥感、资源勘探、侦察、测绘、灾情预报等领域。多普勒(Doppler)中心频率是SAR高精度成像处理所需要的重要参数。一方面， 其用来构造成像所需的方位匹配函数或相位补偿函数；另一方面，通过其可以反推载机的 运动参数，对载机进行运动误差补偿。而多普勒中心频率的误差在成像时将导致图像信号 信噪比下降，方位模糊性增加和成像目标的位置平移。从原理上讲，多普勒中心频率可以由 从全球定位系统(GPS)与惯性导航系统(INQ提供的载机速度和载机速度方向的法线与波 束发射方向的夹角(简称载机斜视角度)信息求出，但由于载机在速度、姿态等方面存在 不确定性，使得计算出的多普勒中心频率与真实值差别较大，精度也通常不能满足要求，因 此，利用SAR的回波数据对多普勒中心频率进行精确测定具有重要的意义。斜视SAR是SAR中的一种，对于斜视SAR而言，由于来自同一散射体(成像目标) 的回波跨越多个距离单元，使得无法直接提取回波的方位向的信号进行多普勒中心频率测 定，并且由于斜视SAR回波的方位向按脉冲重复频率采样，多普勒中心频率存在严重的模 糊。多普勒中心频率由多普勒模糊数与基带多普勒中心频率两部分组成。因此，多普勒模 糊数的测定对于多普勒中心频率的精确测定至关重要。文献《SAR多普勒中心快速解模糊的新方法》(《信号处理》，2005，21 (4A) =503-505) 中，首先估计不同距离单元上的模糊的多普勒中心频率，然后进行线性拟合得到模糊多普 勒中心频率随载机到成像目标的距离(简称载机斜距)变化的斜率，实现多普勒模糊数 测定；但该方法使用传统相关多普勒估计方法获得模糊多普勒中心频率，受部分孔径照射 效应影响大，仅适用于低对比度场景。文献“Adding sensitivity to the MLBF doppler centroid estimator"(IEEE Transactions on Geoscience and Remote Sensing,vol. 45, no. 2，pp. 279492，2007)，采用多视技术进行多普勒模糊数测定，但由于多视处理涉及到额 外的傅里叶变换及逆傅里叶变换产生距离压缩的多视数据，运算量较大。文献《基于Radon 变换的大前斜视SAR多普勒参数估计方法》(信号处理，2009，25 (0 . :210_215。)中，采 用拉登(radon)变换方法进行多普勒中心频率的解模糊，此方法首先对回波进行距离向脉 冲压缩，然后选取高对比度区域数据进行拉登变换，检测出成像目标的距离走动轨迹的倾 角，从而根据多普勒中心频率与成像目标的距离走动轨迹的倾角的关系，测定多普勒中心 频率，但由于多普勒中心频率对成像目标的距离走动轨迹的倾角的检测精度非常敏感，致 使测定的精度差，稳定性差，并且由于拉登变换需在二维平面内沿不同的直线进行积分，同 样存在运算量大的缺点。

发明内容
本发明的目的是为了克服现有技术中存在的计算多普勒中心频率运算量大的不 足，提出了斜视合成孔径雷达多普勒模糊数快速测定方法。本发明的技术方案是斜视合成孔径雷达多普勒模糊数快速测定方法，包括如下 步骤A.获取成像目标的回波数据，并利用常规匹配滤波方法进行距离向脉冲压缩获得 距离压缩后的回波数据；B.针对距离压缩后的回波数据，选取对比度高的的区域，利用Carmy边缘检测算 子进行边缘检测，得到含有线性特征的二值数据；C.根据惯导设备或提供的载机速度、载机斜视角度信息，得到成像目标的距离走 动轨迹的倾角的粗略值，并根据其粗略值确定二值数据的Radon变换区域和Radon变换步 长；D.根据步骤C设定的Radon变换区域和Radon变换步长，对步骤B得到的二值数 据进行Radon变换，得到成像目标的距离走动轨迹的倾角的精确值；E.根据Radon变换得到的成像目标的距离走动轨迹的倾角的精确值，结合斜视 SAR的发射信号波长、距离向采样频率及脉冲重复频率，计算得到多普勒中心频率的粗略 值；F.根据测得的多普勒中心频率的粗略值，结合斜视SAR的脉冲重复频率，计算得 到多普勒模糊数。本发明的有益效果通过引入边缘检测技术，将灰度数据转化为了二值数据，因二 值数据具有稀疏特性，极大的减少了 Radon变换的运算时间；另外通过利用惯导设备提供 的载机速度、载机斜视角度信息，粗略确定Radon变换的区域，避免了 Radon变换在整个二 维平面进行搜索，从而能够快速高效的测定成像目标的距离走动轨迹的倾角，进而测定出 多普勒中心频率的粗略值，并结合脉冲重复频率，快速实现多普勒模糊数测定。与直接采用 拉登变换方法相比，其测定效率(速度)得到了明显提高，节省了大量的处理时间，并且相 比于直接测量多普勒中心频率，对多普勒模糊数的测定更加稳健、可靠。所以本发明具有测 定简捷、准确、稳健，处理时间短、效率高，实时性强等优点。


图1是多普勒中心频率模糊数测定流程图。图2是斜视SAR几何配置示意图。图3为本发明成像目标距离压缩后的回波数据的仿真图。图4为本发明成像目标距离压缩后的回波数据的边缘检测结果仿真图。图5为本发明回波数据的边缘检测结果的拉登变换检测结果仿真图。
具体实施例方式下面结合附图和具体实施例对本发明做进一步的说明。下面结合附图和具体实施 例对本发明做进一步的说明。本实施例中斜视SAR的几何配置如图2所示，在一个xyz空间坐标系中，0表示坐标系的原点，P表示成像目标，载机速度V = 200m/s，载机斜视角度θ =59° ；斜视SAR的 发射信号波长为λ = 0.03125m，发射信号调频斜率为& = 8X1012Hz/s，斜视SAR的成像 目标的回波信号经相干解调后的回波信号表示为^(i,r) = exp|-7^i ·exp|jnKr r-^-R{t) J(1)其中τ表示回波信号的距离向时间，变化范围为[2. 70Χ 10-4 2. 73Χ10-4]秒，t 为回波信号的方位向时间，变化范围为W2]秒，其中t= 1秒为波束中心照射目标的时间， exp W表示以e为底的自然函数，j表示虚部，c表示光速，R(t)为斜视SAR天线到成像目
标的距离，其具体公式为
_ 2 2 β
R{t) = ^J Ιξ + (Vtf - IRfiVt COS0 ^ R0-Vtsme+ Vt2(2)
2R0其中=Rtl表示波束中心照射成像目标时，成像目标与斜视SAR天线之间的距离为 22.824km。斜视合成孔径雷达多普勒模糊数快速测定方法，如图1所示，包括如下步骤A.获取成像目标的回波数据，并利用常规匹配滤波方法进行距离向脉冲压缩获得 距离压缩后的回波数据；首先对回波数据进行距离向FFT(傅立叶变换)，然后通过距离向乘以常规的匹配 滤波器匹配函数H1(L)，获得距离压缩后的回波数据，匹配函数H1(L)的具体公式如下
Hl{fr) = rect{^M-j^^)(3)
BKr其中f；为距离向频率，变化范围为[_4040]MHz，B为斜视合成孔径雷达的发射信 号带宽，本式中取B = 40MHz, rect(*)和exp (*)分别表示矩形函数和以e为底的指数函 数，由于这两个函数是本技术领域中的公知常识，因此不再详细介绍。然后对距离向信号进行IFFT(逆傅立叶变换)，得到距离压缩后的回波数据的信 号时域表达式Sr{t,r) =sinc\^B τ-|(4)式中，sine (*)表示sine函数，由于该函数是本技术领域中的公知常识，因此不再 详细介绍。如图3所示，横坐标表示一定时间范围内载机运动时回波数据的方位向采样点 数，纵坐标表示距离向采样点数，成像目标选取5个点，从图中可以看出，成像目标的每个 点的距离压缩后的回波数据呈现的距离走动轨迹表现为倾斜的直线。B.针对距离压缩后的回波数据，选取对比度高的的区域，利用图像处理中的 Canny边缘检测算子进行边缘检测，得到含有线性特征的二值数据；如图4所示，横坐标表示一定时间范围内载机运动时回波数据的方位向采样点 数，纵坐标表示距离向采样点数，成像目标选取5个点，从图中可以看出，经过边缘检测后 的回波数据，仍然保留了数据中的线性特征，同时距离压缩后的回波灰度数据转化成了二值数据。C.根据惯导设备(如GPS或INQ或提供的载机速度、载机斜视角度θ信息， 得到成像目标的距离走动轨迹的倾角α的粗略值，并根据其粗略值确定二值数据的拉登 (Radon)变换区域和拉登(Radon)变换步长；利用惯导设备提供的载机速度、载机斜视角度信息，本实施例中速度ν = 200m/s, 载机斜视角度θ =59°，波长λ = 0.03125m，依据多普勒中心频率与载机斜视角度的倾 角存在如下关系，
权利要求
1.斜视合成孔径雷达多普勒模糊数快速测定方法，包括如下步骤A.获取成像目标的回波数据，并利用常规匹配滤波方法进行距离向脉冲压缩获得距离 压缩后的回波数据；B.针对距离压缩后的回波数据，选取对比度高的的区域，利用Carmy边缘检测算子进 行边缘检测，得到含有线性特征的二值数据；C.根据惯导设备或提供的载机速度、载机斜视角度信息，得到成像目标的距离走动轨 迹的倾角的粗略值，并根据其粗略值确定二值数据的Radon变换区域和Radon变换步长；D.根据步骤C设定的Radon变换区域和Radon变换步长，对步骤B得到的二值数据进 行Radon变换，得到成像目标的距离走动轨迹的倾角的精确值；E.根据Radon变换得到的成像目标的距离走动轨迹的倾角的精确值，结合斜视SAR的 发射信号波长、距离向采样频率及脉冲重复频率，计算得到多普勒中心频率的粗略值；F.根据测得的多普勒中心频率的粗略值，结合斜视SAR的脉冲重复频率，计算得到多 普勒模糊数。
2.根据权利要求1所述的斜视合成孔径雷达多普勒模糊数快速测定方法，其 特征在于，上述步骤C中成像目标的距离走动轨迹的倾角的粗略值的计算公式为‘_ t (2ν· fs-sm0“arC c-pRF Xl80° if (arctan(^^))<0 πc-PRFa = <ο,2v- />sin<9、π - arctani-—-)r, s - n-c'PRF xl80° if (arctan(2v-/--Sm^))>0、πc · PRF
全文摘要
本发明涉及斜视合成孔径雷达多普勒模糊数快速测定方法，包括步骤A.获取成像目标的回波数据，并利用常规匹配滤波方法进行距离向脉冲压缩获得距离压缩后的回波数据；B.针对距离压缩后的回波数据，选取对比度高的的区域，利用Canny边缘检测算子进行边缘检测，得到含有线性特征的二值数据；C.根据惯导设备或提供的载机速度、载机斜视角度信息，得到成像目标的距离走动轨迹的倾角的粗略值，并根据其粗略值确定二值数据的Radon变换区域和Radon变换步长。本发明的有益效果通过引入边缘检测技术，极大的减少了Radon变换的运算时间，使本发明具有测定简捷、准确、稳健，处理时间短、效率高，实时性强等优点。
文档编号G01S7/41GK102141611SQ201010611139
公开日2011年8月3日 申请日期2010年12月29日 优先权日2010年12月29日
发明者孔令讲, 张晓玲, 李文超, 杨建宇, 杨晓波, 杨海光, 武俊杰, 黄钰林 申请人:电子科技大学