极化SAR目标散射旋转域解译理论与应用

陈思伟; 李永祯; 王雪松; 肖顺平

doi:10.12000/JR17033

极化SAR目标散射旋转域解译理论与应用

DOI: 10.12000/JR17033

国防科技大学电子科学与工程学院电子信息系统复杂电磁环境效应国家重点实验室长沙 410073

基金项目: 国家自然科学基金(41301490, 61490690, 61490692)

详细信息

作者简介:
陈思伟(1984–)，男，四川人，博士，国防科技大学电子科学与工程学院讲师，主要研究方向包括雷达极化信息处理、成像雷达信息处理、目标散射建模与解译、微波遥感大数据、环境遥感与灾害遥感等。E-mail: chenswnudt@163.com

李永祯(1977–)，男，内蒙古人，博士，国防科技大学电子科学与工程学院研究员，电子信息系统复杂电磁环境效应国家重点实验室副主任，主要研究方向为新体制雷达与电子对抗。E-mail: e0061@sina.com

王雪松(1972–)，男，内蒙古人，博士，国防科技大学理学院院长，教授，博士生导师，主要研究方向为雷达极化信息处理、雷达目标识别、新体制雷达技术

肖顺平(1964–)，男，江西人，博士，国防科技大学电子科学与工程学院教授，博士生导师，电子信息系统复杂电磁环境效应国家重点实验室主任，主要研究方向包括雷达极化信息处理、电子信息系统仿真评估技术、雷达目标识别等

通讯作者:
陈思伟 chenswnudt@163.com

中图分类号: TN957.52
计量
- 文章访问数: 4975
- HTML全文浏览量: 1014
- PDF下载量: 733
- 被引次数: 7
出版历程
- 收稿日期: 2017-03-28
- 修回日期: 2017-06-28
- 网络出版日期: 2017-10-28

Polarimetric SAR Target Scattering Interpretation in Rotation Domain: Theory and Application

The State Key Laboratory of Complex Electromagnetic Environment Effects on Electronics and Information System, School of Electronics Science and Engineering, National University of Defense Technology, Changsha 410073, China

Funds: The National Natural Science Foundation of China (41301490, 61490690, 61490692)

摘要

摘要: 雷达目标的后向散射敏感于目标姿态与雷达视线的相对几何关系。雷达目标的这种散射多样性给以极化合成孔径雷达(SAR)为代表的成像雷达信息处理与应用造成诸多不便，是当前目标散射机理精细解译和定量应用面临的主要技术瓶颈之一。该文回顾并介绍一种在绕雷达视线旋转域解译目标散射机理的新思路，主要包括新近提出的统一的极化矩阵旋转理论和极化相干特征旋转域可视化解译理论。目标散射旋转域解译方法的核心思想是将特定几何关系下获得的目标信息拓展到绕雷达视线的旋转域，为目标散射信息深度挖掘和利用奠定基础。该文详细分析上述方法导出的一系列新的极化特征参数集，并开展应用验证。对比研究证实了旋转域解译方法在地物辨识与分类等领域的应用潜力。
- 极化合成孔径雷达 /
- 散射机理 /
- 旋转域 /
- 极化目标分解 /
- 图像解译
Abstract: Backscattering of radar targets is sensitive to the relative geometry between target orientations and the radar line of sight. This scattering diversity makes imaging radar represented by polarimetric Synthetic Aperture Radar (SAR) information processing and applications very difficult. This situation has become one of the main bottlenecks in the interpretation of the target scattering mechanism and quantitative applications. In this work, we review and introduce a new interpretation of the target scattering mechanism in the rotation domain along the radar line of sight. This concept includes the recently established uniform polarimetric matrix rotation theory and polarimetric coherence pattern visualization and interpretation in the rotation domain. The core idea of target scattering interpretation in the rotation domain is to extend the amount of target information acquired at a given geometry to the rotation domain, which then provides fundamentals for the deep mining and utilization of target scattering information. This work mainly focuses on the investigation of derived new polarimetric feature sets and application demonstrations. Comparison study results validate the promising potential for the application of the established interpretation framework in the rotation domain with respect to target discrimination and classification.
- Polarimetric Synthetic Aperture Radar (PolSAR) /
- Scattering mechanism /
- Rotation domain /
- Polarimetric target decomposition /
- Image interpretation

HTML全文

1. 引言

超宽带雷达具有穿透能力强和距离向分辨率高等优点^[1]，它广泛应用于地质勘测、生命救援、穿墙目标跟踪等多种军用和民用领域^[2–4]。传统的超宽带雷达是窄脉冲体制，它实现简单，但抗干扰能力差，在跟踪距离较远、干扰强等复杂的环境中，这种体制的雷达难以满足实际的要求^[5,6]。新体制的超宽带雷达有线性调频体制、频率步进体制和伪随机编码体制，它们的工作方式类似，发射端发射持续时间较长，带宽足够，幅度相对较小的大时间带宽积信号；接收端接收到回波后先进行脉冲压缩或者是傅里叶变换处理，然后才得到雷达的脉冲响应函数。线性调频体制和频率步进体制雷达的信源设计复杂，系统线性度要求高，目前难以在地质勘测、生命救援、穿墙目标跟踪等应用中普及。伪随机编码信号可以由数字芯片直接合成，回波信号经过脉冲压缩而实现高信噪比和强抗干扰特性^[7,8]。上世纪90年代中期，欧洲的火星探测计划促进了伪随机编码体制探地雷达的研制^[9]；德国伊尔梅瑙科技大学的学者研制了发射无载频基带信号的伪随机编码雷达；英国Utsi Electronics公司在2007年推出了中心频率为15 MHz的伪随机编码雷达Groundvue-6^[10]。2012年开始，中科院电子所研制了超深编码探地雷达系统^[11]和伪随机编码体制的穿墙雷达、生命探测雷达^[12]。

近年来，超宽带雷达日益发展成熟，在工程应用中得到普及，但是单发单收雷达日渐显示出自身的不足之处。比如，在地质勘测方面，单发单收雷达观测到的目标信息通常认为是位于测线的正下方；同样地，穿墙目标跟踪和生命探测方面，单发单收雷达探测到的信息通常认为是位于雷达的正前方。虽然合成孔径雷达Synthetic Aperture Radar (SAR)方法能获得目标的距离向和方位向等多维信息^[13]，但是，平台运动的不确定性导致数据处理复杂，另一方面，平台的机械运动使观测信息难以实时融合处理，这使SAR方法难以推广到生命救援雷达和穿墙目标跟踪雷达等超宽带雷达的应用中去。2004年，美国学者Eran Fishler等人讨论总结了多发多收雷达的发展前景，在论文中指出，多发多收系统可以很大程度地提高雷达的探测性能^[14]。多发多收超宽带雷达系统的多个收发通道协同工作，它能实时探测到目标的距离向和方位向信息^[15]。美国的田纳西州大学研制了一款1发8收的多发多收雷达系统，该系统接收通道通过射频开关切换而实现8个通道分时输入^[16,17]。但由于系统基于脉冲体制，脉冲峰值功率和射频开关耐压的矛盾而制约着雷达系统的性能。GSSI公司的多通道雷达系统SIR30可以实现8个通道同时工作，但是该系统结构复杂、成本高、功耗大，在8个通道同时工作时会造成数据量过大而难以实时处理。早在2002年，美国Time Domain公司就成功研制了一款多发多收穿墙雷达样机^[15]，这是第1款伪随机编码体制的多发多收超宽带雷达系统。对于伪随机编码体制的多发多收雷达系统，目前这方面的文献相对较少。

