Coastal wind field retrieval from polarimetric synthetic aperture radar

Coastal winds are strongly influenced by topology and discontinuity between land and sea surfaces. Wind assessment from remote sensing in such a complex area remains a challenge. Space-borne scatterometer does not provide any information about the coastal wind field, as the coarse spatial resolution hampers the radar backscattering. Synthetic aperture radar （SAR） with a high spatial resolution and all-weather observation abilities has become one of the most important tools for ocean wind retrieval, especially in the coastal area. Conventional methods of wind field retrieval from SAR, however, require wind direction as initial information, such as the wind direction from numerical weather prediction models （NWP）, which may not match the time of SAR image acquiring. Fortunately, the polarimetric observations of SAR enable independent wind retrieval from SAR images alone. In order to accurately measure coastal wind fields, this paper proposes a new method of using co-polarization backscattering coefficients from polarimetric SAR observations up to polarimetric correlation backscattering coefficients, which are acquired from the conjugate product of co-polarization backscatter and cross-polarization backscatter. Co-polarization backscattering coefficients and polarimetric correlation backscattering coefficients are obtained form Radarsat-2 single-look complex （SLC） data.The maximum likelihood estimation is used to gain the initial results followed by the coarse spatial filtering and fine spatial filtering. Wind direction accuracy of the final inversion results is 10.67 with a wind speed accuracy of 0.32 m/s. Unlike previous methods, the methods described in this article utilize the SAR data itself to obtain the wind vectors and do not need external wind directional information. High spatial resolution and high accuracy are the most important features of the method described herein since the use of full polarimetric observations contains more information about the space measured.This article is a useful addition to the work of independent SAR wind retrieval. The experimental results herein show that it is feasible to employ the co-polarimetric backscattering coefficients and the polarimetric correlation backscattering coefficients for coastal wind field retrieval.

The polarimetric features of oil spills in full polarimetric synthetic aperture radar images

Compared with single-polarized synthetic aperture radar （SAR） images, full polarimetric SAIl images contain not only geometrical and backward scattering characteristics, but also the polarization features of the scattering targets. Therefore, the polarimetric SAR has more advantages for oil spill detection on the sea surface. As a crucial step in the oil spill detection, a feature extraction directly influences the accuracy of oil spill discrimination. The polarimetric features of sea oil spills, such as polarimetric entropy, average scatter angle, in the full polarimetric SAR images are analyzed firstly. And a new polarimetric parameter P which reflects the proportion between Bragg and specular scattering signals is proposed. In order to investigate the capability of the polarimetric features for observing an oil spill, systematic comparisons and analyses of the multipolarization features are provided on the basis of the full polarimetric SAR images acquired by SIR-C/X-SAR and Radarsat-2. The experiment results show that in C-band SAR images the oil spills can be detected more easily than in L-band SAR images under low to moderate wind speed conditions. Moreover, it also finds that the new polarimetric parameter is sensitive to the sea surface scattering mechanisms. And the experiment results demonstrate that the new polarimetric parameter and pedestal height perform better than other polarimetric parameters for the oil spill detection in the C-band SAR images.

Ground Surface Ruptures and Near-Fault,Large-Scale Displacements Caused by the Wenchuan Ms8.0 Earthquake Derived from Pixel Offset Tracking on Synthetic Aperture Radar Images

The 12 May 2008 Wenchuan Ms8.0 earthquake produced surface displacements along the causative fault, the Yingxiu-Beichuan Fault, which are up to several meters near the fault. Because of the large gradient, satellite synthetic aperture radar （SAR） interferometric data are strongly incoherent; the usual SAR interferometry method does not allow such displacements to be measured. In the present study, we employed another approach, the technique based on pixel offset tracking, to solve this problem. The used image data of six tracks are from the Advanced Land Observing Satellite, Phased Array type L-band Synthetic Aperture Radar （ALOS/PALSAR） dataset of Japan. The results show that the entire surface rupture belt is 238 km long, extending almost linearly in a direction of 42°north-east. It is offset left laterally by a north-west-striking fault at Xiaoyudong, and turns at Gaochuan, where the rupture belt shifts toward the south by 5 km, largely keeping the original trend. In terms of the features of the rupture traces, the rupture belt can be divided into five sections and three types. Among them, the Beichuan-Chaping and Hongkou-Yingxiu sections are relatively complex, with large widths and variable traces along the trend. The Pingtong-Nanba and Qingping-Jingtang sections appear uniform, characterized by straight traces and small widths. West of Yingxiu, the rupture traces are not clear. North of the rupture belt, surface displacements are 2.95 m on average, mostly 2-3.5 m, with 7-9 m the maximum near Beichuan. South of the rupture belt, the average displacement is 1.75 m, dominated by 1-2 m, with 3-4 m at a few sites. In the north, the displacements in the radar line of sight are of subsidence, and in the south, they are uplifted, in accordance with a right-slip motion that moves the northern wall of the fault to the east, and the southern wall to the west, respectively. Along the Guanxian-Jiangyou Fault, there is a uplift zone in the radar line of sight, which is 66 km long, 1.5-6 km wide, and has vertical displacements of approximately 2 m, but no observable rupture traces.

