基于极化相干矩阵的河流水质污染监测初探

利用遥感手段监测水质污染具有监测范围广、实时性强的优点。本研究选择中国南方地区水质变化明显的河流区域作为研究对象,利用C波段星载极化合成孔径雷达(极化SAR)卫星RADARSAT-2提供的数据,对基于极化SAR的河流水质污染监测技术进行了初步研究。首先介绍了电磁波的极化现象以及极化SAR的基本原理;其次结合一次同步测量实验,提取河流区域,并对极化SAR数据与水质常见监测指标的监测数据进行对比分析,通过对10个采样点18组数据的分析发现,极化相干矩阵中的T 22元素与部分水质指标(如五日生化需氧量BOD₅等)具有较强的相关性,从而反映水质的污染状况;并通过实测数据和最小二乘法,拟合得到了利用[T]矩阵元素反演BOD₅的经验公式,拟合系数达到0.82。最后通过对地表散射模型和菲涅尔系数的分析,从理论上探讨了极化相干矩阵中部分元素与水体物理性质存在相关性的原因。初步理论分析和实验数据表明,T 22元素能够反映水体的污染状况。     

DEM Generation Combining SAR Polarimetry and Shape-From-Shading Techniques

Estimation of the polarization orientation angle shifts induced by terrain azimuth slope variations is a recently developed application in radar polarimetry. In general, without any prior knowledge on the terrain, two polarimetric SAR (POLSAR) flight passes are required to derive terrain slopes in perpendicular directions for digital elevation model (DEM) generation. Moreover, we note that SAR intensity is a strong indicator of the range component of the terrain slopes. In this letter, we developed a method for DEM generation requiring only one POLSAR flight pass, by combining orientation angle estimation and a shape-from-shading technique. In particular, when limited POLSAR data are available, this POLSAR technique provides an alternative way for DEM generation. National Aeronautics and Space Administration Jet Propulsion Laboratory (NASA/JPL) AIRSAR L-band POLSAR data over Camp Roberts, California, is used to demonstrate the results of the method proposed in this letter, and a DEM derived from simultaneously measured C-band interferometric SAR from NASA/JPL topographic SAR instrument is selected as the comparative ground truth to validate the effectiveness of this single POLSAR method. Analyses and discussions are also included in this letter.     

Coherence optimization using the polarization state conformation in PolInSAR

A new coherence optimization algorithm is proposed for polarimetric synthetic aperture radar (SAR) interferometry applications by using the polarization state conformation algorithm based on the polarimetric basis transformation along with the polarization signatures. Through application of this algorithm, the resemblance between the scattering mechanisms of the same target in the repeat-pass polarimetric SAR (POLSAR) images is maximized. Then, coherence maps between the repeat-pass POLSAR images, before and after application of the algorithm, are generated. The coherences obtained by this method represent the best coherences or optimized coherences between the POLSAR images. The effects predicted by the theory are confirmed by the POLSAR data acquired by the Jet Propulsion Laboratory Spaceborne Imaging Radar mission.     

Polarization state conformation and its application to change detection in polarimetric SAR data

For polarimetric SAR (POLSAR) images, it is ideal that scattering geometries of the same target should display resemblance between multidate images, which are used in change detection applications, since the scattering mechanisms may change due to the data acquisition geometry. However, sometimes it is difficult to achieve these conditions. An attempt is made to maximize the resemblance between the scattering geometries in multidate images for a specific target. An algorithm is developed based on the polarimetric basis transformation along with the polarization signatures. As a result, the resemblance between the scattering mechanisms of the same target in both images is maximized. The effects predicted by the theory are confirmed by the change detection analysis of POLSAR data acquired by the Jet Propulsion Laboratory Spaceborne Imaging Radar-C mission.     

Polarimetric SAR Data Correction and Terrain Topography Measurement Based on the Radar Target Orientation Angle