本文提出一种基于射频开关切换的伪随机编码超宽带多发多收雷达设计方法。为了降低雷达系统的成本和复杂性，同时为了尽可能地提高扫描率，采用射频开关切换的方法来实现多发多收雷达系统，该系统的最大通道能实现4发16收。选择码元长度为1023的二相调制m序列作为雷达的发射信号，信号中心频率为1.06 GHz，由于射频开关的耐压值有限，发射峰值功率为18 dBm。通过混合采样方法实现超宽带信号的采样接收，由200 MSPS的实时采样率通过64轮混合采样实现12.8 GSPS的混合采样率，并且调用FPGA内部的DSP硬核来进行实时脉冲压缩，以缓解后续数据处理运算量大的问题。实验结果表明，该雷达系统在距离向上有较远的探测距离，并且能实时探测到目标的方位向和距离向信息。

2. 雷达系统关键技术分析

2.1 雷达工作过程

伪随机编码多发多收雷达系统的工作过程如图1所示，图中为m发n收雷达工作过程示意图。图1中T_m 表示第m号发射天线，R_n 表示第n号接收天线；S_nm(t)表示第n号发射天线发射信号、第m号接收天线接收到的回波信号，回波信号已经过脉冲压缩处理。

图 1 雷达工作过程示意图

Figure 1. Schematic diagram of radar operating

下载: 全尺寸图片幻灯片

伪随机编码雷达回波经过脉冲压缩处理后得到窄脉冲信号，如图1所示，因此伪随机编码雷达的探测原理与脉冲体制雷达相同。对于m发n收的超宽带雷达系统，每个发射天线都有n组接收天线，所以m发n收雷达系统总共有m×n组不同的接收数据。由于多个收发天线的空间相对位置不同，对这m×n组数据进行聚集处理便可以得到观测目标的距离向和方位向等信息。

2.2 基于射频开关切换的多发多收雷达技术

多发多收雷达可以通过单发单收雷达扩展得到，事实上，把多个独立的单发单收雷达组合起来工作就可以形成多发多收雷达系统，但这时会涉及到雷达之间的时间同步和数据融合问题；另一方面，多个独立的雷达会受到环境因素等的影响而难以保证各个雷达保持一致性。所以，通常会把多发多收雷达做成一个独立的雷达系统，这样不但解决了同步问题，也使数据融合处理显得简单可行。

在超宽带雷达的应用中，观测目标通常是缓慢变化的，因此可以考虑以开关切换的方法实现多发多收超宽带雷达。对于m发n收雷达，为了得到一帧完整的数据，发射开关要依次切换m轮，同时，发射开关每切换一轮，接收开关要依次切换n轮。即经过m×n轮切换才能得到一帧完整的m×n数据，雷达的成像速率比扫率慢m×n倍，因此要尽可能提高雷达系统的扫描率。在开关切换过程中，只要通道的切换速度足够快，就可以近似认为各个通道同时独立工作，并不丢失目标的实时信息。

目前常用的信号切换开关有机械开关和射频开关两种，机械开关通过机械式接触实现不同通道的切换，它的耐压值高、隔离度好，但是切换速度慢，限制了多发多收雷达系统的切换速度。射频开关由半导体器件和微波电路组成，由电压控制来实现不同通道的切换，它的耐压值有限、隔离度适中，但是切换速度快。对于脉冲超宽带雷达，由于脉冲峰值过高而难以由射频开关切换实现多发多收雷达。伪随机编码体制雷达由于峰值功率小而可以很好地应用于基于射频开关切换的多发多收超宽带雷达。

2.3 伪随机编码信号产生技术

为了便于伪随机编码信号的发射，通常以二相调制的伪随机编码信号作为雷达的发射信号。为了得到二相调制信号，一种方法可以借鉴通信中的调制解调思想，采用混频方法把基本信号调制频带上，但是，由于超宽带信号带宽与中心频率相当，器件的宽带调制实现困难。另一种方法是基于FPGA和高速Digital to Analog Convertor (DAC)器件，通过数字方法直接合成伪随机编码雷达发射信号，研究和实验表明，这种方法能方便、可控地生成所需要的信号。

由于FPGA芯片内部主频一般只能工作在500 MHz以下，为了产生高速的伪随机编码信号，结合FPGA内部的并串转换模块、双倍数据速率Double Data Rate (DDR)接口和多路数据接口的DAC芯片，通过FPGA内部低速时钟产生外部高速伪随机编码信号。图2是高速伪随机编码信号产生原理图，选取四路数据接口的高速DAC，其中f_S是FPGA内部DAC采样数据存储器ROM1～ROM16的工作时钟，f_D是DDR接口的工作时钟，f_DAC是外部高速DAC的工作时钟。根据并串转换、DDR接口和四端口DAC芯片的时钟关系，有f_DAC=8f_D和f_D=2f_S，即FPGA内部只需要低速时钟便能产生高速的伪随机编码信号。DAC产生的信号通过信号频带内的低通或带通滤波器后便能产生平滑的编码信号，编码信号经放大器放大后便是雷达天线的发射信号。

图 2 高速伪随机编码信号产生原理图

Figure 2. Schematic diagram of high-speed pseudo random coded signal

下载: 全尺寸图片幻灯片

2.4 超宽带信号的混合采样技术

对于超宽带信号的采集，目前市场上能满足超宽带信号采样率要求的器件很少见，所以通常会采用混合采样对超宽带信号进行采集^[18]。混合采样的思想是牺牲时间来换取更高的采样率，在多个雷达脉冲重复周期Pulse Repetition Period (PRP)内通过多次发送相同的信号来实现一道完整雷达回波的采集。是混合采样过程的时序图，采样过程中Analog to Digital Converter (ADC)工作时钟f_ADC依次比上一轮的混合采样时钟延迟 $\Delta$ t。则混合采样的采样率为：

${f_{\rm H}} = \frac{1}{{\Delta t}}$

(1)

采集一道完整回波的混合采样次数为：

$n = \frac{1}{{{f_{\rm ADC}} \cdot \Delta t}}$

(2)

图 3 混合采样时序图

Figure 3. Timing diagram of hybrid sampling

下载: 全尺寸图片幻灯片

假设每一轮混合采样的采样点数为N_P，则n轮混合采样得到的数据是一个n×N_P矩阵，如式(3)所示。矩阵中每一轮的采样数据按行存储，而按列读出即是相邻时间间隔为 $\Delta$ t的采样数据。因此，矩阵中的元素按依次从上到下、从左到右读取便可以得到一道顺序排列好的完整回波。