A Review on Applications of Imaging Synthetic Aperture Radar with a Special Focus on Cryospheric Studies

The cryosphere is the frozen part of the Earth’s system. Snow and ice are the main constituents of the cryosphere and may be found in different states, such as snow, freshwater ice, sea ice, perma-frost, and continental ice masses in the form of glaciers and ice sheets. The present review mainly deals with state-of-the-art applications of synthetic aperture radar (SAR) with a special emphasize on cryospheric information extraction. SAR is the most important active microwave remote sensing (RS) instrument for ice monitoring, which provides high-resolution images of the Earth’s surface. SAR is an ideal sensor in RS technology, which works in all-weather and day and night conditions to provide useful unprecedented information, especially in the cryospheric regions which are almost inaccessible areas on Earth. This paper addresses the technological evolution of SAR and its applications in studying the various components of the cryosphere. The arrival of SAR radically changed the capabilities of information extraction related to ice type, new ice formation, and ice thickness. SAR applications can be divided into two broad classes-polarimetric applications and interferometric applications. Polarimetric SAR has been effectively used for mapping calving fronts, crevasses, surface structures, sea ice, detection of icebergs, etc. The paper also summarizes both the operational and climate change research by using SAR for sea ice parameter detection. Digital elevation model (DEM) generation and glacier velocity mapping are the two most important applications used in cryosphere using SAR interferometry or interferometric SAR (InSAR). Space-borne InSAR techniques for measuring ice flow velocity and topography have developed rapidly over the last decade. InSAR is capable of measuring ice motion that has radically changed the science of glaciers and ice sheets. Measurement of temperate glacier velocities and surface characteristics by using airborne and space-borne interferometric satellite images have been the significant application in glaciology and cryospheric studies. Space-borne InSAR has contributed to major evolution in many research areas of glaciological study by measuring ice-stream flow velocity, improving understanding of ice-shelf processes, yielding velocity for flux-gate based mass-balance assessment, and mapping flow of mountain glaciers. The present review summarizes the salient development of SAR applications in cryosphere and glaciology.

Precise Three-Dimensional Deformation Retrieval in Large and Complex Deformation Areas via Integration of Offset-Based Unwrapping and Improved Multiple-Aperture SAR Interferometry:Application to the 2016 Kumamoto Earthquake

Conventional synthetic aperture radar(SAR)interferometry(InSAR)has been successfully used to precisely measure surface deformation in the line-of-sight(LOS)direction,while multiple-aperture SAR interferometry(MAI)has provided precise surface deformation in the along-track(AT)direction.Integration of the InSAR and MAI methods enables precise measurement of the two-dimensional(2D)deformation from an interferometric pair;recently,the integration of ascending and descending pairs has allowed the observation of precise three-dimensional(3D)deformation.Precise 3D deformation measurement has been applied to better understand geological events such as earthquakes and volcanic eruptions.The surface deformation related to the 2016 Kumamoto earthquake was large and complex near the fault line;hence,precise 3D deformation retrieval had not yet been attempted.The objectives of this study were to①perform a feasibility test of precise 3D deformation retrieval in large and complex deformation areas through the integration of offset-based unwrapped and improved multiple-aperture SAR interferograms and②observe the 3D deformation field related to the 2016 Kumamoto earthquake,even near the fault lines.Two ascending pairs and one descending the Advanced Land Observing Satellite-2(ALOS-2)Phased Array-type L-band Synthetic Aperture Radar-2(PALSAR-2)pair were used for the 3D deformation retrieval.Eleven in situ Global Positioning System(GPS)measurements were used to validate the 3D deformation measurement accuracy.The achieved accuracy was approximately 2.96,3.75,and 2.86 cm in the east,north,and up directions,respectively.The results show the feasibility of precise 3D deformation measured through the integration of the improved methods,even in a case of large and complex deformation.