Topographic variations caused by the range and the azimuth terrain slopes induce polarization orientation changes which cause the polarization to rotate about the line of sight. The existence of these variations reduce the accuracy measurement of geophysical parameters from polarimetric synthetic aperture radar (PolSAR) images. For this reason most inversion studies are best done in area of flat earth. In area which has significant terrain variations require compensation for topography. In real situations, terrain slopes rotate the polarization basis of the polarimetric scattering matrices by an orientation angle shift, and induce significant cross-polarization power. In this paper, two methods have been investigated using the polarimetric orientation angle (PAO): the first one involves the rotation of the polarimetric scattering and coherency matrices to achieve maximum azimuthally asymmetry for polarimetric data compensation to ensure accurate estimation of geophysical parameters in rugged terrain areas. The second approach has been developed which measures azimuth and range terrain slopes that are related to shifts in polarization orientation angle. Terrain elevation maps relative to a plane parallel to the radar line of sight can then be generated by integrating these slopes requiring only one PolSAR flight pass by combing orientation angle estimation and a shape-from-shading technique (SFS) which is mostly used by the computer vision community. Experimental results with C-band polarimetric RADARSAT2 data are used evaluate the data compensation algorithm and DEM generation.     

Characterization of forests and deforestation in Cambodia using ALOS/PALSAR observation

In this study, we have demonstrated the capability of full polarimetric ALOS/Phased Array L-band Synthetic Aperture Radar data for the characterization of the forests and deforestation in Cambodia, to support climate change mitigation policies of Reducing Emission from Deforestation and Forest Degradation (REDD). We have observed mean backscattering coefficient (σ°), entropy (H), alpha angle (α), anisotropy (A), pedestal height (PH), Radar Vegetation Index (RVI) and Freeman–Durden three-component decomposition parameters. The observations show that the forest types and deforested area are showing variable polarimetric and backscattering properties because of the structural difference. Evergreen forest is characterized by a high value of σ° HV (?12.96 dB) as compared with the deforested area (σ° HV=?22.2 dB). The value of polarimetric parameters such as entropy (0.93), RVI (0.91), PH (0.41) and Freeman–Durden volume scattering (0.43) is high for evergreen forest, whereas deforested area is characterized by the low values of entropy (0.36) and RVI (0.17). Based on these parameters, it is found that σ° HV, entropy, RVI and PH provide best results among other parameters.     

基于Freeman散射熵和各向异性度的极化SAR影像分类算法研究

极化SAR影像中阴影、水体和裸露的耕地3种地物类型有非常相似的极化散射特性,常规基于非相干分解的分类方法难以将其有效地区分。对此,本文引入基于Freeman分解的散射熵Hf和各向异性度Af两个特征参数,并将其用于极化SAR影像分类。首先利用Hf和Af参数将阴影和水体提取出来,然后将其他地物按散射机制分为3大类,并对每一类再次利用Hf和Af参数进行细分,最后通过基于Wishart分布的聚类和迭代分类,得到最终的分类结果。通过利用Radarsat-2在河南登封获取的全极化SAR数据进行试验,表明该算法执行效率高,能够有效地区分阴影、水体和裸露的耕地,并且对其他地物类型也有很好的分类效果。     

Evaluation of modified four-component scattering power decomposition method over highly rugged glaciated terrain

In recent years, there has been increased utilization of fully polarimetric synthetic aperture radar (POLSAR) data to study glaciated terrain features for glaciological and climate change modelling. This article is concerned with more accurate results and appropriate analysis of POLSAR data over a highly rugged glaciated area in Himalayan region. For this purpose, the modified Yamaguchi four-component scattering power decomposition (4-CSPD) method with a rotation concept of 3 × 3 coherency matrix [T] about line of sight is evaluated. It has been found that the modified Yamaguchi 4-CSPD method significantly improved the decomposition results as compared with the original 4-CSPD by minimizing the cross-polarized Horizontal-Vertical (HV) components. This modified 4-CSPD leads to enhancement in the double bounce scattering and surface scattering components and also avoids the overestimation problem in the volume scattering component as compared with the original 4-CSPD from the sloped terrain. The significant reductions of the negative power occurrence in the surface scattering (3.9%) and the double bounce scattering (19.7%) components have also been noticed as compared with the original 4-CSPD method over the glaciated area in this part of the Indian Himalaya.     

采用螺旋散射的SAR极化定标参考地物提取方法

极化合成孔径雷达观测系统获取的影像需要经过极化定标处理才能进行定量的分析与应用。当前的极化定标方案普遍采用分布式地物解算串扰和交叉极化通道不平衡误差,需要在定标前选取满足一定散射特征的分布式地物作为定标参考样本。利用螺旋散射的极化特征与分布特点,提出了一种新的极化定标参考地物自动提取方法。该方法根据极化目标分解构建螺旋散射比率特征,并采用自适应阈值分割方法自动提取地物样本。采用C波段机载极化雷达影像进行实验,结果表明,所提方法能够保持影像极化定标的准确性,并能提高极化定标的精度。     