${{S}} = \left[ {\begin{array}{*{20}{c}}{{S_1}(0)}&{{S_1}(1)}& \cdots &{{S_1}({N_{\rm P}} - 1)}\\{{S_2}(0)}&{{S_2}(1)}& \cdots &{{S_2}({N_{\rm P}} - 1)}\\ \vdots & \vdots & \ddots & \vdots \\{{S_n}(0)}&{{S_n}(1)}& \cdots &{{S_n}({N_{\rm P}} - 1)}\end{array}} \right]$

(3)

2.5 回波信号的实时脉冲压缩

在伪随机编码体制的超宽带雷达中，原始回波信号的实时脉冲压缩是信号处理的首要和关键步骤。由于超宽带雷达的采样率高、数据量大，而微处理器处理速度和DSP芯片运算能力有限。本文基于时间域的互相关算法，通过调用FPGA中的DSP硬核进行并行计算，并结合流水线模式，实现快速实时脉冲压缩。伪随机编码信号的脉冲压缩公式为：

$R[l] = \sum\limits_{i = 0}^{{N_{\rm r}} - 1} {s[i + l]r[i]} ,\ l = 0, 1, 2, ·\!·\!·, {N_{\rm R}} - 1$

(4)

式中，r是参考码信号，它通常为雷达的发射信号，长度为N_r个采样点；s是雷达回波信号，长度为N_s个采样点。根据式(4)，调用FPGA中的N_r个DSP硬核进行乘法运算，并通过流水线方式进行累加操作，构造图4所示的脉冲压缩并行运算框架。图中D表示1个时钟周期的延时，DSP硬核乘法器的一个输入端分别为各自固定的参考码，另一个输入端同时输入回波信号，压缩结果由最后一级输出端R_sr[i]输出。

图 4 脉冲压缩并行运算框架

Figure 4. Structure of pulse compression computed parallelly

下载: 全尺寸图片幻灯片

由于FPGA中的DSP硬核数量有限，如果参考码长度N_r大于乘加器数量，脉冲压缩并行运算将不能实现。为了解决上面的问题，把长度为N_r的参考码分解成M段，每段长度为N_r0=N_r/M；分段条件满足N_r0小于FPGA中DSP硬核数量，同时通过对参考码补0或者适当选择M值使N_r0为整数。由分段方法，可以得到分段后的第j段参考码与原参考码对应关系为：

${r_j}[i] = r[j \cdot {N_{r0}} + i]$

(5)

式中，j=0, 1, ···,M–1,i=0, 1, ···,N_r0–1。分段后，式(4)可以转化成参考码分段并行运算脉冲压缩，如式(6)：

$R[l] = \sum\limits_{j = 0}^{M - 1} {{R_j}[l]}$

(6)

式中，l=0, 1, ···,N_R–1,R_j 表示第j段参考码的压缩结果，为：

${R_j}[l] = \sum\limits_{i = 0}^{{N_{{\rm r}0}} - 1} {s[j \cdot {N_{{\rm r}0}} + l + i] \cdot {r_j}[i]}$

(7)

即先把各段参考码采用图4的并行运算框架进行脉冲压缩，再把压缩的结果累加，容易证明式(4)和式(6)是等价的。

3. 多发多收雷达系统设计

多发多收雷达系统的原理图如图5所示，主要由4部分组成，从上到下依次为：时钟同步单元、发射机单元、接收机单元和开关切换单元，图中分别用不同的填充色标注。图6是多发多收雷达系统样机，该系统的主要参数如表1所示。

图 5 多发多收雷达系统原理图

Figure 5. Schematic diagram of MIMO radar system

下载: 全尺寸图片幻灯片

图 6 多发多收雷达系统样机

Figure 6. MIMO radar system

下载: 全尺寸图片幻灯片

表 1 多发多收雷达系统主要参数

Table 1. Main parameter of MIMO radar

参数	数值
通道数量	m发n收(最大4发16收)
码元形式	m序列
码元长度	1023
单个码元宽度(ns)	0.9375
发射信号幅值(Vpp)	5
信号中心频率(GHz)	1.06
–10 dB信号带宽(GHz)	0.2～2.1
实时采样率(MSPS)	200
混合采样率(GSPS)	12.8
平均次数(times)	32
接收机灵敏度(dBm)	–97
脉冲重复频率(KHz)	333.3
采样时窗(ns)	300
扫描率(track/s)	164
成像速率(frame/s)	$\frac{{164}}{{m \cdot n}}$
总功耗(W)	20

下载: 导出CSV

| 显示表格

3.1 时钟同步单元的设计

雷达系统需要考虑收发同步问题，特别在混合采样时，对同步精度要求是ns甚至ps级别的，所以雷达系统通常工作在同一时钟基准上。如图5第1部分所示，参考时钟是一个高稳的50 MHz晶振，50 MHz参考时钟经过频率综合器倍频后产生3.2 GHz时钟。3.2 GHz时钟扇出两路，一路作为DAC的数据采样时钟，另一路16分频后生成200 MHz的时钟。200 MHz时钟扇出两路，一路经过可编程延时线，作为ADC的工作时钟，另一路作为FPGA的输入时钟。FPGA的输入时钟由内部锁相环Phase Lock Loop (PLL)生成100 MHz的时钟，100 MHz作为FPGA系统的工作时钟。从上面分析得到，雷达系统的所有时钟都是以50 MHz高稳晶振作为基准，即雷达系统各部分始终是同步的。

3.2 发射机单元的设计

为了提高扫描率同时考虑信号的主旁瓣比，选择码元长度为1023的m序列作为发射信号，单个码元的宽度为0.9375 ns。结合DAC器件的特性同时考虑天线尺寸，选择二相调制信号作为雷达发射信号，信号的载波频率为1.06 GHz。DAC的工作频率为3.2 GHz，对每个码元采样3个有用点，理论上，此时DAC将会产生三角波，但受器件带宽和系统带宽的限制，实际波形会是平滑的双极性脉冲波形；经过低通滤波器后，波形将会变成类似于正弦波的二相调制信号。伪随机编码体制雷达通常发射小功率的信号，因此只需要用低噪声放大器将小信号放大，放大后的信号如图7所示。

图 7 雷达发射信号

Figure 7. Radar transmitting signal

下载: 全尺寸图片幻灯片

3.3 接收机单元的设计

接收天线输出的回波信号先通过低通滤波器，把不需要的高频成份滤去；由于回波信号通常很微弱，为了达到ADC芯片最低检测电平，低通滤波器后级联低噪声放大器把回波信号放大18 dB。ADC的工作时钟为200 MHz，由64轮混合采样得到采样率为12.8 GSPS，此时通过FPGA控制可编程延时线使ADC工作时钟每一轮依次比上一轮延迟78 ps。为了进一步提高回波信号质量，对每个通道的回波信号各自进行32次累加平均，即平均次数32次。

