利用全极化SAR(Polarimetric SAR, POLSAR)进行GMTI的研究分为单通道与多通道2类方法。在单通道POLSAR GMTI领域，加拿大****Liu[11-12]提出利用当前分时极化测量体制POLSAR所获取的交叉极化项进行干涉处理以提高运动目标检测性能的方法，并采用CV-580 SAR实测全极化数据证明了交叉极化进行GMTI的潜力，但本质上这种利用了HV与VH极化的互易性的方法仍可看作单极化ATI方法的拓展。对于多通道POLSAR GMTI，文献[13]指出利用单航过极化干涉SAR进行GMTI是未来极化干涉技术的重要应用，并以可测速区间为优化准则对全极化顺轨干涉SAR(AT-POLINSAR，添加前缀“AT-”用以区分极化交轨干涉SAR(POLINSAR))系统进行了优化设计，但该研究立足于现有的分时极化测量SAR，并不适用于未来同时极化测量体制且没有对新系统进行具体的检测器设计。文献[7-10]指出具有精确极化散射特性测量能力的同时极化测量技术是未来乃至相当长一段时期内极化测量体制发展的主流趋势，因而开展同时测量AT-POLINSAR慢动目标自动检测算法的研究对POLSAR GMTI理论体系的建立以及AT-POLINSAR系统设计都具有重要的意义。当前可直接用于AT-POLINSAR进行GMTI的算法主要有2类：一是文献[1]提出的基于多频段多孔径极化SAR(MFMA POLSAR)的方法，即通过最优极化干涉理论[14-16]对MFMA POLSAR回波降维后采用杂波对消进行运动目标成像；二是文献[17]提出的利用对数似然比(Likelihood-Ratio Test, LRT)[18]对目标极化干涉协方差矩阵进行检测的方法。但这2种方法都难以推导检测量的统计分布，且依赖于对目标先验信息的获取，不宜扩展成恒虚警(Constant False Alarm Rate, CFAR)检测器。综上所述，为设计实用的AT-POLINSAR慢动目标自动检测算法，主要有2方面亟待研究的问题：一是建立POLSAR顺轨干涉的统计建模与参数估计理论，与传统单极化ATI不同的是，全极化ATI除考虑同一极化通道内的干涉相干度外，还要额外考虑不同极化通道之间的相干度以充分挖掘运动目标与杂波背景之间的差异；二是构造性能较优且解析门限易于求解的检测量，单极化条件下利用干涉图幅度和相位联合CFAR检测量的检测性能较优[2-4]，那么如何将这种方法扩展至全极化或者针对极化矢量干涉回波重新构造新的检测量将是实现慢动目标CFAR检测的关键。
针对上述问题，本文首先对单基线AT-POL-INSAR回波进行建模分析；然后针对其6维极化矢量干涉数据，提出基于背景杂波平均相干度最优的降维新方法，从高维干涉矢量中合成出次优极化标量干涉数据，从而可以将单极化顺轨干涉SAR(AT-INSAR)慢动目标检测算法直接应用到全极化系统进行目标检测；最后，分别采用蒙特卡罗仿真实验与场景检测实验对次优极化、单极化的慢动目标CFAR检测进行性能对比。
1 全极化顺轨干涉SAR系统 本节基于MFMA POLSAR[1]与分时极化AT-INSAR[13]的设计思路，提出一种设备复杂度相对较低、可实现度高的单基线AT-POLINSAR的新系统。基于对系统工作方式、天线收发模式与信号形式的设计，建立了极化矢量干涉模型, 将SAR单极化顺轨干涉扩展成全极化矢量干涉。本节的AT-POLINSAR的信号模型可扩展到MFMA POLINSAR，同时，矢量干涉模型也适用于分时极化测量顺轨干涉SAR。
1.1 AT-POLINSAR系统信号与回波建模 图 1为单基线AT-POLINSAR系统的成像几何示意图，设雷达平台在高度为H的直线轨迹以速度vp匀速飞行，沿雷达平台运动方向配置2副相同的全极化天线分别作为前、后向干涉通道m和s，并在不同的方位采样时刻同时发射h极化与v极化状态的正交波形，从而在h、v极化天线之间形成了共同的物理干涉基线2Lb。关于天线工作方式，本文假设2副天线按照“乒乓”模式工作，即2副天线轮流发射全极化信号并同时接收极化回波，相较于全孔径发射子孔径接收的方案，这种方式在基线配置、系统实现等方面具有更高的灵活性[19]。
 |