基于四分量散射模型的多极化SAR图像分类

基于四分量散射模型提出了一种多极化SAR(synthetic aperture radar)图像非监督分类算法。与Freeman三分量散射模型不同,四分量散射模型在Freeman三分量的基础上增加了螺旋散射分量(helix),该分量反映了复杂地貌和不规则城市建筑的散射机理,可以用来处理复杂的场景图像。算法强调了初始分类的重要性,在初始分类中考虑了混合散射机制像素的存在,从而提高了分类结果的精确度。聚类过程中,采用由四个散射分量组成的特征向量进行迭代聚类。为了实现算法的完全非监督,利用特征向量给出了一种新的聚类终止准则。NASA/JPL实验室AIRSAR全极化数据分类实验结果表明,该算法具有较好的分类效果,并获得了较高的分类精度。     

简缩极化特征值分析的溢油检测

为有效利用简缩极化SAR进行海洋溢油检测,本文基于简缩极化特征值分析,提出了3个用于简缩极化溢油检测的参数,引入了基于简缩极化特征值分解的简缩极化熵Hc(Compact Polarization Entropy)、简缩极化比参数PFc(Compact Polarization Fraction)、简缩极化基准高度PHc(Compact Polarization Pedestal Height)特征进行海洋溢油检测。海表的散射类型主要为低熵散射(小粗糙面发生的Bragg散射),为弱去极化、弱散射过程随机性状态,由于溢油会阻尼海水的Bragg散射,使其熵值变高、呈去极化、强散射过程随机性状态,故简缩极化熵、简缩极化比参数和简缩极化基准高度可以用来检测海洋溢油。本文采用C波段的Radarsat-2、SIR-C/X-SAR数据进行了实验,结果表明:简缩极化熵、简缩极化比参数和简缩极化基准高度能够有效抑制疑似溢油,使海水与疑似溢油差异变小;突出溢油区域,使海水与溢油的可区分性变大。     

极化干涉SAR数据地表土地类型分类

基于新疆和田地区1994年10月9日和10日SIR-C-L波段全极化雷达数据。首先对极化干涉测量的基本原理和数据处理流程进行了详细的阐述，接着，用Cloude相干最优算法得到了与3种地物散射机制相对应的3个最优相干图。并且就地物相干性对极化的强烈依赖和3种散射机制中地物的最优相干特性进行了分析，具有最高相干值的相位图在提取DEM方面较有利，具有最低相干值的相干图在地物识别方面较有利。最后，在对最优相干系数。后向散射系数和熵进行数据相关性分析基础上，利用得到的最优相干系数，熵和后向散射系数数据进行了土地类型的识别和分类，得到了很好的效果。     

An improved building detection approach using L-band POLSAR two-dimensional time-frequency decomposition over oriented built-up areas

The objective of this study is to efficiently extract detailed information about various man-made targets in oriented built-up areas using polarimetric synthetic aperture radar (POLSAR) images. This paper develops an improved approach for building detection by utilizing Two-Dimensional Time-Frequency (2-D TF) decomposition. This method performs outstandingly in distinguishing between man-made and natural targets based on the isotropic behaviors, frequency-sensitive responses, and scattering mechanisms of objects. The proposed method can preserve the spatial resolution and exploit the advantages of TF decomposition; specifically, the exact outlines of buildings can be effectively located, and more types of features (e.g., flat roofs, roads, and walls that are oblique to the radar illumination) can be distinguished from forests in complex built-up areas by 2-D TF decomposition. The coarser-resolution subaperture images that are produced in the azimuth direction, which correspond to different looking angles, are beneficial for detecting man-made structures with main scattering centers oriented at oblique angles with respect to the radar illumination. In the range direction, the obtained subaperture images, which correspond to various observation frequencies, can be helpful in distinguishing flat roofs and roads from forests. This method was successfully implemented to analyze both NASA/JPL L-band AIRSAR and L-band EMISAR data sets. The building detection results of the proposed method exhibit a significant improvement over those of other methods and reach an overall accuracy over 80%, with approximately 20% higher than the accuracies of K-means clustering and the entropy/alpha-Wishart classifier and approximately 10% higher than the accuracy of the support vector machine method. Moreover, building details can be precisely detected, obliquely oriented buildings can be identified, and the distinction between buildings and forests is significantly improved, as both visually and statistically indicated. This method is highly adaptable and has substantial application value.     

