CN106019277B - 星载sar方位成像中的电离层传播效应影响判决方法 - Google Patents
星载sar方位成像中的电离层传播效应影响判决方法 Download PDFInfo
- Publication number
- CN106019277B CN106019277B CN201610300209.9A CN201610300209A CN106019277B CN 106019277 B CN106019277 B CN 106019277B CN 201610300209 A CN201610300209 A CN 201610300209A CN 106019277 B CN106019277 B CN 106019277B
- Authority
- CN
- China
- Prior art keywords
- msub
- mrow
- satellite
- mfrac
- msup
- 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
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO 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/00—Systems 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/88—Radar or analogous systems specially adapted for specific applications
- G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
- G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO 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/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/40—Means for monitoring or calibrating
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
本发明提供一种星载SAR方位成像中的电离层传播效应影响判决方法。技术方案为:以星载SAR系统参数为输入,第一步:针对观测场景中任意目标点P,计算得到合成孔径中心时刻t0的电离层穿刺点坐标第二步:针对穿刺点坐标利用在线IRI模型获取t0时刻的第三步:判断与门限值的大小,得到判决结果。本发明针对星载SAR方位成像处理,提出了电离层传播效应的影响判决方法,给出了判决门限和判决方法,在星载SAR电离层影响分析与校正处理中有广泛应用。
Description
技术领域
本发明属于航天与微波遥感结合的交叉技术领域,特别涉及一种星载SAR(Synthetic Aperture Radar,合成孔径雷达)方位成像中的电离层传播效应的影响判决方法。
背景技术
随着航天技术及微波探测技术的日益发展,星载SAR作为空间探测的重要手段,发挥着越来越重要的作用。当前,星载SAR的频段范围和分辨能力进一步扩展,功能进一步增强,其中工作在低频段的高分辨星载SAR是主要的发展方向之一。低频段高分辨的星载SAR具有良好的穿透性能,既能够发现隐蔽的军事目标,还广泛应用于森林生物量反演。然而,随着工作频率的降低和分辨率的提高,电离层传播效应对其成像性能的影响愈加严重。
在现有的星载SAR电离层效应研究中,关注的重点在于电离层色散效应对距离向成像的影响。然而,随着系统分辨率的不断提高,电离层对星载SAR方位向成像也会产生一定的影响。在星载SAR方位成像的过程中,电磁波在电离层内的传播路径不断变化,导致电离层STEC(Slant Total Electron Content,斜距向电子总量)不断变化,进而引入沿方位向变化的相位误差,可能造成方位向成像质量的下降。
电离层传播效应对星载SAR方位成像的影响集中反映在STEC(Slant TotalElectron Content,斜距向电子总量)在合成孔径时间内的变化上,目前尚未查到相应的影响判决方法。
发明内容
本发明的目的是:提出一种电离层传播效应对星载SAR方位成像的影响判决方法,可以应用于星载SAR电离层影响分析与校正处理。
本发明的技术方案是:
已知星载SAR的系统参数包括:载频fc,合成孔径时间Ts,方位向分辨率ρa,观测场景中任意目标点P在东北天坐标系下的坐标为(x0,y0,z0),该点对应的合成孔径中心时刻为t0,星载SAR在方位向慢时间tn时刻的东北天坐标系轨道坐标数据为n的取值为任意整数并且|tn-t0|≤Ts/2。电离层的高度为zi。
第一步,针对观测场景中任意目标点P,利用空间直线方程,计算得到其合成孔径中心时刻t0的星载SAR波束在电离层上的穿刺点(Ionospheric Penetration Point,IPP)坐标
第二步,针对穿刺点坐标利用在线IRI(International ReferenceIonosphere,国际参考电离层)模型获取合成孔径中心时刻t0的电离层垂直向TEC(TotalElectron Content,电子总量)值再利用下式得到电离层STEC值:
其中,几何变换因子γ(t0)为:
第三步,如果下述不等式成立,
则对观测场景中任意目标点P的方位向成像时忽略电离层传播效应的影响,其中c=3×108m/s为真空中光速,K=40.28m3/s2;否则,必须考虑其影响。
采用本发明可取得以下技术效果:
本发明针对星载SAR方位成像处理,提出了电离层传播效应的影响判决方法,给出了判决门限和判决方法,在星载SAR电离层影响分析与校正处理中有广泛应用。
附图说明
图1是本发明提供的星载SAR方位成像中的电离层传播效应影响判决流程图;
图2是星载SAR方位成像与电离层传播效应的几何示意图;
图3是电离层传播效应超过判决门限时的星载SAR方位成像结果;
图4是电离层传播效应未超过判决门限时的星载SAR方位成像结果。
具体实施方式
下面将结合附图对本发明提供的方法进行详细说明。
图1为本发明提供的星载SAR方位成像中的电离层传播效应影响判决流程图,总共分为三步。以星载SAR系统参数为输入,第一步:针对观测场景中任意目标点P,计算得到合成孔径中心时刻t0的电离层穿刺点坐标第二步:针对穿刺点坐标利用在线IRI模型获取t0时刻的第三步:判断与门限值的大小,得到判决结果。
图2为星载SAR方位成像与电离层传播效应的几何关系示意图。其中,O为东北天坐标系原点,P为观测场景中的任意目标点,星载SAR的运行轨迹近似为直线(图中细实线所示),粗虚线代表抽象为薄屏的电离层,实心圆点IPP代表t0时刻的星载SAR到目标点P的雷达波束在电离层上的穿刺点,Zi为电离层高度,取值范围通常为350~400km。穿刺点到目标点P的距离代表了电磁波在电离层内的传播路径,也表征了STEC的大小。可以看出,随着星载SAR的运动,电磁波在电离层内的传播路径不断变化,即STEC不断变化。整个合成孔径内信号传播路径变化引入的方位向QPE(Quadratic Phase Error,二次相位误差)为:
其中,合成角ψsyn由方位向分辨率决定:
通常方位向QPE≤π/4时可以忽略传播路径变化的影响,进而可以得到的判决门限为:
图3为电离层传播效应超过判决门限时的星载SAR方位成像仿真结果。仿真参数设置为:fc=470MHz,ρa=4.45m,因此判决门限为43.3TECU,利用本发明进行判决的结果是需要考虑传播效应的影响。为了便于对比,图3同时给出了未包含电离层传播效应(理想情况,虚线表示)和包含电离层传播效应(实际情况,实线表示)两种情况下的方位向成像结果。由图3可以发现,由于超过了判决门限,相比于理想情况,星载SAR方位向成像质量明显下降,主要表现在主瓣稍显展宽,副瓣明显抬高。
图4是电离层传播效应未超过判决门限时的星载SAR方位成像结果。仿真参数设置为:fc=470MHz,ρa=4.45m,因此判决门限为43.3TECU,利用本发明进行判决的结果是不需要考虑传播效应的影响。类似地,图4同时给出了未包含电离层传播效应(理想情况,虚线表示)和包含电离层传播效应(实际情况,实线表示)两种情况下的方位向成像结果。由图4可以发现,由于低于判决门限,星载SAR方位向成像质量相比于理想情况几乎没有下降。
Claims (1)
1.一种星载SAR方位成像中的电离层传播效应影响判决方法,SAR是指合成孔径雷达,已知星载SAR的系统参数包括:载频fc,合成孔径时间Ts,方位向分辨率ρa;观测场景中任意目标点P在东北天坐标系下的坐标为(x0,y0,z0),该点对应的合成孔径中心时刻为t0;星载SAR在方位向慢时间tn时刻的东北天坐标系轨道坐标数据为n的取值为任意整数并且|tn-t0|≤Ts/2;电离层的高度为zi,其特征在于,包括下述步骤:
第一步,针对观测场景中任意目标点P,利用空间直线方程,计算得到其合成孔径中心时刻t0的星载SAR波束在电离层上的穿刺点坐标
<mrow>
<msub>
<mi>x</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<msub>
<mi>t</mi>
<mn>0</mn>
</msub>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>z</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<msub>
<mi>z</mi>
<mn>0</mn>
</msub>
</mrow>
<mrow>
<msub>
<mi>z</mi>
<msub>
<mi>t</mi>
<mn>0</mn>
</msub>
</msub>
<mo>-</mo>
<msub>
<mi>z</mi>
<mn>0</mn>
</msub>
</mrow>
</mfrac>
<mo>&CenterDot;</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<msub>
<mi>t</mi>
<mn>0</mn>
</msub>
</msub>
<mo>-</mo>
<msub>
<mi>x</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>x</mi>
<mn>0</mn>
</msub>
</mrow>
<mrow>
<msub>
<mi>y</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<msub>
<mi>t</mi>
<mn>0</mn>
</msub>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>z</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<msub>
<mi>z</mi>
<mn>0</mn>
</msub>
</mrow>
<mrow>
<msub>
<mi>z</mi>
<msub>
<mi>t</mi>
<mn>0</mn>
</msub>
</msub>
<mo>-</mo>
<msub>
<mi>z</mi>
<mn>0</mn>
</msub>
</mrow>
</mfrac>
<mo>&CenterDot;</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>y</mi>
<msub>
<mi>t</mi>
<mn>0</mn>
</msub>
</msub>
<mo>-</mo>
<msub>
<mi>y</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>y</mi>
<mn>0</mn>
</msub>
</mrow>
第二步,针对穿刺点坐标利用在线国际参考电离层模型获取合成孔径中心时刻t0的电离层垂直向电子总量的值再利用下式得到电离层值:
<mrow>
<mi>STEC</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<msub>
<mi>t</mi>
<mn>0</mn>
</msub>
</mrow>
</msub>
<mo>,</mo>
<msub>
<mi>y</mi>
<msub>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
<mn>0</mn>
</msub>
</msub>
<mo>,</mo>
<msub>
<mi>z</mi>
<mi>i</mi>
</msub>
<mo>,</mo>
<msub>
<mi>t</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>TEC</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<msub>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
<mn>0</mn>
</msub>
</msub>
<mo>,</mo>
<msub>
<mi>y</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<msub>
<mi>t</mi>
<mn>0</mn>
</msub>
</mrow>
</msub>
<mo>,</mo>
<msub>
<mi>z</mi>
<mi>i</mi>
</msub>
<mo>,</mo>
<msub>
<mi>t</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>&CenterDot;</mo>
<mi>&gamma;</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>t</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
</mrow>
其中,几何变换因子γ(t0)为:
<mrow>
<mi>&gamma;</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>t</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<msqrt>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<msub>
<mi>t</mi>
<mn>0</mn>
</msub>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>x</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>y</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<msub>
<mi>t</mi>
<mn>0</mn>
</msub>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>y</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>z</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<msub>
<mi>z</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
<mrow>
<msub>
<mi>z</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<msub>
<mi>z</mi>
<mn>0</mn>
</msub>
</mrow>
</mfrac>
</mrow>
第三步,如果下述不等式成立,
<mrow>
<mi>S</mi>
<mi>T</mi>
<mi>E</mi>
<mi>C</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<msub>
<mi>t</mi>
<mn>0</mn>
</msub>
</mrow>
</msub>
<mo>,</mo>
<msub>
<mi>y</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<msub>
<mi>t</mi>
<mn>0</mn>
</msub>
</mrow>
</msub>
<mo>,</mo>
<msub>
<mi>z</mi>
<mi>i</mi>
</msub>
<mo>,</mo>
<msub>
<mi>t</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>&le;</mo>
<mfrac>
<mrow>
<msub>
<mi>cf</mi>
<mi>c</mi>
</msub>
</mrow>
<mrow>
<mn>8</mn>
<mi>K</mi>
<mi> </mi>
<msup>
<mi>tan</mi>
<mn>2</mn>
</msup>
<mrow>
<mo>(</mo>
<mn>0.886</mn>
<mo>&times;</mo>
<mfrac>
<mi>c</mi>
<mrow>
<mn>4</mn>
<msub>
<mi>f</mi>
<mi>c</mi>
</msub>
<msub>
<mi>&rho;</mi>
<mi>a</mi>
</msub>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
</mrow>
则对观测场景中任意目标点P的方位向成像时忽略电离层传播效应的影响,其中c=3×108m/s为真空中光速,K=40.28m3/s2;否则,必须考虑其影响。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610300209.9A CN106019277B (zh) | 2016-05-09 | 2016-05-09 | 星载sar方位成像中的电离层传播效应影响判决方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610300209.9A CN106019277B (zh) | 2016-05-09 | 2016-05-09 | 星载sar方位成像中的电离层传播效应影响判决方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106019277A CN106019277A (zh) | 2016-10-12 |
CN106019277B true CN106019277B (zh) | 2018-03-20 |
Family
ID=57099131
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610300209.9A Active CN106019277B (zh) | 2016-05-09 | 2016-05-09 | 星载sar方位成像中的电离层传播效应影响判决方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106019277B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108008367B (zh) * | 2017-12-01 | 2019-09-13 | 中国人民解放军国防科技大学 | 星载单航过InSAR系统电离层误差校正方法 |
CN108845293A (zh) * | 2018-04-18 | 2018-11-20 | 中国人民解放军国防科技大学 | 一种星载低波段全极化sar电离层fra估计方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6919839B1 (en) * | 2004-11-09 | 2005-07-19 | Harris Corporation | Synthetic aperture radar (SAR) compensating for ionospheric distortion based upon measurement of the group delay, and associated methods |
CN103675775A (zh) * | 2013-12-12 | 2014-03-26 | 北京理工大学 | 背景电离层对geo sar成像影响的分析方法 |
CN103969660A (zh) * | 2014-05-16 | 2014-08-06 | 中国科学院光电研究院 | 电离层误差修正方法 |
CN104793191A (zh) * | 2015-02-03 | 2015-07-22 | 北京理工大学 | 一种背景电离层对geo sar成像影响分析方法及其验证方法 |
-
2016
- 2016-05-09 CN CN201610300209.9A patent/CN106019277B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6919839B1 (en) * | 2004-11-09 | 2005-07-19 | Harris Corporation | Synthetic aperture radar (SAR) compensating for ionospheric distortion based upon measurement of the group delay, and associated methods |
CN103675775A (zh) * | 2013-12-12 | 2014-03-26 | 北京理工大学 | 背景电离层对geo sar成像影响的分析方法 |
CN103969660A (zh) * | 2014-05-16 | 2014-08-06 | 中国科学院光电研究院 | 电离层误差修正方法 |
CN104793191A (zh) * | 2015-02-03 | 2015-07-22 | 北京理工大学 | 一种背景电离层对geo sar成像影响分析方法及其验证方法 |
Also Published As
Publication number | Publication date |
---|---|
CN106019277A (zh) | 2016-10-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102346249B (zh) | 合成孔径雷达宽测绘带对地观测步进扫描模式的实现方法 | |
CN106249237B (zh) | 一种曲线轨迹下大斜视sar频域成像方法 | |
CN104316924B (zh) | 机载超高分辨率sar反投影图像的自聚焦运动补偿方法 | |
Carrano et al. | Simulating the impacts of ionospheric scintillation on L band SAR image formation | |
CN105785327B (zh) | 频率分集阵列合成孔径雷达高分辨宽测绘带成像方法 | |
CN102680974B (zh) | 一种星载滑动聚束合成孔径雷达的信号处理方法 | |
CN106093932B (zh) | 一种扫描波束高分辨率的雷达散射计 | |
CN105487074B (zh) | 一种双基地合成孔径雷达数值距离多普勒成像方法 | |
CN103235305B (zh) | 星载超高分辨率滑动聚束sar成像方法 | |
CN105372657A (zh) | 基于回波数据的视频合成孔径雷达运动补偿成像方法 | |
CN104407349B (zh) | 一站固定式双站低频超宽带sar的频域成像方法 | |
CN104280566A (zh) | 基于空时幅相估计的低空风切变风速估计方法 | |
CN103487798A (zh) | 一种相控阵雷达测高方法 | |
CN106019277B (zh) | 星载sar方位成像中的电离层传播效应影响判决方法 | |
Mohanty et al. | Ionospheric scintillation observation using space-borne synthetic aperture radar data | |
Ji et al. | Spaceborne P-band SAR imaging degradation by anisotropic ionospheric irregularities: A comprehensive numerical study | |
CN105785336B (zh) | Geosar方位成像中的电离层时变效应影响判决方法 | |
CN107271995A (zh) | 基于波束指向调整的系统灵敏度优化设计方法 | |
Ji et al. | Extended scintillation phase gradient autofocus in future spaceborne P-band SAR mission | |
CN104020472A (zh) | 一种便于实时处理的方位ncs大斜视sar成像方法 | |
CN101846739A (zh) | Sar扩展场景原始数据混合域仿真方法 | |
Fan et al. | High frame-rate and low-latency video SAR based on robust Doppler parameters estimation in the terahertz regime | |
CN106872949B (zh) | 一种基于自适应平衡加载的杂波谱配准补偿方法 | |
CN105044700A (zh) | 基于星载高分辨率sar图像聚焦性能的地面高程反演方法 | |
CN107271996A (zh) | 一种机载cssar地面运动目标成像方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |