作者中心
专家审稿
责编办公
编辑办公
-
摘要: 地基雷达是近20几年逐渐发展成熟的微波遥感成像技术，目前已广泛应用于滑坡、崩塌等地质灾害的监测中。地基雷达通过干涉测量原理可以监测到目标区域发生的微小形变，然而受人为因素、地质因素、气象因素等影响，导致雷达图像失相干严重，给长期定量化监测带来较大的难度。因此，迫切需要在定量监测的基础上，进一步开展变化检测方面的应用，为长期全面了解监测区域的动态变化提供有效信息。针对上述问题，该文提出了一种基于改进的模糊C均值聚类(FCM)算法对地基雷达图像进行无监督变化检测，该方法首次利用相干系数图和均值对数比值图进行非下采样轮廓波变换(NSCT)和局部能量法得到合成差异图，然后利用主成分分析(PCA)提取合成差异图中每个像素的特征向量，根据地基雷达图像特点对FCM进行改进，通过改进的FCM对每个像素的特征向量进行聚类得到最终的变化检测结果。利用地基雷达LSA对中国西南某省出现的堰塞体的治理过程进行监测，获取监测区域的地基雷达图像，监测过程中受降水等影响监测体出现滑坡，使用该文方法对其进行变化检测，结果表明该文方法更容易进行聚类分割，变化检测结果在保留变化区域的同时噪声点明显减少。
-
关键词:
- 地基雷达图像 /
- 变化检测 /
- 无监督 /
- 相干系数 /
- 改进的模糊C均值聚类
Abstract: Ground-based radar is a microwave remote sensing imaging technology that has been gradually developed throughout the past 20 years so that it has become mature. At present, it has been widely used in monitoring geological disasters such as landslides and collapses. Ground-based radars can detect micro-variations in target areas through the principle of interferometry. However, due to human factors, geological factors, and meteorological factors, the radar image of the monitored area is incoherent, which makes long-term quantitative monitoring difficult. Therefore, further developing the application of change detection while considering quantitative monitoring is urgent, to provide effective information on long-term changes and comprehensively understand the dynamic changes in the monitored area. To solve the above problems, an unsupervised change detection method using ground-based radar images and based on an improved Fuzzy C-Means clustering (FCM) algorithm is proposed in this paper. In this method, for the first time, the Nonsubsampled Contourlet Transform (NSCT) is performed on the coherence coefficient map and the mean log ratio map to obtain the fusion difference map. Then, principal component analysis is used to extract the feature vectors of each pixel in the fusion difference image. The FCM is improved according to the characteristics of the ground-based radar images. The improved FCM is used to cluster the feature vectors of each pixel to obtain the change detection result. A ground-based radar LSA was used to monitor the treatment process of a dam in southwest China. During the monitoring process, landslides occurred in the monitored area affected by precipitation and other factors. This method is used to detect the change of the radar image before and after the landslide. The results show that the proposed method allows for easier clustering and segmenting, and the change detection results can significantly reduce the noise points while retaining the change area. -
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中Tm 表示第m号发射天线，Rn 表示第n号接收天线；Snm(t)表示第n号发射天线发射信号、第m号接收天线接收到的回波信号，回波信号已经过脉冲压缩处理。
伪随机编码雷达回波经过脉冲压缩处理后得到窄脉冲信号，如图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，其中fS是FPGA内部DAC采样数据存储器ROM1～ROM16的工作时钟，fD是DDR接口的工作时钟，fDAC是外部高速DAC的工作时钟。根据并串转换、DDR接口和四端口DAC芯片的时钟关系，有fDAC=8fD和fD=2fS，即FPGA内部只需要低速时钟便能产生高速的伪随机编码信号。DAC产生的信号通过信号频带内的低通或带通滤波器后便能产生平滑的编码信号，编码信号经放大器放大后便是雷达天线的发射信号。
图 2 高速伪随机编码信号产生原理图Figure 2. Schematic diagram of high-speed pseudo random coded signal2.4 超宽带信号的混合采样技术
对于超宽带信号的采集，目前市场上能满足超宽带信号采样率要求的器件很少见，所以通常会采用混合采样对超宽带信号进行采集[18]。混合采样的思想是牺牲时间来换取更高的采样率，在多个雷达脉冲重复周期Pulse Repetition Period (PRP)内通过多次发送相同的信号来实现一道完整雷达回波的采集。图3是混合采样过程的时序图，采样过程中Analog to Digital Converter (ADC)工作时钟fADC依次比上一轮的混合采样时钟延迟 Δt。则混合采样的采样率为：
fH=1Δt (1) 采集一道完整回波的混合采样次数为：
n=1fADC⋅Δt (2) 假设每一轮混合采样的采样点数为NP，则n轮混合采样得到的数据是一个n×NP矩阵，如式(3)所示。矩阵中每一轮的采样数据按行存储，而按列读出即是相邻时间间隔为 Δt的采样数据。因此，矩阵中的元素按依次从上到下、从左到右读取便可以得到一道顺序排列好的完整回波。
S=[S1(0)S1(1)⋯S1(NP−1)S2(0)S2(1)⋯S2(NP−1)⋮⋮⋱⋮Sn(0)Sn(1)⋯Sn(NP−1)] (3) 2.5 回波信号的实时脉冲压缩
在伪随机编码体制的超宽带雷达中，原始回波信号的实时脉冲压缩是信号处理的首要和关键步骤。由于超宽带雷达的采样率高、数据量大，而微处理器处理速度和DSP芯片运算能力有限。本文基于时间域的互相关算法，通过调用FPGA中的DSP硬核进行并行计算，并结合流水线模式，实现快速实时脉冲压缩。伪随机编码信号的脉冲压缩公式为：
R[l]=Nr−1∑i=0s[i+l]r[i], l=0,1,2,···,NR−1 (4) 式中，r是参考码信号，它通常为雷达的发射信号，长度为Nr个采样点；s是雷达回波信号，长度为Ns个采样点。根据式(4)，调用FPGA中的Nr个DSP硬核进行乘法运算，并通过流水线方式进行累加操作，构造图4所示的脉冲压缩并行运算框架。图中D表示1个时钟周期的延时，DSP硬核乘法器的一个输入端分别为各自固定的参考码，另一个输入端同时输入回波信号，压缩结果由最后一级输出端Rsr[i]输出。
由于FPGA中的DSP硬核数量有限，如果参考码长度Nr大于乘加器数量，脉冲压缩并行运算将不能实现。为了解决上面的问题，把长度为Nr的参考码分解成M段，每段长度为Nr0=Nr/M；分段条件满足Nr0小于FPGA中DSP硬核数量，同时通过对参考码补0或者适当选择M值使Nr0为整数。由分段方法，可以得到分段后的第j段参考码与原参考码对应关系为：
rj[i]=r[j⋅Nr0+i] (5) 式中，j=0, 1, ···,M–1,i=0, 1, ···,Nr0–1。分段后，式(4)可以转化成参考码分段并行运算脉冲压缩，如式(6)：
R[l]=M−1∑j=0Rj[l] (6) 式中，l=0, 1, ···,NR–1,Rj 表示第j段参考码的压缩结果，为：
Rj[l]=Nr0−1∑i=0s[j⋅Nr0+l+i]⋅rj[i] (7) 即先把各段参考码采用图4的并行运算框架进行脉冲压缩，再把压缩的结果累加，容易证明式(4)和式(6)是等价的。
3. 多发多收雷达系统设计
多发多收雷达系统的原理图如图5所示，主要由4部分组成，从上到下依次为：时钟同步单元、发射机单元、接收机单元和开关切换单元，图中分别用不同的填充色标注。图6是多发多收雷达系统样机，该系统的主要参数如表1所示。
表 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) 164m⋅n 总功耗(W) 20 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所示。
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序列的特性一致，所以接收机单元满足雷达系统设计的需求。
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 为了测试4发16收各个通道的一致性，对雷达系统4发16收开关的任意收发组合进行闭环测试，最终闭环测试得到的64组数据结果如图10所示。从图10可以知道，64组数据的脉冲压缩波形近似相同、幅度近似相等，除了在回波延时存在一定的差别外，其他特性基本相一致。而对于回波的延时差异，通过雷达系统校准可以在数据后处理时补偿回来。
对系统的接收灵敏度进行测试，任意选择收发组合通道，发射信号和接收信号都经过低噪声放大器，中间级联可调衰减器后闭环。通过调节衰减器的大小，测试最小能检测到的信号电平。不同衰减值时的测试结果如图11所示，当衰减120 dB时虽然看到峰值点，事实上难以分辨。所以可以认为衰减115 dB时得到系统接收机的最小检测电平，由发射处输出电平为5 Vpp，通过计算得到接收机的灵敏度为–97 dBm。
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以内。脉冲雷达随着探测目标距离增大，信噪比迅速下降；而伪随机编码雷达可以通过脉冲压缩而进一步提高系统的信噪比。
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是由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测试墙体为一般楼房的砖墙，厚度为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 地基雷达图像无监督变化检测流程图
Figure 1. Flow chart for unsupervised change detection based on ground-based radar image
图 8 变化检测结果与地形3维数据融合图
Figure 8. Effect of fusion of change detection results and UAV tilt photography
表 1 改进的FCM算法计算步骤
Table 1. Improved FCM algorithm calculation steps
算法：改进的FCM算法计算步骤 输入：融合差异图每个像素的特征向量 Yη (1) 初始化隶属度 u ，满足式(10)；(2) 通过式(12)计算聚类中 Ci ；(3) 根据式(17)计算目标函数 J ；(4) 迭代计算：当迭代次数 ≤q 时，根据式(18)更新隶属度矩阵，返回步骤(2)继续执行 当迭代次数 >q 时，停止迭代； 输出：隶属度矩阵 uη1 和uη2 
下载: 导出CSV
表 2 地基雷达LSA系统参数
Table 2. Ground-based radar LSA system parameters
参数 数值 参数 数值 频率范围 16.5～17.5 GHz 监测范围(方位向) 60°× 30° 形变监测精度 0.1 mm 脉冲重复频率 优于500 Hz 作用距离 5 km 位移观测周期 240～600 s 图像分辨率 0.3 m × 0.0054 rad 供电电源 AC 220 V/DC 24 V 设备总重 ≤150kg 工作温度 –30～ 50 ℃
下载: 导出CSV
-
[1] WANG Yan, DU Lan, and DAI Hui. Unsupervised SAR image change detection based on SIFT keypoints and region information[J]. IEEE Geoscience and Remote Sensing Letters, 2016, 13(7): 931–935. doi: 10.1109/LGRS.2016.2554606 [2] GONG Maoguo, SU Linzhi, JIA Meng, et al. Fuzzy clustering with a modified MRF energy function for change detection in synthetic aperture radar images[J]. IEEE Transactions on Fuzzy Systems, 2014, 22(1): 98–109. doi: 10.1109/TFUZZ.2013.2249072 [3] GAO Feng, DONG Junyu, LI Bo, et al. Automatic change detection in synthetic aperture radar images based on PCANet[J]. IEEE Geoscience and Remote Sensing Letters, 2016, 13(12): 1792–1796. doi: 10.1109/LGRS.2016.2611001 [4] GONG Maoguo, CAO Yu, and WU Qiaodi. A neighborhood-based ratio approach for change detection in SAR images[J]. IEEE Geoscience and Remote Sensing Letters, 2012, 9(2): 307–311. doi: 10.1109/LGRS.2011.2167211 [5] GONG Maoguo, ZHOU Zhiqiang, and MA Jingjing. Change detection in synthetic aperture radar images based on image fusion and fuzzy clustering[J]. IEEE Transactions on Image Processing, 2012, 21(4): 2141–2151. doi: 10.1109/TIP.2011.2170702 [6] SUMAIYA M N and KUMARI R S S. Logarithmic mean-based thresholding for SAR image change detection[J]. IEEE Geoscience and Remote Sensing Letters, 2016, 13(11): 1726–1728. doi: 10.1109/LGRS.2016.2606119 [7] 杨祥立, 徐德伟, 黄平平, 等. 融合相干/非相干信息的高分辨率SAR图像变化检测[J]. 雷达学报, 2015, 4(5): 582–590. doi: 10.12000/JR15073YANG Xiangli, XU Dewei, HUANG Pingping, et al. Change detection of high resolution SAR images by the fusion of coherent/incoherent information[J]. Journal of Radars, 2015, 4(5): 582–590. doi: 10.12000/JR15073 [8] ELGUEBALY T and BOUGUILA N. A Bayesian approach for SAR images segmentation and changes detection[C]. The 25th Biennial Symposium on Communications, Kingston, Canada, 2010: 24–27. doi: 10.1109/BSC.2010.5473011. [9] 邵宁远, 邹焕新, 陈诚, 等. 面向变化检测的SAR图像超像素协同分割算法[J]. 系统工程与电子技术, 2019, 41(7): 1496–1503. doi: 10.3969/j.issn.1001-506X.2019.07.09SHAO Ningyuan, ZOU Huanxin, CHEN Cheng, et al. Change detection oriented superpixel cosegmentation algorithm for SAR images[J]. Systems Engineering and Electronics, 2019, 41(7): 1496–1503. doi: 10.3969/j.issn.1001-506X.2019.07.09 [10] CARINCOTTE C, DERRODE S, and BOURENNANE S. Unsupervised change detection on SAR images using fuzzy hidden Markov chains[J]. IEEE Transactions on Geoscience and Remote Sensing, 2006, 44(2): 432–441. doi: 10.1109/TGRS.2005.861007 [11] BEZDEK J C, EHRLICH R, and FULL W. FCM: The fuzzy c-means clustering algorithm[J]. Computers & Geosciences, 1984, 10(2/3): 191–203. doi: 10.1016/0098-3004(84)90020-7 [12] PHAM D L. Spatial models for fuzzy clustering[J]. Computer Vision and Image Understanding, 2001, 84(2): 285–297. doi: 10.1006/cviu.2001.0951 [13] TRAN T N, WEHRENS R, and BUYDENS L M C. Clustering multispectral images: A tutorial[J]. Chemometrics and Intelligent Laboratory Systems, 2005, 77(1/2): 3–17. doi: 10.1016/j.chemolab.2004.07.011 [14] LI Zhenfang, BAO Zheng, LI Hai, et al. Image autocoregistration and InSAR interferogram estimation using joint subspace projection[J]. IEEE Transactions on Geoscience and Remote Sensing, 2006, 44(2): 288–297. doi: 10.1109/TGRS.2005.860984 [15] HYDE P, DUBAYAH R, WALKER W, et al. Mapping forest structure for wildlife habitat analysis using multi-sensor (LiDAR, SAR/InSAR, ETM+, Quickbird) synergy[J]. Remote Sensing of Environment, 2006, 102(1/2): 63–73. doi: 10.1016/j.rse.2006.01.021 [16] ADU J, GAN Jianhong, WANG Yan, et al. Image fusion based on nonsubsampled contourlet transform for infrared and visible light image[J]. Infrared Physics & Technology, 2013, 61: 94–100. doi: 10.1016/j.infrared.2013.07.010 [17] WU Chao and WU Yiquan. Multitemporal images change detection using nonsubsampled contourlet transform and kernel fuzzy C-means clustering[C]. The 2nd International Symposium on Intelligence Information Processing and Trusted Computing, Wuhan, China, 2011: 22–23. doi: 10.1109/IPTC.2011.31. [18] 李青松, 覃锡忠, 贾振红, 等. 非下采样Contourlet域融合和参数化内核图割的SAR图像无监督水灾变化检测[J]. 中国图象图形学报, 2014, 19(6): 971–978. doi: 10.11834/jig.20140619LI Qingsong, QIN Xizhong, JIA Zhenhong, et al. Unsupervised detection of flood changes with SAR images combining nonsubsampled Contourlet domain fusion and parametric kernel graph cuts[J]. Journal of Image and Graphics, 2014, 19(6): 971–978. doi: 10.11834/jig.20140619 [19] 刘陆洋, 贾振红, 杨杰, 等. 利用双差异图和PCA的SAR图像变化检测[J]. 计算机工程与设计, 2019, 40(7): 2002–2006.LIU Luyang, JIA Zhenhong, YANG Jie, et al. SAR image chang detective using double difference images and PCA algorithm[J]. Computer Engineering and Design, 2019, 40(7): 2002–2006. [20] 慕彩红, 霍利利, 刘逸, 等. 基于小波融合和PCA-核模糊聚类的遥感图像变化检测[J]. 电子学报, 2015, 43(7): 1375–1381. doi: 10.3969/j.issn.0372-2112.2015.07.019MU Caihong, HUO Lili, LIU Yi, et al. Change detection for remote sensing images based on wavelet fusion and PCA-Kernel fuzzy clustering[J]. Acta Electronica Sinica, 2015, 43(7): 1375–1381. doi: 10.3969/j.issn.0372-2112.2015.07.019 [21] EL-DAWLATLY S, OSMAN H, and SHAHEIN H I. New spatial FCM approach with application to SAR target clustering[C]. The 8th International Conference on Signal Processing, Beijing, China, 2006. doi: 10.1109/ICOSP.2006.345896. 期刊类型引用(4)
1. 张晶，王博，孟国祥，陈双艳. 频控阵雷达对雷达侦察测向的欺骗机理研究. 航空兵器. 2024(01): 71-76 . 
百度学术2. 冯一伦，蒋彦雯，曾桂兰，范红旗. 干涉仪对频分正交MIMO雷达的测向性能分析. 航空兵器. 2024(01): 103-110 . 百度学术3. 唐硕，臧伟旺，李震宇，高兵，张怀根. 机载雷达无源探测技术研究. 电子技术与软件工程. 2022(11): 180-183 . 百度学术4. 王文钦，张顺生. 频控阵雷达技术研究进展综述. 雷达学报. 2022(05): 830-849 . 
本站查看其他类型引用(3)
-
下载:
百度学术
计量
- 文章访问数: 3403
- HTML全文浏览量: 1301
- PDF下载量: 255
- 被引次数: 7