ADAPTIVE MODEL-BASED SCATTERING DECOMPOSITION OF POLARIMETRIC SAR INTERFEROMETRY

In this paper,a new decomposition method is proposed to solve the problems that vegetation component is overestimated and is not sensitive to directional scattering features with traditional polarimetric Synthetic Aperture Radar(SAR)decomposition.It uses a Polarimetric Interferometric Similarity Parameter(PISP)calculated from Polarimetric SAR Interferometry(PolInSAR)datasets to the scattering decomposition.The PISP is proposed to reveal the geometric sensitivity of SAR interferometry.It is defined by three optimized mechanisms obtained from PolInSAR datasets,therefore,it not only relates to the coherent scattering mechanism closely,but also sufficiently uses the phase and amplitude information.The PISP of building is high,and forest’s PISP is low.The proposed method uses the PISP as a judge condition to select different vegetation model adaptively.The decomposition results show the proposed method can effectively solve the vegetation ingredients overestimation problem.In addition,it is sensitive to the directional scattering.

利用InSAR形变反演滑坡滑面运动特征——以汶川县布瓦村滑坡为例

当前滑坡滑面运动特征获取多采用接触式手段,费时费力、成本高昂,且仅能获取稀疏点位滑面运动状态。因此,本文提出一种非接触式滑坡滑面运动特征获取方法,其基本思路是:基于位错理论建立滑坡表面形变与滑面运动间的关系模型,同时引入滑面运动平滑约束条件,最终以InSAR等技术获取的地表形变为约束,反演获得滑面运动特征。本文选择四川省汶川县布瓦村滑坡作为研究对象,以时序InSAR技术获取的滑坡坡面升降轨形变数据为约束,基于发展的理论方法反演获得布瓦村滑坡滑面几何与运动学参数,并结合野外勘查、Lidar观测以及模型残差较好地验证了反演结果的可靠性,其中反演升降轨模型残差仅3.1 mm·a^(-1)和3.4 mm·a^(-1)。反演结果显示,布瓦村滑坡滑面位于坡表以下平均~15.0 m深度处,滑面整体以沿坡向运动为主,最大运动量级达~80 mm·a^(-1),同时在坡体西侧中下部发现微量的垂直坡向运动。本文提出方法可方便、快速的获取空间连续的滑坡滑面运动特征,为滑坡灾害预警和失稳风险评估提供更为直接可靠的科学参考。

基于极化SAR梯度和复Wishart分类器的舰船检测

舰船检测是极化SAR系统的重要应用之一。现有的舰船检测方法容易受到旁瓣泄露的干扰,使得舰船目标的形态难以提取,导致检测结果不符合真实情况。此外,在舰船过于密集、尺度不一致的情况下,相邻舰船由于旁瓣的影响有时会被认为是单个目标,从而造成漏检。针对这些问题,该文提出一种基于极化SAR梯度和复Wishart分类器的舰船检测方法。首先,将似然比检验(LRT)梯度引入对数比值梯度框架,使其适用于极化SAR数据;基于LRT梯度图进行恒虚警(CFAR)检测,提取舰船的边缘信息,消除伪影的同时抑制强旁瓣对舰船精细轮廓提取的影响。其次,利用复Wishart迭代分类器对舰船强散射部分进行检测,可排除大部分的杂波干扰且保持舰船形态细节。最后,将二者信息融合,从而可以保持舰船形态细节的同时克服旁瓣和伪信号的虚警。该文在3幅来自ALOS-2卫星的极化SAR图像上进行了对比实验,实验表明与其他方法相比,该文所提算法具有更少的虚警和漏检,且能够有效克服旁瓣泄露,保持舰船形态细节。

SBAS-InSAR技术融合CNN-LSTM模型的矿区开采沉陷监测与预测

