CN101937077B - 双基地前视、斜视合成孔径雷达多普勒中心频率的测定方法 - Google Patents
双基地前视、斜视合成孔径雷达多普勒中心频率的测定方法 Download PDFInfo
- Publication number
- CN101937077B CN101937077B CN2010102367820A CN201010236782A CN101937077B CN 101937077 B CN101937077 B CN 101937077B CN 2010102367820 A CN2010102367820 A CN 2010102367820A CN 201010236782 A CN201010236782 A CN 201010236782A CN 101937077 B CN101937077 B CN 101937077B
- Authority
- CN
- China
- Prior art keywords
- entropy
- slope
- waveform entropy
- distance
- waveform
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Abstract
该发明属于合成孔径雷达(SAR)无模糊多普勒中心频率测定方法,包括:距离向脉冲压缩,获取初始斜率,频域、时域转换并获取序列的波形熵,确定最小波形熵,确定无模糊多普勒中心频率。该发明利用合成孔径雷达回波在距离时间、方位时间域的几何特征及移动平台本身的信息,对距离走动轨迹进行反复校正确定最小熵、以衡量其轨迹是否得到很好校正,进而通过所对应轨迹的斜率与所发射信号的波长之比得到无模糊多普勒中心频率,其测定效率较拉登变换提高了五倍以上;因而具有对双基地前视和斜视合成孔径雷达(SAR)多普勒中心频率的测定简捷、准确,处理时间短、效率高,实时性强,可为后续高精度成像、运动补偿提供准确、可靠的数据等特点。
Description
技术领域
本发明属于合成孔径雷达(SAR)参数的测定方法,特别是涉及双基地前视及斜视SAR无模糊多普勒中心频率的测定方法,该方法也适用于单基地斜视合成孔径雷达的无模糊多普勒中心频率的测定。
背景技术
与光学传感器相比,合成孔径雷达(SAR)具有穿透性强,能全天时、全天候工作的独特优点,目前已被广泛应用在地球遥感、资源勘探、侦察、测绘、灾情预报等领域。而与常规侧视SAR相比,前视及斜视SAR在战机对地侦察打击、自主着陆、物资空投、导弹精确末制导等方面具有独特的优势,目前得到了广泛的研究。
多普勒中心频率是SAR高精度处理所需要的重要参数。一方面,其用来构造方位匹配函数或相位补偿函数;另一方面,通过其可以反推运动参数,进行运动误差补偿。而多普勒中心频率误差加大将导致图像信号的信噪比下降,方位模糊性增加,并产生目标位置的平移。从原理上讲,多普勒中心频率可以由载机速度及下视(斜视)角度信息求出,但由于载机(移动)平台往往在速度、姿态(角度)等方面存在不确定性,使得计算出的多普勒参数与实际值差别较大;而从全球定位系统(GPS)与惯性导航系统(INS)得到的航迹数据来计算多普勒中心频率,精度也不能满足要求,因此,利用雷达回波信号来对多普勒中心频率测进行精确测定就具有重大的意义。
相比于常规单基地SAR,双基地前视及斜视SAR多普勒中心频率测定面临更多的难题,比如:天线的收、发分别设置使得相位中心存在同步误差,造成回波方位功率谱不再对称,因而不能利用方位功率谱对称原理对其进行多普勒中心频率测定;前、斜配置造成来自同一散射体的回波跨越多个距离单元,又使得无法直接提取SAR方位向信号进行多普勒中心频率测定,此外还存在严重的多普勒中心频率模糊等等问题。
目前针对单基地SAR,应用比较广泛的方法是相关多普勒中心频率估计器,其本质是利用方位功率谱对称原理进行多普勒中心频率的测定,并且只能实现基带多普勒中心频率的测定,因此无法直接应用于双基地前视及斜视SAR;而在《基于最小熵DCFT的双基地SAR多普勒参数估计》(《现代雷达》2007,29(10):P53-56)的文献中、提出采用基于离散余弦傅立叶变换(DCFT)方法进行双基地SAR多普勒参数测定,由于是直接提取方位向信号进行估计,无法实现多普勒中心频率解模糊,不适用于存在大距离走动(运动)的前视及斜视SAR;在“Comparison of Doppler centroid estimators in bistatic airborneSAR”(International Geoscience and Remote Sensing Symposium(IGARSS),2005,pp.1963-1966,Korea)及“基于Radon变换的双站SAR多普勒参数估计”(《电子与信息学报》,2008,30(6):1331-1335)文献所记载的技术中,采用拉登(radon)变换方法进行双基地SAR多普勒中心频率的测定,此类方法首先对回波进行距离向脉冲压缩,然后选取高对比度区域的信息进行拉登变换,检测出距离走动轨迹的倾角,从而根据多普勒中心频率与距离走动轨迹倾角的关系、测定多普勒中心频率;但此类方法,由于拉登变换涉及到需在二维平面内沿不同的直线进行积分,并且需要确定变换的区域,因而存在测定方法复杂,处理时间长、实时性差,变换区域的过程需人工干预、自动调节能力差等缺陷。
发明内容
本发明的目的是针对背景技术存在的缺陷,研究设计一种双基地前视、斜视合成孔径雷达多普勒中心频率的测定方法,以克服常规多普勒中心频率测定方法无法用于双基地前视和斜视SAR,以及处理时间长、效率低等问题;达到简化双基地前视和斜视合成孔径雷达(SAR)多普勒中心频率的测定、缩短处理时间,有效提高测定的效率和测定的实时性,以及为后续高精度成像、运动补偿提供准确、可靠的数据等目的。
本发明的解决方案是根据双基地前视及斜视SAR具有大距离范围内的走动(运动)轨迹为直线、小距离轨迹为曲线的特点,且其点目标响应的走动轨迹在距离时间、方位时间域主要表现为一条斜直线。根据该几何特点,本发明将移动平台本身的速度、倾斜角度作为基本信息,通过循环(不断)校正距离范围内走动的轨迹、更新其走动轨迹的斜率,直至得到最小熵,并根据最小熵所对应走动轨迹的斜率与发射信号的波长之比、即得到该距离范围(合成孔径时间)内的无模糊多普勒中心频率。因此,本发明方法包括:
A.距离向脉冲压缩:对所获取发射信号的二维回波数据,利用常规匹配滤波方法进行距离向脉冲压缩;
B.获取初始斜率:根据移动平台本身的速度、倾斜角度信息,确定距离走动轨迹的理想斜率,并将该斜率作为初始斜率;
C.频域、时域转换并获取序列的波形熵:根据信号的时移与频移的对应关系及经步骤A压缩后的数据,利用步骤B所得初始斜率或经步骤D更新后的斜率对距离走动轨迹进行校正,然后将经轨迹校正后的数据通过逆傅立叶变换(IFFT)为距离时域、并将表示该距离时域的直线沿方位向投影,最后经计算获得沿方位向投影所得序列的波形熵,并存储所得波形熵;
D.确定最小波形熵:首先判断步骤C所得波形熵是否为第一次获得的波形熵,若是、则将该第一次获得的波形熵所对应的斜率直接更新步骤B所得的初始斜率;否则、与前一循环所存储的波形熵进行大小比较,若小于前一循环所存储的波形熵、则存储该波形熵,同时返回步骤C进行下一迭代循环,直至所得波形熵大于前一循环所得的波形熵止,此时前一循环所得波形熵即为(所测距离范围内的)最小熵;
E.确定无模糊多普勒中心频率:步骤D获得的最小熵所对应走动轨迹的斜率与所发射信号的波长之比,即为该距离范围(合成孔径时间)内的无模糊多普勒中心频率。
所述获得沿方位向投影所得序列的熵,其熵H(x)通过下式获得:
其中:x(n)为沿方位向投影序列,n为距离向采样点数、即n=0,1,…,N-1,
而所述无模糊多普勒中心频率,其中心频率为:
其中:fdc为无模糊多普勒中心频率,λ为载波波长,R(t)为发射站到目标再到接收站的瞬时距离,t为方位时间,t=0为波束中心通过目标的时刻,K为斜率。
发明效果
本发明利用合成孔径雷达回波在距离时间、方位时间域的几何特征及移动平台本身的信息,对距离走动轨迹进行反复校正确定最小熵、以衡量距离走动轨迹是否得到很好校正,进而通过该项最小熵所对应走动轨迹的斜率与所发射信号的波长之比,即得拟测无模糊多普勒中心频率;本发明方法不但可对双基地前视及斜视合成孔径雷达无模糊多普勒中心频率进行测定,且由于利用了移动平台的信息,与拉登变换相比,其测定效率(速度)提高了五倍以上,节省了大量的处理时间;同时,因无须人工确定变换区域,自动处理能力更强;因而,本发明具有对双基地前视和斜视合成孔径雷达(SAR)多普勒中心频率的测定简捷、准确,处理时间短、效率高,实时性强,可为后续高精度成像、运动补偿提供准确、可靠的数据等特点。
附图说明
图1为本发明方法流程示意图(方框图);
图2为本实施方式双基地前视SAR几何配置示意图;
图3为本实施方式取得最小波形熵时所对应的走动轨迹校正结果示意图(坐标图);
图4为走动轨迹校正后往距离向投影所得序列波形示意图(坐标图)。
具体实施方式
图2为本实施方式双基地前视SAR几何配置示意图,
发射机固定,通过雷达发射机发射信号波长为λ=0.03125m的线性调频信号,调频斜率为Kr=8×1012Hz/s;接收机速度为v=100m/s、本实施方式仿真中速度误差设为+0.5m/s,下视角为双基地前视SAR点目标回波信号经相干解调后,表示为:
其中:τ表示距离快时间,变化范围由发射机到目标再到接收机的距离和决定,t为方位时间,变化范围为[0 2]秒,其中t=1秒为波束中心照射目标的时间,R(t)为发射机到目标再到接收机的距离,有:
其中:Rt=201km、Rr=20km分别为合成孔径中心时刻,发射机、接收机到目标的距离;
本实施方式的具体测定方法为:
A.首先对回波进行距离向FFT,然后通过距离向乘以匹配滤波器匹配函数H1(fr),实现距离向脉冲压缩,
其中:fr为距离向频率,变化范围为[-40 40]MHz,B=40MHz为发射信号带宽;
B.根据移动平台本身的速度、倾斜角信息,确定距离走动轨迹的理想斜率K0,并将该斜率作为初始斜率:
代入相应参数后,本实施方式的初始斜率为:K0=-86.60m/s;
C.然后针对距离向脉冲压缩后的数据,在距离频率、方位时间域乘以一个利用距离走动轨迹斜率表示的线性相位因子H2(fr,t),进行距离走动轨迹校正,
其中:c为光速,K为距离走动轨迹斜率(初始斜率为-86.60m/s);
熵:指的是体系的混乱程度,它在控制论、概率论、数论、天体物理、生命科学等领域都有重要应用。本发明用其衡量距离走动轨迹是否得到很好的校正。
将距离走动轨迹校正后的数据沿方位向投影得到序列x(n),n=0,1,…,1023计算序列的波形熵,并存储波形熵:
其中
D.判断是否为第一次获得的波形熵,若是、则将该波形熵所对应的斜率直接更新步骤B所得初始斜率,本实施方式中更新步长为0.01m/s;否则、与所存储的波形熵进行大小比较,若小于所存储的波形熵,则存储该波形熵、更新距离走动斜率,并返回步骤C进行下一迭代循环,本实施方式经过43次迭代,即获得最小波形熵;图3所示,即为得到该最小波形熵时所对应的距离走动轨迹校正结果示意图;图4为距离走动轨迹校正后往距离向的投影所得序列波形示意图;
E.根据步骤D所得最小波形熵所对应走动轨迹的斜率K=-87.03m/s与所发射信号的波长λ=0.03125m之比:
即为本实施方式在该距离范围(合成孔径时间)内的无模糊多普勒中心频率;其误差在理论范围内。
在相同条件下采用本实施方式与拉登变换处理方法相比:本实施方式对无模糊多普勒中心频率的测定仅需1.42秒,而采用拉登变换处理本实施方式则需7.89秒;本实施方式相对于拉登变换方法,效率提高了5.6倍。
Claims (3)
1.一种双基地前视、斜视合成孔径雷达多普勒中心频率的测定方法,包括:
步骤A.距离向脉冲压缩:对所获取发射信号的二维回波数据,利用常规匹配滤波方法进行距离向脉冲压缩;
步骤B.获取初始斜率:根据移动平台本身的速度、倾斜角度信息,确定距离走动轨迹的理想斜率,并将该斜率作为初始斜率;
步骤C.频域、时域转换并获取序列的波形熵:根据信号的时移与频移的对应关系及经步骤A压缩后的数据,利用步骤B所得初始斜率或经步骤D更新后的斜率对距离走动轨迹进行校正,即针对距离向脉冲压缩后的数据,在距离频率、方位时间域乘以一个利用距离走动轨迹斜率表示的线性相位因子H2(fr,t),进行距离走动轨迹校正,
其中:c为光速,K为距离走动轨迹斜率;
然后将经轨迹校正后的数据通过逆傅立叶变换为距离时域、并将表示该距离时域的直线沿方位向投影,最后经计算获得沿方位向投影所得序列的波形熵,并存储所得波形熵;
步骤D.确定最小波形熵:首先判断步骤C所得波形熵是否为第一次获得的波形熵,若是、则将该第一次获得的波形熵所对应的斜率直接更新步骤B所得的初始斜率;否则、与前一循环所存储的波形熵进行大小比较,若小于前一循环所存储的波形熵、则存储该波形熵,同时返回步骤C进行下一迭代循环,直至所得波形熵大于前一循环所得的波形熵止,此时前一循环所得波形熵即为最小熵;
步骤E.确定无模糊多普勒中心频率:步骤D获得的最小熵所对应走动轨迹的斜率与所发射信号的波长之比,即为该距离范围内的无模糊多普勒中心频率。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2010102367820A CN101937077B (zh) | 2010-07-23 | 2010-07-23 | 双基地前视、斜视合成孔径雷达多普勒中心频率的测定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2010102367820A CN101937077B (zh) | 2010-07-23 | 2010-07-23 | 双基地前视、斜视合成孔径雷达多普勒中心频率的测定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101937077A CN101937077A (zh) | 2011-01-05 |
CN101937077B true CN101937077B (zh) | 2012-11-07 |
Family
ID=43390507
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2010102367820A Active CN101937077B (zh) | 2010-07-23 | 2010-07-23 | 双基地前视、斜视合成孔径雷达多普勒中心频率的测定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101937077B (zh) |
Families Citing this family (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102798842B (zh) * | 2012-07-16 | 2014-03-12 | 西安电子科技大学 | 前视阵列雷达的俯仰-慢时间空时二维自适应处理方法 |
CN102937709B (zh) * | 2012-10-16 | 2014-05-07 | 北京理工大学 | 基于直达波的双基地sar无同步回波聚焦成像方法 |
CN102967858B (zh) * | 2012-11-14 | 2014-03-05 | 电子科技大学 | 雷达前视超分辨成像方法 |
CN104459628A (zh) * | 2013-09-21 | 2015-03-25 | 杨瑞娟 | 准正交频分复用多载波线性调频雷达信号设计及处理方法 |
CN103823210B (zh) * | 2014-03-12 | 2016-08-31 | 北京理工大学 | 一种非合作式星地双基地sar时频同步方法 |
CN103869318B (zh) * | 2014-04-08 | 2016-02-17 | 电子科技大学 | 双基地前视合成孔径雷达动目标速度误差函数构造方法 |
CN105158745A (zh) * | 2015-08-05 | 2015-12-16 | 电子科技大学 | 移变双基前视合成孔径雷达距离徙动校正方法 |
CN107589414B (zh) * | 2017-09-07 | 2021-01-26 | 电子科技大学 | 基于相位中心点追踪的斜前视sar多普勒质心估计方法 |
CN109507667B (zh) * | 2018-12-21 | 2021-12-28 | 曲卫 | 基于导航卫星信号的双站sar接收机运动误差补偿方法 |
CN109655795B (zh) * | 2019-01-24 | 2020-06-30 | 南京莱斯电子设备有限公司 | 基于波形熵的米波雷达同频窄脉冲干扰抑制方法和系统 |
CN110488271A (zh) * | 2019-08-01 | 2019-11-22 | 中国人民解放军陆军工程大学 | 一种一维距离像速度补偿方法 |
-
2010
- 2010-07-23 CN CN2010102367820A patent/CN101937077B/zh active Active
Non-Patent Citations (6)
Title |
---|
卢晶琦 等.发射机固定的斜视双基地SAR的NLCS算法.《雷达科学与技术》.2006,第4卷(第1期),27-30. * |
吕睿 等.基于最小熵DCFT的双基地SAR多普勒参数估计.《现代雷达》.2007,第29卷(第10期),53-56. * |
朱振波 等.基于Radon变换的双站SAR多普勒参数估计.《电子与信息学报》.2008,第30卷(第6期),1331-1335. * |
武俊杰 等.双基地斜视SAR距离徙动成像算法.《电子科技大学学报》.2007,第36卷(第5期),1106-1109. * |
王霖郁 等.改进的机载SAR多普勒中心估计方法.《测绘学报》.2010,第39卷(第3期),271-275. * |
陈云 等.大斜视角下双基地合成孔径雷达成像算法研究.《信息通信》.2008,(第6期),11-13. * |
Also Published As
Publication number | Publication date |
---|---|
CN101937077A (zh) | 2011-01-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101937077B (zh) | 双基地前视、斜视合成孔径雷达多普勒中心频率的测定方法 | |
US7298325B2 (en) | Technique for accurate estimate of large antenna inertial two dimensional orientation using relative GPS spatial phase | |
US7038612B2 (en) | Method for SAR processing without INS data | |
CN105759263B (zh) | 一种高分辨率大场景下的星载斜视sar雷达成像方法 | |
CN102707269B (zh) | 一种机载雷达距离走动校正方法 | |
US7301495B2 (en) | Interrupt SAR implementation for range migration (RMA) processing | |
CN101887122B (zh) | 可消除地面高程误差的星载sar图像目标定位方法 | |
CN110488292B (zh) | 一种基于卫星编队的遥感系统 | |
CN113439220A (zh) | 一种基于卫星编队的遥感系统及星座系统 | |
CN102288962B (zh) | 一种超宽带合成孔径雷达实时运动补偿方法 | |
CN102778681A (zh) | 固定发射站双基地前视合成孔径雷达成像方法 | |
CN110456350B (zh) | 一种星载sar星座系统 | |
CN101915919A (zh) | 利用北斗卫星信号的双基sar成像系统 | |
CN110389366B (zh) | 一种基于多源sar卫星的海上目标运动估计方法 | |
CN113050134A (zh) | 基于卫星导航信息的海面风场反演观测方法 | |
Gray et al. | A revised calibration of the interferometric mode of the CryoSat-2 radar altimeter improves ice height and height change measurements in western Greenland | |
JP3987451B2 (ja) | 合成開口レーダを用いる地表面変動量計測方法 | |
Baumgartner et al. | Large along-track baseline SAR-GMTI: First results with the TerraSAR-X/TanDEM-X satellite constellation | |
CN102141611B (zh) | 斜视合成孔径雷达多普勒模糊数快速测定方法 | |
CN108828585B (zh) | 一种空间平台目标无伺服跟踪测角方法 | |
Zhang et al. | Error analysis of dynamic sea surface height measurement by near-nadir interferometric SAR | |
Rouveure et al. | Terrestrial mobile mapping based on a microwave radar sensor. Application to the localization of mobile robots | |
Mao et al. | Doppler centroid estimation for Doppler beam sharpening imaging based on the morphological edge detection method | |
Zhang et al. | Study of accurate ocean-altimetry with GNSS-R | |
Jiang et al. | Terrain observation using ground-based interferometric millimeter-wave sar imaging |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |