作者中心
专家审稿
责编办公
编辑办公
Research on the Technology of Airborne Multi-channel Wide Angle Staring SAR Ground Moving Target Indication
-
摘要: 机载广角凝视合成孔径雷达(WasSAR)是一种可对观测区域实施多角度长时间凝视成像探测的新兴SAR成像技术。将机载WasSAR成像与地面运动目标指示(GMTI)技术相结合，则可对重点区域内出现的地面运动目标实施持续成像跟踪监视，从而获取准确的动态感知信息。该文首先建立了机载多通道WasSAR动目标回波模型，分析了WasSAR动目标特性；然后，通过采用偏移相位校正和改进二维自适应校正方法，消除了载机姿态误差与通道非均衡的影响；在此基础上，提出了机载多通道WasSAR动目标检测跟踪算法，实现了复杂路况上行驶动目标的准确检测与跟踪；最后，提出了机载多通道WasSAR动目标行驶轨迹重构算法，实现了起伏路面下的动目标行驶轨迹精确重构。此外，该文中给出了作者团队利用自主研制机载多通道WasSAR-GMTI系统开展的外场飞行试验和实测数据处理结果，验证了地面运动目标持续跟踪监视的有效性和实用性，为后续开展更加深入的研究提供基础。Abstract: Airborne Wide angle staring Synthetic Aperture Radar (WasSAR) is a novel SAR imaging technique that enables multiangle and long-time staring imaging of an observation region. The combination of WasSAR and Ground Moving Target Indication (GMTI) facilitates the long-time tracking of moving targets in key areas, thereby obtaining dynamic sensing information. Herein, we initially constructed a moving target echo model for airborne multi-channel WasSAR, followed by an analysis of the characteristics of a moving target during WasSAR imaging. Further, group phase shift calibration and modified A2DC methods are proposed to mitigate the influence of platform attitude errors and channel imbalances. In accordance with this, an extended airborne multi-channel WasSAR ground moving target detection and tracking method is proposed, which achieves accurate detection and tracking of targets moving on complex roads. Finally, a trajectory reconstruction method for moving targets in airborne multi-channel WasSAR-GMTI is presented, demonstrating accurate trajectory reconstruction for targets on complex roads. Moreover, a flight experiment using our independently developed airborne multi-channel WasSAR-GMTI system, along with the associated real data processing results, is presented. These results confirm the validity and practicability of airborne WasSAR-GMTI for moving target dynamic surveillance and serve as a foundation for future research.
-
1. 引言
潮间带区域是一种位于陆地与海洋之间的特殊海岸区域，在低潮时露出水面而在高潮时被水淹没[1]，其构成了一个独特的湿地生态系统，提供了多种植被必需的生长环境，是水产养殖、海洋环境观测、海洋经济开发以及海岸防御的重点区域。为了更好地开发利用潮间带区域，利用遥感对潮间带进行观测和研究的工作已得到开展[2–5]。如何基于遥感图像自动准确地将潮间带区域与水域及其他陆地区域分开，并能够准确区分潮间带区域不同植被或养殖类型，是潮间带区域监测所关心的问题。合成孔径雷达(Synthetic Aperture Radar, SAR)尤其是全极化SAR，由于其全天时全天候工作的优势和分类能力强的特点，在该方面的应用得到了关注[6–8]。
针对极化SAR数据，通过特征值分解以及散射分析等手段，可以从中提取出大量的多极化特征来描述地物的散射信息[9,10]。这些特征为极化SAR数据的应用提供了新的视角。近年来，基于极化特征的研究有了突破性的进展，广泛地应用在地物分类、典型目标检测提取等领域[11–13]。相比于其他特征，极化熵被证实可以更好地描述地物的极化相关的纹理信息(polarization-dependent variation of texture)[14]和表面粗糙度[15]，并且有较好的目标识别能力[16]。因此，极化熵广泛地应用在地物分类、海岸带典型目标检测[3,17]等领域。比较典型的有Cloude和Pottier1997年提出的基于极化熵的地物分类框架[16]以及Eun-Sung等人利用极化熵基于GEV分布实现对水下养殖区域的检测等[3]。
另外，从对极化SAR数据的统计特性的分析出发，已有研究将有限混合模型(Finite Mixture Models, FMM)等理论引入到基于极化SAR地物分类应用中[18]。从传统的混合Wishart模型[18]到混合U分布模型[19]，基于FMM的极化SAR地物分类方法在城区、农田以及海洋等区域达到较好的分类效果。
在潮间带区域，其在退潮后露出水面的部分主要由淤泥滩、浅留水体以及水产养殖等部分组成，而上述部分的电磁散射强度与水面差别不大，因此在极化SAR图像中表现为类似的散射强度和散射机制，不易区分。现有针对极化SAR数据进行地物分类的算法在潮间带这种特殊的地物场景中并不能得到很好的分类结果。现有基于SAR数据对于潮间带区域也有较多的研究，主要针对有潮海域的地貌变动监控[7,20–22]、土壤湿度分析[1,8]以及特定地物的散射特性分析[3,17,23,24]等。比较典型的有Li Zhen等人针对1999–2006年波罗的海海域的潮间带地貌演变研究[7,20–22]以及Hoonyol Lee等人对于韩国南部海域潮间带的土壤湿度变化研究[1]等。
本文针对潮间带的特点，基于极化特征和FMM提出了一种新的地物分类方法。首先，本文针对4种典型的极化特征进行分析和筛选，得到一组最适合描述潮间带区域的多极化特征：极化熵和反熵。在此基础上，基于对极化熵图像的统计特性分析和极值理论，本文优选GEV分布对潮间带区域的极化熵分布进行描述[20,25]。随后构建了一种基于GEV混合模型(GEV Mixture Model, GEVMM)在极化熵特征中对潮间带地物进行标记。最后结合反熵的信息实现潮间带地物分类。本文利用RadarSAT-2获取的上海崇明岛东滩区域数据进行实验验证。与经典多极化SAR分类方法结果及实地考察真值的对比结果表明，本文方法能够有效地将潮间带区域与水域进行区分，并对潮间带内部不同场景类型进行分类，效果优于Wishart-H/A/α等经典方法。
本文的结构如下：第2节详细介绍本文提出的分类方法，第3节基于实测数据对方法进行验证，并对结果进行分析讨论，第4节总结全文。
2. 极化熵结合GEV模型的分类方法
2.1 信号建模与极化特征
在全极化SAR模式中，SAR系统以水平极化和垂直极化两种模式的组合来发射和接收电磁波[26]。极化散射矩阵为：
S=[SHHSHVSVHSVV] (1) 一般来说， SHV=SVH，因此Pauli散射矢量定义为：
k=1√2[SHH+SVVSHH−SVV2SHV]T (2) 其中，T表示转置。基于k可以得到极化相干矩阵T3：
T3=1LL∑n=1kk∗T (3) 其中，L为视数，上标 ∗代表了复共轭。从极化散射矩阵S和极化相干矩阵T3中可以提取出来很多极化特征，其中4个最典型的多极化特征(Span, H, A, α)如下。
极化总功率Span是一个典型的极化特征，其定义可以从极化散射矩阵S得到：
Span=|SHH|2+|SHV|2+|SVH|2+|SVV|2 (4) Span影像为4个极化通道图像的功率之和，其中包含了各通道图像的信息。
从极化相干矩阵T3中获取极化特征一般通过特征值分解，其中比较典型的为 H/A/α分解。 H/A/α已经在极化图像分类领域广泛使用，其定义为：
T3=3∑i=1λiT3i=3∑i=1λiuiui∗T (5) 其中， λi是T3矩阵的特征值，代表着第i个归一化分量 T3i的比重，其概率可以用Pi来表示：
Pi=λi/λi3∑i=1λi3∑i=1λi (6) 极化熵和反熵是 H/A/α分解后得到的两个重要特征，可以通过对 λi进行计算得到：
H=−3∑i=1PilogPi (7) A=λ2−λ3λ2+λ3 (8) 其中，H是极化熵，A是反熵，计算时特征值按照 λ1>λ2>λ3来排序。其中极化熵H描述了散射的随机性，其取值范围为0到1。当极化熵值较低时，可以认为只有一种主导散射机制；当极化熵值较高时，此时表现出较为强烈的去极化特性，散射机制趋向随机；当极化熵接近于1时，散射机制接近于噪声。反熵A是对极化熵H的一个有效补充，是对第2个和第3个特征值的差异的描述，并且对极化熵值大于0.7的高熵区域有着较好的识别性能[26]。一般来说，潮间带的大多数区域对于反熵的响应较弱，只有部分高熵区域，例如露出水面的金属隔离网，才有较高的反熵值。
另外一个重要的 H/A/α分解特征是平均散射角，可以通过特征向量计算得到：
α=3∑i=1Picos−1(|ui(1)|) (9) 其中，α为平均散射角，其值域范围为0° and 90°[26]。
2.2 极化特征比较与选择
在潮间带区域，退潮残留的水体、淤泥滩和水下农场等构成了复杂的地貌。残留水体和淤泥滩均有相对光滑的表面，且绝大部分的水下农场养殖区域可以被称为“浅层水域”，仅有少量的隔离网和部分植被的枝叶露出水面。因此，该区域在各个通道的雷达功率图像中均会表现出和海水相似的图像灰度值，如图1(a)所示。所以Span图虽然综合了全极化SAR的4个通道功率图像的幅度信息，但是这种简单的功率和并不能增强潮间带区域和海水的区分能力。
潮间带混杂的地物也导致了其散射机制的随机性变化较大，而极化熵H和反熵A均是对散射随机性的描述，因此这两个特征均有可能用来描述在潮间带区域。潮间带区域的极化熵和反熵示例分别如图1(b)和图1(c)所示。可以看出，相比于Span，极化熵可以较好地区分海洋和潮间带区域，并且可以区分潮间带内部地物类型。而反熵则可以明显将高熵和低熵区域区分开来。因此对于潮间带区域，理论上来说极化熵是一种更普适的特征，而反熵更适合用来区分一些高熵区域。
平均散射角α则是对潮间带区域的散射机制的直接描述。较为光滑的表面会产生较低的α值，随着α值的增大，散射机制依次从表面散射，体散射到二次散射过渡。简而言之，对于潮间带区域，平均散射角α可以认为是一种对表面粗糙度的衡量。淤泥滩和露出水面的水下植被与水面相比拥有不同粗糙度的表面，从而在α特征图中会有不同的表现，如图1(d)所示。但是，由于在潮汐作用下潮间带区域的地物表面相对平滑，其散射机制大多为不同程度的布拉格表面散射，因此这种粗糙度带来的差异性不足以用来对潮间带区域的地物分类。
基于上述分析，本文优选极化熵H和反熵A来进行潮间带的分类，其中H为主要特征，A则作为补充。
2.3 极化熵的统计特性分析与混合GEV模型
极化熵是一种对极化散射机制随机性的衡量。然而，对于中等分辨率(如Radarsat-2全极化模式分辨率为8 m)图像，一个像素点中可能包含多种地物，即有着较强的散射随机性，意味着比较高的极化熵值，当散射随机性很强时，极化熵值趋近于1，接近噪声水平。在潮间带中，虽然总体上来说其极化熵值不是很高，但是部分露出水面的植被以及金属隔离网区域有着较强的散射随机性会产生接近于1的极高熵值。这些高熵值在极化熵图像的统计直方图中导致一种long-tails现象。long-tails现象是风险预测或者灾害预警等极端事件研究中的重点，包括“左拖尾(厚尾)”和“右拖尾(薄尾)”[27]。对于这两种现象，传统的统计分布模型例如高斯分布模型或者Gamma分布模型难以给出好的拟合效果[25]。而常用的一些描述long-tails现象的模型，包括广义Gamma分布和k分布，均对“右拖尾”有较好的描述能力，但是不能很好地适应“左拖尾”。而且，k分布是基于乘性噪声模型提出，仅适用于均匀区域和弱纹理区域，对于极化熵图像不太适用[26]。极值理论方面的研究表明，GEV分布可以很好地描述有极端事件发生的随机过程。在近年来的研究中GEV分布已经被大量地用在SAR图像对特殊目标的分析中[3,20,25]。其中Eun-Sung的论文中对比了广义Gamma分布族中的Weibull分布和log-normal分布，证明了GEV分布能够更好地描述水下养殖等类似区域[3]。因此本文采用GEV来描述潮间带区域的极化熵图像。GEV分布是一族连续的极值概率分布，并且具有有限分布和稳定的特性。GEV分布的概率密度函数如式(10)所示：
gev(x,μ,σ,ξ)={e−(ξx−μσ+1)−1/ξ(ξx−μσ+1)−1−1/ξσ,ξx−μσ+1>0, σ>0, ξ≠0e−(x−μσ)−e−(x−μσ)σ,σ>0, ξ=0 (10) GEV分布的参数有3个：即形状参数 ξ，尺度参数σ和位置参数μ。图2中给出了对于不同 ξ下的GEV分布的形态，其中 ξ分别为–0.5, 0, 0.5, σ和μ均为0.1和0.5。当 ξ等于0时为Gumbel分布，当 ξ大于0时为Weibull分布，而当 ξ小于0时为Frechet分布。其中Weibull描述的是一种“右拖尾分布”，Frechet描述的是一种“左拖尾分布”[27]，而Gumbel分布为两者之间的过渡。因此，GEV分布对long-tails有着广泛的适应性以及较好的描述能力。
通过上面的分析可以得出，GEV分布能够较好地描述潮间带区域的极化熵图像。由于单个概率分布只有一个峰值，而对于潮间带区域，多种地物必然会导致多个峰值的出现。因此基于GEV分布，可以通过建立一个有限混合模型(Finite Mixture Model, FMM)来实现对潮间带的地物分类。有限混合GEV模型(GEV Mixture Models, GEVMM)n个GEV分布的线性叠加，其定义如式(11)所示：
f(x;ξ,σ,μ)=n∑i=1aifi(x;ξi,σi,μi) (11) 其中，n为GEV分布的个数，在潮间带地物分类中代表着潮间带的地物种类，fi(·)代表着第i个种地物所服从的GEV分布。
2.4 基于GEV模型的分类方法
对于潮间带区域的极化熵图像H的地物分类可以分为两个步骤：首先是GEV混合模型的计算，然后将模型映射到图像空间中实现分类。其中最重要的一步是GEVMM的计算。本文以EM算法为基础，来实现对GEVMM的构建。
EM算法是一种极大似然估计算法，以迭代的方式来实现对模型的估计。算法的每一次迭代都是由E步(求期望)和M步(将期望最大化)构成[28]。在本文中，首先给定类别个数和参数的初始值，然后在E步，依照当前模型参数，计算分模型i对图像中每个点xj的响应度 ˆγji，其计算公式为：
ˆγji=aifi(xj;ξi,σi,μi|θi)n∑iaifi(xj;ξi,σi,μi|θi) (12) 其中，j为样本标记，i为类别标记，式中θi指在第i个分模型。在M步中，我们先根据响应度大小对图像像素进行归类，然后用最小二乘法估计出每个分模型的参数[29]。实验中最小二乘法估计通过MATLAB函数gevfit实现。如此迭代，直到参数收敛。
另一方面，传统的FMM进行模型估计时分模型个数是固定的，难以适应不同的地物类型的特征值分布。因此，本文采用了自适应的的GEVMM估计方法。在迭代的过程中，如果对某一个类别产生最大响应的点的数目少于某个阈值时，该类别即认为无意义并删除。在本文的实验中，这个阈值被设为50。迭代收敛后得到GEVMM以及其中的各个分量，然后通过计算图像中的每个点对模型各分量的响应度，得到一个分类标记图。
反熵作为极化熵的补充，其在高于0.7时采用较好的识别能力[26]。因此可以通过0.7这个阈值对反熵图像进行阈值分割，实现二值分类。基于此，本文通过将两个特征得到的结果进行融合，更好地实现潮间带的地物分类。假设极化熵得到的M类集合为 {A1,A2,⋅⋅⋅,AM}，反熵得到的两类为 {B1,B2}，融合算法要处理的集合为各类的交集，即：
C={Ai∩Bj,Ai∩Bj≠∅,i∈[1,M],j∈[1,2]} (13) 对于C中的每一个像素x，分别计算其对Ai与Bj的从属度 mx(Ai)和 mx(Bj)：
mx(Ai)=|H(x)−xAi| (14) mx(Bj)=|A(x)−xBj| (15) 其中，H(x)为极化熵在x位置的值， xAi为极化熵图像中 Ai类所有像素的均值，即类别中心。同理， A(x)为反熵在x位置的值， xBj为反熵图像中 Bj类所有像素的均值。本文中采用如下准则来对交集C进行判决：如果 ∑x∈Cmx(Ai)>∑x∈Cmx(Bj)，那么C属于Ai类，反之C属于Bj类。这样就实现了分类结果的融合。
GEVMM的构建流程以及图像分类算法流程如图3所示。
3. 实验与分析
基于本文的方法，利用Radarsat-2卫星的全极化数据进行了实验验证，实验数据的分辨率为8 m。研究区域位于上海崇明岛东滩潮间带，其卫星数据获取时间为2015年5月16日，此时的潮水处于低潮状态。图4描述了实验区域，其中图4(a)给出了整幅数据的Pauli伪彩色图像，其中的黄色矩形框代表着我们选定的实验区域，大小为20.0 km×5.6 km。数据获取时的实验区域的实时状态如图4(b)所示，从近到远依次是低矮植被区域，淤泥滩以及露出水面的水产养殖场的隔离网。
实验中，本文首先基于实测数据对潮间带区域的极化特征进行目视解译与定量衡量来佐证前面的分析；然后基于所示的流程对潮间带进行地物分类；最后利用PolSAR Pro软件对研究区域进行Wishart-H/A/α分类，并和本文方法的实验结果进行对比和分析。
3.1 基于实测数据的极化特征分析
图5给出了图4(a)黄框所示的研究区域的极化特征。基于上面的分析，下面结合实地考察结果以及图4(a)中所示的Pauli伪彩色图对4个特征进行目视解译。在图5(a)所示的Span影像中，我们可以识别出潮间带的大致轮廓以及金属隔离网等，但是潮间带内部细节的识别能力较差，并且潮间带和海洋的边界比较模糊。图5(b)给出了极化熵图像，从图中我们可以识别出海面，以及潮间带中的淤泥滩，草地以及暗色的残留水体等，并且和海面的区分较为显著，但是对金属隔离网的识别能力较差。图5(c)中的反熵可以很好地识别金属隔离网，但是无法识别出其他地物。而对于图5(d)，平均散射角可以对潮间带区域中的地物进行一定的识别，但是不能很好地区分海面和潮间带中的水产养殖。对比Span、极化熵、反熵和平均散射角4张图，可以发现极化熵可以对除去金属隔离网之外的地物有着最好的视觉判别而反熵作为补充可以很好地识别出金属隔离网区域。因此，极化熵结合反熵，可以给出最好的视觉判别效果。
为了进一步地佐证，本文采用Michelson对比度准则来衡量上述图像的类间对比度(between-region contrast)[30]。 Michelson准则是一种经典的图像衡量标准，并且广泛地用在SAR图像评价中[3,31]。其定义如下：
C=Fmax−FminFmax+Fmin (16) 其中，Fmax和Fmin分别为图像的最大值和最小值。对比度值C位于[0, 1]区间。一般来说，C越高，对比度越强，区分能力也越强。表1给出了4个极化特征的类间对比度，从表中可以看出，反熵拥有最高的对比度值，极化熵的对比度紧随其后且高于其他两个特征。从前面的分析中可以得知，反熵对高熵部分响应强烈并且抑制了低熵部分的值，导致了较高的对比度值，也确保了其对潮间带的高熵部分的识别能力。从图5(c)可以看出，反熵对于具有高熵值的金属隔离网有着很好的区分能力，但是对其他低熵区域，区分能力很差。极化熵有着第二高的对比度值，从图5(b)中也可以看出极化熵对整个潮间带区域的细节有着很好的描述以及较好的区分能力。平均散射角α的对比度低于极化熵，从图5(d)中也可以看出其对潮间带的细节有一定的描述能力，但是无论是边缘的清晰度还是Michelson对比度都明显弱于极化熵。而Span图像在视觉判别能力和对比度方面表现最差。综上所述，极化熵和反熵的组合可以对潮间带区域地物进行很好的区分。
表 1 各极化特征的Michelson类间对比度Table 1. Michelson between-region contrast of different features极化特征 类间对比度 Span 0.4092 Entropy 0.7703 Anisotropy 0.9959 α 0.6757 3.2 基于GEVMM的极化熵统计建模
基于图5(b)所示的极化熵图像，我们按照图3所示的流程进行GEVMM建模和地物分类。在GEVMM构建中，根据实地考察结果结合极化熵图像我们初始指定模型的分量数目为8，经过迭代收敛后得到的模型中有5个分量，图6(f)给出了最终的GEVMM和极化熵直方图以及模型中各分量的图示。为了确保GEVMM的优越性，本文用Gamma分布和log-normal分布分别来拟合GEVMM得到的标记图中对应的每种类别的归一化直方图，与GEVMM的各个分量的拟合效果进行对比，如图6(a)–图6(e)所示。图6(a)–图6(e)展示了拟合对比结果，其中蓝色部分为各类别的直方图，绿线代表GEVMM的对应分量拟合结果，黑线代表Gamma模型拟合结果，红线代表了log-normal模型拟合结果。可以看出，GEV分布可以很好地拟合图6(d)和图6(e)两图的“左拖尾现象”，明显优于其他两种模型。同样，图6(b)中GEV分布的拟合效果也明显优于其他两种分布。
图 6 GEVMM及其各分量与Gamma分布和log-normal分布的对比：(a)–(e)分别为5个分量与Gamma分布，log-normal分布以及对应标记区域的直方图的对比，其中蓝色区域为归一化直方图，绿线是GEV拟合结果，黑线是Gamma拟合结果，红线是log-normal拟合结果，(f)给出了GEVMM及其各个分量与研究区域直方图的对比结果，其中蓝线为归一化直方图，红线为GEVMM，绿线为GEVMM的各个分量Figure 6. Fitness comparison among GEV distribution and Gamma distribution and Log-normal distribution of each component in GEVMM: (a)–(e) represent the five components of the GEVMM and the fitting results by the Gamma distribution and log-normal distribution for the histograms, which are marked as blue, the green lines represent the GEV fitting results, the black lines represent the most fitted Gamma distribution and the red lines represent the most fitted log-normal distribution, (f) shows the five components of GEVMN as green lines and the respective histograms as blue lines, the red line represents the final model为了进一步确认GEVMM的准确性，本文采用赤池信息量准则(Akaike Information Criterion, AIC)对各个模型进行衡量[32]。AIC准则的定义如下：
AIC=2k−2Ln (17) 其中，k为参数个数，n为样本个数，L为对数相似度(log likelihood)。对数相似度的定义如下：
L=−n2ln(2π)−n2ln(ssen)−n2 (18) 其中，sse为直方图和拟合结果之间的残差和。对于AIC准则来说，k越小或者L越大，AIC的值越小，模型的正确性也越高。其中k越小，模型越简洁，而L越大，模型越准确。表2给出了3种分布在每个类别对应的AIC值，可以看出GEV分布明显优于其他两种分布。结合图6和表2，可以得出GEV分布模型有着最好的拟合能力。另外，最终得到的GEVMM和统计直方图之间的AIC值为4.0661，也说明了模型的有效性。
表 2 GEV分布，Gamma分布和log-normal分布在每种类别中的拟合结果的AIC值Table 2. The AIC values of the fitting results between the GEV distribution, the Gamma distribution and log-normal distributionAIC 1 2 3 4 5 GEV 6.3988 4.0095 4.1105 4.7009 6.2093 Gamma 11.3878 9.3573 7.6928 9.0210 8.3989 log-normal 11.3878 9.3573 7.6928 9.0213 8.3997 3.3 地物分类结果与分析
图7展示了地物分类实验结果，其中图7(a)给出了研究区域的Pauli伪彩色图。基于极化熵的GEVMM得到的图像分类标记图如图7(b)所示，通过融合图7(c)给出的反熵的分割结果，得到的最终分类结果如图7(d)所示。图7(e)给出了Wishart-H/A/α方法的分类结果。研究区域真值图如图7(f)所示，其中蓝色为海洋区域，浅蓝色和浅绿色区域分别为潮间带中的水下养殖与残留水体，黄色和棕黄色代表植被散射区域和淤泥滩，深红色代表金属隔离网。
对比发现极化熵得到的结果可以很好地识别出潮间带区域的大多地物类型，与真值较为一致，但是该图中地物只是分为5类，不能很好地识别出图7(f)中深红色的金属隔离网区域。而在图7(c)的反熵的阈值分割结果可以较好地提取出金属隔离网区域。图7(d)中融合后得到的分类结果和真值图相比有较为一致分类效果，并优于图7(b)和图7(c)的结果。而图7(e)所示的Wishart-H/A/α的分类结果和真值图有着较大的不同。在研究区域，图7(e)只能将金属隔离网区域和其他地物区分，对于潮间带的其他5类，均不具备较好的区分度，并且潮间带区域和海洋区域未能进行很好地区分。
针对本文提出的方法和Wishart-H/A/α分类方法，实验中采用Kappa系数来衡量两种方法的分类效果。结合真值数据，两种方法的Kappa系数分别为0.9019和0.3662，总体分类精度分别为0.8912和0.3712，也验证了本文方法的准确性。可见，在潮间带区域的地物分类中，本文的方法要优于Wishart-H/A/α分类。
4. 结束语
针对潮间带区域的地物分类问题，本文提出了结合极化熵和反熵的潮间带地物分类方法。基于极化熵图像，本文对潮间带区域建立GEVMM。在此基础上通过和反熵的信息融合得到最终的分类结果。实验中利用C波段Radarsat-2全极化数据实现了对上海崇明岛东滩潮间带区域的地物分类。对比真值图，本文方法可以给出较好目视结果并优于Wishart-H/A/α方法。最后的定量分析证明本文方法具有较好的准确性。
-
图 5 当
vx=0.6m/mss 时的运动目标成像Figure 5. The imaging results of moving target (
vx=0.6m/mss )
图 6 当
vy=0.2m/mss 时的运动目标成像Figure 6. The imaging results of moving target (
vy=0.2m/mss )
图 13 通道误差校正后的RD域干涉相位结果
Figure 13. The results of interferometric phase after channel error calibration
图 14 通道误差校正后的图像域DPCA检测结果与ATI检测结果
Figure 14. The detection results of DPCA and ATI after channel error calibration
图 15 机载多通道WasSAR动目标检测与跟踪处理流程图
Figure 15. The processing flow of moving target detection and tracking in WasSAR
图 16 机载WasSAR子孔径回波、CSI及RELAX检测结果
Figure 16. The sub-aperture echo of airborne WasSAR, CSI and RELAX detection results
图 17 不同跟踪算法下的RD域多目标轨迹对比图
Figure 17. The comparison results of multiple target trajectory in RD domain by different tracking algorithms
图 21 合作车辆目标的多子孔径轨迹重构结果
Figure 21. The multi channel trajectory reconstruction of cooperative targets
图 23 起伏路面下的运动目标轨迹重构处理流程图
Figure 23. The processing flow of 3-D moving target trajectory reconstruction
图 24 起伏路面下运动目标的原始观测几何与等效观测几何
Figure 24. The real geometry and equivalence of 3-D moving target
图 26 观测场景中的拱桥与合作车辆
Figure 26. The pictures of arch bridge and cooperative pickup truck of observation scene
图 27 起伏路面动目标三维速度估计结果与三维位置重构结果
Figure 27. The three-dimensional velocity estimation and trajectory reconstruction of moving target in the three-dimensional field
表 1 机载WasSAR系统仿真参数
Table 1. The simulated parameters of airborne WasSAR system
参数 数值 参数 数值 飞行半径 ra 2000 m 工作频段 Ku波段 飞行高度H 2000 m 距离分辨率 ρr 0.167 m 飞机速度 va 180 km/h 脉冲重复频率 PRF 3125 Hz 
下载: 导出CSV
-
[1] CARRARA W G, GOODMAN R S, and MAJEWSKI R M. Spotlight Synthetic Aperture Radar: Signal Processing Algorithms[M]. Boston: Artech House, 1995: 13–78. [2] 徐丰, 王海鹏, 金亚秋. 深度学习在SAR目标识别与地物分类中的应用[J]. 雷达学报, 2017, 6(2): 136–148. doi: 10.12000/JR16130XU Feng, WANG Haipeng, and JIN Yaqiu. Deep learning as applied in SAR target recognition and terrain classification[J]. Journal of Radars, 2017, 6(2): 136–148. doi: 10.12000/JR16130 [3] 王广学. UWB SAR叶簇隐蔽目标变化检测技术研究[D]. [博士论文], 国防科学技术大学, 2011.WANG Guangxue. Foliage-concealed target change detection for UWB SAR[D]. [Ph.D. dissertation], National University of Defense Technology, 2011. [4] 李田, 程晓, 关真富, 等. 基于SAR数据的南极冰山分布监测[J]. 南京信息工程大学学报: 自然科学版, 2020, 12(2): 231–235. doi: 10.13878/j.cnki.jnuist.2020.02.011LI Tian, CHENG Xiao, GUAN Zhenfu, et al. Investigation of Antarctic iceberg distribution based on SAR images[J]. Journal of Nanjing University of Information Science & Technology: Natural Science Edition, 2020, 12(2): 231–235. doi: 10.13878/j.cnki.jnuist.2020.02.011 [5] XIANG Deliang, TANG Tao, BAN Yifang, et al. Man-made target detection from polarimetric SAR data via nonstationarity and asymmetry[J]. IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing, 2016, 9(4): 1459–1469. doi: 10.1109/jstars.2016.2520518 [6] Nepal earthquake displacement[EB/OL]. http://www.esa.int/ESA_Multimedia/Images/2015/04/Nepal_earthquake_displacement, 2015. [7] HUANG Yan, LIAO Guisheng, XU Jingwei, et al. GMTI and parameter estimation via time-Doppler chirp-varying approach for single-channel airborne SAR system[J]. IEEE Transactions on Geoscience and Remote Sensing, 2017, 55(8): 4367–4383. doi: 10.1109/tgrs.2017.2691742 [8] DA SILVA A B C, JOSHI S K, BAUMGARTNER S V, et al. Phase correction for accurate DOA angle and position estimation of ground-moving targets using multi-channel airborne radar[J]. IEEE Geoscience and Remote Sensing Letters, 2022, 19: 4021605. doi: 10.1109/LGRS.2022.3144735 [9] KIRSCHT M. Detection and imaging of arbitrarily moving targets with single-channel SAR[J]. IEE Proceedings - Radar, Sonar and Navigation, 2003, 150(1): 7–11. doi: 10.1049/ip-rsn:20030076 [10] SUWA K, YAMAMOTO K, TSUCHIDA M, et al. Image-based target detection and radial velocity estimation methods for multichannel SAR-GMTI[J]. IEEE Transactions on Geoscience and Remote Sensing, 2017, 55(3): 1325–1338. doi: 10.1109/tgrs.2016.2622712 [11] 贺雄鹏. 阵列雷达宽幅地面运动目标检测方法研究[D]. [博士论文], 西安电子科技大学, 2020.HE Xiongpeng. An investigation of ground moving target detection approaches for wide-swath array radar[D]. [Ph.D. dissertation], Xidian University, 2020. [12] SCARBOROUGH S M, CASTEEL JR C H, GORHAM L, et al. A challenge problem for SAR-based GMTI in urban environments[C]. SPIE Algorithms for Synthetic Aperture Radar Imagery XVI, Orlando, USA, 2009: 73370G. [13] SOUMEKH M. Synthetic Aperture Radar Signal Processing with MATLAB Algorithms[M]. New York, USA: Wiley, 1999: 486–551. [14] CHEN Jingwei, AN Daoxiang, WANG Wu, et al. Extended polar format algorithm for large-scene high-resolution was-SAR imaging[J]. IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing, 2021, 14: 5326–5338. doi: 10.1109/JSTARS.2021.3081515 [15] JIA Gaowei, BUCHROITHNER M F, CHANG Wenge, et al. Fourier-based 2-D imaging algorithm for circular synthetic aperture radar: Analysis and application[J]. IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing, 2016, 9(1): 475–489. doi: 10.1109/JSTARS.2015.2502430 [16] HORN R, NOTTENSTEINER A, REIGBER A, et al. F-SAR—DLR’s new multifrequency polarimetric airborne SAR[C]. 2009 IEEE International Geoscience and Remote Sensing Symposium, Cape Town, South Africa, 2009: II-902–II-905. [17] PINHEIRO M, PRATS P, SCHEIBER R, et al. Tomographic 3D reconstruction from airborne circular SAR[C]. 2009 IEEE International Geoscience and Remote Sensing Symposium, Cape Town, South Africa, 2009: III-21–III-24. [18] PONCE O, PRATS-IRAOLA P, PINHEIRO M, et al. Fully polarimetric high-resolution 3-D imaging with circular SAR at L-band[J]. IEEE Transactions on Geoscience and Remote Sensing, 2014, 52(6): 3074–3090. doi: 10.1109/tgrs.2013.2269194 [19] 安道祥, 陈乐平, 冯东, 等. 机载圆周SAR成像技术研究[J]. 雷达学报, 2020, 9(2): 221–242. doi: 10.12000/JR20026AN Daoxiang, CHEN Leping, FENG Dong, et al. Study of the airborne circular synthetic aperture radar imaging technology[J]. Journal of Radars, 2020, 9(2): 221–242. doi: 10.12000/JR20026 [20] 丁金闪. 视频SAR成像与动目标阴影检测技术[J]. 雷达学报, 2020, 9(2): 321–334. doi: 10.12000/JR20018DING Jinshan. Focusing algorithms and moving target detection based on video SAR[J]. Journal of Radars, 2020, 9(2): 321–334. doi: 10.12000/JR20018 [21] 林赟, 洪文. 圆迹合成孔径雷达成像与应用[M]. 北京: 电子工业出版社, 2020: 1–30.LIN Yun and HONG Wen. Circular Synthetic Aperture Radar Imaging and Its Applications[M]. Beijing: Publishing House of Electronics Industry, 2020: 1–30. [22] 林赟, 张琳, 韦立登, 等. 无先验模型复杂结构设施SAR全方位三维成像方法研究[J]. 雷达学报, 2022, 11(5): 909–919. doi: 10.12000/JR22148LIN Yun, ZHANG Lin, WEI Lideng, et al. Research on full-aspect three-dimensional SAR imaging method for complex structural facilities without prior model[J]. Journal of Radars, 2022, 11(5): 909–919. doi: 10.12000/JR22148 [23] LIN Yun, HONG Wen, TAN Weixian, et al. Interferometric circular SAR method for three-dimensional imaging[J]. IEEE Geoscience and Remote Sensing Letters, 2011, 8(6): 1026–1030. doi: 10.1109/lgrs.2011.2150732 [24] 王武. 机载圆周SAR成像及地面动目标指示技术研究[D]. [博士论文], 国防科技大学, 2019.WANG Wu. Study on circular SAR imaging and ground moving target indication[D]. [Ph.D. dissertation], National University of Defense Technology, 2019. [25] UYSAL F, MURTHY V, and SCARBOROUGH S M. Blind phase calibration for along-track interferometry: Application to gotcha data set[C]. SPIE Algorithms for Synthetic Aperture Radar Imagery XXI, Baltimore, USA, 2014: 90930O. [26] PAGE D, OWIRKA G, NICHOLS H, et al. Detection and tracking of moving vehicles with gotcha radar systems[J]. IEEE Aerospace and Electronic Systems Magazine, 2014, 29(1): 50–60. doi: 10.1109/MAES.2014.130075 [27] HERSEY R K and CULPEPPER E. Radar processing architecture for simultaneous SAR, GMTI, ATR, and tracking[C]. 2016 IEEE Radar Conference, Philadelphia, USA, 2016: 1–5. [28] DEMING R W, MACINTOSH S, and BEST M. Three-channel processing for improved geo-location performance in SAR-based GMTI interferometry[C]. SPIE Algorithms for Synthetic Aperture Radar Imagery XIX, Baltimore, USA, 2012: 83940F. [29] DEMING R, BEST M, and FARRELL S. Simultaneous SAR and GMTI using ATI/DPCA[C]. SPIE Algorithms for Synthetic Aperture Radar Imagery XXI, Baltimore, USA, 2014: 90930U. [30] PILLAI U, LI Keyong, and SCARBOROUGH S M. Geolocation of moving targets in Gotcha data using multimode processing[C]. SPIE Algorithms for Synthetic Aperture Radar Imagery XXII, Baltimore, USA, 2015: 947508. [31] PILLAI U, LI Keyong, and SCARBOROUGH S M. Target geolocation in Gotcha data using panoramic processing[C]. 2015 IEEE Radar Conference, Arlington, USA, 2015: 21–26. [32] RAYNAL A M, BICKEL D L, and DOERRY A W. Stationary and moving target shadow characteristics in synthetic aperture radar[C]. SPIE Radar Sensor Technology XVIII, Baltimore, USA, 2014: 90771B. [33] WELLS L, SORENSEN K, DOERRY A, et al. Developments in SAR and IFSAR systems and technologies at Sandia National Laboratories[C]. 2003 IEEE Aerospace Conference Proceeding, Big Sky, USA, 2003: 2_1085–2_1095. [34] Sandia national laboratories. Pathfinder radar ISR & SAR systems[EB/OL]. https://www.sandia.gov/radar/pathfinder-radar-isr-and-synthetic-aperture-radar-systems/video/. [35] DAMINI A, BALAJI B, PARRY C, et al. A videoSAR mode for the x-band wideband experimental airborne radar[C]. SPIE Algorithms for Synthetic Aperture Radar Imagery XVII, Orlando, USA, 2010: 76990E. [36] CHEN Jingwei, AN Daoxiang, WANG Wu, et al. A novel generation method of high quality video image for high resolution airborne ViSAR[J]. Remote Sensing, 2021, 13(18): 3706. doi: 10.3390/rs13183706 [37] PALM S, SOMMER R, JANSSEN D, et al. Airborne circular W-band SAR for multiple aspect urban site monitoring[J]. IEEE Transactions on Geoscience and Remote Sensing, 2019, 57(9): 6996–7016. doi: 10.1109/tgrs.2019.2909949 [38] HENKE D and MEIER E. Tracking and refocussing of moving targets in multichannel SAR data[C]. 2015 IEEE International Geoscience and Remote Sensing Symposium (IGARSS), Milan, Italy, 2015: 3735–3738. [39] SHEN Wenjie, LIN Yun, YU Lingjuan, et al. Single channel circular SAR moving target detection based on logarithm background subtraction algorithm[J]. Remote Sensing, 2018, 10(5): 742. doi: 10.3390/rs10050742 [40] 聊蕾, 左潇丽, 云涛, 等. 基于图像序列的VideoSAR动目标检测方法[J]. 雷达科学与技术, 2016, 14(6): 563–567, 573. doi: 10.3969/j.issn.1672-2337.2016.06.001LIAO Lei, ZUO Xiaoli, YUN Tao, et al. An approach to detect moving target in VideoSAR imagery sequence[J]. Radar Science and Technology, 2016, 14(6): 563–567, 573. doi: 10.3969/j.issn.1672-2337.2016.06.001 [41] 张营, 朱岱寅, 俞翔, 等. 一种VideoSAR动目标阴影检测方法[J]. 电子与信息学报, 2017, 39(9): 2197–2202. doi: 10.11999/JEIT161394ZHANG Ying, ZHU Daiyin, YU Xiang, et al. Approach to moving targets shadow detection for VideoSAR[J]. Journal of Electronics & Information Technology, 2017, 39(9): 2197–2202. doi: 10.11999/JEIT161394 [42] ZHONG Chao, DING Jinshan, and ZHANG Yuhong. Joint tracking of moving target in single-channel video SAR[J]. IEEE Transactions on Geoscience and Remote Sensing, 2022, 60: 5212718. doi: 10.1109/TGRS.2021.3115491 [43] QIN Siqi, DING Jinshan, WEN Liwu, et al. Joint track-before-detect algorithm for high-maneuvering target indication in video SAR[J]. IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing, 2021, 14: 8236–8248. doi: 10.1109/jstars.2021.3104603 [44] LIU Zhongkang, AN Daoxiang, and HUANG Xiaotao. Moving target shadow detection and global background reconstruction for VideoSAR based on single-frame imagery[J]. IEEE Access, 2019, 7: 42418–42425. doi: 10.1109/ACCESS.2019.2907146 [45] CHEN Leping, AN Daoxiang, and HUANG Xiaotao. A backprojection-based imaging for circular synthetic aperture radar[J]. IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing, 2017, 10(8): 3547–3555. doi: 10.1109/JSTARS.2017.2683497 [46] CHEN Leping, AN Daoxiang, and HUANG Xiaotao. Resolution analysis of circular synthetic aperture radar noncoherent imaging[J]. IEEE Transactions on Instrumentation and Measurement, 2020, 69(1): 231–240. doi: 10.1109/TIM.2019.2890932 [47] WANG Wu, AN Daoxiang, LUO Yuxiao, et al. The fundamental trajectory reconstruction results of ground moving target from single-channel CSAR geometry[J]. IEEE Transactions on Geoscience and Remote Sensing, 2018, 56(10): 5647–5657. doi: 10.1109/TGRS.2018.2823310 [48] LUO Yuxiao, AN Daoxiang, WANG Wu, et al. Local road area extraction in CSAR imagery exploiting improved curvilinear structure detector[J]. IEEE Transactions on Geoscience and Remote Sensing, 2022, 60: 5227615. doi: 10.1109/TGRS.2022.3172227 [49] LI Jianpeng, AN Daoxiang, WANG Wu, et al. A novel method for single-channel CSAR ground moving target imaging[J]. IEEE Sensors Journal, 2019, 19(19): 8642–8649. doi: 10.1109/JSEN.2019.2912863 [50] AN Daoxiang, WANG Wu, and ZHOU Zhimin. Refocusing of ground moving target in circular synthetic aperture radar[J]. IEEE Sensors Journal, 2019, 19(19): 8668–8674. doi: 10.1109/JSEN.2019.2922649 [51] GE Beibei, FAN Chongyi, AN Daoxiang, et al. A novel phase calibration method for dual-channel CSAR-GMTI processing[J]. IEEE Geoscience and Remote Sensing Letters, 2020, 17(4): 636–640. doi: 10.1109/lgrs.2019.2926436 [52] GE Beibei, AN Daoxiang, CHEN Leping, et al. Ground moving target detection and trajectory reconstruction methods for multichannel airborne circular SAR[J]. IEEE Transactions on Aerospace and Electronic Systems, 2022, 58(4): 2900–2915. doi: 10.1109/TAES.2022.3141332 [53] GE Beibei, AN Daoxiang, LIU Jinyuan, et al. Modified adaptive 2-D calibration algorithm for airborne multichannel SAR-GMTI[J]. IEEE Geoscience and Remote Sensing Letters, 2023, 20: 4004805. doi: 10.1109/LGRS.2023.3267148 [54] 葛蓓蓓. 机载多通道WasSAR运动目标检测与轨迹重构技术研究[D]. [博士论文], 国防科技大学, 2023.GE Beibei. Study on moving target detection and trajectory reconstruction methods for airborne multi-channel WasSAR[D]. [Ph.D. dissertation], National University of Defense Technology, 2023. [55] TIAN Xiaoqing, LIU Jing, MALLICK M, et al. Simultaneous detection and tracking of moving-target shadows in ViSAR imagery[J]. IEEE Transactions on Geoscience and Remote Sensing, 2021, 59(2): 1182–1199. doi: 10.1109/TGRS.2020.2998782 [56] HENKE D, MAGNARD C, FRIOUD M, et al. Moving-target tracking in single-channel wide-beam SAR[J]. IEEE Transactions on Geoscience and Remote Sensing, 2012, 50(11): 4735–4747. doi: 10.1109/tgrs.2012.2191561 [57] HENKE D, DOMINGUEZ E M, SMALL D, et al. Moving target tracking in single- and multichannel SAR[J]. IEEE Transactions on Geoscience and Remote Sensing, 2015, 53(6): 3146–3159. doi: 10.1109/tgrs.2014.2369060 [58] ZHANG Yun, MU Huilin, JIANG Yicheng, et al. Moving target tracking based on improved GMPHD filter in circular SAR system[J]. IEEE Geoscience and Remote Sensing Letters, 2019, 16(4): 559–563. doi: 10.1109/LGRS.2018.2878467 [59] 穆慧琳. 多通道SAR地面运动目标检测与成像研究[D]. [博士论文], 哈尔滨工业大学, 2021.MU Huilin. Research on ground moving target detection and imaging in multichannel SAR system[D]. [Ph.D. dissertation], Harbin Institute of Technology, 2021. [60] DEMING R W. Along-track interferometry for simultaneous SAR and GMTI: Application to Gotcha challenge data[C]. SPIE Algorithms for Synthetic Aperture Radar Imagery XVIII, Orlando, USA, 2011: 80510P. [61] 刘向阳. 机载多通道SAR-GMTI误差分析与补偿方法研究[D]. [博士论文], 西安电子科技大学, 2010.LIU Xiangyang. Study on error analysis and compensation for multi-channel airborne SAR-GMTI systems[D]. [Ph.D. dissertation], Xidian University, 2010. [62] ENDER J H G. The airborne experimental multi-channel SAR-system AER-II[C]. European SAR Conference, Konigswinter, Germany, 1996: 49–52. [63] GIERULL C H. Digital channel balancing of along-track interferometric SAR data[R]. Technical Report DRDC TM 2003-024, 2003. [64] SONG Yongping, LOU Jun, and JIN Tian. A novel II-CFAR detector for ROI extraction in SAR image[C]. 2013 IEEE International Conference on Signal Processing, Communication and Computing (ICSPCC 2013), Kunming, China, 2013: 1–4. [65] PALM S, MARESCH A, and STILLA U. Investigation on circular mapping by FMCW-SAR on small airplanes[J]. The International Archives of the Photogrammetry, Remote Sensing and Spatial Information Sciences, 2013, 40-1(W-1): 281–286. doi: 10.5194/isprsarchives-XL-1-W1-281-2013 期刊类型引用(0)
其他类型引用(2)
-
下载:
计量
- 文章访问数: 736
- HTML全文浏览量: 190
- PDF下载量: 195
- 被引次数: 2