针对传统矿区开采沉陷监测方法耗费人力财力和预测预警模型较少的问题,研究提出一种基于短基线集合成孔径雷达干涉测量(Small Baseline Subset-Interferometry Synthetic Aperture Radar,SBAS-InSAR)技术和卷积神经网络(Convolutional Neural Networks,CNN)与长短期记忆网络(Long Short-Term Memory,LSTM)相结合的矿区开采沉陷监测预测方法。首先,利用SBAS-InSAR技术对建新煤矿进行矿区开采沉陷监测,获取了该矿区的年平均沉降速率和累计沉降值。用GNSS监测数据与SBAS-InSAR结果进行对比验证,其拟合效果较好。其次,在此基础上利用CNN-LSTM模型预测后6期沉降数据,其结果与CNN和LSTM预测结果进行对比。研究显示,CNN-LSTM模型的平均绝对误差(S_(MAE))和均方根误差(S_(RMSE))比单一的CNN和LSTM分别至少降低了44.8%和40.6%,其决定系数均高于98%。最后,进一步预测前6期和中6期沉降数据,验证了CNN-LSTM预测模型在时间上的一致性。因此,SBAS-InSAR融合CNN-LSTM模型在类似矿山开采沉陷监测和预测中有较好的应用前景。

基于最优PS点获取方法的矿山工业广场沉降监测

为保障矿山地面和井下工作人员、工业广场周边居民以及运煤系统运行安全,实时监测地表变形情况,提出永久散射体(PS)点的优化获取方法,以提升PS-合成孔径雷达干涉测量(InSAR)技术在沉降监测中的适用性。首先,选择淮北市郊煤矿工业广场煤柱Ⅱ513工作面开采沉陷期间2020年10月31日—2022年8月22日共26景SAR卫星影像,运用PS-InSAR技术确定最优的PS点目标选取方法;然后,基于该方法获取该工作面回采过程中的开采损害保护区地面沉降速率、累计沉降量;最后,基于地面水准点实测数据验证PS-InSAR监测精度并分析研究区地表及建(构)筑物动态沉降情况。结果表明:相干系数和振幅离差指数双阈值法较适应于研究区沉降监测;研究区内可探测到的最大下沉速率为-26.5 mm/a,最大累计下沉值为-53.7 mm;同时,探测到研究区西北部存在岩溶塌陷地质灾害因素;利用水准数据验证精度,均方根误差(RMSE)仅为3.8 mm,决定系数达到0.91。

2023年积石山Ms6.2级地震InSAR同震形变探测与断层滑动分布反演

2023年12月18日,甘肃省临夏回族自治州积石山县发生了Ms6.2级地震,地震灾害损失显著高于同震级地震,查明其同震形变场与断层滑动分布特征,有助于揭示此次地震灾害损失严重的原因。利用Sentinel-1卫星升、降轨雷达影像,采用合成孔径雷达干涉技术获取了本次地震的同震形变场,进而基于Okada弹性位错模型,确定了本次地震的震源参数,并基于分布式滑动模型反演了本次地震断层面上的滑动分布。结果表明,积石山地震为逆冲型地震,升、降轨同震形变场沿视线向最大形变量分别为7.1 cm和7.8 cm,断层最大滑动量为0.31 m,主要集中在地下0~8 km。与2019年6月17日四川长宁Ms 6.0级地震对比发现,长宁地震地表最大形变量为8 cm,断层最大滑动量为0.38 m,均大于积石山地震,推测积石山地震灾害损失严重的主要原因为震区居民房屋相对集中、房屋抗震性能差和黄土区场地地震波放大效应等。

面向目视解译的全极化SAR船只精细化特征表征方法

随着卫星技术的发展,极化合成孔径雷达(PolSAR)数据的分辨率和数据质量得到大幅提升,为人造目标的精细化目视解译提供了良好的数据条件。目前主要采用多分量分解的方法,但是易造成像素错分问题,为此,该文结合Yamaguchi极化分解和极化熵提出了一种非固定阈值划分的方法用于实现全极化SAR图像船只结构精细化特征表征。Yamaguchi极化分解能够识别基本散射机制,其修正后的体散射模型更符合实测数据,可有效对人造目标进行表征。极化熵H在弱去极化状态下可以看成某一指定等效点的目标散射机制,能够有效突出船只主散射特征。因此,该文通过将Yamaguchi极化分解算法的非固定三分量与极化熵的低中高熵内嵌,将其分为非固定阈值的九分类成分,从而降低硬阈值处理在阈值边界处受噪声影响产生的类别随机性。并且将二次散射和单次散射均显著的区域称为混合散射(MSM),以更好匹配实验中船只典型结构的散射类型。在此基础上,利用广义相似性参数进一步缩短类内距离,采用改进后的GSP-Wishart分类器进行迭代聚类,旨在通过提高二次散射和混合散射机制以提高不同类型船只可区分度。最后,该文采用中国上海某港口的高分三号全极化SAR数据进行实验,为了验证每艘船只特征表征正确性,通过船舶自动识别系统(AIS)收集并筛选了该港口船只信息及光学数据,并与极化SAR数据中每艘船只进行匹配。实验结果表明该方法可有效区分散货船、集装箱船和油轮3种类型船只。

基于时序InSAR技术的中贵天然气管道天水市段沿线滑坡隐患识别与形变分析

中贵天然气管道天水市段沿线地处黄土高原,地貌类型多样、构造活跃、复杂地质条件下发育多种具有隐蔽性和潜伏性的地质灾害,其中滑坡灾害对中贵天然气管道的安全危害极大,因此,对中贵天然气管道天水段沿线滑坡隐患进行有效识别与分析具有重要意义。采用Sentinel-1A升降轨卫星数据,基于时序合成孔径雷达干涉测量(synthetic aperture radar interferometry,InSAR)对中贵天然气天水段沿线滑坡隐患进行了解译识别、现场复核、发育特征与典型滑坡形变分析。研究区共识别17处滑坡隐患点,现场复核最终确定13处,其余4处为人类工程活动区,其中常沟村和磨峪沟村滑坡形变受降雨影响。统计发现研究区管道沿线滑坡多发生在坡度40°~50°,高差400~600 m,东北坡向和距管道100~150 m范围内,且岩土体强度较低,地层以第四系全新统冲洪积层(Q_(4)^(a1+p1))为主。结果表明:联合升降轨的时序InSAR技术可以有效识别管道沿线滑坡隐患,为油气管线的安全运营及今后油气管道选线提供重要的科学依据和参考。

基于RIME-CNN-SVR模型的麦田土壤水分反演

土壤水分监测对于农业生产和作物产量预估具有重要意义。近年来深度学习技术在土壤水分反演领域得到广泛应用,但大多侧重于模型结构增强和优化,对模型超参数优化研究探索不足。该研究提出了一种基于霜冰优化算法(rime optimization algorithm,RIME)的卷积神经网络(convolutional neural network,CNN)超参数优化模型,结合极化分解技术来校正植被对土壤水分反演精度的影响,以提高冬小麦农田土壤水分反演性能。首先利用RIME优化CNN超参数以构建RIME-CNN模型,然后使用RIME-CNN模型对特征参数进行自适应提取和挖掘,之后对这些特征参数进行正则化处理并输入到支持向量回归(support vector regression,SVR)模型,构建RIME-CNN-SVR模型进行土壤水分估算。为验证所建RIME-CNN-SVR模型的有效性,利用合成孔径雷达(synthetic aperture radar,SAR)数据结合光学遥感数据,在河南省开封市冬小麦农田区进行试验验证和精度分析。结果表明,该方法在不增加模型结构复杂性和可学习参数的前提下,显著提升了模型的预测性能,决定系数可达0.72,均方根误差为2.78%,平均绝对误差为2.20%。该研究可为农业生产提供一个更为准确、可靠的土壤水分监测手段。

基于多干扰机协同的极化SAR干扰方法

针对单极化干扰机生成的干扰会被极化合成孔径雷达(SAR)鉴别和对消的问题,本文提出了一种基于多干扰机协同的极化SAR干扰方法。该方法在不改变现有干扰机单极化发射结构的基础上,通过时序协同不同位置的多个单极化干扰机,实现对极化SAR的干扰。具体来说,根据极化SAR抗压制和欺骗干扰处理方式的不同,该方法利用不同位置的多个单极化干扰机同时发射,使极化SAR无法准确估计出干扰机发射天线极化比,从而无法对消压制干扰,实现压制干扰的目的;使用不同位置的多个单极化干扰机交替发射,有效破坏相邻两个接收脉冲之间的相关性,使极化SAR无法实现假目标对消,从而实现对极化SAR欺骗干扰的目的,并且多部单极化干扰机交替发射也能实现方位向间歇采样的干扰效果,增加欺骗干扰生成虚假目标的复杂性。由于实际应用中不需要改变现有单极化干扰机发射结构,该文所提方法为实现对极化SAR干扰带来了便利。并且该方法在生成压制和欺骗干扰时均使用模板调制干扰方式,通过简单的模板修改就能实现可控的压制和欺骗干扰效果。文中首先分析了极化SAR抗压制和欺骗干扰的方法,然后分别给出了多干扰机生成压制和欺骗干扰的时序协同方式,最后以模板调制干扰为基础具体分析了多干扰机协同生成干扰的方法,并分析了多干扰机的布局和干扰机数量选取准则。理论推导和仿真实验证明了本文方法能够实现对极化SAR的干扰。