| 图 1 单基线AT-POLINSAR成像几何示意图 Fig. 1 Geometry schematic diagram of single-baseline AT-POLINSAR acquisition |
| 图选项 |
图 1中，β为合成孔径雷达下视角，Ra为雷达成像斜距，θ为斜视角，Lb为基线长度，vp为平台运动速度，设雷达平台运动方向(方位向)为x轴正方向，设地距方向为y轴正方向。在t=0时刻，前向干涉通道m的天线相位中心的坐标为(0, 0, H)，后向干涉通道s的天线相位中心的坐标为(－Lb, 0, H)，设目标以初始坐标T0(x0, y0, 0)作匀速直线运动，其方位向与距离向的速度分量分别为vx、vy。m通道天线发射全极化线性调频(Linear Frequency Modulation, LFM)信号s=[sh(τ), sv(τ)]T[9-10]，sh(τ)与sv(τ)分别为水平、垂直极化，两者相互正交，且具有如下的共同形式：
 | (1) |
式中：τ为合成孔径处理的快时间；fc为发射信号载频；Kr为线性调频率；Tp为脉冲宽度；|·|表示取绝对值算；Re{·}表示取实部运算。为讨论方便，下文SAR信号处理时不区分极化状态。发射信号s(τ)目标T0散射后分别由m、s通道的全极化天线进行接收。
 | (2) |
其中：sxyi(t, τ)表示i通道x极化天线接收到的y极化状态发射信号的回波，t为慢时间，x, y∈{h, v}，i∈{m, s}，这里h、v表示极化状态，m、s则表示干涉通道类型；Axy为目标散射系数；τdi为i通道接收到回波的时延，若t时刻目标与m通道和i通道天线相位中心的瞬时斜距分别为RTm(t)与RTi(t)，则
 | (3) |
其中：sinc(·)表示冲激响应函数；Tr为脉冲重复时间；λ为载波波长，对瞬时斜距进行抛物线近似[20]，i通道回波的相位即为
 | (4) |
式中：R0为最短距离[20]，mi的下标i=m或s，mm=0，ms=1。对式(4)的信号进行方位向压缩，并忽略幅度系数，可以得到目标在目标存在的某一固定距离单元内的复图像为
 | (5) |
其中：u为方位向像素；ux0为静止目标方位向像素位置；Δu=ΔX/ΔRx为目标偏移像素数，ΔRx为方位向分辨率，ΔX为方位向偏移距离，且ΔX=
 | (6) |
式中：
 | (7) |
式中：vr为径向速度。以上是m通道天线发射、s通道天线接收情形下天线接收回波信号情况，若s通道天线发射、m通道天线接收，则其干涉基线将变为2Lb，φ′i=2φi。相对于天线全孔径发射、子孔径接收的方案而言，采用“乒乓”工作模式可使信号处理变得更加灵活，且根据实际应用选取基线的长短类型。上述工作模式不仅适用于同时极化AT-INSAR，也适用于分时极化AT-INSAR系统；在当前技术条件下，若忽略分时极化体制对散射矩阵测量结果的影响[7-8]，则无需重新设计AT-POLINSAR即可利用两极化SAR沿航向配置并交替发射全极化信号进行ATI，这种方案尤其适用于具有固定轨道的星载SAR进行顺轨干涉。
1.2 AT-POLINSAR系统极化矢量干涉模型 由上述分析，2种极化状态的回波经过前后干涉通道的h极化与v极化天线接收后，可分别在前后干涉通道内测量到2个2×2极化散射矩阵[S]m与[S]s。为从极化散射矩阵获取干涉相位信息[14]并获取更高阶的统计量，通常需要先将散射矩阵向量化，这里采用字典式基ΨL对S进行基变换，得到前后干涉通道所获取的极化矢量分别为
 | (8) |
式中：上标H表示复共轭转置；|·|为秩运算；q=6，为散射矢量维数；
 | (9) |
其中：〈·〉表示统计平均，
 | (10) |
其中：
 | (11) |
采用R进行表示的优势在于能够同时分析极化通道与干涉通道之间的相干度，且二者都能用于区分运动目标与静止目标。因而，对于单极化SAR ATI来说，以hh极化回波(h极化发射、h极化接收)为例，只需要考虑c11、c33即可。在式(9)的复高斯分布假设下，C将服从复Wishart分布[5-6]：
 | (12) |
式中：Γ(·)为不完全伽马函数；n为多视处理时的名义视图数；
2 极化矢量干涉降维 由第1节的信号模型可知，相较于传统的极化SAR或干涉SAR，AT-POLINSAR信号处理维度急剧增加，为从其回波中检测出运动目标并估计运动状态，需要对高维极化干涉矢量数据进行降维处理，本节在回顾经典的最优极化干涉降维方法的基础上，提出了一种新的次优极化干涉降维方法。
2.1 经典最优极化干涉降维方法 对极化矢量进行降维处理并定义一种“新的单极化”状态[14-15]，可得到前后干涉通道的标量图像分别为
 | (13) |
式(13)即为全极化矢量干涉的一般表达式，通过其可以充分利用极化矢量提高干涉数据的质量。从式(13)中还可以提取出极化干涉图为I=
 | (14) |
对于式(14)需要说明的是，相干度的任意损失都将引起干涉相位差估计精度的下降，因而“信号s1之间的与s2相干度ρ越大，干涉相位φi的估计就越准确”[1]。为获得最大相干度ρmax，Cloude和Papathanassiou[14]首次提出通过式(15)求解最优权矢量wm与ws的方法：
 | (15) |
通过式(15)，可以求得km与ks的最优加权矢量分别为wm, opt与ws, opt，利用最优权矢量加权得到的新的合成极化状态μm, opt与μs, opt称为最优极化，同样还可得到最优干涉矩阵Jopt，最优干涉图Iopt以及最优干涉相位φi, opt。
2.2 基于背景杂波平均相干度最优的次优极化干涉降维新方法 传统的干涉图模型之所以不再适用于最优极化状态，是因为极化矢量km和ks最优加权矢量wm, opt与ws, opt是从每一像素单元逐一计算得到的[16]，而这种非线性的加权系数会严重影响最优干涉图的分布类型和分布函数的复杂度。但如果对整副全极化图像采用统一的线性加权矢量进行降维，分布推导的问题将可得到简化，但线性加权矢量的选取并不是任意的，其必须满足一定的准则。本文为对极化数据线性降维处理并使其降维后的杂波背景区域仍然能够像最优干涉那样达到较高的相干度，提出了针对背景杂波的次优干涉相干度准则(Suboptimal Interferometric Coherence Principle, SICP)。
假设待检测场景具有均匀的杂波背景，AT-POLINSAR能先后从前、后向干涉通道中获取该区域的两幅极化SAR图像单视(SLC)数据。令两极化干涉通道的线性加权矢量分别为wm, subopt和ws, subopt，可得加权后前、后向干涉通道的标量图像分别为
 | (16) |
式中：
 | (17) |
同式(9)，
 | (18) |
对其求偏导，得
 | (19) |
式中：λ1与λ2分别为引入的拉格朗日乘数。令式(19)中Lc的偏导为0求解
 | (20) |
根据式(20)的特征值解λ1和λ2，即可定义公共特征值Λ=λ1λ2*，那么
下面讨论次优极化的统计特性。加权后杂波次优极化干涉矢量
 | (21) |
式中：上标q=2，表示干涉矢量的维数；
 | (22) |
 | (23) |
式(22)中：I0(·)与Kn－1(·)分别为第1类和第2类修正的贝塞尔函数，二者的阶数分别为0阶与n－1阶[22]。式(23)中：－π < ψ≤π；2F1为广义超几何函数，分布函数的最大取值是ψ=θ，对于静止杂波背景而言θ=0 rad。对AT-POLINSAR分布式杂波进行仿真，可以得到式(22)和式(23)的理论PDF曲线对实验数据统计直方图，拟合效果如图 2(a)和(b)所示，由图可见单极化干涉图统计分布理论对次优极化的适用性。文献[23]对基于干涉图的自动检测算法进行了研究，提出了检测性能较优的利用幅度相位进行联合检测的M&P检测量并推导了检测量所服从的S分布[见文献[23]式(6)~式(13)]。本文将M&P检测量应用于次优极化干涉图幅度、相位并将其命名为次优极化检测量(SICD)，采用仿真实验可以得到S分布对次优极化SICD检测量的统计分布拟合效果如图 2(c)所示。
 |
| 图 2 次优极化统计分布 Fig. 2 Statistical distribution of suboptimal polarization |
| 图选项 |
图 2(a)~(c)所示的拟合效果充分验证了式(21)结论的正确性，表明次优极化干涉与单极化干涉具有相同的统计分布类型，因而可直接将单极化多通道SAR慢动目标检测方法应用到次优极化进行慢动目标CFAR检测。图 2(d)为仿真得到的各单极化状态、最优极化状态、次优极化状态的相干度统计分布，由图可知，次优极化同样能够达到较高的相干度，从而可以提高分布式杂波背景的相位估计精度。
图 3为hh单极化与次优极化的干涉图在极坐标系下的二维分布散点图(图中径向表示相干度幅度，紫色数字标注)，从图中也看出上述结论。对比图 3(a)、(b)可以发现, 次优干涉图的杂波分布比较集中，与目标之间的区分度明显提高，这也是次有干涉准则能够利用全极化优势以提高检测性能的原因。但需要注意的是，由于目标采用了与杂波相同的加权系数，其相干度往往不能达到最优，干涉相位中心也会发生随机的偏移，目标干涉相干度与相估计中心偏移的方向、大小与杂波特性密切相关，但总的来说，采用这种降维方式可较充分地利用全极化信息以提高运动目标检测时的信杂比。
 |
| 图 3 干涉图二维散布特性 Fig. 3 scatter characteristics of interferograms |
| 图选项 |
3 AT-POLINSAR慢小目标CFAR检测 3.1 基于SICP的慢小目标CFAR检测流程 第2节讨论了次优极化的统计特性及干涉图散布特性，结果表明现有单极化AT-INSAR慢动目标CFAR检测方法亦适用于次优极化，因而采用SICP准则对AT-POLINSAR进行慢动目标检测的流程可归纳为图 4所示。
 |
| 图 4 采用SICP进行运动目标检测流程图 Fig. 4 Flowchart of moving target detection using SICP |
| 图选项 |
如图 4所示，检测时首先针对前后干涉通道的SLC数据划分保护区域和协方差矩阵估计区域，并从中估计出前后干涉通道的杂波的多视协方差矩阵
1) 利用干涉协方差矩阵估计结果
2) 将wm, subopt与ws, subopt分别应用到前后通道的全极化图像，矢像加权结果作为“次优极化”输出，即μm, subopt与μs, subopt。
3) 对前后干涉通道应用慢动目标检测算法进行慢动目标检测。
3.2 检测实验与结果
3.2.1 蒙特卡罗仿真实验 为检验次优极化的检测性能，首先将次优检测与单极化检测、全极化LRT方法[17]进行性能对比。采用蒙特卡罗实验，对文献[23]中单极化干涉相位检测器、天线相位中心偏置(DPCA)检测器、双曲线检测器以及M&P检测器应用到次优极化时的检测性能接收机工作特性(ROC)曲线进行仿真分析。实验前首先需要对AT-POLINSAR进行数据仿真，本文借鉴了文献[24]的数据生成方法对AT-POLINSAR分布式杂波与目标数据进行仿真：即产生6维相互独立的复高斯分布数据后对其进行二维相关滤波得到服从6维联合复高斯分布的干涉矢量回波序列。但与文献[24]不同的是，本文采用特征值分解方法代替Choleskey分解方法进行6维复高斯分布序列的产生，从而降低了对干涉协方差矩阵正定性的要求，只需要满足半正定即可对分布式目标或杂波的回波数据进行仿真。为使比较更有针对性，数据仿真的关键参数参考文献[17]进行设置，信杂比(Target-to-Clutter Ratio, TCR)、目标运动速度等参数如表 1所示，其中雷达平台运动速度为126 m/s，重复观测时间为1.6 ms。
表 1 分布式杂波、目标仿真关键参数 Table 1 Key parameters of distributed clutter and target simulation
| 项目 | σxy/σhh | TCR | SNR | ρxy-hh | ρxy-hv | |
| 杂波 | hh | 1 | 1 | 1 | 0.2 | |
| hv | 0.25 | 0.5 | 0.2 | 1 | ||
| vv | 2 | 2 | 0.7 | 0.2 | ||
| 目标1 (0.8 m/s) | hh | 1 | 0.9 | 4 | 1 | 0.5 |
| hv | 0.25 | 0.9 | 1 | 0.5 | 1 | |
| vv | 4 | 0.9 | 8 | 0.7 | 0.5 | |
| 目标2 (1.6 m/s) | hh | 1 | 1.1 | 4 | 1 | 0.5 |
| hv | 0.25 | 1.1 | 1 | 0.5 | 1 | |
| vv | 4 | 1.1 | 8 | 0.7 | 0.5 |
表选项
表 1中，TCR为hh极化状态下的目标与杂波的RCS之比；信噪比(Signal-to-Noise Ratio，SNR)参数用于计算各单极化状态下顺轨干涉的相干度，即还需手动设置杂波与目标极化相干矩阵RC, RT(但不是任意地选取)，本文假设目标交叉极化之间的相干度比杂波低，将RC与RT设为
 | (24) |