发射信号的时间为960 ns，采样时窗为300 ns，实际采样时间为1260 ns。对于12.8 GSPS的采样率，参考码采样点数为12288点，回波采样点数为16128点，把参考码分成154段调用80个DSP硬核来完成实时脉冲压缩。图8是雷达系统闭环测试得到的脉冲压缩结果图，由于发射单元和接收单元都级联低噪放，为了防止在测试时接收机饱和，测试时在发射和接收机之间级联了40 dB的衰减器。脉冲压缩的旁瓣始终在–30 dB附近，它与m序列的特性一致，所以接收机单元满足雷达系统设计的需求。

图 8 脉冲压缩结果

Figure 8. Results of pulse compression

下载: 全尺寸图片幻灯片

3.4 开关切换单元的设计

选择Mini-Circuits公司的一切四射频开关模块，射频开关的详细参数如表2所示。为了实现最大阵列为4发16收的的多发多收雷达系统，发射通道采用1个上述的一切四射频开关，接收通道采用5个一切四的开关由树型结构级联成一切十六射频开关，图9是收发开关的级联结构图。对于以上6个开关的控制，采用独立的开关控制板，开关控制板与雷达系统通信，由开关控制板实现任意收发组合的阵列形式。

表 2 射频开关详细参数

Table 2. Main parameter of MIMO radar

参数	数值
带宽	0～3 GHz
插损	1.5 dB@0.5～2 GHz
1 dB压缩点	25 dBm@0.5～2 GHz
隔离度	37 dB@0.5～2 GHz
切换速度	25 ns/Turn on; 45 ns/Turn off

下载: 导出CSV

| 显示表格

图 9 收发开关级联结构图

Figure 9. Structure of radio frequency switches

下载: 全尺寸图片幻灯片

为了测试4发16收各个通道的一致性，对雷达系统4发16收开关的任意收发组合进行闭环测试，最终闭环测试得到的64组数据结果如图10所示。从图10可以知道，64组数据的脉冲压缩波形近似相同、幅度近似相等，除了在回波延时存在一定的差别外，其他特性基本相一致。而对于回波的延时差异，通过雷达系统校准可以在数据后处理时补偿回来。

图 10 各个通道一致性测试结果

Figure 10. Results of consistency of all channels

下载: 全尺寸图片幻灯片

对系统的接收灵敏度进行测试，任意选择收发组合通道，发射信号和接收信号都经过低噪声放大器，中间级联可调衰减器后闭环。通过调节衰减器的大小，测试最小能检测到的信号电平。不同衰减值时的测试结果如图11所示，当衰减120 dB时虽然看到峰值点，事实上难以分辨。所以可以认为衰减115 dB时得到系统接收机的最小检测电平，由发射处输出电平为5 Vpp，通过计算得到接收机的灵敏度为–97 dBm。

图 11 不同衰减值的闭环测试结果

Figure 11. Results of different attenuation

下载: 全尺寸图片幻灯片

4. 实验结果

4.1 单发单收测试对比结果

为了验证雷达系统性能，选择单发单收模式下进行穿墙动目标测试，测试墙体为一般楼房的砖墙，厚度为25 cm，并同实验室研制的同频段脉冲雷达进行对比实验。同频段脉冲雷达的中心频率是1.2 GHz, –10 dB带宽为0.4～2.0 GHz，脉冲幅值为36 Vpp，采样率为13.6 GHz，线性平均次数为64次。对比实验测试结果如图12所示，从图中可以知道，伪随机编码雷达在15 m内跟踪结果非常清晰，而同频段的脉冲雷达的跟踪范围在10 m以内。脉冲雷达随着探测目标距离增大，信噪比迅速下降；而伪随机编码雷达可以通过脉冲压缩而进一步提高系统的信噪比。

图 12 不同体制雷达的穿墙测试结果

Figure 12. Results of different radar systems though the wall

下载: 全尺寸图片幻灯片

4.2 2发4收空气中测试结果

为了方便验证多发多收雷达的性能，选用2发4收系统作测试，此时天线口径为1.8 m的线阵，其中2个发射天线在两侧，4个接收天线在中间；两侧发射天线和接收天线的间距是22.5 cm，中间相邻接收天线的间距为45 cm。首先在空气中对2发4收系统进行测试，主要是对空气中目标的距离向和方位向信息进行识别。图13是同时对空气中两个直径为7 cm的可比拟金属杯的测试场景，图14是两个目标中心点在不同间距时的测试结果图。

图 13 2发4收对空气中两个目标测试场景

Figure 13. Experimental scene using 4-input 2-output radar to detect two targets in the air

下载: 全尺寸图片幻灯片

图 14 两个目标在不同间距时的测试结果图

Figure 14. Results of different distance of two targets

下载: 全尺寸图片幻灯片

图14是由BP算法处理得到的结果，从图中可以得到，两个目标的距离向固定在2 m处。当两个目标相距35 cm以上时，可以很容易地把它们区分开来，如图14(a)和图14(b)；当目标相距25 cm时，测试结果中两个目标开始重合，但仍然能区分开来，如图14(c)；当目标相距20 cm时，两个目标的测试结果已经重合在一起，这时并不能把它们区分开来。

4.3 2发4收穿墙测试结果

为了进一步验证多发多收雷达的性能，仍然用以上空气中的2发4收阵列做穿墙实验，探测墙体后目标的距离向和方位向信息。图15是穿墙测试场景，此时直接对墙体后的人进行识别，图16是对墙体后目标的测试结果。测试实验增加了对背景数据(即无目标时)的记录，并将其与有目标时数据相减得到相关结果。

图 15 2发4收穿墙测试场景

Figure 15. Experimental scene using 4-input 2-output radar to detect human through the wall

下载: 全尺寸图片幻灯片

图 16 墙体后目标的测试结果图

Figure 16. Results of human behind the wall

下载: 全尺寸图片幻灯片

测试墙体为一般楼房的砖墙，厚度为25 cm，图16是同样由BP算法处理得到的结果。图16(a)是20 cm×20 cm的金属板在墙后1 m处的测试结果，主要用来对系统进行标定；图16(b)是对墙后正前方人体测试的结果，目标在墙后4.2 m处；图16(c)是对墙后右前方的人体测试的结果，目标在墙后距离向4.2 m，方位向–1.7 m处；图16(d)是对墙后右前方的人体测试的结果，目标在墙后距离向3.8 m，方位向1.3 m处。从图中可以看出，2发4收雷达能对墙体后的目标进行方位向和距离向的区分。

5. 结束语

本文提出一种基于射频开关切换的伪随机编码超宽带多发多收雷达设计方法。采用射频开关切换的方法来实现多发多收雷达系统，基于FPGA和DAC产生伪随机编码信号，由混合采样方法实现高速超宽带信号的接收采样，并直接在FPGA中进行实时脉冲压缩。实验结果表明，该雷达系统能实时探测到目标的方位向和距离向信息，在地质勘测、生命救援、穿墙目标跟踪等领域具有重要的应用价值。

图 1 研究区域荷兰Flevoland

