作者中心
专家审稿
责编办公
编辑办公
Atmospheric Effects on the Performance of Geosynchronous Orbit SAR Systems
-
摘要: 地球同步轨道合成孔径雷达(GEO SAR)合成孔径时间长、观测范围大，易受到大气层效应时空变化的影响，使成像的聚焦质量和差分干涉处理精度严重下降。该文针对常规对流层和背景电离层等大气层缓变的干扰部分，建立了高精度时频混合GEO SAR信号模型，分析了不同大气层参数的时间变化率对成像质量和差分干涉处理精度的影响。针对大气层干扰中对流层湍流和电离层闪烁等随机扰动造成的影响，基于幂律功率谱模型，建立了大气扰动参数和成像质量的定量化分析模型，获得了大气层随机干扰的强度与成像评估指标间的关系。最后，通过仿真验证了模型的有效性，并分析了长孔径时间内缓变和随机扰动的大气层误差对成像和差分干涉处理质量的影响，仿真结果表明：L波段GEO SAR成像和差分干涉处理受时空变电离层干扰的影响十分严重，必须予以补偿；对流层干扰对其影响较小，仅当积累时间达到数百秒时才需要考虑它对成像性能的恶化。
-
关键词:
- 地球同步轨道合成孔径雷达 /
- 合成孔径雷达成像 /
- 对流层干扰影响 /
- 电离层干扰影响
Abstract: GEOsynchronous orbit Synthetic Aperture Radar (GEO SAR) has a long synthetic aperture and large observation region; therefore, it is easily affected by atmospheric spatial-temporal changes, which results in a serious degradation in the focusing quality and performance of differential interference processing. In this paper, a high-precision spatial-temporal hybrid GEO SAR signal model is established for the slowly disturbed parts of the atmosphere such as the background troposphere and the ionosphere. The effects of the time rate of different atmospheric parameters on the image quality and accuracy of differential interference processing are analyzed. Considering the influence of random disturbances such as tropospheric turbulence and ionospheric scintillation in atmosphere, a quantitative analysis model is established based on a power law power spectrum model using the atmospheric disturbance and imaging quality parameters. The relationships of the random disturbances intensity in the atmosphere and the imaging evaluation index are obtained. Finally, the model is verified by simulations, and the effects of slow-varying atmospheric errors in a long aperture time on the quality of imaging and differential interference processing are analyzed. Simulation results show that spatial-temporal variable ionospheric disturbances seriously affect L-band GEO SAR imaging and differential interference processing, which must be compensated, and they are slightly affected by tropospheric disturbances. Moreover, it is necessary to consider the tropospheric effects on imaging performance only when the integration time reaches several hundred seconds. -
1. 引言
当目标沿雷达视线(Light Of Sight, LOS)方向运动时，其回波信号的载频将发生偏移，即产生多普勒现象。除目标整体运动外，若目标或目标上的某些结构还存在独立的振动或旋转，则称其为微动。微动会在目标主体运动对应的主多普勒谱周围产生边带，即产生微多普勒效应[1]。直升机、飞机旋转叶片、小型卫星和空间碎片等航空航天目标的典型微动包括自旋、进动和章动[2]等。
对于空间微动目标，其高分辨雷达回波[3,4]蕴含着散射中心2维或3维分布等结构信息，同时包含着自旋频率、进动频率及进动角等运动信息，上述信息为准确的目标分类、识别提供了重要支撑。目前，典型空间微动目标的高分辨雷达成像与微动参数估计方法研究[5–9]已受到雷达成像与雷达自动目标识别领域的广泛关注。
空间微动目标的高分辨成像方法包括参数化方法[10–14]与非参数化[15–21]方法两类。其中，参数化成像方法首先建立各种微动形式的参数化模型，进而采用基于模型的参数估计方法实现高分辨成像。主要包括基于固定散射中心模型的成像方法[10,22–24]及基于滑动散射中心模型[11,13]的进动目标成像方法。对于章动等复杂微动形式，需要建立非常复杂的参数化模型，并实现大量未知参数的准确求解，由于目标的散射中心坐标与微动参数耦合，因此求解运算量很大。非参数化成像方法则主要包括自适应时频分析[25]与散射中心航迹关联[16,17,26–28]成像两类。与参数化成像方法相比，非参数化成像方法具有各种微动形式具有鲁棒性，能够避免由于模型失配而引起的较大误差，计算效率较高。
对于非参数化方法，基于航迹矩阵分解的成像方法[17]可实现自旋、进动、章动等微动目标的高分辨成像。该类方法的关键步骤之一是在距离-慢时间域实现散射中心航迹的精确估计和关联。现有方法采用卡尔曼滤波器和最小欧氏距离准则，实现基于1维斜距信息的航迹关联[17,29]，当散射中心回波包络交叉点较多、相距较近时容易产生较大的关联误差。此外，获取高质量的高分辨距离像(High Resolution Range Profile, HRRP)也是非参数化成像的关键步骤。
为了解决上述问题，本文提出一种基于距离-瞬时多普勒(Range-Instantaneous Doppler, RID)像序列的微动目标高分辨3维成像新方法。该方法充分利用散射中心在距离-瞬时多普勒域2维分布比距离-慢时间域1维分布可分性更强等特性，提出基于RID像序列的散射中心航迹关联方法，提高了航迹交叉点散射中心的可分性。进而通过带约束条件的矩阵分解求得散射中心3维分布和等效雷达视线矩阵，实现空间微动目标高分辨3维成像。最后，仿真数据证明了算法的有效性。
本文结构如下：第2节介绍了RID序列的生成方法；第3节研究了基于RID序列的航迹矩阵关联方法，以及基于现代谱估计的航迹矩阵精估计方法；第4节研究了基于航迹矩阵分解的微动目标高分辨3维成像方法；第5节以锥体章动目标为例，给出目标航迹关联及3维成像结果；最后一节进行了总结。
2. 距离-瞬时多普勒序列
对于信号
s(t) ，其短时傅里叶变换(STFT, Short-Time Fourier Transform)满足[15,30]：STFT(τ,ω)=∫s(t)w(t−τ)exp{−jωt}dt (1) 其中，
ω 表示角频率，τ 表示时延，w(⋅) 为窗函数。为实现散射中心航迹的准确关联，需要获得其距离-瞬时多普勒像序列。假设雷达发射大时宽-带宽积脉冲信号，距离脉压后回波共包含
Nr 个距离单元，则对存在回波的距离单元rn (n∈[N1,N2] ,[N1,N2] 为存在回波的距离单元区间)分别做STFT以得到其时频图In(fd,tm) ，其中fd 表示多普勒，tm 表示慢时间。随后，将时频图堆成3维矩阵Q(rn,fd,tm) 。最后，沿时间轴tm=ti (i∈[1,Na] ,Na 为方位单元数)取出2维矩阵切片，即得到ti 时刻的RID像。连续变换ti 即可获得RID像序列。该过程示意图如图1所示。3. 航迹矩阵关联及航迹矩阵精估计
3.1 基于分水岭法的RID像散射中心提取
为了利用RID像序列实现散射中心航迹关联，需要提取每幅图像中散射中心的2维坐标。分水岭(watershed)算法[31]能够精确定位图像中的微弱边缘，并获得封闭且连续的分割曲线，因此适用于提取RID像中的散射中心支撑域。此外，相比于基于统计学的图像分割算法，该算法计算量小且分割较为准确，适用于图像数据的实时处理。因此，本文首先使用采用分水岭方法对RID像进行图像分割以获得每个散射中心对应的支撑域，然后计算每个支撑域对应的散射中心质心，并将此质心作为散射中心2维坐标的粗估计。基于watershed方法的RID像分割过程实现方法如下：
Step1：将原始图像归一化后，通过设定门限值将其转化为二值图像；
Step2：计算二值图像中每个像素点到其最近非零点的距离(如果像素本身非零，则其本身为最近的非零点，因此距离为0)，用于替代该像素点的像素值，得到矩阵
D ；Step3：令
D=−D ，得到梯度图像；Step4：采用分水岭方法对Step3中得到的梯度图像进行分割[31]。
散射中心分割完成后，提取分割后每个散射中心的支撑域。具体步骤为：首先将分割后图像的1值和0值点赋为0，并将其他点赋为255；然后求二值图像的连通域；最后取其边界得到微动目标每个散射中心的支撑域。最后，将每个散射中心对应支持区的质心作为RID图像中每个散射中心2维坐标的估计。其中，质心计算方法如下：
xc=∑u∑vuf(u,v)∑u∑vf(u,v),yc=∑u∑vvf(u,v)∑u∑vf(u,v) (2) 其中，
f(u,v) 表示(u,v) 点处的像素值，u, v分别表示像素点的横坐标和纵坐标。将每一时刻的RID像都做上述处理，则可获得散射中心在各个时刻对应的坐标。
3.2 航迹矩阵关联
由于微动目标具有惯性，因此认为相邻两幅距离-瞬时多普勒图像中同一散射中心的坐标连续变化，从而基于最近邻法实现航迹关联。设第
i (i∈[1,Na] )幅RID图像中的第j (j∈[1,P] ,P 为散射中心个数)个散射中心的坐标向量为aij ，计算该点与第i+1 幅RID图像中各个散射中心坐标向量的欧氏距离，选取与其欧氏距离最小的散射中心作为与该点关联的散射中心，即计算式(3)：minj‖aij−ai+1,j‖2 (3) 依次计算RID序列中相邻两幅图像中各散射中心的关联点，从而实现RID图像中各散射中心的关联，并得到微动目标航迹矩阵
W 的粗估计。该矩阵的每一列对应一个散射中心在观测时间内的瞬时斜距。具体而言，基于RID序列的航迹关联实现方法如下：Step1：初始化航迹矩阵
WNa×P ，令所有元素都为0；Step2：将
P 个散射中心的初始时刻瞬时斜距写入WNa×P 的第1行；Step3：令
i=1 ,j=1 ，计算第i 幅RID像中第j 个散射中心与第i+1 幅RID像中所有散射中心的欧氏距离，根据式(3)将最小欧氏距离对应的散射中心瞬时斜距写入WNa×P(i+1,j) ；Step4：令
j=j+1 ，重复step3直到j=P ，实现第i 幅RID像与第i+1 幅RID像的2维航迹关联；Step5：令
i=i+1 ，重复step3—step4直到i=Na−1 ，获得矩阵WNa×P 。3.3 航迹矩阵精估计
由于RID像的距离分辨率为
ρr=c/2B ，对于X波段雷达而言通常为10–2 m量级，精度较低；而利用散射中心支撑域的质心对其2维坐标进行近似也会导致较大误差，从而使航迹矩阵产生抖动，影响微动目标3维成像的精度。针对该问题，在获得散射中心2维关联结果的基础上，可以进一步采用Root-MUSIC等谱估计方法[32]对散射中心的瞬时斜距进行精估计，并对航迹矩阵进行修正，从而提高对微动目标散射中心3维坐标估计的准确性。微动目标经运动补偿后的回波信号可表示为：
\begin{aligned} {s_0}\left( {f,{t_{\rm{m}}}} \right) =& \sum\limits_p {A_{{p}}}{\rm{rect}}\left( {\frac{f}{B}} \right)\\ {\rm{}}& \cdot \exp \left( {{\rm j}\frac{{4{{π}} }}{c}\left( {{f\!_{\rm{c}}} + f} \right)\Delta {R_{{p}}}\left( {{t_{\rm{m}}}} \right)} \right) \end{aligned} (4) 其中，
p \in \left[ {1,{P}} \right] 表示散射中心序号，{A_{{p}}} 表示其幅度，B 为带宽，c 为光速，{f\!_{\rm{c}}} 表示载频，{R_{{p}}} 表示第p 个散射中心与参考点之间的瞬时斜距。若忽略距离窗，则式(4)可被改写为{s_1}\left( {n,{t_{\rm{m}}}} \right) = \sum\limits_p {{{A'}\!\!_{{p}}}\exp \left( {{\rm j}{\omega _{{p}}}n} \right)} (5) 其中，
{A'\!\!_{{p}}} = {A_{{p}}}\exp \Bigr( {{\rm j}4{{π}} \left( {{f\!_{\rm{c}}} - B/2} \right)\Delta {R_{{p}}}\left( {{t_{\rm{m}}}} \right)/c} \Bigr) ，散射中心对应的角频率为{\omega _{{p}}} = 4{{π}} \Delta f\Delta {R_{{p}}}\left( {{t_{\rm{m}}}} \right)/c ,\Delta f = B/{N_{\rm{r}}} ,{N_{\rm{r}}} 为距离单元数，n \in \left[ {1,{N_{\rm{r}}}} \right] 。接下来，对每次距离向回波精确估计{\omega _{{p}}} 以求出\Delta {R_{{p}}}\left( {{t_{\rm{m}}}} \right) ，从而得到抑制旁瓣和噪声后的高质量高分辨1维距离像(HRRP)。为了采用Root-MUSIC方法，首先构造距离回波的协方差矩阵：
\hat {{R}} = \frac{1}{{{N_{\rm{r}}} - m}}\sum\limits_{n = m}^{{N_{\rm{r}}}} {{{\tilde {{S}}}_{\rm{r}}}\left( {n,{t_{\rm m}}} \right)\tilde {{S}}_1^*\left( {n,{t_{\rm m}}} \right)} (6) 其中，
m 表示窗长，且\begin{align} \tilde {{S}}_{\rm{r}}\left( {n,{t_{\rm{m}}}} \right) =& \Bigr[ {{s_1}\left( {n,{t_{\rm{m}}}} \right)}\ {{s_1}\left( {n - 1,{t_{\rm{m}}}} \right)}\ ·\!·· \\ {\rm{}}& \quad \ {s_1}\left( {n - m + 1,{t_{\rm{m}}}} \right) \Bigr]^{\rm{T}} \end{align} (7) 通过Z变换找到与单位元距离最近的P个根可以求得角频率
{\omega _{{p}}} 。随后，由\Delta {R_{{p}}}\left( {{t_{\rm{m}}}} \right) = {\omega _{{p}}}c/4{{π}} \Delta f 得到精估计的瞬时斜距\Delta {R_{{p}}}\left( {{t_{\rm{m}}}} \right) 。最后，通过最小欧氏距离准则将{t_{\rm{m}}} 时的瞬时斜距写入{{{W}}\!_{{{{N}}_{\rm{a}}} \times {{P\,}}}} 的相应行中，即可得到精估计的航迹矩阵{{{W}}'\!\!_{{{{N}}_{\rm{a}}} \times {{P}}}} 。此外，当回波的信噪比较低时，可以通过构造观测字典，采用噪声稳健的稀疏信号重构方法[33,34]获得HRRP，并实现航迹矩阵的精估计。4. 基于航迹矩阵分解的微动目标高分辨3维成像
根据运动的相对性，对于微动目标上的固定散射中心，其在距离-慢时间域的航迹矩阵可以表示为：
{{{W}}\!_{{{{N}}_{\rm{a}}} \times {{P}}}} = {{{R}}_{{{{N}}_{\rm{a}}} \times 3}}{{{S}}_{3 \times {{P}}}} (8) 其中，
P 为散射中心个数，{N_{\rm{a}}} 为方位单元数，矩阵{{R}} 表示不同时刻的等效雷达视线矩阵，{{S}} 表示目标的散射中心坐标矩阵。根据式(8)可知，从{{W}} 中重构矩阵{{S}} 则可得到目标3维散射中心坐标。本文采用基于矩阵奇异值分解的方法重构矩阵{{S}} [17,21]。利用矩阵奇异值分解法，航迹矩阵可以分解为
{{{W}}\!_{1{{{N}}_{\rm{a}}} \times {{K}}}} = {{{U}}\!_{{{{N}}_{\rm{a}}} \times {{K}}}}{{{Σ}} _{{{K}} \times {{K}}}}{{V}}_{{{K}} \times {{K}}}^{\rm{T}} 。对于3维微动，根据矩阵秩的特性，{{{Σ}} _{{{K}} \times {{K}}}} 的前3个奇异值较大，而其余奇异值趋近于零。因此可做如下近似：\begin{array}{l} {{{W}}\!_1} = \left[ {{{\left( {{{{U}}\!_1}} \right)}_{{{{N}}_{\rm{a}}} \times 3}},{{\left( {{{{U}}\!_2}} \right)}_{{{{N}}_{\rm{a}}} \times \left( {{{K}} - 3} \right)}}} \right]\left[ {\begin{array}{*{20}{l}} {{{\left( {{{{Σ}} _1}} \right)}_{3 \times 3}}} & 0\\ \quad \ 0 & 0 \end{array}} \right]\\ \quad\quad\quad \cdot\left[ \begin{array}{l} {\quad \left( {{{{V}}\!_1}} \right)_{3 \times {{K}}}}\\ {\left( {{{{V}}\!_2}} \right)_{\left( {{{K}} - 3} \right) \times {{K}}}} \end{array} \right] \approx {{{U}}\!_1}\left( {{{{Σ}} _1}{{{V}}\!_1}} \right)\\ \quad\quad= {{R}}'{{S}}' \end{array} (9) 其中，近似后得到
{{R}}' = {{{U}}\!_1} ,{{S}}' = {{{Σ}} _1}{{{V}}\!_1} 。并且对于任意可逆矩阵{{{A}}_{3 \times 3}} ,{{R}}'{{S}}' = \left( {{{R}}'{{A}}} \right)\left( {{{{A}}^{ - 1}}{{S}}'} \right) 成立。根据
{{R}} 的定义可知{{R}} 各行构成的行向量的模为1。将{{R}} 用行向量的形式表示为{{R}} = {\left[ {{{{l}}_1}\;{{{l}}_2}\; ·\!·\!· \;{{{l}}_{{{{N}}_{\rm{a}}}}}} \right]^{\rm{T}}} ，则下列等式成立：{{{l}}_n}{{A}}{{{A}}^{\rm{T}}}{{l}}_n^{\rm{T}} = {{I}},\;\;\;\;n \in \left[ {1,{N_{\rm{a}}}} \right] (10) 其中，
{{I}} 是单位矩阵。估计值{\hat {{A}\,}} 为式(10)的最小均方解，则{{\hat{{A}\,}}^{ - 1}}{{S}}' 相当于{{S}} 的等距映射。对于任意的正交矩阵
{{{A}}_1} ，满足下列关系：{{R}}{{S}} = \left(\! {{{R}}'{\hat{{A}\,}}{{{A}}_1}} \!\right)\left(\! {{{A}}_1^{\rm{T}}{{{\hat{{A}\,}}}^{ - 1}}{{S}}'} \!\right),\;\;{\rm{s}}.{\rm{t}}.\; {{{A}}_{1}}{{A}}_1^{\rm{T}} = {{I}} (11) 其中，
{{{R}}'{\hat{{A}\,}}} 与{{{A}}_1} 相乘相当于旋转雷达视线，{{\hat{{A}\,}}^{ - 1}}{{S}}' 与{{A}}_1^{\rm{T}} 相乘相当于散射中心关于原点旋转。由于满足{{{A}}_1}{{A}}_1^{\rm{T}} = {{I}} ，根据式(12)计算矩阵{{{A}}_1} ：{\hat{{l}}}{{{A}}_1} = {{\hat{{l}}}_0} (12) 其中，
{{\hat{{l}}}_0} 是初始时刻雷达视线方向矢量，令{\hat{{A}\,}}{{\hat{{A}\,}}^{\rm{T}}} = \left( {\begin{array}{*{20}{c}} {{{{l}}_{{1}}}}&{{{{l}}_{{2}}}} \\ {{{{l}}_{{2}}}}&{{{{l}}_{{3}}}} \end{array}} \right) ，则{\hat{{l}}} = \left[ {\begin{array}{*{20}{c}} {{{l}}_{{1}}} \\ {{{l}}_{{2}}} \\ {{{l}}_{{3}}} \\ \end{array}} \right] 。结合矩阵奇异值分解所得的{{R}}' ,{{S}}' 以及估计出的矩阵{\hat{{A}\,}} ，可以得到等效雷达视线矩阵为{{R}} = {{{R}}'{\hat{{A}\,}}}{{{A}}_1} ，散射中心3维坐标矩阵为{{S}} = {{A}}_1^{\rm{T}}{{\hat{{A}\,}}^{ - 1}}{{S}}' 。通过上述航迹矩阵分解方法可以获得微动目标3维散射中心分布，进而实现空间微动目标高分辨3维成像。整体算法流程图如图2所示。
图 2 基于航迹矩阵分解的微动目标高分辨成像算法流程图Figure 2. The flow chart for high-resolution imaging of micro-motion targets based on trajectory matrix decomposition5. 实验与分析
本节采用仿真数据对所提算法进行验证。微动目标散射中心分布如图3(a)所示，该目标由9个散射中心组成。仿真参数为：带宽2 GHz，载频10 GHz，脉冲重复频率PRF=2000 Hz，观测时间为1 s。章动目标自旋角频率为1 Hz，锥旋角频率为0.4 Hz，摆动角频率为0.1 Hz，摆动幅度为5°。回波信号的信噪比为20 dB。
距离脉压后的目标回波如图3(b)所示，其中最底部曲线对应锥顶散射中心。9个散射中心航迹交叉点较多，基于1维距离像关联难度较大。采用watershed方法从图3(c)所示RID像中提取散射中心支撑域的结果如图3(d)所示，进而从中计算出各散射中心坐标，如图3(e)所示，其中蓝色圆圈表示散射中心支撑域轮廓，红色标记表示通过计算得到的散射中心坐标。由图可知，散射中心轮廓清晰，分割效果良好。基于RID像序列的距离-多普勒-慢时间3维关联结果如图3(f)所示，在距离-时间维的关联结果如图3(g)所示，其中不同颜色代表不同散射中心的航迹。由该图可知，该方法能够有效避免交叉点处关联错误等问题，获得准确的散射中心航迹关联结果。
利用Root-MUSIC的谱估计方法对航迹矩阵进行精估计，结果如图3(h)所示。最后，采用航迹矩阵分解法获得微动目标3维散射中心分布的结果如图4(a)所示，其中红色星号表示估计值，蓝色圆圈表示真实值。可以看出，成像结果与真实散射中心分布一致，从而证明了本文所提算法的有效性。等效雷达视线矩阵估计结果如图4(b)所示。
为测试所提成像方法的抗噪性能，在保持其他参数不变的条件下，给目标回波中分别加入信噪比为0 dB, 5 dB, 10 dB, 15 dB, 20 dB的高斯白噪声。在每个信噪比下做50次蒙特卡洛实验，并按照式(13)计算均方根误差(Root Mean Square Error, RMSE)：
{\rm{RMSE}} = \sqrt {\frac{{\displaystyle\sum\limits_{n = 1}^{{N_{\rm{m}}}} {{{\sum\limits_{p = 1}^{P} {\left[ {{{\left( {S_{{{x}},p}^n - T_{{{x}},p}^n} \right)}^2} + {{\left( {S_{{{y}},p}^n - T_{{{y}},p}^n} \right)}^2} + {{\left( {S_{{{z}},p}^n - T_{{{z}},p}^n} \right)}^2}} \right]} }\biggr/ {P}}} }}{{{N_{\rm{m}}}}}} (13) 其中，
n \in \left[ {1,{N_{\rm{m}}}} \right] ,{N_{\rm{m}}} 为蒙特卡洛实验次数，p \!\in\! \left[ {1,{P}}\, \right] ,{P} 为散射中心个数，\left[ {T_{{{x}},p}^n,T_{{{y}},p}^n,T_{{{z}},p}^n} \right] 和\left[ {S_{{{x}},p}^n,S_{{{y}},p}^n,S_{{{z}},p}^n} \right] 分别表示第n 次蒙特卡洛实验中第p 个散射点的真实坐标和估计坐标。最终，不同信噪比下的RMSE曲线如图5所示，可以看出，RMSE随着SNR的增加而降低。6. 结束语
针对传统参数化成像方法对复杂微动目标建模困难，未知参数求解运算量大等问题，本文提出一种基于RID图像序列的微动目标非参数化高分辨3维成像方法。该方法首先基于watershed法对RID图像进行分割提取散射中心，进而基于最近邻准则对散射中心航迹进行关联，接着通过Root-MUSIC方法实现航迹矩阵的精估计。最终，通过航迹矩阵分解实现微动目标的高分辨3维成像。该方法有效避免了参数化成像方法未知参数求解困难，易产生模型失配等不足。同时，2维关联方法克服了散射中心航迹交叉严重时，传统1维关联方法引起的关联误差，实现了复杂微动目标的高分辨3维成像。
在未来工作中，将研究低信噪比环境下的散射中心关联方法及非参数化微动目标高分辨3维成像方法，并进一步研究基于高分辨图像及等效雷达视线矩阵的微动目标特征提取及识别方法。
-
图 2 基于大气折射率廓线数据获得的对流层信号延迟
Figure 2. Tropospheric signal delays based on atmospheric refractive index profile data
图 3 不同波段、不同积累时间下对流层对GEO SAR成像影响的评估结果
Figure 3. Results of tropospheric GEO SAR imaging at different wave bands and different integrate times
图 6 电离层闪烁影响下成像指标参数分析
Figure 6. Analysis of imaging quality parameters in the presence of ionospheric scintillation
图 7 时空变电离层引入的干涉相位屏误差和评估
Figure 7. The interference phase screen errors and evaluation introduced by TSV ionosphere
表 1 GEO SAR对流层延迟的各阶时间变化情况
Table 1. Temporal variability of each order of GEO SAR tropospheric delay
项目 数值 \Delta \,{r_0} (m) 2.21 {q_1} (m/s) 2.52×10–4 {q_2} (m/s2) 2.71×10–7 {q_3} (m/s3) 1.64×10–13 
下载: 导出CSV
表 2 电离层TEC随时间变化情况
Table 2. Rate of ionospheric TEC changing with time
项目 数值 TEC0(TECU) 68.3 k1 (TECU/s) 0.0068 k2 (TECU/s2) 7.32×10–6 k3 (TECU/s3) 4.42×10–12
下载: 导出CSV
表 3 背景电离层对GEO SAR成像影响评估结果
Table 3. Results of GEO SAR imaging effected by background ionosphere
合成孔径时间(s) 距离向PSLR (dB) 方位向PSLR (dB) 方位向移位(m) 100 –13.23 –13.21 4.3 300 –13.25 –12.11 4.3 500 –13.22 –10.39 4.3
下载: 导出CSV
-
[1] Tomiyasu K. Synthetic aperture radar in geosynchronous orbit[C]. Proceedings of 1978 Antennas and Propagation Society International Symposium, Washington, DC, USA, 1978: 42–45 [2] Tomiyasu K and Pacelli J L. Synthetic aperture radar imaging from an inclined geosynchronous orbit[J]. IEEE Transactions on Geoscience and Remote Sensing, 1983, 21(3): 324–329. DOI: 10.1109/TGRS.1983.350561 [3] NASA and JPL. Global earthquake satellite system: A 20-year plan to enable earthquake prediction[EB/OL]. http://solidearth.jpl.nasa.gov/GESS/3123_GESS_Rep_2003.pdf, 2015, 9 [4] Edelstein W N, Madsen S N, Moussessian A, et al.. Concepts and technologies for synthetic aperture radar from MEO and geosynchronous orbits[C]. Proceedings of SPIE 5659, Enabling Sensor and Platform Technologies for Spaceborne Remote Sensing, Honolulu, Hawaii, 2005: 195–203 [5] Laurence Gray A, Mattar K E, and Sofko G. Influence of ionospheric electron density fluctuations on satellite radar interferometry[J]. Geophysical Research Letters, 2000, 27(10): 1451–1454. DOI: 10.1029/2000GL000016 [6] Meyer F, Bamler R, Jakowski N, et al. The potential of low-frequency SAR systems for mapping ionospheric TEC distributions[J]. IEEE Geoscience and Remote Sensing Letters, 2006, 3(4): 560–564. DOI: 10.1109/LGRS.2006.882148 [7] Meyer F J and Nicoll J. The impact of the ionosphere on interferometric SAR processing[C]. Proceedings of 2008 IEEE International Geoscience and Remote Sensing Symposium, Boston, MA, USA, 2008: II-391–II-394 [8] Meyer F J. Performance requirements for ionospheric correction of low-frequency SAR data[J]. IEEE Transactions on Geoscience and Remote Sensing, 2011, 49(10): 3694–3702. DOI: 10.1109/TGRS.2011.2146786 [9] Sun J P, Bi Y K, Wang Y P, et al.. High resolution SAR performance limitation by the change of tropospheric refractivity[C]. Proceedings of 2011 IEEE CIE International Conference on Radar, Chengdu, China, 2011 [10] Zhang F, Li G J, Li W, et al. Multiband microwave imaging analysis of ionosphere and troposphere refraction for spaceborne SAR[J]. International Journal of Antennas and Propagation, 2014, 2014: 913056. [11] Hanssen R F. Radar Interferometry: Data Interpretation and Error Analysis[M]. Dordrecht: Kluwer Academic Publishers, 2001 [12] Danklmayer A, Doring B J, Schwerdt M, et al. Assessment of atmospheric propagation effects in SAR images[J]. IEEE Transactions on Geoscience and Remote Sensing, 2009, 47(10): 3507–3518. DOI: 10.1109/TGRS.2009.2022271 [13] Hobbs S, Mitchell C, Forte B, et al. System design for geosynchronous synthetic aperture radar missions[J]. IEEE Transactions on Geoscience and Remote Sensing, 2014, 52(12): 7750–7763. DOI: 10.1109/TGRS.2014.2318171 [14] Bruno D and Hobbs S E. Radar imaging from geosynchronous orbit: Temporal decorrelation aspects[J]. IEEE Transactions on Geoscience and Remote Sensing, 2010, 48(7): 2924–2929. DOI: 10.1109/TGRS.2010.2042062 [15] Bruno D, Hobbs S E, and Ottavianelli G. Geosynchronous synthetic aperture radar: Concept design, properties and possible applications[J]. Acta Astronautica, 2006, 59(1–5): 149–156. DOI: 10.1016/j.actaastro.2006.02.005 [16] Ruiz Rodon J, Broquetas A, Guarnieri A M, et al. Geosynchronous SAR focusing with atmospheric phase screen retrieval and compensation[J]. IEEE Transactions on Geoscience and Remote Sensing, 2013, 51(8): 4397–4404. DOI: 10.1109/TGRS.2013.2242202 [17] Ruiz-Rodon J, Broquetas A, Makhoul E, et al. Nearly zero inclination geosynchronous SAR mission analysis with long integration time for earth observation[J]. IEEE Transactions on Geoscience and Remote Sensing, 2014, 52(10): 6379–6391. DOI: 10.1109/TGRS.2013.2296357 [18] Guarnieri A M, Rocca F, and Ibars A B. Impact of atmospheric water vapor on the design of a Ku band geosynchronous SAR system[C]. Proceedings of 2009 IEEE International Geoscience and Remote Sensing Symposium, Cape Town, South Africa, 2009 [19] Guarnieri A M, Tebaldini S, Rocca F, et al.. GEMINI: Geosynchronous SAR for earth monitoring by interferometry and imaging[C]. Proceedings of 2012 IEEE International Geoscience and Remote Sensing Symposium, Munich, Germany, 2012 [20] Monti Guarnieri A, Broquetas A, Recchia A, et al. Advanced radar geosynchronous observation system: ARGOS[J]. IEEE Geoscience and Remote Sensing Letters, 2015, 12(7): 1406–1410. DOI: 10.1109/LGRS.2015.2404214 [21] Monti Guarnieri A, Leanza A, Recchia A, et al. Atmospheric phase screen in GEO-SAR: Estimation and compensation[J]. IEEE Transactions on Geoscience and Remote Sensing, 2018, 56(3): 1668–1679. DOI: 10.1109/TGRS.2017.2766084 [22] Hu C, Long T, Zeng T, et al. The accurate focusing and resolution analysis method in geosynchronous SAR[J]. IEEE Transactions on Geoscience and Remote Sensing, 2011, 49(10): 3548–3563. DOI: 10.1109/TGRS.2011.2160402 [23] Hu C, Tian Y, Yang X P, et al. Background ionosphere effects on geosynchronous SAR focusing: Theoretical analysis and verification based on the BeiDou Navigation Satellite System (BDS)[J]. IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing, 2016, 9(3): 1143–1162. DOI: 10.1109/JSTARS.2015.2475283 [24] Dong X C, Hu C, Tian Y, et al. Experimental study of ionospheric impacts on geosynchronous SAR using GPS signals[J]. IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing, 2016, 9(6): 2171–2183. DOI: 10.1109/JSTARS.2016.2537401 [25] Tian Y, Hu C, Dong X C, et al.. Analysis of effects of time variant troposphere on geosynchronous SAR imaging[C]. Proceedings of 2016 IEEE International Geoscience and Remote Sensing Symposium, Beijing, China, 2016: 5051–5054 [26] 田野, 董锡超, 胡程. 对流层对地球同步轨道SAR成像的影响研究[J]. 信号处理, 2015, 31(12): 1562–1567. DOI: 10.3969/j.issn.1003-0530.2015.12.003Tian Ye, Dong Xi-chao, and Hu Cheng. Analysis of troposphere impacts on geosynchronous SAR imaging[J]. Journal of Signal Processing, 2015, 31(12): 1562–1567. DOI: 10.3969/j.issn.1003-0530.2015.12.003 [27] Tian Y, Hu C, Dong X C, et al. Theoretical analysis and verification of time variation of background ionosphere on geosynchronous SAR imaging[J]. IEEE Geoscience and Remote Sensing Letters, 2015, 12(4): 721–725. DOI: 10.1109/LGRS.2014.2360235 [28] Wang R, et al. Joint amplitude-phase compensation for ionospheric scintillation in GEO SAR imaging[J]. IEEE Transactions on Geoscience and Remote Sensing, 2017, 55(6): 3454–3465. [29] Hu C, Li Y H, Dong X C, et al. Performance analysis of L-band geosynchronous SAR imaging in the presence of ionospheric scintillation[J]. IEEE Transactions on Geoscience and Remote Sensing, 2017, 55(1): 3454–3465. DOI: 10.1109/TGRS.2016.2602939 [30] Li Y H, Hu C, Dong X C, et al. Impacts of ionospheric scintillation on geosynchronous SAR focusing: Preliminary experiments and analysis[J]. Science China Information Sciences, 2015, 58(10): 1–3. [31] Hu C, Li Y H, Dong X C, et al. Impacts of temporal-spatial variant background ionosphere on repeat-track GEO D-InSAR system[J]. Remote Sensing, 2016, 8(11): 916. DOI: 10.3390/rs8110916 [32] Kou L L, Wang X Q, Xiang M S, et al. Effect of orbital errors on the geosynchronous circular synthetic aperture radar imaging and interferometric processing[J]. Journal of Zhejiang University Science C, 2011, 12(5): 404–416. DOI: 10.1631/jzus.C1000170 [33] Kou L L, Xiang M S, Wang X Q, et al. Tropospheric effects on L-band geosynchronous circular SAR imaging[J]. IET Radar,Sonar&Navigation, 2013, 7(6): 693–701. [34] Ji Y F, Zhang Q L, Zhang Y S, et al. L-band geosynchronous SAR imaging degradations imposed by ionospheric irregularities[J]. Science China Information Sciences, 2017, 60(6): 060308. DOI: 10.1007/s11432-016-9064-1 [35] Li D X, Rodriguez-Cassola M, Prats-Iraola P, et al. Modelling of tropospheric delays in geosynchronous synthetic aperture radar[J]. Science China Information Sciences, 2017, 60(6): 060307. DOI: 10.1007/s11432-016-9065-1 [36] Tofsted D H. Turbulence simulation: Outer scale effects on the refractive index spectrum[R]. Technical Report ARL-TR-548. US Army Research Lab. NM, 2000 [37] Von Kármán T. Progress in the statistical theory of turbulence[J]. Proceedings of the National Academy of Sciences of the United States of America, 1948, 34(11): 530–539. DOI: 10.1073/pnas.34.11.530 [38] Tunick A D. The Refractive Index Structure Parameter/Atmospheric Optical Turbulence Model: CN2[M]. U.S. Adelphi: A.R. Laboratory, 1998 [39] Liu Z P, Hu C, Zeng T, et al. Improved secondary range compression focusing method in GEO SAR[C]. 2011 IEEE International Conference on Acoustics, Speech and Signal Processing (ICASSP), Prague, Czech Republic, 2011: 1373–1376 [40] Long T, Hu C, Ding Z, et al.. Geosynchronous SAR: System and Signal Processing[M]. Singapore: Springer, 2018 [41] Zhang D D, Chen Z Y, Dong X C, et al.. Simulating the impacts of ionospheric scintillation on geosynchronous SAR[C]. Proceedings of the 18th International Radar Symposium (IRS), Prague, Czech Republic, 2017 -
下载:
计量
- 文章访问数:
- HTML全文浏览量:
- PDF下载量:
- 被引次数: 0
