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

一种星载sar电离层色散效应的判决与处理方法
技术领域
1.本发明涉及星载合成孔径雷达（synthetic aperture radar, sar）领域，尤其是一种星载sar电离层色散效应的判决与处理方法。


背景技术：

2.合成孔径雷达（synthetic aperture radar, sar）是一种主动式微波成像雷达，具有全天时、全天候的工作特点，在环境保护、灾害监测、海洋观测、资源勘查和军事侦察等方面的应用上具有独特的优势。
3.对于许多的军事应用或民事应用，例如生物量反演、灾害预警、环境监测等等，通常会使用具有更好穿透性能的低频雷达信号（通常小于2 ghz）。但是，由于电离层的存在，包括色散效应等问题会对信号造成影响，使得成像结果出现散焦和几何偏移等问题。针对电离层色散效应的问题，传统的校正思路是在干涉处理前进行校正。但是当频段较低（例如l波段和p波段）和太阳活动强烈时，将导致成像结果出现较严重的散焦现象。因此，有必要在成像处理前加入电离层色散效应的校正步骤以获得精聚焦的成像结果。


技术实现要素：

4.有鉴于此，本发明的主要目的在于提出星载sar电离层色散效应的判决与处理方法，从而建立电离层色散效应校正的处理流程。
5.本发明利用国际参考电离层（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电离层色散效应校正的处理流程。
6.为达到上述目的，本发明的技术方案为：一种星载sar电离层色散效应的判决与处理方法，包括如下步骤：1、利用iri模型获得tec参数，并结合星载sar系统参数计算电离层色散效应的特征参数，即aio、qpe和cpe；2、判断aio、qpe和cpe是否超过各自的门限值，并获得各自的判决标志位；3、根据判决标志位决定是否进行电离层色散效应的各项校正步骤。
7.有益效果：本发明将对星载sar电离层色散效应中的空变性、成像影响和干涉影响进行判决，根据判决结果决定是否执行相应的电离层色散效应校正步骤。当特征参数小于门限值时，将不执行相应的校正步骤以提高星载sar的信号处理效率；而当特征参数大于门限值时，将执行相应的校正步骤，以抑制电离层色散效应。本发明提出的判决方法和标准化电离层色
散效应校正的处理流程，兼顾了信号处理效率与高精度星载sar产品的需求，对于未来低频段星载sar地面处理系统的研制具有重要的参考意义。
附图说明
8.图1为结合所提的判决方法建立的电离层色散效应校正的处理流程；图2表示某星载l波段sar的系统参数下，方位向图像偏移随方位向tec梯度值的变化情况；图3表示某星载l波段sar的系统参数下，带宽为150 m时qpe随tec的变化情况；图4表示某星载l波段sar的系统参数下，带宽为150 m时cpe随tec的变化情况；图5表示某星载l波段sar的系统参数下，在
∆
teca取值分别为0 tecu、0.005 tecu和0.5 tecu时的星载sar方位向成像结果。
9.图6表示某星载l波段sar的系统参数下，在tec取值分别为10 tecu、40 tecu和80 tecu时的星载sar距离向成像结果。
具体实施方式
10.本发明提出了星载sar电离层色散效应的判决方法与处理流程，完整流程如图1所示。
11.本发明考虑到不是每一景星载sar回波数据都会面临严重的电离层色散效应（电离层色散效应的程度与数据获取的时间有很大的关系，一般来说，在中午时分的太阳活动强烈，电离层色散效应较严重），同时还需要兼顾星载sar的信号处理效率，有必要对电离层色散效应引起的主要影响进行判决从而决定是否需要进行相应的校正步骤。方位向偏移量（azimuth image offset, aio）、二次相位误差（quadratic phase error, qpe）和三次相位误差（cubical phase error, cpe）是三个典型的电离层色散效应影响的特征参数，利用这三个特征参数对电离层色散效应的成像影响、干涉影响和空变性进行衡量，从而建立一个标准化的电离层色散效应校正的处理方法。
12.该发明的星载sar电离层色散效应的判决与处理方法包括如下步骤：步骤101：利用iri模型获得tec参数，并结合星载sar系统参数计算电离层色散效应的特征参数，即aio、qpe和cpe：aio的计算方式如下：
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(1)其中，f0表示载频；k表示标准电子参数，取值通常为40.28m3/s2；r0表示卫星到场景中心点的最近斜距；
∆
teca表示沿方位向的tec梯度值，用以表征电离层色散效应的空变性。
13.qpe的计算方式如下：
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(2)其中，br表示带宽；c表示光速；tec表示信号传播路径上的总电子含量。
14.cpe的计算方式如下所示：
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(3)式(2)和式(3)中的tec可以从iri模型中直接获得，而式(1)中
∆
teca的则可以利用iri模型和空间几何关系，通过计算获得。设星载sar卫星在方位η时刻的地球固定坐标系下的坐标为[xs(η),ys(η),zs(η)]，电离层高度是zi，场景中某一点的坐标为(x0,y0,z0)，利用空间直线方程，雷达波束与电离层在方位η时刻的交点p坐标的可以由下式计算：
ꢀꢀꢀ
(4)进一步地，方位向尺寸变换因子da可以计算得到：
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(5)利用iri模型获得交点p位置处的tec值为tec(p;η)，进而获得沿方位向的tec分量teca(p;η)：teca(p;η)=tec(p;η)
·
daꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(6)沿方位向的tec梯度值
∆
teca可以计算得到：(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。
[0015]
步骤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，此时需要在干涉处理前进行电离层色散效应校正。
[0016]
下面结合具体实施例对本发明的技术方案作进一步详细说明。
[0017]
实施例1本处采用实际l波段星载sar的系统参数进行仿真，展示判决门限的设置对于电离层色散效应的影响是否合理。图2给出了某星载l波段sar的系统参数下，方位向图像偏移随方位向tec梯度值的变化情况。假设距离向分辨率为3 m，可见当方位向tec梯度大于0.0095 tecu/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的系统参数下，在
∆
teca取值分别为0 tecu/km、0.005 tecu/km和0.5 tecu/km时的星载sar方位向成像结果。可见，其中
∆
teca的取值是参考图2的结果，0.005 tecu/km使得aio小于门限值；0.05 tecu/km使得aio大于门限值。可见，当
∆
teca取值为0.005 tecu/km时，方位向几乎不存在偏移；而当
∆
teca取值为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时，距离向完全散焦，分辨率严重恶化。上述仿真结果说明本发明所提的判决方法是合理的。
[0018]
以上所述，仅为本发明的一个实施例而已，并非用于限定本发明的保护范围。