Figure 1. Study area of Flevoland, the Netherlands

下载: 全尺寸图片幻灯片

图 2 导出的极化角参数

Figure 2. Derived polarimetric angle parameters

下载: 全尺寸图片幻灯片

图 3 旋转域极化相干特征的可视化图示例

Figure 3. Illustration of polarimetric coherence pattern example

下载: 全尺寸图片幻灯片

图 4 极化相干特征旋转域优化前后对比图

Figure 4. Comparison of polarimetric coherence without and with optimization in rotation domain

下载: 全尺寸图片幻灯片

图 5 极化相干特征旋转域优化前后的统计分布图

Figure 5. Histograms of polarimetric coherence for full scene AIRSAR data

下载: 全尺寸图片幻灯片

AIRSAR数据中7种已知地物的旋转域极化相干特征的可视化图。(a1)–(a4) 茎豆，(b1)–(b4) 油菜，(c1)–(c4)豌豆，(d1)–(d4)土豆，(e1)–(e4)紫苜蓿，(f1)–(f4)小麦，(g1)–(g4)甜菜。其中，1–4分别代表 $\left| {{\gamma _{\left( {{\rm{HH}} + {\rm{VV}}} \right) - \left( {{\rm{HH}} - {\rm{VV}}} \right)}}\left( \theta \right)} \right|$ , $\left| {{\gamma _{\left( {{\rm{HH}} - {\rm{VV}}} \right) - \left( {{\rm{HV}}} \right)}}\left( \theta \right)} \right|$ , $\left| {{\gamma _{{\rm{HH}} - {\rm{VV}}}}\left( \theta \right)} \right|$ 和 $\left| {{\gamma _{{\rm{HH}} - {\rm{HV}}}}\left( \theta \right)} \right|$

Polarimetric coherence patterns of the seven crop types from AIRSAR data. (a1)–(a4)stembeans, (b1)–(b4)rapeseed, (c1)–(c4)peas, (d1)–(d4)potatoes, (e1)–(e4)lucerne, (f1)–(f4)wheat and (g1)–(g4)beet. The numbers 1–4 indicate $\left| {{\gamma _{\left( {{\rm{HH}} + {\rm{VV}}} \right) - \left( {{\rm{HH}} - {\rm{VV}}} \right)}}\left( \theta \right)} \right|$ , $\left| {{\gamma _{\left( {{\rm{HH}} - {\rm{VV}}} \right) - \left( {{\rm{HV}}} \right)}}\left( \theta \right)} \right|$ , $\left| {{\gamma _{{\rm{HH}} - {\rm{VV}}}}\left( \theta \right)} \right|$ and $\left| {{\gamma _{{\rm{HH}} - {\rm{HV}}}}\left( \theta \right)} \right|$ , respectively

下载: 全尺寸图片幻灯片

图 7 旋转域极化相干特征典型值对比图

Figure 7. Errorbar plots of typical polarimetric coherence parameters

下载: 全尺寸图片幻灯片

图 8 旋转域极化相干度对比图

Figure 8. Errorbar plots for polarimetric coherence degree in rotation domain

下载: 全尺寸图片幻灯片

图 9 旋转域极化相干起伏度对比图

Figure 9. Errorbar plots of polarimetric coherence fluctuation in rotation domain

下载: 全尺寸图片幻灯片

图 10 旋转域极化相干对比度对比图

Figure 10. Errorbar plots of polarimetric coherence contrast in rotation domain

下载: 全尺寸图片幻灯片

图 11 旋转域极化相干宽度对比图

Errorbar plots of polarimetric coherence beamwidth ${\rm{B}}{{\rm{W}}_{0.95}}$ in rotation domain

下载: 全尺寸图片幻灯片

图 12 旋转域极化相干特征最大化角对比图

Errorbar plots of maximum rotation angles ${\theta _{\gamma {\rm{ - }}\max }}$ (in deg) which produce maximized coherence in rotation domain

下载: 全尺寸图片幻灯片

图 13 旋转域极化相干特征最小化角对比图

Errorbar plots of minimum rotation angles ${\theta _{\gamma {\rm{ - }}\min }}$ (in deg) which produce minimized coherence in rotation domain

下载: 全尺寸图片幻灯片

图 14 基于极化特征组合的农作物辨识结果

Figure 14. Crops discrimination results based on combinations of derived polarimetric parameters

下载: 全尺寸图片幻灯片

表 1 旋转域极化相干矩阵的振荡参数集^[40]

Table 1. Oscillation parameter set of polarimetric coherence matrix in rotation domain

