CN108983188B - 一种基于干涉成像高度计的雷达回波的搜索跟踪方法 - Google Patents

一种基于干涉成像高度计的雷达回波的搜索跟踪方法 Download PDF

Info

Publication number
CN108983188B
CN108983188B CN201810778659.8A CN201810778659A CN108983188B CN 108983188 B CN108983188 B CN 108983188B CN 201810778659 A CN201810778659 A CN 201810778659A CN 108983188 B CN108983188 B CN 108983188B
Authority
CN
China
Prior art keywords
echo
search
interference imaging
noise
subpattern
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.)
Expired - Fee Related
Application number
CN201810778659.8A
Other languages
English (en)
Other versions
CN108983188A (zh
Inventor
石晓进
翟文帅
杨杰芳
张云华
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
National Space Science Center of CAS
Original Assignee
National Space Science Center of CAS
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by National Space Science Center of CAS filed Critical National Space Science Center of CAS
Priority to CN201810778659.8A priority Critical patent/CN108983188B/zh
Publication of CN108983188A publication Critical patent/CN108983188A/zh
Application granted granted Critical
Publication of CN108983188B publication Critical patent/CN108983188B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/41Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S13/00Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/88Radar or analogous systems specially adapted for specific applications
    • G01S13/882Radar or analogous systems specially adapted for specific applications for altimeters
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S13/00Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/88Radar or analogous systems specially adapted for specific applications
    • G01S13/89Radar or analogous systems specially adapted for specific applications for mapping or imaging
    • G01S13/90Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
    • G01S13/9021SAR image post-processing techniques
    • G01S13/9023SAR image post-processing techniques combined with interferometric techniques

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Physics & Mathematics (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • General Physics & Mathematics (AREA)
  • Electromagnetism (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明公开了一种基于干涉成像高度计的雷达回波的搜索跟踪方法,所述方法包括:步骤1)干涉成像高度计进入噪声测量子模式:在不同的AGC下采集***噪声,计算所有AGC下的***噪声等效幅度;步骤2)干涉成像高度计进入回波搜索子模式:设定采样范围,采样窗在采样范围内移动进行回波搜索;步骤3)在搜索到雷达回波后,干涉成像高度计转入脉冲簇测量子模式:采用高脉冲重复频率对采样窗进行调整,然后进行回波跟踪。本发明的方法通过设计三个子工作模式,实现了干涉成像高度计***噪声测量、大范围回波搜索和回波精细跟踪的功能,满足了干涉成像高度计在高脉冲重复频率下准确搜索跟踪雷达回波信号的要求。

Description