Radar Polarization Orientation Shifts in Built-Up Areas

Polarization orientation angle shifts can be seen not only in rugged terrain areas but also in urban areas. The latter is explained by backscatter from a wall of a building or house, which is equivalent to a tilted ground-surface patch. From the scattering model of built-up areas, the polarization orientation angle shift in the built-up areas is given as , where is the wall or street orientation angle, and is the radar incidence angle. Japan Aerospace Exploration Agency Pi-SAR L-band polarimetric data of Gifu, Japan, show a good agreement with the theory. The phase difference between VH and HH polarizations is used to demonstrate the contribution of double-bounce scattering ground-wall and wall-ground over a wide range of wall orientation angles.     

On the degree of polarization for SAR sea oil slick observation

A polarimetric model to relate the degree of polarization, DoP, to the sea surface scattering with and without oil slicks, under low-to-moderate wind conditions, is proposed. DoP, measured directly from the Mueller scattering matrix, is shown to be a reliable measure of the departure from Bragg scattering; a phenomenon that, under low-to-moderate wind conditions, occurs when an oil slick is present. Following this theoretical rationale, a simple filter is developed to observe oil slicks in quad-polarimetric full-resolution Synthetic Aperture Radar (SAR) data. Experiments, undertaken on a meaningful set of quad-polarization Single Look Complex (SLC) C-band RADARSAT-2 SAR data, where both well-known oil slicks and a weak-damping look-alike are in place, demonstrate the soundness of the model and its effectiveness from an operational viewpoint.     

北京地区冬小麦冠层光谱数据与叶面积指数统计关系研究

以北京地区冬小麦为研究对象,利用TM传感器的光谱响应函数处理地面测量获得的冬小麦冠层光谱数据,得到对应于TM传感 器红光波段和近红外波段的反射率,进而计算出冬小麦冠层的归一化植被指数NDVI。建立了LAI与NDVI之间的不同经验关 系模型,对实验结果进行分析后得出,LAI与NDVI之间具有高度的指数相关性。     

Development of an inversion algorithm for dry snow density estimation and its application with ENVISAT-ASAR dual co-polarization data

Radar remote sensing has great potential to determine the extent and properties of snow cover. Availability of space-borne sensor dual-polarization C-band data of environmental satellite- (ENVISAT-) advanced synthetic aperture radar (ASAR) can enhance the accuracy in measurement of snow physical parameters as compared with single polarization data measurement. This study shows the capability of C-band synthetic aperture radar (SAR) data for estimating dry snow density over snow covered rugged terrain in Himalayan region. The snow density is an important parameter for the snow hydrology and avalanche forecasting related studies. An algorithm has been developed for estimating snow density, based on snow volume scattering and snow-ground scattering components. The radar backscattering coefficients of both horizontal–horizontal (hh) and vertical–vertical (vv) polarization and incidence angle are used as inputs in the algorithm to provide the snow dielectric constant which can be used to derive snow density using Looyenga's semi-empirical formula. Comparison was made between snow density estimated from algorithm using ENVISAT-ASAR hh and vv polarization data and the measured field value. The mean absolute error between estimated and measured snow density was found to be 0.024 g/cm³.     

The potential of linear discriminative Laplacian eigenmaps dimensionality reduction in polarimetric SAR classification for agricultural areas

In this paper, the linear discriminative Laplacian eigenmaps (LDLE) dimensionality reduction (DR) algorithm is introduced to C-band polarimetric synthetic aperture radar (PolSAR) agricultural classification. A collection of homogenous areas of the same crop class usually presents physical parameter variation, such as the biomass and soil moisture. Furthermore, the local incidence angle also impacts a lot on the same crop category when the vegetation layer is penetrable with C-band radar. We name this phenomenon as the “observed variation of the same category” (OVSC). The most common PolSAR features, e.g., the Freeman–Durden and Cloude–Pottier decompositions, show an inadequate performance with OVSC. In our research, more than 40 coherent and incoherent PolSAR decomposition models are stacked into the high-dimensionality feature cube to describe the various physical parameters. The LDLE algorithm is then performed on the observed feature cube, with the aim of simultaneously pushing the local samples of the same category closer to each other, as well as maximizing the distance between local samples of different categories in the learnt subspace. Finally, the classification result is obtained by nearest neighbor (NN) or Wishart classification in the reduced feature space. In the simulation experiment, eight crop blocks are picked to generate a test patch from the 1991 Airborne Synthetic Aperture Radar (AIRSAR) C-band fully polarimetric data from of Flevoland test site. Locality preserving projections (LPP) and principal component analysis (PCA) are then utilized to evaluate the DR results of the proposed method. The classification results show that LDLE can distinguish the influence of the physical parameters and achieve a 99% overall accuracy, which is better than LPP (97%), PCA (88%), NN (89%), and Wishart (88%). In the real data experiment, the Chinese Hailaer nationalized farm RadarSat2 PolSAR test set is used, and the classification accuracy is around 94%, which is again better than LPP (90%), PCA (88%), NN (89%), and Wishart (85%). Both experiments suggest that the LDLE algorithm is an effective way of relieving the OVSC phenomenon.     

Full year crop monitoring and separability assessment with fully-polarimetric L-band UAVSAR: A case study in the Sacramento Valley,California

Spatial and temporal information on plant and soil conditions is needed urgently for monitoring of crop productivity. Remote sensing has been considered as an effective means for crop growth monitoring due to its timely updating and complete coverage. In this paper, we explored the potential of L-band fully-polarimetric Uninhabited Aerial Vehicle Synthetic Aperture Radar (UAVSAR) data for crop monitoring and classification. The study site was located in the Sacramento Valley, in California where the cropping system is relatively diverse. Full season polarimetric signatures, as well as scattering mechanisms, for several crops, including almond, walnut, alfalfa, winter wheat, corn, sunflower, and tomato, were analyzed with linear polarizations (HH, HV, and VV) and polarimetric decomposition (Cloude–Pottier and Freeman–Durden) parameters, respectively. The separability amongst crop types was assessed across a full calendar year based on both linear polarizations and decomposition parameters. The unique structure-related polarimetric signature of each crop was provided by multitemporal UAVSAR data with a fine temporal resolution. Permanent tree crops (almond and walnut) and alfalfa demonstrated stable radar backscattering values across the growing season, whereas winter wheat and summer crops (corn, sunflower, and tomato) presented drastically different patterns, with rapid increase from the emergence stage to the peak biomass stage, followed by a significant decrease during the senescence stage. In general, the polarimetric signature was heterogeneous during June and October, while homogeneous during March-to-May and July-to-August. The scattering mechanisms depend heavily upon crop type and phenological stage. The primary scattering mechanism for tree crops was volume scattering (>40%), while surface scattering (>40%) dominated for alfalfa and winter wheat, although double-bounce scattering (>30%) was notable for alfalfa during March-to-September. Surface scattering was also dominant (>40%) for summer crops across the growing season except for sunflower and tomato during June and corn during July-to-October when volume scattering (>40%) was the primary scattering mechanism. Crops were better discriminated with decomposition parameters than with linear polarizations, and the greatest separability occurred during the peak biomass stage (July-August). All crop types were completely separable from the others when simultaneously using UAVSAR data spanning the whole growing season. The results demonstrate the feasibility of L-band SAR for crop monitoring and classification, without the need for optical data, and should serve as a guideline for future research.     

LiDAR不同强度校正法对樟子松叶面积指数估测的影响

机载激光雷达（LiDAR）强度数据在获取过程中受多种因素影响,各因素影响的有效量化及校正对机载LiDAR强度校正及应用具有重要意义。本文以雷达方程为基础,分别采用距离、入射角及距离和入射角对LiDAR点云强度进行校正,从中提取冠层总强度和强度比值两类参数,用于估测森林叶面积指数（LAI）,以期量化各影响因素强度校正对不同类型参数估测森林LAI的影响。结果表明：强度经距离校正能够提高森林LAI的估测精度,而强度经数字高程模型衍生入射角校正非但没能提高估测精度,反而降低了估测精度。强度经距离和入射角综合校正虽能提高森林LAI的估测精度,但结果却低于距离单独校正的结果。与此同时,对冠层总强度参数而言,强度校正前后森林LAI估测结果的差异较为明显,而对强度比值参数而言,强度校正前后森林LAI估测结果差异不大。综上可知,不同因素强度校正对森林LAI估测的影响不同,且影响程度与所用参数变量类型密切相关。因此,在未来强度应用研究中,应根据变量参数类型选择合适的校正方式,以避免不恰当校正造成的成本浪费及精度降低。