元素项	$A = \sqrt \cdot$	B	$\omega$	${\theta _0} = \displaystyle \frac{1}{\omega }{\rm{Angle}}\left\{ \cdot \right\}$
${\mathop{\rm Re}\nolimits} \left[ {{T_{12}}\left( \theta \right)} \right]$	${{\mathop{\rm Re}\nolimits} ^2}\left[ {{T_{12}}} \right] + {{\mathop{\rm Re}\nolimits} ^2}\left[ {{T_{13}}} \right]$	0	2	${\mathop{\rm Re}\nolimits} \left[ {{T_{13}}} \right] + j{\mathop{\rm Re}\nolimits} \left[ {{T_{12}}} \right]$
${\mathop{\rm Re}\nolimits} \left[ {{T_{13}}\left( \theta \right)} \right]$	${{\mathop{\rm Re}\nolimits} ^2}\left[ {{T_{12}}} \right] + {{\mathop{\rm Re}\nolimits} ^2}\left[ {{T_{13}}} \right]$	0	2	$- {\mathop{\rm Re}\nolimits} \left[ {{T_{12}}} \right] + j{\mathop{\rm Re}\nolimits} \left[ {{T_{13}}} \right]$
${\mathop{\rm Im}\nolimits} \left[ {{T_{12}}\left( \theta \right)} \right]$	${{\mathop{\rm Im}\nolimits} ^2}\left[ {{T_{12}}} \right] + {{\mathop{\rm Im}\nolimits} ^2}\left[ {{T_{13}}} \right]$	0	2	${\mathop{\rm Im}\nolimits} \left[ {{T_{13}}} \right] + j{\mathop{\rm Im}\nolimits} \left[ {{T_{12}}} \right]$
${\mathop{\rm Im}\nolimits} \left[ {{T_{13}}\left( \theta \right)} \right]$	${{\mathop{\rm Im}\nolimits} ^2}\left[ {{T_{12}}} \right] + {{\mathop{\rm Im}\nolimits} ^2}\left[ {{T_{13}}} \right]$	0	2	$- {\mathop{\rm Im}\nolimits} \left[ {{T_{12}}} \right] + j{\mathop{\rm Im}\nolimits} \left[ {{T_{13}}} \right]$
${\mathop{\rm Re}\nolimits} \left[ {{T_{23}}\left( \theta \right)} \right]$	$\frac{1}{4}{\left( {{T_{33}} - {T_{22}}} \right)^2} + {{\mathop{\rm Re}\nolimits} ^2}\left[ {{T_{23}}} \right]$	0	4	$\frac{1}{2}\left( {{T_{33}} - {T_{22}}} \right) + j{\mathop{\rm Re}\nolimits} \left[ {{T_{23}}} \right]$
${T_{22}}\left( \theta \right)$	$\frac{1}{4}{\left( {{T_{33}} - {T_{22}}} \right)^2} + {{\mathop{\rm Re}\nolimits} ^2}\left[ {{T_{23}}} \right]$	$\frac{1}{2}\left( {{T_{22}} + {T_{33}}} \right)$	4	${\mathop{\rm Re}\nolimits} \left[ {{T_{23}}} \right] + j\frac{1}{2}\left( {{T_{22}} - {T_{33}}} \right)$
${T_{33}}\left( \theta \right)$	$\frac{1}{4}{\left( {{T_{33}} - {T_{22}}} \right)^2} + {{\mathop{\rm Re}\nolimits} ^2}\left[ {{T_{23}}} \right]$	$\frac{1}{2}\left( {{T_{22}} + {T_{33}}} \right)$	4	$- {\mathop{\rm Re}\nolimits} \left[ {{T_{23}}} \right] + j\frac{1}{2}\left( {{T_{33}} - {T_{22}}} \right)$
${\left\| {{T_{12}}\left( \theta \right)} \right\|^2}$	${{\mathop{\rm Re}\nolimits} ^2}\left[ {{T_{12}}T_{13}^ * } \right] + \frac{1}{4}{\left( {{{\left\| {{T_{13}}} \right\|}^2} - {{\left\| {{T_{12}}} \right\|}^2}} \right)^2}$	$\frac{1}{2}\left( {{{\left\| {{T_{12}}} \right\|}^2} + {{\left\| {{T_{13}}} \right\|}^2}} \right)$	4	${\mathop{\rm Re}\nolimits} \left[ {{T_{12}}T_{13}^ * } \right] + j\frac{1}{2}\left( {{{\left\| {{T_{12}}} \right\|}^2} - {{\left\| {{T_{13}}} \right\|}^2}} \right)$
${\left\| {{T_{13}}\left( \theta \right)} \right\|^2}$	${{\mathop{\rm Re}\nolimits} ^2}\left[ {{T_{12}}T_{13}^ * } \right] + \frac{1}{4}{\left( {{{\left\| {{T_{13}}} \right\|}^2} - {{\left\| {{T_{12}}} \right\|}^2}} \right)^2}$	$\frac{1}{2}\left( {{{\left\| {{T_{12}}} \right\|}^2} + {{\left\| {{T_{13}}} \right\|}^2}} \right)$	4	$- {\mathop{\rm Re}\nolimits} \left[ {{T_{12}}T_{13}^ * } \right] + j\frac{1}{2}\left( {{{\left\| {{T_{13}}} \right\|}^2} - {{\left\| {{T_{12}}} \right\|}^2}} \right)$
${\left\| {{T_{23}}\left( \theta \right)} \right\|^2}$	$\frac{1}{4}{\left\{ {\frac{1}{4}{{\left( {{T_{33}} - {T_{22}}} \right)}^2} + {{{\mathop{\rm Re}\nolimits} }^2}\left[ {{T_{23}}} \right]} \right\}^2}$	$\begin{array}{l}\frac{1}{2}\left\{ {\frac{1}{4}{{\left( {{T_{33}} - {T_{22}}} \right)}^2} + {{{\mathop{\rm Re}\nolimits} }^2}\left[ {{T_{23}}} \right]} \right\}\\\quad + {{\mathop{\rm Im}\nolimits} ^2}\left[ {{T_{23}}} \right]\end{array}$	8	$\begin{array}{l}\frac{1}{2}\left( {{T_{33}} - {T_{22}}} \right){\mathop{\rm Re}\nolimits} \left[ {{T_{23}}} \right] \\\quad+ j\frac{1}{2}\left[ {{{{\mathop{\rm Re}\nolimits} }^2}\left[ {{T_{23}}} \right] - \frac{1}{4}{{\left( {{T_{33}} - {T_{22}}} \right)}^2}} \right]\end{array}$

下载: 导出CSV

表 2 AIRSAR数据极化相干特征旋转域优化前后对比结果

Table 2. Comparison of polarimetric coherence without and with optimization in rotation domain for AIRSAR data

极化相干特征	原始值均值	最优值均值	增强百分比(%)
$\left\| {{\gamma _{\left( {{\rm{HH}} + {\rm{VV}}} \right) - \left( {{\rm{HH}} - {\rm{VV}}} \right)}}\left( \theta \right)} \right\|$	0.30	0.33	10.00
$\left\| {{\gamma _{\left( {{\rm{HH}} - {\rm{VV}}} \right) - {\rm{HV}}}}\left( \theta \right)} \right\|$	0.11	0.48	336.36
$\left\| {{\gamma _{{\rm{HH}} - {\rm{VV}}}}\left( \theta \right)} \right\|$	0.35	0.64	82.86
$\left\| {{\gamma _{{\rm{HH}} - {\rm{HV}}}}\left( \theta \right)} \right\|$	0.13	0.45	246.15

下载: 导出CSV

参考文献(49)