一种基于干涉成像高度计的雷达回波的搜索跟踪方法
技术领域
本发明涉及雷达高度计技术领域,尤其涉及一种基于干涉成像高度计的雷达回波的搜索跟踪方法。
背景技术
干涉成像高度计是一种新型雷达高度计,采用了短基线、小角度干涉测高的新技术体制。相对于传统雷达高度计,干涉成像高度计采用合成孔径技术和新型的高度跟踪技术实现海面高度的高精度测量,并具有宽刈幅和可成像的优点。但是由于干涉成像高度计需要进行方位向合成孔径,所以干涉成像高度计的脉冲重复频率远高于传统雷达高度计。为了满足干涉成像高度计在高脉冲重复频率下正常工作和稳定搜索回波信号的需求,需要根据其具体特点进行工作模式设计。目前,天宫二号所搭载的三维成像微波高度计在国际上首次验证了干涉成像高度计的工作机理。
发明内容
本发明的目的在于,针对干涉成像高度计的工作要求和特点,提供一种基于干涉成像高度计的雷达回波的搜索跟踪方法,可以满足干涉成像高度计在高脉冲重复频率下捕获跟踪雷达回波的要求。
为了实现上述目的,本发明提供了一种基于干涉成像高度计的雷达回波的搜索跟踪方法,所述方法包括:
步骤1)干涉成像高度计进入噪声测量子模式:在不同的AGC下采集***噪声,计算所有AGC下的***噪声等效幅度;
步骤2)干涉成像高度计进入回波搜索子模式:设定采样范围,采样窗在采样范围内移动进行回波搜索;
步骤3)在搜索到雷达回波后,干涉成像高度计转入脉冲簇测量子模式:采用高脉冲重复频率对采样窗进行调整,然后进行回波跟踪。
作为上述方法的一种改进,所述步骤1)具体为:所述干涉成像高度计不发射雷达信号,在不同的AGC下采集***噪声,则噪声的等效幅度Anoise为:
其中,Anoise为噪声等效幅度,N为采样窗内采样点数,xn为第n个采样点的幅度值;记录所有AGC下的***噪声等效幅度。
作为上述方法的一种改进,所述步骤2)进一步包括:
步骤2-1)根据卫星位置计算雷达回波时延初始搜索值Td
其中,X,Y,Z为卫星位置,c为光速,RE为地球半径;
步骤2-2)设定回波搜索范围ΔT,采样窗在(Td-ΔT,Td+ΔT)范围内移动并搜索回波;
步骤2-3)计算采样窗内回波数据的等效幅度A,并与相同AGC下的Anoise进行比较,若A>S·Anoise,S为自定义因子,则判断搜索到回波,并记录此时的雷达回波时延T;
步骤2-4)当连续2次搜索到回波,转入步骤3),否则,转入步骤2-2)。
作为上述方法的一种改进,所述回波搜索范围ΔT满足:ΔT<Ts/2;其中,Ts为回波搜索子模式脉冲重复周期。
作为上述方法的一种改进,所述步骤3)进一步包括:
步骤3-1)在搜索到雷达回波后,干涉成像高度计转入脉冲簇测量子模式;计算雷达发射与接收信号之间的间隔脉冲数Nburst和脉冲簇测量子模式子脉冲重复周期PRTb为:
Nburst=floor(T/Tbmin) (3)
PRTb=T/Nburst (4)
其中,floor代表向下取整;T为步骤2)中搜索到回波时的雷达回波时延,Tbmin为脉冲簇测量子模式下最小脉冲簇重复周期;
步骤3-2)根据计算的Nburst和PRTb参数,调整干涉成像高度计在脉冲簇测量子模式下的脉冲重复频率和发射接收信号间延时,对回波进行搜索;
步骤3-3)计算采样窗内回波重心G,调整采样窗位置,使回波处于采样窗重心处,对回波进行跟踪;
步骤3-4)计算采样窗内回波数据的等效幅度A,并与相同AGC下的Anoise进行比较,若A<S·Anoise,则判断未搜索到回波;
步骤3-5)连续2次未搜索到回波,则判断回波位置不正确,转入步骤2)。
本发明的优点在于:
本发明通过设计三个子工作模式:噪声测量子模式,回波搜索子模式和脉冲簇测量子模式,实现了干涉成像高度计***噪声测量、大范围回波搜索和回波精细跟踪的功能,满足了干涉成像高度计在高脉冲重复频率下准确搜索跟踪雷达回波信号的要求。
附图说明
图1为本发明的基于干涉成像高度计的雷达回波的搜索跟踪方法的流程图;
图2为本发明的噪声测量子模式时序示意图;
图3为本发明的回波搜索子模式时序示意图;
图4为本发明的脉冲簇测量子模式时序示意图。
具体实施方式
现结合附图对本发明作进一步的描述。
如图1所示,本发明提供了一种基于干涉成像高度计的雷达回波的搜索跟踪方法。
本发明所设计的干涉成像高度计工作模式可描述为由噪声测量子模式,回波搜索子模式,脉冲簇测量子模式组成的工作模式,具体包括:
步骤1):干涉成像高度计首先进入采噪声子模式,不发射雷达信号,只采集***噪声,如图2所示。由于***噪声幅度与接收机AGC(程控增益)有关,所以需要计算并记录所有AGC下的噪声等效幅度Anoise。等效幅度计算公式如下式所示。
其中A为等效幅度,N为采样窗内采样点数,xn为第n个采样点的幅度值。
步骤2):在记录完所有AGC下的***噪声等效幅度之后,干涉成像高度计进入回波搜索子模式,如图3所示。首先需要根据GPS参数中的卫星位置计算雷达回波时延Td,表达式为:
其中X,Y,Z为GPS获得的卫星位置,c为光速,RE为地球半径。
步骤3):设定回波搜索范围ΔT,采样窗在(Td-ΔT,Td+ΔT)范围内移动并搜索回波,若回波搜索子模式脉冲重复周期为Ts,则要求Ts>2ΔT。
步骤4):计算采样窗内回波数据的等效幅度A,并与相同AGC下的Anoise进行比较,若A>S·Anoise,则判断搜索到回波。其中S为自定义因子,代表对***信噪比的预期,即采样窗内回波等效幅度大于***噪声等效幅度S倍时,判断此时的回波为有效信号。
步骤5):由于雷达回波信号强度会受到地表地形,云雨等外界条件的影响,回波信号不一定保持稳定。若出现在(Td-ΔT,Td+ΔT)范围内未搜索到回波的异常情况,应重新进行搜索,只有连续2次搜索到回波,才能判断稳定搜索到回波并转入脉冲簇测量子模式。
步骤6):在转入脉冲簇测量子模式后,需要根据雷达回波时延,计算脉冲簇测量子模式的间隔脉冲数和子脉冲重复周期,如图4所示。设回波搜索子模式下,稳定搜索到回波时,雷达回波时延为T,脉冲簇测量子模式下最小脉冲簇重复周期为Tbmin,则雷达发射与接收信号之间的间隔脉冲数Nburst和脉冲簇测量子模式子脉冲重复周期PRTb可以表示为:
Nburst=floor(T/Tbmin) (3)
PRTb=T/Nburst (4)
其中floor代表向下取整。
根据计算的Nburst和PRTb参数,调整干涉成像高度计在脉冲簇测量子模式下的脉冲重复频率和发射接收信号间延时,保证以高脉冲重复频率搜索到雷达回波;
步骤7):计算采样窗内回波重心G,微调采样窗位置,使回波处于采样窗中心处。回波重心G的公式如下:
其中N为采样窗内采样点数,xn为第n个采样点的幅度值。
步骤8):计算采样窗内回波数据的等效幅度A,并与相同AGC下的Anoise进行比较,若A<S·Anoise,则判断未搜索到回波,若出现异常,连续2次未搜索到回波,则判断回波位置不正确,应转入回波搜索子模式,重新搜索回波。
本发明已成功应用在天宫二号干涉成像高度计项目中。在该项目中,噪声测量子模式对***128个AGC值进行逐一测量,并记录各AGC下的噪声等效幅度。回波搜索子模式脉冲重复频率为1000Hz-1300Hz,采样窗可以在500us的大范围内进行回波搜索。脉冲簇子模式脉冲重复频率为15500Hz-16000Hz,采样窗在2us的范围内进行微调,对回波进行跟踪。在出现异常情况,比如目标后向散射系数较低,回波信噪比低,导致未能正常搜索或跟踪上回波时,***可以正确进行子工作模式切换,重新进入回波搜索子模式搜索回波。
最后所应说明的是,以上实施例仅用以说明本发明的技术方案而非限制。尽管参照实施例对本发明进行了详细说明,本领域的普通技术人员应当理解,对本发明的技术方案进行修改或者等同替换,都不脱离本发明技术方案的精神和范围,其均应涵盖在本发明的权利要求范围当中。

Claims (5)

1.一种基于干涉成像高度计的雷达回波的搜索跟踪方法,所述方法包括:
步骤1)干涉成像高度计进入噪声测量子模式:在不同的AGC下采集***噪声,计算所有AGC下的***噪声等效幅度;
步骤2)干涉成像高度计进入回波搜索子模式:设定采样范围,采样窗在采样范围内移动进行回波搜索;
步骤3)在搜索到雷达回波后,干涉成像高度计转入脉冲簇测量子模式:采用高脉冲重复频率对采样窗进行调整,然后进行回波跟踪。
2.根据权利要求1所述的基于干涉成像高度计的雷达回波的搜索跟踪方法,其特征在于,所述步骤1)具体为:所述干涉成像高度计不发射雷达信号,在不同的AGC下采集***噪声,则噪声的等效幅度Anoise为:
其中,Anoise为噪声等效幅度,N为采样窗内采样点数,xn为第n个采样点的幅度值;记录所有AGC下的***噪声等效幅度。
3.根据权利要求2所述的基于干涉成像高度计的雷达回波的搜索跟踪方法,其特征在于,所述步骤2)进一步包括:
步骤2-1)根据卫星位置计算雷达回波时延初始搜索值Td
其中,(X,Y,Z)为卫星位置,c为光速,RE为地球半径;
步骤2-2)设定回波搜索范围ΔT,采样窗在(Td-ΔT,Td+ΔT)范围内移动并搜索回波;
步骤2-3)计算采样窗内回波数据的等效幅度A,并与相同AGC下的Anoise进行比较,若A>S·Anoise,S为自定义因子,则判断搜索到回波,并记录此时的雷达回波时延T;
步骤2-4)当连续2次搜索到回波,转入步骤3),否则,转入步骤2-2)。
4.根据权利要求3所述的基于干涉成像高度计的雷达回波的搜索跟踪方法,其特征在于,所述回波搜索范围ΔT满足:ΔT<Ts/2;其中,Ts为回波搜索子模式脉冲重复周期。
5.根据权利要求4所述的基于干涉成像高度计的雷达回波的搜索跟踪方法,其特征在于,所述步骤3)进一步包括:
步骤3-1)在搜索到雷达回波后,干涉成像高度计转入脉冲簇测量子模式;计算雷达发射与接收信号之间的间隔脉冲数Nburst和脉冲簇测量子模式子脉冲重复周期PRTb为:
Nburst=floor(T/Tbmin) (3)
PRTb=T/Nburst (4)
其中,floor代表向下取整;T为步骤2)中搜索到回波时的雷达回波时延,Tbmin为脉冲簇测量子模式下最小脉冲簇重复周期;
步骤3-2)根据计算的Nburst和PRTb参数,调整干涉成像高度计在脉冲簇测量子模式下的脉冲重复频率和发射接收信号间延时,对回波进行搜索;
步骤3-3)计算采样窗内回波重心G,调整采样窗位置,使回波处于采样窗重心处,对回波进行跟踪;
步骤3-4)计算采样窗内回波数据的等效幅度A,并与相同AGC下的Anoise进行比较,若A<S·Anoise,则判断未搜索到回波;
步骤3-5)连续2次未搜索到回波,则判断回波位置不正确,转入步骤2)。
CN201810778659.8A 2018-07-16 2018-07-16 一种基于干涉成像高度计的雷达回波的搜索跟踪方法 Expired - Fee Related CN108983188B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810778659.8A CN108983188B (zh) 2018-07-16 2018-07-16 一种基于干涉成像高度计的雷达回波的搜索跟踪方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810778659.8A CN108983188B (zh) 2018-07-16 2018-07-16 一种基于干涉成像高度计的雷达回波的搜索跟踪方法

