CN114578310B - 一种星载sar电离层色散效应的判决与处理方法 - Google Patents

一种星载sar电离层色散效应的判决与处理方法 Download PDF

Info

Publication number
CN114578310B
CN114578310B CN202210478091.4A CN202210478091A CN114578310B CN 114578310 B CN114578310 B CN 114578310B CN 202210478091 A CN202210478091 A CN 202210478091A CN 114578310 B CN114578310 B CN 114578310B
Authority
CN
China
Prior art keywords
dispersion effect
qpe
cpe
aio
ionospheric
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
Application number
CN202210478091.4A
Other languages
English (en)
Other versions
CN114578310A (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.)
Aerospace Information Research Institute of CAS
Original Assignee
Aerospace Information Research Institute 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 Aerospace Information Research Institute of CAS filed Critical Aerospace Information Research Institute of CAS
Priority to CN202210478091.4A priority Critical patent/CN114578310B/zh
Publication of CN114578310A publication Critical patent/CN114578310A/zh
Application granted granted Critical
Publication of CN114578310B publication Critical patent/CN114578310B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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/02Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
    • G01S13/0218Very long range radars, e.g. surface wave radar, over-the-horizon or ionospheric propagation systems
    • 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)
  • Position Fixing By Use Of Radio Waves (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明公开一种星载SAR电离层色散效应的判决与处理方法,以实现星载SAR数据产品的高精度和信号处理的高效率。结合所提的判决方法,建立的处理流程包括:利用IRI模型获得TEC参数,并结合星载SAR系统参数计算电离层色散效应的特征参数,即AIO、QPE和CPE;判断AIO、QPE和CPE是否超过各自的门限值,并获得相应的判决标志位;根据判决标志位决定是否进行电离层色散效应的各项校正步骤。通过本发明,能够建立起高效的星载SAR电离层色散效应校正的处理流程。

Description

一种星载SAR电离层色散效应的判决与处理方法
技术领域
本发明涉及星载合成孔径雷达(Synthetic Aperture Radar, SAR)领域,尤其是一种星载SAR电离层色散效应的判决与处理方法。
背景技术
合成孔径雷达(Synthetic Aperture Radar, SAR)是一种主动式微波成像雷达,具有全天时、全天候的工作特点,在环境保护、灾害监测、海洋观测、资源勘查和军事侦察等方面的应用上具有独特的优势。
对于许多的军事应用或民事应用,例如生物量反演、灾害预警、环境监测等等,通常会使用具有更好穿透性能的低频雷达信号(通常小于2 GHz)。但是,由于电离层的存在,包括色散效应等问题会对信号造成影响,使得成像结果出现散焦和几何偏移等问题。针对电离层色散效应的问题,传统的校正思路是在干涉处理前进行校正。但是当频段较低(例如L波段和P波段)和太阳活动强烈时,将导致成像结果出现较严重的散焦现象。因此,有必要在成像处理前加入电离层色散效应的校正步骤以获得精聚焦的成像结果。
发明内容
有鉴于此,本发明的主要目的在于提出星载SAR电离层色散效应的判决与处理方法,从而建立电离层色散效应校正的处理流程。
本发明利用国际参考电离层(International Reference Ionosphere, IRI)模型获得总电子含量(Total Electron Content, TEC)参数,并结合星载SAR系统参数计算电离层色散效应的特征参数,即方位向偏移量(Azimuth Image Offset, AIO)、二次相位误差(Quadratic Phase Error, QPE)和三次相位误差(Cubical Phase Error, CPE)。判断AIO、QPE和CPE是否超过各自的门限值,并获得相应的判决标志位。根据判决标志位决定是否进行电离层色散效应的各项校正步骤。此外,结合上述判决方法,本发明将建立星载SAR电离层色散效应校正的处理流程。
为达到上述目的,本发明的技术方案为:一种星载SAR电离层色散效应的判决与处理方法,包括如下步骤:
1、利用IRI模型获得TEC参数,并结合星载SAR系统参数计算电离层色散效应的特征参数,即AIO、QPE和CPE;
2、判断AIO、QPE和CPE是否超过各自的门限值,并获得各自的判决标志位;
3、根据判决标志位决定是否进行电离层色散效应的各项校正步骤。
有益效果:
本发明将对星载SAR电离层色散效应中的空变性、成像影响和干涉影响进行判决,根据判决结果决定是否执行相应的电离层色散效应校正步骤。当特征参数小于门限值时,将不执行相应的校正步骤以提高星载SAR的信号处理效率;而当特征参数大于门限值时,将执行相应的校正步骤,以抑制电离层色散效应。本发明提出的判决方法和标准化电离层色散效应校正的处理流程,兼顾了信号处理效率与高精度星载SAR产品的需求,对于未来低频段星载SAR地面处理系统的研制具有重要的参考意义。
附图说明
图1为结合所提的判决方法建立的电离层色散效应校正的处理流程;
图2表示某星载L波段SAR的系统参数下,方位向图像偏移随方位向TEC梯度值的变化情况;
图3表示某星载L波段SAR的系统参数下,带宽为150 M时QPE随TEC的变化情况;
图4表示某星载L波段SAR的系统参数下,带宽为150 M时CPE随TEC的变化情况;
图5表示某星载L波段SAR的系统参数下,在∆TEC a 取值分别为0 TECU、0.005 TECU和0.5 TECU时的星载SAR方位向成像结果。
图6表示某星载L波段SAR的系统参数下,在TEC取值分别为10 TECU、40 TECU和80TECU时的星载SAR距离向成像结果。
具体实施方式
本发明提出了星载SAR电离层色散效应的判决方法与处理流程,完整流程如图1所示。
本发明考虑到不是每一景星载SAR回波数据都会面临严重的电离层色散效应(电离层色散效应的程度与数据获取的时间有很大的关系,一般来说,在中午时分的太阳活动强烈,电离层色散效应较严重),同时还需要兼顾星载SAR的信号处理效率,有必要对电离层色散效应引起的主要影响进行判决从而决定是否需要进行相应的校正步骤。方位向偏移量(Azimuth Image Offset, AIO)、二次相位误差(Quadratic Phase Error, QPE)和三次相位误差(Cubical Phase Error, CPE)是三个典型的电离层色散效应影响的特征参数,利用这三个特征参数对电离层色散效应的成像影响、干涉影响和空变性进行衡量,从而建立一个标准化的电离层色散效应校正的处理方法。
该发明的星载SAR电离层色散效应的判决与处理方法包括如下步骤:
步骤101:利用IRI模型获得TEC参数,并结合星载SAR系统参数计算电离层色散效应的特征参数,即AIO、QPE和CPE:
AIO的计算方式如下:
Figure 362077DEST_PATH_IMAGE001
(1)
其中,f 0表示载频;K表示标准电子参数,取值通常为40.28m3/s2R 0表示卫星到场景中心点的最近斜距;∆TEC a 表示沿方位向的TEC梯度值,用以表征电离层色散效应的空变性。
QPE的计算方式如下:
Figure 89730DEST_PATH_IMAGE002
(2)
其中,B r 表示带宽;c表示光速;TEC表示信号传播路径上的总电子含量。
CPE的计算方式如下所示:
Figure 733201DEST_PATH_IMAGE003
(3)
式(2)和式(3)中的TEC可以从IRI模型中直接获得,而式(1)中∆TEC a 的则可以利用IRI模型和空间几何关系,通过计算获得。设星载SAR卫星在方位η时刻的地球固定坐标系下的坐标为[x s (η),y s (η),z s (η)],电离层高度是z i ,场景中某一点的坐标为(x 0,y 0,z 0),利用空间直线方程,雷达波束与电离层在方位η时刻的交点P坐标的
Figure 333947DEST_PATH_IMAGE004
可以由下式计算:
Figure 62868DEST_PATH_IMAGE005
(4)
进一步地,方位向尺寸变换因子D a 可以计算得到:
Figure 344945DEST_PATH_IMAGE006
(5)
利用IRI模型获得交点P位置处的TEC值为TEC(P;η),进而获得沿方位向的TEC分量TEC a (P;η):
TEC a (P;η)=TEC(P;ηD a (6)
沿方位向的TEC梯度值∆TEC a 可以计算得到:
Figure 842923DEST_PATH_IMAGE007
(7)
步骤102:判断AIO、QPE和CPE是否超过各自的门限值,并获得各自的判决标志位:
对于特征参数AIO,取判决门限值为ρ a /2,其中ρ a 为方位向分辨率。如果计算出的AIO大于ρ a /2,则令AIO的判决标志位F1为TRUE;否则F1为FALSE。对于特征参数QPE,取判决门限值为π/4。如果计算出的QPE大于π/4,则令QPE的判决标志位F2为TRUE;否则F2为FALSE。对于特征参数CPE,取判决门限值为π/8。如果计算出的CPE大于π/8,则令CPE的判决标志位F3为TRUE;否则F3为FALSE。
步骤103:根据判决标志位决定是否进行电离层色散效应的各项校正步骤:
首先,对AIO进行判决。当AIO小于判决门限时,令AIO的判决标志位F1为FALSE,此时不需要进行电离层空变性校正;当AIO大于判决门限时,令AIO的判决标志位F1为TRUE,此时需要进行电离层空变性校正。其次,以QPE的判决结果为依据分成两种方式。方式一:当QPE小于判决门限时,令QPE的判决标志位F2为FALSE,此时不需要执行电离层色散效应校正操作,直接进行成像处理和干涉处理即可。方式二:当QPE大于判决门限时,令QPE的判决标志位F2为TRUE。此时需要先在成像处理前执行电离层色散效应校正,再对CPE的标志位进行判决。当CPE小于判决门限时,令CPE的判决标志位F3为FALSE,此时不需要在干涉处理前进行电离层色散效应校正;当CPE大于判决门限时,令CPE的判决标志位F3为TRUE,此时需要在干涉处理前进行电离层色散效应校正。
下面结合具体实施例对本发明的技术方案作进一步详细说明。
实施例1
本处采用实际L波段星载SAR的系统参数进行仿真,展示判决门限的设置对于电离层色散效应的影响是否合理。图2给出了某星载L波段SAR的系统参数下,方位向图像偏移随方位向TEC梯度值的变化情况。假设距离向分辨率为3 m,可见当方位向TEC梯度大于0.0095TECU/km时,方位向偏移将超过距离向分辨率的一半,此时偏移是不可忽略的,需要进行电离层空变性校正。图3给出了某星载L波段SAR的系统参数下,带宽为150 M时QPE随TEC的变化情况。可见,当TEC大于16 TECU时,QPE将大于门限值π/4,此时色散效应将对距离向成像造成不可忽略的影响,需要在成像处理前进行电离层色散效应校正。图4给出了星载L波段SAR的系统参数下,带宽为150M时CPE随TEC的变化情况。可见,当TEC大于70 TECU时,QPE将大于门限值π/8,此时色散效应对于成像性能的恶化将尤为严重,因此不仅需要在成像处理前进行电离层色散效应校正,还需要在干涉处理前进行电离层色散效应校正。图5给出了某星载L波段SAR的系统参数下,在∆TEC a 取值分别为0 TECU/km、0.005 TECU/km和0.5 TECU/km时的星载SAR方位向成像结果。可见,其中∆TEC a 的取值是参考图2的结果,0.005 TECU/km使得AIO小于门限值;0.05 TECU/km使得AIO大于门限值。可见,当∆TEC a 取值为0.005TECU/km时,方位向几乎不存在偏移;而当∆TEC a 取值为0.05 TECU/km时,方位向存在不可忽略的偏移。图6给出了某星载L波段SAR的系统参数下,在TEC取值分别为10 TECU、40 TECU和80 TECU时的星载SAR距离向成像结果。其中TEC的取值是参考图2和图3的结果,10 TECU使得QPE和CPE均小于门限值;40 TECU使得QPE大于门限值,CPE小于门限值;80 TECU使得QPE和CPE均大于门限值。可见,当TEC取值为10 TECU时,距离向聚焦情况几乎不受影响;当TEC取值为40 TECU时,旁瓣提升严重,峰值旁瓣比恶化,距离向出现了很明显的散焦现象;当TEC取值为80 TECU时,距离向完全散焦,分辨率严重恶化。上述仿真结果说明本发明所提的判决方法是合理的。
以上所述,仅为本发明的一个实施例而已,并非用于限定本发明的保护范围。

Claims (3)

1.一种星载SAR电离层色散效应的判决与处理方法,其特征在于,包括如下步骤:
步骤S101、利用IRI模型获得TEC参数,并结合星载SAR系统参数计算电离层色散效应的特征参数,即AIO、QPE和CPE;利用IRI模型获得TEC参数,并结合星载SAR系统参数计算电离层色散效应的特征参数,即方位向偏移量AIO、二次相位误差QPE和三次相位误差CPE,具体包括:
首先,利用IRI模型获取随空间位置变化的电离层TEC参数;
其次,利用TEC参数和星载SAR系统参数计算电离层色散效应的特征参数;电离层色散效应的影响包括距离向成像质量恶化和图像偏移;用QPE表征是否需要在成像前进行电离层色散效应校正,用 CPE表征是否需要在干涉前进行电离层色散效应校正,用AIO表征是否需要进行电离层色散效应的空变性校正;
步骤S102、判断AIO、QPE和CPE是否超过各自的门限值,并获得各自的判决标志位:
步骤103:根据判决标志位决定是否进行电离层色散效应的各项校正,基于决定的结果执行校正处理流程。
2.根据权利要求1所述的星载SAR电离层色散效应的判决与处理方法,其特征在于,判断AIO、QPE和CPE是否超过各自的门限值,并获得各自的判决标志位:
将计算出来的AIO、QPE和CPE与各自的门限值对比,判断是否需要考虑相应的影响,并更新判决标志位;如果AIO超过门限值,表示电离层色散效应存在空变性;如果QPE超过门限值,表示需要在成像处理前进行电离层色散效应校正;如果CPE超过门限值,表示不仅需要在成像处理前进行电离层色散效应校正,还需要在干涉处理前进行电离层色散效应校正;AIO、QPE和CPE对应的标志位分别为F1、F2和F3,对应的门限值分别为方位向分辨率的一半、π/4和π/8。
3.根据权利要求1所述的星载SAR电离层色散效应的判决与处理方法,其特征在于,根据判决标志位决定是否进行电离层色散效应的各项校正步骤具体包括:
首先,对AIO进行判决;当AIO小于判决门限时,令AIO的判决标志位F1为FALSE,此时不需要进行电离层空变性校正;当AIO大于判决门限时,令AIO的判决标志位F1为TRUE,此时需要进行电离层空变性校正;
其次,以QPE的判决结果为依据分成两种处理方式,方式一:当QPE小于判决门限时,令QPE的判决标志位F2为FALSE,此时不需要执行电离层色散效应校正操作,直接进行成像处理和干涉处理即可;方式二:当QPE大于判决门限时,令QPE的判决标志位F2为TRUE,此时需要先在成像处理前执行电离层色散效应校正,再对CPE的标志位进行判决,当CPE小于判决门限时,令CPE的判决标志位F3为FALSE,此时不需要在干涉处理前进行电离层色散效应校正;当CPE大于判决门限时,令CPE的判决标志位F3为TRUE,此时需要在干涉处理前进行电离层色散效应校正。
CN202210478091.4A 2022-05-05 2022-05-05 一种星载sar电离层色散效应的判决与处理方法 Active CN114578310B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210478091.4A CN114578310B (zh) 2022-05-05 2022-05-05 一种星载sar电离层色散效应的判决与处理方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210478091.4A CN114578310B (zh) 2022-05-05 2022-05-05 一种星载sar电离层色散效应的判决与处理方法

Publications (2)

Publication Number Publication Date
CN114578310A CN114578310A (zh) 2022-06-03
CN114578310B true CN114578310B (zh) 2022-07-29

Family

ID=81778853

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210478091.4A Active CN114578310B (zh) 2022-05-05 2022-05-05 一种星载sar电离层色散效应的判决与处理方法

Country Status (1)

Country Link
CN (1) CN114578310B (zh)

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106019279B (zh) * 2016-05-09 2018-04-06 中国人民解放军国防科学技术大学 星载sar方位成像中的电离层空变效应影响判决方法
CN105785336B (zh) * 2016-05-09 2018-02-16 中国人民解放军国防科学技术大学 Geosar方位成像中的电离层时变效应影响判决方法
CN114019509A (zh) * 2021-11-22 2022-02-08 中国地质大学(北京) 星载全极化p波段宽带sar电离层色散误差补偿方法

Also Published As

Publication number Publication date
CN114578310A (zh) 2022-06-03

Similar Documents

Publication Publication Date Title
CN100588988C (zh) 缺乏高精度运动补偿系统下实现机载雷达重聚焦成像方法
CN102288962B (zh) 一种超宽带合成孔径雷达实时运动补偿方法
CN111175749B (zh) 一种星载sar成像处理方法
CN111342882B (zh) Ka频段对星快速校相方法、存储介质和电子设备
US8212714B1 (en) Using doppler radar images to estimate aircraft navigational heading error
CN107783156B (zh) 卫星跟踪方法
CN111948464A (zh) 偏馈无线闭环自跟踪校相系统
CN113078939A (zh) 一种通信卫星多波束指向标校测量误差估计方法
Park et al. Effect of outlier removal from temporal ASF corrections on multichain Loran positioning accuracy
CN116772838B (zh) 一种用于机械相控阵天线的惯导误差补偿的方法
CN114076945A (zh) 一种nlfm信号的电离层色散效应解析模型的建立方法
CN114578310B (zh) 一种星载sar电离层色散效应的判决与处理方法
CN111736120B (zh) 一种基于天波传播校正源信号的阵列误差校正方法
US20220394661A1 (en) Estimating geolocation of a user terminal
CN111929669B (zh) 基于距离补偿的无线电侦测测距的修正方法及系统
CN111025358B (zh) 基于导航卫星信号单短基线的定向方法
Tian et al. Experimental research on an FM broadcast based T N R passive radar system
CN112147587B (zh) 一种雷达波束方位中心海上标定方法
Pi et al. Measuring Small-and Medium-Scale TEC Spatial Variations and Irregularities from Ground-Based GNSS Observations
Meller et al. On Radar DoA Estimation and Tilted Rotating Electronically Scanned Arrays
CN113595657B (zh) 基于太阳噪声的无线电测控设备校相方法及装置
Jun et al. A novel Ka band phase calibration method for satellite ground TT&C antenna based on UAV
Liu et al. An Imaging Compensation Scheme for Correcting Ionospheric Effect on High-Resolution Spaceborne P-Band SAR
Li et al. An Autofocus Algorithm Based on Topography Searching for Airborne SAR
Chen et al. Inclined MIMO GEO SAR System Design and Performance Analysis

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