[1]	Lee J S and Pottier E. Polarimetric Radar Imaging: From Basics to Applications[M]. Boca Raton: CRC Press, 2009.
[2]	Cloude S R. Polarisation Application in Remote Sensing[M]. Oxford: Oxford University Press, 2009.
[3]	金亚秋, 徐丰. 极化散射与SAR遥感信息理论与方法[M]. 北京: 科学出版社, 2008. Jin Ya-qiu and Xu Feng. Theory and Approach for Polarimetric Scattering and Information Retrieval of SAR Remote Sensing[M]. Beijing: Science Press, 2008.
[4]	张红, 王超, 刘萌, 等. 极化SAR理论、方法与应用[M]. 北京: 科学出版社, 2015. Zhang Hong, Wang Chao, Liu Meng, et al.. Theory, Approach and Application of Polarimetric SAR[M]. Beijing: Science Press, 2015.
[5]	Dell’Acqua F and Gamba P. Remote sensing and earthquake damage assessment: Experiences, limits, and perspectives[J]. Proceedings of the IEEE, 2012, 100(10): 2876–2890. doi: 10.1109/JPROC.2012.2196404
[6]	Sato M, Chen S W, and Satake M. Polarimetric SAR analysis of Tsunami damage following the March 11, 2011 East Japan Earthquake[J]. Proceedings of the IEEE, 2012, 100(10): 2861–2875. doi: 10.1109/JPROC.2012.2200649
[7]	Chen S W and Sato M. Tsunami damage investigation of built-up areas using multitemporal spaceborne full polarimetric SAR images[J]. IEEE Transactions on Geoscience and Remote Sensing, 2013, 51(4): 1985–1997. doi: 10.1109/TGRS.2012.2210050
[8]	Chen S W, Wang X S, and Sato M. Urban damage level mapping based on scattering mechanism investigation using fully polarimetric SAR data for the 3.11 East Japan earthquake[J]. IEEE Transactions on Geoscience and Remote Sensing, 2016, 54(12): 6919–6929. doi: 10.1109/TGRS.2016.2588325
[9]	吴一戎. 多维度合成孔径雷达成像概念[J]. 雷达学报, 2013, 2(2): 135–142. http://radars.ie.ac.cn/CN/abstract/abstract93.shtml Wu Yi-rong. Concept of multidimensional space joint-observation SAR[J]. Journal of Radars, 2013, 2(2): 135–142. http://radars.ie.ac.cn/CN/abstract/abstract93.shtml
[10]	洪文. 圆迹SAR成像技术研究进展[J]. 雷达学报, 2012, 1(2): 124–135. http://radars.ie.ac.cn/CN/abstract/abstract29.shtml Hong Wen. Progress in circular SAR imaging technique[J]. Journal of Radars, 2012, 1(2): 124–135. http://radars.ie.ac.cn/CN/abstract/abstract29.shtml
[11]	张杰, 张晰, 范陈清, 等. 极化SAR在海洋探测中的应用与探讨[J]. 雷达学报, 2016, 5(6): 596–606. http://radars.ie.ac.cn/CN/abstract/abstract390.shtml Zhang Jie, Zhang Xi, Fan Chen-qing, et al.. Discussion on application of polarimetric synthetic aperture radar in marine surveillance[J]. Journal of Radars, 2016, 5(6): 596–606. http://radars.ie.ac.cn/CN/abstract/abstract390.shtml
[12]	许成斌, 周伟, 丛瑜, 等. 基于峰值区域的高分辨率极化SAR舰船目标特征分析与鉴别[J]. 雷达学报, 2015, 4(3): 367–373. http://radars.ie.ac.cn/CN/abstract/abstract213.shtml Xu Cheng-bin, Zhou Wei, Cong Yu, et al.. Ship analysis and detection in high-resolution Pol-SAR imagery based on peak zone[J]. Journal of Radars, 2015, 4(3): 367–373. http://radars.ie.ac.cn/CN/abstract/abstract213.shtml
[13]	王雪松. 雷达极化技术研究现状与展望[J]. 雷达学报, 2016, 5(2): 119–131. http://radars.ie.ac.cn/CN/abstract/abstract346.shtml Wang Xue-song. Status and prospects of radar polarimetry techniques[J]. Journal of Radars, 2016, 5(2): 119–131. http://radars.ie.ac.cn/CN/abstract/abstract346.shtml
[14]	杨汝良, 戴博伟, 李海英. 极化合成孔径雷达极化层次和系统工作方式[J]. 雷达学报, 2016, 5(2): 132–142. http://radars.ie.ac.cn/CN/abstract/abstract337.shtml Yang Ru-liang, Dai Bo-wei, and Li Hai-ying. Polarization hierarchy and system operating architecture for polarimetric synthetic aperture radar[J]. Journal of Radars, 2016, 5(2): 132–142. http://radars.ie.ac.cn/CN/abstract/abstract337.shtml
[15]	代大海, 廖斌, 肖顺平, 等. 雷达极化信息获取与处理的研究进展[J]. 雷达学报, 2016, 5(2): 143–155. http://radars.ie.ac.cn/CN/abstract/abstract316.shtml Dai Da-hai, Liao Bin, Xiao Shun-ping, et al.. Advancements on radar polarization information acquisition and processing[J]. Journal of Radars, 2016, 5(2): 143–155. http://radars.ie.ac.cn/CN/abstract/abstract316.shtml
[16]	赵春雷, 王亚梁, 阳云龙, 等. 雷达极化信息获取及极化信号处理技术研究综述[J]. 雷达学报, 2016, 5(6): 620–638. http://radars.ie.ac.cn/CN/abstract/abstract389.shtml Zhao Chun-lei, Wang Ya-liang, Yang Yun-long, et al.. Review of radar polarization information acquisition and polarimetric signal processing techniques[J]. Journal of Radars, 2016, 5(6): 620–638. http://radars.ie.ac.cn/CN/abstract/abstract389.shtml
[17]	Huynen J R. Phenomenological Theory of Radar Targets[D]. [Ph.D. dissertation], Delft University of Technology, 1970.
[18]	Cloude S R and Pottier E. An entropy based classification scheme for land applications of polarimetric SAR[J]. IEEE Transactions on Geoscience and Remote Sensing, 1997, 35(1): 68–78. doi: 10.1109/36.551935
[19]	Touzi R. Target scattering decomposition in terms of roll-invariant target parameters[J]. IEEE Transactions on Geoscience and Remote Sensing, 2007, 45(1): 73–84. doi: 10.1109/TGRS.2006.886176
[20]	Paladini R, Martorella M, and Berizzi F. Classification of man-made targets via invariant coherency-matrix eigenvector decomposition of polarimetric SAR/ISAR images[J]. IEEE Transactions on Geoscience and Remote Sensing, 2011, 49(8): 3022–3034. doi: 10.1109/TGRS.2011.2116121
[21]	Chen S W, Li Y Z, Wang X S, et al.. Modeling and interpretation of scattering mechanisms in polarimetric synthetic aperture radar: Advances and perspectives[J]. IEEE Signal Processing Magazine, 2014, 31(7): 79–89.
[22]	Xu F and Jin Y Q. Deorientation theory of polarimetric scattering targets and application to terrain surface classification[J]. IEEE Transactions on Geoscience and Remote Sensing, 2005, 43(10): 2351–2364. doi: 10.1109/TGRS.2005.855064
[23]	An W T, Cui Y, and Yang J. Three-component model-based decomposition for polarimetric SAR data[J]. IEEE Transactions on Geoscience and Remote Sensing, 2010, 48(6): 2732–2739. doi: 10.1109/TGRS.2010.2041242
[24]	Yamaguchi Y, Sato A, Boerner W M, et al.. Four-component scattering power decomposition with rotation of coherency matrix[J]. IEEE Transactions on Geoscience and Remote Sensing, 2011, 49(6): 2251–2258. doi: 10.1109/TGRS.2010.2099124
[25]	Van Zyl J J, Arii M, and Kim Y. Model-based decomposition of polarimetric SAR covariance matrices constrained for nonnegative eigenvalues[J]. IEEE Transactions on Geoscience and Remote Sensing, 2011, 49(9): 3452–3459. doi: 10.1109/TGRS.2011.2128325
[26]	Arii M, Van Zyl J J, and Kim Y. Adaptive model-based decomposition of polarimetric SAR covariance matrices[J]. IEEE Transactions on Geoscience and Remote Sensing, 2011, 49(3): 1104–1113. doi: 10.1109/TGRS.2010.2076285
[27]	Neumann M, Ferro-Famil L, and Reigber A. Estimation of forest structure, ground, and canopy layer characteristics from multibaseline polarimetric interferometric SAR data[J]. IEEE Transactions on Geoscience and Remote Sensing, 2010, 48(3): 1086–1104. doi: 10.1109/TGRS.2009.2031101
[28]	Chen S W, Wang X S, Xiao S P, et al.. General polarimetric model-based decomposition for coherency matrix[J]. IEEE Transactions on Geoscience and Remote Sensing, 2014, 52(3): 1843–1855. doi: 10.1109/TGRS.2013.2255615
[29]	Ballester-Berman J D and Lopez-Sanchez J M. Applying the Freeman-Durden decomposition concept to polarimetric SAR interferometry[J]. IEEE Transactions on Geoscience and Remote Sensing, 2010, 48(1): 466–479. doi: 10.1109/TGRS.2009.2024304
[30]	Chen S W, Wang X S, Li Y Z, et al.. Adaptive model-based polarimetric decomposition using PolInSAR coherence[J]. IEEE Transactions on Geoscience and Remote Sensing, 2014, 52(3): 1705–1718. doi: 10.1109/TGRS.2013.2253780
[31]	Wang C L, Yu W D, Wang R, et al.. Comparison of nonnegative eigenvalue decompositions with and without reflection symmetry assumptions[J]. IEEE Transactions on Geoscience and Remote Sensing, 2014, 52(4): 2278–2287. doi: 10.1109/TGRS.2013.2259177
[32]	Li H Z, Wang C, Zhang H, et al.. A unified three-component scattering model for polarimetric coherent target decomposition[J]. International Journal of Remote Sensing, 2012, 33(9): 2868–2891. doi: 10.1080/01431161.2011.622728
[33]	Zhu F Y, Zhang Y H, and Li D. An extension of a complete model-based decomposition of polarimetric SAR data[J]. IEEE Geoscience and Remote Sensing Letters, 2016, 13(2): 287–291. doi: 10.1109/LGRS.2015.2511076
[34]	Deng L and Yan Y N. Improving the Yamaguchi4 decomposition method using selective polarization orientation compensation[J]. Canadian Journal of Remote Sensing, 2016, 42(2): 125–135. doi: 10.1080/07038992.2016.1160774
[35]	Zou B, Zhang Y, Cao N, et al.. A four-component decomposition model for PolSAR data using asymmetric scattering component[J]. IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing, 2015, 8(3): 1051–1061. doi: 10.1109/JSTARS.2014.2380151
[36]	Wang Y H, Liu H W, and Jiu B. PolSAR coherency matrix decomposition based on constrained sparse representation[J]. IEEE Transactions on Geoscience and Remote Sensing, 2014, 52(9): 5906–5922. doi: 10.1109/TGRS.2013.2293663
[37]	Cloude S R and Pottier E. A review of target decomposition theorems in radar polarimetry[J]. IEEE Transactions on Geoscience and Remote Sensing, 1996, 34(3): 498–518.
[38]	张腊梅, 段宝龙, 邹斌. 极化SAR图像目标分解方法的研究进展[J]. 电子与信息学报, 2016, 38(12): 3289–3297. http://youxian.cnki.com.cn/yxdetail.aspx?filename=LDAX20170807001&dbname=CAPJ2015 Zhang Lamei, Duan Baolong, and Zou Bin. Research development on target decomposition method of polarimetric SAR image[J]. Journal of Electronics & Information Technology, 2016, 38(12): 3289–3297. http://youxian.cnki.com.cn/yxdetail.aspx?filename=LDAX20170807001&dbname=CAPJ2015
[39]	Chen S W, Ohki M, Shimada M, et al.. Deorientation effect investigation for model-based decomposition over oriented built-up areas[J]. IEEE Geoscience and Remote Sensing Letters, 2013, 10(3): 273–277.
[40]	Chen S W, Wang X S, and Sato M. Uniform polarimetric matrix rotation theory and its applications[J]. IEEE Transactions on Geoscience and Remote Sensing, 2014, 52(8): 4756–4770. doi: 10.1109/TGRS.2013.2284359
[41]	Chen S W, Li Y Z, and Wang X S. A visualization tool for polarimetric SAR data investigation[C]. The 11th European Conference on Synthetic Aperture Radar, Hamburg, Germany, 2016: 579–582.
[42]	Chen S W and Wang X S. Polarimetric coherence pattern: A visualization tool for PolSAR data investigation[C]. The IEEE International Geoscience and Remote Sensing Symposium (IGARSS), Beijing, China, 2016: 7509–7512.
[43]	Xiao S P, Chen S W, Chang Y L, et al.. Polarimetric coherence optimization and its application for manmade target extraction in PolSAR data[J]. IEICE Transactions on Electronics, 2014, 97(6): 566–574.
[44]	Chen S W, Li Y Z, and Wang X S. Crop discrimination based on polarimetric correlation coefficients optimization for PolSAR data[J]. International Journal of Remote Sensing, 2015, 36(16): 4233–4249. doi: 10.1080/01431161.2015.1079345
[45]	Lee J S, Schuler D L, Ainsworth T L, et al.. On the estimation of radar polarization orientation shifts induced by terrain slopes[J]. IEEE Transactions on Geoscience and Remote Sensing, 2002, 40(1): 30–41. doi: 10.1109/36.981347
[46]	Chen S W, Wang X S, and Sato M. PolInSAR complex coherence estimation based on covariance matrix similarity test[J]. IEEE Transactions on Geoscience and Remote Sensing, 2012, 50(11): 4699–4710. doi: 10.1109/TGRS.2012.2192937
[47]	陶臣嵩, 陈思伟, 李永祯, 等. 结合旋转域极化特征的极化SAR地物分类[J]. 雷达学报, 录用待刊. doi: 10.12000/JR16131 Tao Chensong, Chen Siwei, Li Yongzhen, et al. Polarimetric SAR terrain classification using polarimetric feature derived from rotation domain[J]. Journal of Radars. doi: 10.12000/JR16131
[48]	Ainsworth T L, Schuler D L, and Lee J S. Polarimetric SAR characterization of man-made structures in urban areas using normalized circular-pol correlation coefficients[J]. Remote Sensing of Environment, 2008, 112(6): 2876–2885. doi: 10.1016/j.rse.2008.02.005
[49]	Yamaguchi Y, Yamamoto Y, Yamada H, et al.. Classification of terrain by implementing the correlation coefficient in the circular polarization basis using X-band POLSAR data[J]. IEICE Transactions on Communications, 2008, E91B(1): 297–301.

施引文献

期刊类型引用(4)

1.	张晶，王博，孟国祥，陈双艳. 频控阵雷达对雷达侦察测向的欺骗机理研究. 航空兵器. 2024(01): 71-76 . 百度学术
2.	冯一伦，蒋彦雯，曾桂兰，范红旗. 干涉仪对频分正交MIMO雷达的测向性能分析. 航空兵器. 2024(01): 103-110 . 百度学术
3.	唐硕，臧伟旺，李震宇，高兵，张怀根. 机载雷达无源探测技术研究. 电子技术与软件工程. 2022(11): 180-183 . 百度学术
4.	王文钦，张顺生. 频控阵雷达技术研究进展综述. 雷达学报. 2022(05): 830-849 . 本站查看