Publications (2)

Publication Number Publication Date
CN108983188A CN108983188A (zh) 2018-12-11
CN108983188B true CN108983188B (zh) 2019-09-24

Family

ID=64549404

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810778659.8A Expired - Fee Related CN108983188B (zh) 2018-07-16 2018-07-16 一种基于干涉成像高度计的雷达回波的搜索跟踪方法

Country Status (1)

Country Link
CN (1) CN108983188B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112799056B (zh) * 2020-12-28 2022-12-13 上海卫星工程研究所 一种星载雷达高度计***及方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR2886022A1 (fr) * 2005-05-18 2006-11-24 Agence Spatiale Europeenne Procede pour l'etablissement d'images cartographiques de vecteurs de vitesses de courants marins de surface et systeme de radar altimetrique mettant en oeuvre le procede
EP1631833B1 (en) * 2003-06-11 2009-04-29 Honeywell International Inc. Methods and apparatus for calibrating the antenna-orientation in an interferometric radar altimeter
EP3056922A2 (en) * 2015-02-11 2016-08-17 Honeywell International Inc. Velocity and attitude estimation using an interferometric radar altimeter
CN106483524A (zh) * 2016-09-21 2017-03-08 中国科学院国家空间科学中心 基于全极化干涉合成孔径雷达数据的建筑物提取方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1631833B1 (en) * 2003-06-11 2009-04-29 Honeywell International Inc. Methods and apparatus for calibrating the antenna-orientation in an interferometric radar altimeter
FR2886022A1 (fr) * 2005-05-18 2006-11-24 Agence Spatiale Europeenne Procede pour l'etablissement d'images cartographiques de vecteurs de vitesses de courants marins de surface et systeme de radar altimetrique mettant en oeuvre le procede
EP3056922A2 (en) * 2015-02-11 2016-08-17 Honeywell International Inc. Velocity and attitude estimation using an interferometric radar altimeter
CN106483524A (zh) * 2016-09-21 2017-03-08 中国科学院国家空间科学中心 基于全极化干涉合成孔径雷达数据的建筑物提取方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
Demonstration of ocean target detection by Tiangong-2 interferometric imaging radar altimeter;Yunhua Zhang et al.;《2018 22nd International Microwave and Radar Conference (MIKON)》;20180517;第261-264页 *
Evaluation of BAQ on Tiangong-2 interferometric imaging radar altimeter data compression;Xiaojin Shi et al.;《2018 22nd International Microwave and Radar Conference (MIKON)》;20180517;第623-624页 *
HY-2A卫星雷达高度计设计及其;许可 等;《中国工程科学》;20130715;第15卷(第7期);第25-32页 *