由上述参数即可对AT-POLINSAR在低TCR条件下的慢动目标和杂波图像进行仿真，对仿真数据进行4视处理，分别采用次优极化、单极化对目标1、2进行107次的蒙特卡罗检测试验，可以得到如图 5所示的ROC曲线。文献[20]指出，实际应用中对于ROC曲线只需要关注其低虚警区间，如(10-5, 10-3)。因而由图 5可得，对于当前2类试验场景，次优极化检测性能排序依次为LRT>SICD>DPCA>双曲线检测器测器。另外，还能发现全极化信息对每一类检测器所带来的检测性能增益是不同的，但总的来说，各检测器均能通过次优极化提高检测概率。
 |
| 图 5 ROC曲线仿真结果 Fig. 5 Simulated curves of ROC |
| 图选项 |
3.2.2 场景检测实验 下面进行场景检测实验，实验前首先采用如图 6(a)所示的AT-POLINSAR场景图像仿真流程进行数据仿真，为显示需要，首先设TCR=2，令目标散射矩阵为
 |
| 图 6 场景图像仿真流程和仿真结果 Fig. 6 Generated flow and result of simulated scene image |
| 图选项 |
表 2 AT-POLINSAR场景图像仿真参数
| 参数 | 数值 |
| 雷达工作距离/m | 4 000.0 |
| 载机飞行速度/(m·s-1) | 276.0 |
| 波长/m | 0.1 |
| 脉冲重复频率/Hz | 662.0 |
| 基线长度/m | 3.34 |
| 合成孔径时间/s | 0.869 6 |
| 重复观测时间/s | 0.012 1 |
| 发射脉冲宽度/μs | 1.0 |
| 发射信号带宽/MHz | 216.0 |
| 发射信号中心频率/GHz | 3.0 |
| 过采样率 | 1.2 |
| 斜视角/(°) | 0 |
| 分辨率(方位/m×距离/m) | 1.0×1.0 |
| 分辨率加权因子 | 1.2 |
| 目标1坐标(距离向/m，方位向/m) | (25, 0) |
| 目标2坐标(距离向/m，方位向/m) | (75, 0) |
| 目标3坐标(距离向/m，方位向/m) | (125, 0) |
| 目标4坐标(距离向/m，方位向/m) | (175, 0) |
| 目标1径向速度/(m·s-1) | 0 |
| 目标2径向速度/(m·s-1) | 0.9 |
| 目标3径向速度/(m·s-1) | -1.2 |
| 目标4径向速度/(m·s-1) | 0 |
| ??注：成像算法为距离多普勒；杂波类型为复Wishart。 | |
表选项
下面进行场景图像检测实验。重新对表 2中的4个点目标进行信号级成像仿真，得如图 7(a)所示点目标响应图像，然后设置TCR=1，得到图 7(e)所示的带有杂波背景的AT-POLINSAR场景图像，由图知，此时所有目标均淹没于杂波背景中，故传统单通道幅度检测已无法对目标进行提取。现首先对AT-POLINSAR场景图像的hh极化通道应用M&P检测器，得到图 6(b)所示的检测量输出；然后利用SICP方法获得求解次优极化，并对次优极化应用DPCA、SICD，可以得到如图 7(c)、(d)所示的SICD输出。经计算，此时M&P(hh极化)、DPCA(次优极化)、SICD对目标2的TCR输出分别为4.77、13.11、17.16 dB，对目标3的TCR输出分别为11.25、13.52、17.39 dB。
 |
