-
摘要: 地基雷达是近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]。太赫兹辐射的光子能量低,对穿透物不会造成损伤,并且可以穿过大多数介电物质,实现无损检测。太赫兹波具有穿透性,能够实现对隐蔽物体的有效检测,可应用于安检相关的领域。太赫兹频段相比于微波频段频率更高,更容易发射大带宽信号,具有更高的分辨率,具有海量的频谱资源,可用于超宽带超高速无线通信。太赫兹波段目标表面的细微结构、粗糙度等细节会显著影响其后向散射特性,实现更小尺寸目标的探测、更精确目标的运动与物理参数反演[3]。太赫兹(terahertz, THz)波段位于微波与红外波之间,其频率范围为0.1~10 THz (1 THz=1012 Hz),对应的波长为30 μm~3 mm。太赫兹频段目标散射特性是太赫兹雷达探测和成像应用的物理基础[4,5],同时也是太赫兹雷达系统进行链路设计、特征提取以及成像算法的重要依据。国内首都师范大学太赫兹实验室研制了太赫兹数字全息成像系统,对太赫兹电磁波的振幅、相位、频率及偏振等全部光学信息的3维空间分布进行精确测量[6]。针对太赫兹波段目标的散射特性,美国麻省LOWELL大学毫米波实验室利用1.56 THz源在紧缩场中对粗糙表面圆柱体的目标散射特性进行了研究[4]。天津大学太赫兹研究中心搭建了以0.2 THz返波管振荡器源、热释电探测器、小型自动旋转光学平台等组成的太赫兹波目标散射特性实验测试系统,并对粗糙铜面的散射特性等进行了研究[7,8]。对于介质[9]和涂覆目标的太赫兹散射,北京航空航天大学江月松等人考虑粗糙度修正表面的散射系数研究了基于经验公式的涂覆目标的太赫兹散射特性[10]。
本文区别于以往采用经验公式[10]以粗糙度修正散射系数的研究方法,把随机粗糙面的建模理念应用到太赫兹波段表面粗糙目标的建模中。首先模拟生成了分形粗糙面近似代替实际复杂的粗糙面,对生成的分形粗糙面进行坐标变换导入计算机辅助设计(Computer Aided Design, CAD)建模软件建立具有粗糙表面的目标模型;然后对表面粗糙目标按照入射波的频率以满足物理光学近似的要求进行剖分。根据菲涅尔反射系数求得表面电流进而计算涂覆粗糙目标的雷达散射截面(Radar Cross Section, RCS)。并针对不同频率以及不同涂覆厚度的表面粗糙涂覆目标,分别进行了仿真分析。
2. 表面粗糙目标模型
2.1 分形粗糙面
自1982年Mandelbrot首次提出“分形”的概念[11],指的是组成部分与整体以某种方式相形似,分形理论就在很多领域中得到应用。“分形”不同于通常意义上的长度、面积、体积等几何概念,分形内部的任何一个相对独立的部分,在一定程度上都应该是整体的再现和缩影,分形几何体内部存在无穷层次、具有见微知著、由点及面的自相似结构,即具有自相似性。由于粗糙面一般具有非线性的几何结构,因此采用非线性的方法模拟粗糙面更能反映其物理本质。自然界的许多物体,如地、海表面、植被和森林等都在一定尺度范围内存在统计意义上的自相似性,由此很多学者将分形理论应用于电磁散射领域中,用于粗糙面的模拟[12,13]。
1维带限Weierstrass-Mandelbrot分形函数的表达式为:
f(x)=√2δ[1−b(2D−4)]1/2[b(2D−4)N1−b(2D−4)(N2+1)]1/2⋅N2∑n=N1b(D−2)ncos(2πsbnx+φn) (1) 其中,
δ 为高度的均方根,b是空间基频,D为分形维数(1<D<2),s为标度因子(s=K0/2π , K0为空间波数),φn 为(0,2π) 上均匀分布的随机相位,该函数具有零均值。一般取b>1,b为有理数时,f(x)表现为周期函数;b为无理数时,f(x)为准周期函数。标度因子s决定频谱的位置,f(x)的无标度区间一般取(sbN1)−1 和(sbN2)−1 ,N=N2−N1+1 ,随着N的增加,越来越多的频率分量加到准周期。图1给出了1维分形粗糙面模型,当分维数D增加时,高频分量比重加大,低频分量作用减小,分形粗糙面的粗糙程度增大。根据瑞利判据,粗糙面相对于入射波的粗糙程度,除与粗糙面的高度函数有关还和入射波的频率有关。如普通的目标表面对于微波段来说是光滑的,但相对于太赫兹频段的波来说却是粗糙的。本文主要研究太赫兹波段下目标表面的微粗糙对其散射特性的影响。2.2 分形表面粗糙目标模型
目标表面粗糙度引起的表面起伏一般在其对应的光滑表面的法线方向[14]。因此,对于轴对称旋转目标而言,其表面的粗糙度可近似考虑为对应母线的起伏。将生成的1维分形粗糙面叠加到光滑目标模型对应的母线进行坐标变换,建立具有分形粗糙表面的目标模型。
对于如图2(a)所示的顶部为半球的粗糙钝锥模型,其母线可以表示为:
x={(r1+f(x))cosα,r1+Δhtanβ+f(x)cosβ,y>0y<0 (2) y={(r1+f(x))sinα,Δh+f(x)sinβ,y>0y<0 (3) 其中,r1为顶部半球半径,r2为底面半径,h为下部锥台高度,
β=atan((r2−r1)/h) 。将生成的圆锥母线导入CAD建模软件,对其绕Y轴旋转并进行坐标变换生成如图2(b)所示的具有分形粗糙表面的钝锥模型。2.3 算法理论
由Stratton-Chu积分公式,目标远区散射场利用物理光学可表示为[15]:
Es=−jk4πexp(−jkr)r∫sˆs×[ˆn×E−Z0ˆs×(ˆn×H)]exp(jkr⋅ˆs)ds (4) 其中,k和Z0分别为自由空间的波数和本征阻抗,
ˆs 为散射波的单位矢量,r为表面上一点的位置矢量,ˆn 为目标表面向外单位法矢量,E和H分别为边界上总的电场和总的磁场。涂覆介质表面的散射示意图如图3所示。其中
θi 为入射角,ˆi 和ˆs 分别为入射波和散射波的单位矢量,矢量ˆei∥ 和ˆer∥ 分别为入射电场、反射电场平行入射面的极化方向,矢量ˆe⊥ 为入射电场和反射电场垂直入射面的极化方向。Ei=E⊥ˆe⊥+E∥ˆe∥i,Es=R⊥E⊥ˆe⊥+R∥E∥ˆe∥r (5) 其中,
Ei 为边界上入射电场,Es 为边界上散射电场,E⊥ 和E∥ 分别为入射电场在ˆe⊥ 和ˆei∥ 方向的场分量,R⊥ 和R∥ 分别为涂覆介质表面在垂直极化和水平极化时的反射系数[16]。涂覆目标雷达散射截面的计算公式为:
σ=limR→∞4πR2|Es|2|Ei|2 (6) 3. 数值结果及讨论
3.1 验证算例
为了验证算法的正确性,先通过下面的模型算例加以说明。图4给出了3 GHz平面波TM极化入射下涂覆半球的双站雷达散射截面,其中半球的半径为0.5 m,涂覆厚度为d=2 cm,涂层介质相对介电常数为
εr=36.0 ,相对磁导率为μr=1.0 。RCS结果曲线可以看出物理光学法和多层快速多极子方法(MLFMA)吻合良好,验证了程序的正确性。图5给出了频率为3 THz的平面波入射下导体立方体的单站雷达散射截面,结果与文献[3]中采用多层快速多极子方法结果一致,可以看出物理光学方法用于计算THz频段目标散射的有效性。
3.2 数值结果
对于图2(b)所示的具有分形粗糙表面的钝锥模型,其顶部半球半径r1=1 mm,底面半径r2=3 mm,锥台高度h=12 mm,分形粗糙面的分维数D=1.5,b=1.5,均方根高度
δ=0.02mm 。涂覆材料相对介电常数εr=(4.0,−1.5) ,相对磁导率μr=(2.0,−1.0) ,涂覆层厚度d=0.03 mm。首先对钝锥导入CAD建模软件进行满足物理光学近似的网格剖分,根据菲涅尔反射系数得出钝锥表面电流分布进而计算其散射场。从图6中结果可以看出,对于模型尺寸相同的光滑钝锥与表面粗糙钝锥的单站雷达散射截面曲线走势基本一致,随着入射角的增大,RCS增大,垂直于锥面照射时达到最大峰值。图6(a)入射频率为30 GHz的情况下光滑钝锥与分形粗糙钝锥的RCS除了小角度基本上重合,可以看出在微波频段目标表面的微粗糙度对RCS的影响很小,可以忽略。图6(b)、图6(c)表明太赫兹波段下光滑钝锥和分形粗糙钝锥目标雷达散射截面出现差异,表面的分形粗糙度引起目标RCS曲线震荡起伏,且频率越高起伏越明显,曲线波动越大。因此在太赫兹波段,目标表面的粗糙度对其散射特性的影响需要考虑。
图7给出了入射波频率为3 THz的不同涂层厚度的粗糙表面目标的后向RCS。可以看出相对于表面为导体的情况,涂覆介质以后,钝锥目标的雷达散射截面几乎在所有角度都有明显减小,并且随着涂层厚度的增大,雷达散射截面持续减小。涂覆介质层对雷达散射截面的缩减有明显的作用,在一定范围内随着涂层厚度的增大,涂覆介质对电磁波的吸收增加表面粗糙钝锥的后向RCS减小。
图8给出了不同入射频率下钝锥单站RCS。随着频率的升高,表面粗糙钝锥的后向RCS多数角度下降,且频率越高RCS值下降得越多。随着频率的增大,入射波的波长变小,目标表面的粗糙度与入射波长的比值增大,粗糙度引起的漫散射效应增大,目标RCS受到表面粗糙度的影响,曲线峰值变得不明显。
图9给出了不同表面粗糙度的圆柱模型单站雷达散射截面,其半径为r=16.25 mm,高度为h=102 mm,入射波频率为0.3 THz。
图10给出了不同表面粗糙度的锥柱模型单站雷达散射截面,半径r=16.25 mm,顶部圆锥高度h1=48.5 mm,底部圆柱部分高度h2=102 mm,入射波频率为0.3 THz。从图9和图10给出的结果可以看出,随着均方根高度的增加,目标表面的粗糙度变大,相对于0.3 THz的入射波其波长仅有1 mm,目标更加粗糙,粗糙度对目标的散射结果影响增大。当粗糙度较小时,RCS曲线可以看作是在光滑模型散射结果叠加小起伏震荡;粗糙度增大以后由表面粗糙度引起的RCS起伏甚至在某些角度可以改变光滑模型的散射曲线。
4. 结论
本文参考分形粗糙面模拟随机环境的方法来建立具有分形粗糙表面目标,采用基于基尔霍夫近似的物理光学方法研究了涂覆目标的太赫兹散射特性。分析了不同的入射波频率以及不同涂层厚度的分形粗糙表面模型在太赫兹波段的散射特性。相对于微波频段波长远大于目标表面微米量级的粗糙度,粗糙度的影响可以不考虑,而在太赫兹波段,波长与粗糙度处于等量级,必须考虑到粗糙度对于目标散射结果的影响。目标表面有涂覆介质材料时,目标的雷达散射截面小于导体情况下的结果,且在一定的范围内涂覆层越厚,目标雷达散射截面吸收越明显。
-
表 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 表 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 ℃ -
[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. 杨帆,何嘉岳,杨瑶佳,金一飞,许慎恒,李懋坤. 界面电磁学的理论与应用. 微波学报. 2023(05): 52-61 . 百度学术
2. 王禄炀,兰峰,宋天阳,何贵举,潘一博,张雅鑫,陈智,杨梓强. 多功能动态波束调控的太赫兹编码超表面. 无线电通信技术. 2022(02): 247-252 . 百度学术
3. 周嵩林,唐隽文,刘罗颢,吴优,刘长昊,金一飞,杨帆,许慎恒,李懋坤. 基于电磁表面的阵列天线及应用概述. 通信学报. 2022(12): 13-23 . 百度学术
4. 李国英,嵇成高,于刚刚,关浩. 相控雷达成像测井仪器中收发天线系统设计. 测井技术. 2022(06): 696-700+706 . 百度学术
其他类型引用(1)
-