Also Published As

Publication number Publication date
CN108983188A (zh) 2018-12-11

Similar Documents

Publication Publication Date Title
US6188348B1 (en) Radar ice sounder with parallel doppler processing
Gogineni et al. An improved coherent radar depth sounder
US5113193A (en) Autonomous synchronization of a bistatic synthetic aperture radar (SAR) system
US7417578B1 (en) Removal of spurious aircraft detections on weather radar
EP3432026A1 (en) Systems and methods for measuring velocity and acceleration with a radar altimeter
JP5502262B2 (ja) 測地距離データを得る方法及び装置
CN108535730A (zh) 一种多普勒气象雷达解速度模糊方法和***
CN109581362A (zh) 合成孔径雷达高度计在可变脉冲簇模式下的信号处理方法
CN108415013A (zh) 基于调频斜率精搜索的低信噪比信号参数提取方法
CN108983188B (zh) 一种基于干涉成像高度计的雷达回波的搜索跟踪方法
Roca et al. The EnviSat RA-2 instrument design and tracking performance
Christensen et al. A low-cost glacier-mapping system
CN100356191C (zh) 一种适用于低对比度场景合成孔径雷达成像的自聚焦方法
Kaimal et al. An acoustic Doppler sounder for measuring wind profiles in the lower boundary layer
US7149148B2 (en) Localization of high speed vehicles using continuous transmit waves
RU2334244C1 (ru) Способ определения местоположения источника радиоизлучения
CN108845312A (zh) 基于脉冲体制无线电高度表的测高方法
CN115616629B (zh) 一种基于天基外辐射源信号的运动目标探测补偿方法
WO2005050240A2 (en) Method and system for determining the range and velocity of a moving object
GB2298538A (en) Dual Mode Radar System
US8138962B2 (en) Method for processing measured vertical profiles of the power of the echoes returned following a transmission of radar signals
Zirizzotti et al. Glacio RADAR system and results
RU2697509C2 (ru) Способ обнаружения, измерения дальности и скорости низколетящей малоскоростной цели в импульсно-доплеровских радиолокационных станциях при высокой частоте повторения импульсов и инвертируемой линейной частотной модуляции
US6204805B1 (en) Dual target tracking altimeter
KR101945795B1 (ko) 확장 및 휴대 가능한 계측 레이더 시스템

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20190924