| 图 7 AT-POLINSAR信号级仿真图像和不同检测算法的检测结果 Fig. 7 Signal level simulation image AT-POLINSAR and its detection results by different algorithms |
| 图选项 |
最后，设置虚警率Pfa=10－4，对图 7(b)~(d)所示检测量输出进行CFAR检测，可以得到如图 7(f)~(h)所示的检测结果，对比可见，2次优极化检测器对目标区域的检测结果较为完整，且无运动目标漏警发生，而单极化检测器中检测性能较优的M&P检测量则产生了一次漏警，未能检测出目标2。通过次优极化，AT-POLINSAR不仅能够提高对低信杂比条件下慢动的目标检测概率，还可利用全极化信息对检测出的运动目标进行最优相位估计，对慢动目标的径向速度进行高精度的估计，从而计算并校正方位向的成像位置偏移，即运动目标成像，如图 7(h)所示。
3.2.3 AT-POLINSAR速度估计精度 采用次优极化方法对目标的径向运动速度进行估计，需要先对图像进行GMTI，然后再对检测结果中心区域重新进行最优干涉处理和径向速估计，这是由于SICP对均匀图像背景应用统一的权矢量，会导致目标像素位置处的干涉相位具有一定的偏差。虽然采用这种二次估计的速度估计流程不如单极化情形速度提取的方便，但却能充分利用全极化信息，图 8为目标径向速度估计相对误差曲线，从图中可见，在各个测试速度点上，利用AT-POLINSAR全极化信息进行相位估计，均能达到最优的径向速度估计精度(亦可验证文献[14]中“干涉相干度越高，相位估计越准确”的结论)。
 |
| 图 8 目标径向速度估计相对误差 Fig. 8 Relative estimation error of target radial velocity |
| 图选项 |
3.2.4 极化干涉协方差矩阵估计点数的选取 背景杂波极化干涉协方差矩阵的估计是SICP处理的重要步骤，协方差矩阵的估计精度会直接影响次优极化干涉相干度及慢动目标CFAR检测器的检测性能，如图 9所示。从图 9的不同估计点数下相干度分布和SICD检测性能ROC曲线结果可以发现，对于常用的4视处理极化顺轨干涉SAR图像，估计点数选择64即可获取较高的杂波相干度和慢动目标检测性能，因而3.2节开展实验时，SICP处理选用8×8窗口对协方差矩阵进行估计。
 |
| 图 9 估计点数对SICP性能的影响 Fig. 9 Performance of SICP vs. estimation point number |
| 图选项 |
4 结论 本文对如何利用AT-POLINSAR进行慢动目标CFAR检测的问题进行了研究。在对单基线AT-POLINSAR进行回波建模的基础上，针对其高维矢量干涉数据提出了基于次优干涉相干度准则。实验结果表明：
1) 次优极化能使杂波背景与目标均获取较高的相干度，但目标的相干度取值与目标特性有关。
2) 次优极化因其具有线性加权矢量, 故统计分布类型与单极化相同，从而单极化AT-INSAR慢动目标检测理论能够通过其扩展至全极化并实现AT-POLINSAR的慢动目标CFAR检测。
3) 将单极化检测器应用到次优极化，检测器的性能排序与单极化情形基本一致，但通过全极化方法能显著提高干涉系统对慢小目标的检测概率。
AT-POLINSAR在慢小目标CFAR检测方面具有绝对的优势，但需要指出的是，目前仍然没有实际的AT-POLINSAR系统，为快速获取AT-POLINSAR半实测等效数据进行检测算法研究，并将极化干涉方法应用于实际POLSAR系统的GMTI，在下一步的研究中可参考文献[25]的单通道ATI方法展开单通道POLSAR ATI的研究。
参考文献
| [1] | WANG G Y, XIA X G, CHEN V C. Radar imaging of moving targets in foliage using multifrequency multiaperture polarimetric SAR[J]. IEEE Transactions on Geoscience and Remote Sensing, 2003, 41(8): 1755-1764. DOI:10.1109/TGRS.2003.813501 |
| [2] | GIERULL C H. Statistical analysis of multilook SAR interferograms for CFAR detection of ground moving targets[J]. IEEE Transactions on Geoscience and Remote Sensing, 2004, 42(4): 691-701. DOI:10.1109/TGRS.2003.821886 |
| [3] | GAO G, WANG X, LAI T. Detection of moving ships based on a combination of magnitude and phase in along-track interferometric SAR-Part Ⅱ:Statistical modelling and CFAR detection[J]. IEEE Transactions on Geoscience and Remote Sensing, 2015, 53(7): 3582-3599. DOI:10.1109/TGRS.2014.2379351 |
| [4] | GIERULL C H, SIKANETA I, CERUTTI-MAORI D.Space-based SAR ground moving target indication[M]//KLEMM R, NICKEL U, GIERULL C, et al.Novel radar techniques and applications volume 1: Real aperture array radar, imaging radar, and passive and multistatic radar.London: SciTech Publishing, 2015: 313-360. |
| [5] | 李仲森, 埃里克·鲍狄埃.极化雷达成像基础与应用[M].洪文, 李洋, 尹婧, 等, 译.北京: 电子工业出版社, 2013: 82-83. LEE J S, POTTIER E.Polarimetric radar imaging: From basics to applications[M].HONG W, LI Y, YIN Q, et al.translated.Beijing: Publishing House of Electronics Industry, 2013: 82-83(in Chinese). |
| [6] | 匡纲要, 陈强, 陆军. 极化合成孔径雷达基础理论及其应用[M]. 长沙: 国防科技大学出版社, 2011: 1-20. KUANG G Y, CHEN Q, LU J. Polarimetric synthetic aperture radar basic principles and its application[M]. Changsha: NUDT Publish House, 2011: 1-20. (in Chinese) |
| [7] | WANG X S, LI Y Z, DAI H Y, et al. Research on instantaneous polarization radar system and external experiment[J]. Chinese Science Bulletin, 2010, 55(15): 1560-1567. DOI:10.1007/s11434-010-3102-y |
| [8] | 李永祯, 李棉全, 程旭, 等. 雷达极化测量体制研究综述[J]. 系统工程与电子技术, 2013, 35(9): 1873-1877. LI Y Z, LI M Q, CHENG X, et al. Summarization of radar polarization measurement modes[J]. System Engineering and Electronic Technology, 2013, 35(9): 1873-1877. DOI:10.3969/j.issn.1001-506X.2013.09.12 (in Chinese) |
| [9] | 王雪宾.同时极化雷达测量波形优化及检测技术研究[D].哈尔滨: 哈尔滨工业大学, 2016. WANG X B.Research on waveform optimization and detection technology of simultaneous polarimetric radar[D].Harbin: Harbin Institute of Technology, 2016(in Chinese). http://www.wanfangdata.com.cn/details/detail.do?_type=degree&id=D01099124 |
| [10] | 王福来, 庞晨, 李永祯, 等. 一种同时全极化雷达正交多相编码波形设计方法[J]. 雷达学报, 2017, 6(4): 340-348. WANG F L, PANG C, LI Y Z, et al. Orthogonal polyphase coded waveform design method for simultaneous fully polarimetric radar'[J]. Journal of Radar, 2017, 6(4): 340-348. (in Chinese) |
| [11] | LIU C. Effects of target motion on polarimetric SAR images[J]. Canadian Journal of Remote Sensing, 2006, 32(2): 51-56. DOI:10.5589/m06-013 |
| [12] | LIU C. Time-frequency analysis of POLSAR moving target data[J]. Canadian Journal of Remote Sensing, 2007, 33(4): 237-249. DOI:10.5589/m07-019 |
| [13] | 李廷伟, 杜湘瑜, 黄海风, 等. 全极化沿航向干涉SAR系统参数优化设计[J]. 系统工程与电子技术, 2011, 33(1): 69-73. LI T W, DU X Y, HUANG H F, et al. Optimization design for the fully polarimetric along-track interferometric SAR system parameters[J]. Systems Engineering and Electronics, 2011, 33(1): 69-73. DOI:10.3969/j.issn.1001-506X.2011.01.14 (in Chinese) |
| [14] | CLOUDE S R, PAPATHANASSIOU K P. Polarimetric SAR interferometry[J]. Remote Sensing Technology and Application, 1999, 36(5): 1551-1565. |
| [15] | PAPATHANASSIOU K P, CLOUDE S R. Single-baseline polarimetric SAR interferometry[J]. IEEE Transactions on Geoscience and Remote Sensing, 2001, 39(11): 2352-2363. DOI:10.1109/36.964971 |
| [16] | 刘忠训, 代大海, 王雪松, 等. 极化SAR动目标定位中的残差图像干涉相干度分析[J]. 电子与信息学报, 2009, 31(3): 606-609. LIU Z X, DAI D H, WANG X S, et al. Analysis of the interferometric coherence between residual images in location of moving target with polarimetric SAR[J]. Journal of Electronics and Information Technology, 2009, 31(3): 606-609. (in Chinese) |
| [17] | STACY N, PREISS M.Polarimetric ATI slow target detection in a log likelihood framework[C]//IEEE Geoscience and Remote Sensing Symposium.Piscataway, NJ: IEEE Press, 2014: 3163-3166. |
| [18] | STACY N, PREISS M.ATI slow target detection in a log likelihood framework[C]//European Conference on Synthetic Aperture Radar(EUSAR), 2008. |
| [19] | 郑春明.合成孔径雷达地面动目标检测技术研究[D].北京: 中国科学院电子学研究所, 2007: 25-65. ZHENG C M.Research on synthetic aperture radar ground moving target detection technology[D].Beijing: Institute of Electronics of Chinese Academy of Sciences, 2007: 25-65(in Chinese). https://www.ixueshu.com/document/f0f0cd7e7ff773b7318947a18e7f9386.html |
| [20] | ZHANG P, ZHANG J F, MA D, et al.Contrastive analysis of the performances of radar multi-look polarimetric detectors[C]//IEEE Radar Conference.Piscataway, NJ: IEEE Press, 2017: 1373-1378. |
| [21] | 罗鹏飞, 张文明. 随机信号分析与处理[M]. 北京: 清华大学出版社, 2012: 20-21. LUO P F, ZHANG W M. Analyses and processes to random signals[M]. Beijing: Tsinghua University Press, 2012: 20-21. (in Chinese) |
| [22] | 《数学手册》编写组. 数学手册[M]. 北京: 高等教育出版社, 2010: 627-628. Writing group of "Mathematical Manual". Mathematical manual[M]. Beijing: Higher Education Press, 2010: 627-628. (in Chinese) |
| [23] | 时公涛.基于干涉图的多通道SAR地面慢动目标自动检测技术研究[D].长沙: 国防科学技术大学, 2009: 67-166. SHI G T.The research on multi-channel SAR automatic detection of slow ground moving targets based on the interferogram[D].Changsha: National University of Defence Technology, 2009: 67-166(in Chinese). http://cdmd.cnki.com.cn/Article/CDMD-90002-2010165389.htm |
| [24] | PIPIA, L, FàBREGAS X.Generation of Pol-SAR and POL-in-SAR data for homogeneous distributed targets simulation[C]//Proceedings of the POLINSAR 2005 Workshop, 2005: 1-4. |
| [25] | 高飞, 孙进平, 袁运能, 等. 单通道SAR实现ATI动目标检测的新方法[J]. 北京航空航天大学学报, 2007, 33(8): 890-893. GAO F, SUN J P, YUAN Y N, et al. New method of moving target detection using single channel SAR ATI[J]. Journal of Beijing University of Aeronautics and Astronautics, 2007, 33(8): 890-893. DOI:10.3969/j.issn.1001-5965.2007.08.006 (in Chinese) |