基于融合距离的极化SAR图像非局部均值滤波

在极化合成孔径雷达(synthetic aperture radar,SAR)图像降噪领域,常见的非局部均值滤波仅依靠像素间的统计距离进行相似性度量,忽略了像素点的空间信息。本文结合极化SAR数据统计特性和图像空间特征作为像素间的相似性度量,提出了一种利用融合距离来计算相邻窗口权重的方法——基于融合距离的非局部均值滤波器。融合距离的引入使得滤波器能够更全面的评估像素间的相似性,从而得到更合适的像素权重。此外,本方法还引进变异系数对邻域窗口的权重进行评估,通过该参数可以控制滤波的程度。在多幅极化SAR图像上的实验结果表明,所提出的滤波器能够在有效抑制斑点噪声的同时保留较为完整的图像边缘信息和极化散射特性。

非邻域相干系数估计的PolInSAR树高反演算法

为了提高全极化干涉合成孔径雷达树高反演精度,提出了一种基于非邻域窗口相干系数估计的Pol In SAR树木高度反演算法。该算法首先求取4对全极化干涉SAR图像的最优干涉相位图,再以局部窗口内像素的幅度和相位联合条件概率密度为准则选取与当前像素满足独立同分布的样本,进而获得相干系数的准确估计值,最后采用三段法进行树木高度反演。对仿真数据的处理结果表明,所提算法能够提高观测场景中树木高度的反演精度,并且能够很好地保持图像中的边界细节信息。

基于多源遥感图像多级协同融合的舰船识别算法

针对极化合成孔径雷达(polarimetric synthetic aperture radar,PolSAR)图像存在斑点噪声严重、可视性差、直接影响目标识别精度的问题,提出一种基于多源遥感图像多级协同融合的舰船识别算法。通过采用多级协同融合方式,丰富图像的特征量,提高舰船识别精度。所提方法首先进行多源遥感数据的像素级融合,然后在上一步基础上进行特征级融合,最终得到新的目标特征。所提方法充分发挥了不同频段的PolSAR与多光谱图像的信息互补优势,不仅保留了多频段PolSAR对目标的极化散射特征,也保留了多光谱数据的空-谱信息。所提方法在可视性与检测精度上表现都较为出色,与传统的单一遥感数据相比,识别精度至少提高了5.12%。

基于GAMMA软件二次开发的地表形变监测系统

InSAR技术具有大范围、高精度、全天时、全天候的优点,已被广泛地应用于地表形变监测与高程获取。GAMMA软件具有支持数据源丰富、处理功能全面、结果精度可靠和效率高等优势,是当前国际上常用的InSAR数据处理软件。但该软件采用命令行交互,操作命令繁杂、对专业要求程度高,不易于非专业用户使用。针对以上问题,本文利用界面友好、开源的QGIS对GAMMA进行二次开发。结合GAMMA提供的py_gamma库和QGIS提供的py_QGIS API接口,采用模块化设计、算法并行优化等技术,开发了一套界面化InSAR形变监测系统,并以西藏沃达地区106景哨兵1号数据进行处理试验,从软件性能和处理结果方面证明了基于QGIS和GAMMA二次开发的可行性和可靠性。相关内容可为命令行用户界面软件转图形化用户界面提供参考。

基于模糊聚类的PolInSAR数据非监督分类方法

提出了一种结合Freeman分解和模糊聚类的极化干涉合成孔径雷达(polarimetric interferometric synthetic aperture radar,PolInSAR)数据非监督分类方法。针对基于Freeman分解的极化SAR图像分类方法中提取的3种散射机理:表面散射、体散射、偶次散射占主导的区域之间存在模糊的缺点,利用PolInSAR处理中的最优干涉相干系数引入的参数——最优相干熵HInt和最优相干各项异性度AInt,将每种散射机理主导区域划分为单个散射机制或多个散射机制共同作用的区域。并将模糊理论引入到HInt/AInt平面的区域边界划分,得到初始分割图像。对初始分割图像进行合并,模糊聚类等操作,得到最终分类结果。采用ESAR Oberpfaffenhofen地区PolInSAR数据实验的结果验证了本文方法的有效性。