An extended multiple-component scattering model for PolSAR images

An extended multiple-component scattering model (MCSM) is proposed for polarimetric synthetic aperture radar (PolSAR) image decomposition. The MCSM is an extension of the three-component scattering model (TCSM), and it describes single-bounce, double-bounce, volume, helix and wire scattering as elementary scattering mechanisms in the analysis of PolSAR images. The proposed MCSM is demonstrated with German Aerospace Centre (DLR) experimental SAR (ESAR) L-band fully polarized images of the Oberpfaffenhofen Test Site Area (DE), Germany. Double-bounce, helix and wire scattering are found to be predominant in urban areas and the results confirm that the MCSM is effective for analysis of buildings in urban areas. A comparison of the TCSM and its extended models is also implemented.     

Mapping of flooded vegetation by means of polarimetric Sentinel-1 and ALOS-2/PALSAR-2 imagery

This article presents for the first time the combination of dual-polarimetric C-band Sentinel-1 synthetic aperture radar (SAR) data and quad-polarimetric L-band ALOS-2/PALSAR-2 imagery for mapping of flooded areas with a special focus on flooded vegetation. L-band SAR data is well suited for mapping of flooded vegetation, while C-band enables an accurate extraction open water areas. Polarimetric decomposition-based unsupervised Wishart classification is combined with object-based post-classification refinement and the integration of spatial contextual information and global auxiliary data. In eight different scenarios, focusing on single datasets or fusion of classification results of several ones, respectively, different polarimetric decomposition and classification principles, including the entropy/anisotropy/alpha and the Freeman–Durden–Wishart classification, were investigated. The helix scattering component of the Yamaguchi decomposition, derived from ALOS-2 imagery, showed high suitability to refine the Sentinel-1-based detection of flooded vegetation. A test site at the Evros River (Greek/Turkish border region) was chosen, which was affected by a flooding event that occurred in spring 2015. The validation was based on high spatial resolution optical WorldView-2 imagery acquired with short temporal delay to the SAR data.     

Unsupervised classification using polarimetric interferometric decomposition and Wishart classifier

ABSTRACT

In this paper, a decomposition scheme of the coherency matrix is presented to parse the information of polarimetric interferometric synthetic aperture radar (PolInSAR) images in detail. First, the decomposition method is improved by the polarimetric interferometric similarity parameter (PISP) to relief the overestimation occurred in the traditional four-component decomposition method. Second, after using the improved four-component decomposition results as the original inputs, the decomposition method is applied to retrieve scattering mechanisms or identify scatters, with the image separated into seven subsections. Finally, based on the modified decomposition results, the basic classification results are regarded as the feature training sets, and the Wishart classifier is then used as an optimized classification process. The applications of the decomposition and classification scheme are shown with typical representative L-band E-SAR images, which are used to show the robustness of the method, as well as with the first published airborne C-band PolInSAR data collected by the Institute of Electronics, Chinese Academy of Sciences, in November 2017. Experimental results demonstrate that the obtained decomposition and classification results are in good agreement with the actual physical scattering mechanisms.     

Decomposition of polarimetric synthetic aperture radar backscatter from upland and flooded forests

The goal of this research was to decompose polarimetric Synthetic Aperture Radar (SAR) imagery of upland and flooded forests into three backscatter types: single reflection, double reflection, and cross-polarized backscatter. We used a decomposition method that exploits the covariance matrix of backscatter terms. First we applied this method to SAR imagery of dihedral and trihedral corner reflectors positioned on a smooth, dry lake bed, and verified that it accurately isolated the different backscatter types. We then applied the method to decompose multi-frequency Jet Propulsion Laboratory (JPL) airborne SAR (AIRSAR) backscatter from upland and flooded forests to explain scattering components in SAR imagery from forested surfaces. For upland ponderosa pine forest in California, as SAR wavelength increased from C-band to P-band, scattering with an odd number of reflections decreased and scattering with an even number of reflections increased. There was no obvious trend with wavelength for cross-polarized scattering. For a bald cypress-tupelo floodplain forest in Georgia, scattering with an odd number of reflections dominated at C-band. Scattering power with an even number of reflections from the flooded forest was strong at L-band and strongest at P-band. Cross-polarized scattering may not be a major component of total backscatter at all three wavelengths. Various forest structural classes and land cover types were readily distinguishable in the imagery derived by the decomposition method. More importantly, the decomposition method provided a means of unraveling complex interactions between radar signals and vegetated surfaces in terms of scattering mechanisms from targets. The decomposed scattering components were additions to the traditional HH and V V backscatter. One cautionary note: the method was not well suited to targets with low backscatter and a low signal-to-noise ratio.     

Similarity-enhanced target detection algorithm using polarimetric SAR images

Target detection and analysis using polarimetric synthetic aperture radar (PolSAR) images are currently of great interest in synthetic aperture radar (SAR) applications. For a complex target, the scattering characteristics are determined by different independent sub-scatterers and their interaction; therefore, the scattering characteristics should be described by a statistical method due to randomness and depolarization. Furthermore, the inherent speckle in SAR data must be reduced by spatial averaging at the expense of loss of spatial resolution. The polarimetric similarity parameter (PSP) is an effective parameter to analyse target characteristics. In order to describe a complex distributed target, two new methods for calculating PSP are proposed, namely Stokes matrix-based PSP (S-PSP) and multiple PolSAR similarity parameter (MPSP). The characteristics of a target can be described and extracted on the basis of the polarimetric similarity, and then the similarity-enhanced target detection methods using S-PSP and MPSP are implemented and demonstrated with German Aerospace Centre (DLR) experimental SAR L-band multiple temporal PolSAR images of Oberpfaffenhofen test site (DE), Germany. The results confirmed that the proposed methods are effective for detection and analysis of buildings in urban areas.     

Model-based target decomposition with the π/4 mode compact polarimetry dataπ/4模式简缩极化基于模型的目标分解

In this paper, model-based (surface, dihedral and volume scattering) target decomposition technique is proposed to decompose the π/4 mode compact polarimetric radar data. A general relationship between fully polarimetric coherence matrix and the Stokes vector of the π/4 mode compact polarimetric data is first established. Based on the Stokes vector, a proposed algorithm to retrieve the power of three scattering mechanisms is given in details. We validate this algorithm with L-band AIRSAR, San Francisco Bay, and results of decomposition are discussed and assessed in detail by being compared with the quad-pol Freeman-Durden decomposition results. Finally, the π/4 mode decomposition is compared with the CTLR (circular transmitting and linear reeving) mode, and with the π/4 mode m ? δ targets decomposition. The comparison results are analyzed and discussed in detail.     

Radar observations of basaltic lava flows,Craters of the Moon,Idaho

Abstract

Radar images were assessed to determine the backscatter characteristics of basaltic lava flows of predominantly pahoehoe textures and the ability to detect fissure vents. The images were obtained from synthetic aperture side-looking airborne radar systems—X-band HH, X-band HV, L-band HH and L-band H V. Smooth, collapsed blisters of shelly pahoehoe have weak returns in all four radar images. These returns are identical to those from pahoehoe surfaces covered with smooth mantles of windblown sediments. Hummocky pahoehoe flows have strong backscatter in all four images, most likely due to the large range in surface roughness causing multiple scattering at both radar wavelengths. Aa lava flows show the greatest variation in backscatter intensities—strong XHH, weak XHV, strong LHH and very strong LHV returns. This variation is due to an increase in multiple scattering at the L-band scale. Although smooth and rough surface textures can be differentiated in the radar images, there are constraints in tracing textural changes back to a particular fissure vent, in part because near-vent flows do not have unique radar signatures. Eruptive fissures are detectable in the radar images by virtue of associated parallel spatter ramparts which have diagnostic, strong backscatter in the X-band images that are in contrast to the weak backscatter of the surrounding shelly pahoehoe lava. However, spatter ramparts are not delineated in the L-band images. The centimetre-scale relief of the agglutinate spatter may cause scattering of the X-band energy more than the L-band energy. Although the structures are several metres high, the look directions for both imaging systems are approximately parallel to the trend of the ramparts. The rampart walls do not serve as reflectors. Such findings emphasize the importance of look direction in the use of radar images to characterize terrains.     

Integration of radar and Landsat-derived foliage projected cover for woody regrowth mapping, Queensland, Australia

In Queensland, Australia, forest areas are discriminated from non-forest by applying a threshold (∼ 12%) to Landsat-derived Foliage Projected Cover (FPC) layers (equating to ∼ 20% canopy cover), which are produced routinely for the State. However, separation of woody regrowth following agricultural clearing cannot be undertaken with confidence, and is therefore not mapped routinely by State Agencies. Using fully polarimetric C-, L- and P-band NASA AIRSAR and Landsat FPC data for forests and agricultural land near Injune, central Queensland, we corroborate that woody regrowth dominated by Brigalow (Acacia harpophylla) cannot be discriminated using either FPC or indeed C-band data alone, because the rapid attainment of a canopy cover leads to similarities in both reflectance and backscatter with remnant forest. We also show that regrowth cannot be discriminated from non-forest areas using either L-band or P-band data alone. However, mapping can be achieved by thresholding and intersecting these layers, as regrowth is unique in supporting both a high FPC (> ∼ 12%) and C-band SAR backscatter (> ~ − 18 dB at HV polarisation) and low L-band and P-band SAR backscatter (e.g. < =∼ 14 dB at L-band HH polarisation). To provide a theoretical explanation, a wave scattering model based on that of Durden et al. [Durden, S.L., Van Zyl, J.J. & Zebker, H.A. (1989). Modelling and observation of radar polarization signature of forested areas. IEEE Trans. Geoscience and Remote Sensing, 27, 290-301.] was used to demonstrate that volume scattering from leaves and small branches in the upper canopy leads to increases in C-band backscattering (particularly HV polarisations) from regrowth, which increases proportionally with FPC. By contrast, low L-band and P-band backscatter occurs because of the lack of double bounce interactions at co-polarisations (particularly HH) and volume scattering at HV polarisation from the stems and branches, respectively, when their dimensions are smaller than the wavelength. Regrowth maps generated by applying simple thresholds to both FPC and AIRSAR L-band data showed a very close correspondence with those mapped using same-date 2.5 m Hymap data and an average 73.7% overlap with those mapped through time-series comparison of Landsat-derived land cover classifications. Regrowth mapped using Landsat-derived FPC from 1995 and JER-1 SAR data from 1994-1995 also corresponded with areas identified within the time-series classification and true colour stereo photographs for the same period. The integration of Landsat FPC and L-band SAR data is therefore expected to facilitate regrowth mapping across Queensland and other regions of Australia, particularly as Japan's Advanced Land Observing System (ALOS) Phase Arrayed L-band SAR (PALSAR), to be launched in 2006, will observe at both L-band HH and HV polarisations.     

基于Yamaguchi分解模型的全极化SAR图像分类

针对利用Yamaguchi分解模型的四个散射分量直接进行类别归属判断精度不高并且所分类别有限的问题,结合模糊C均值的理论,提出了一种基于Yamaguchi分解模型的全极化SAR分类算法,把四个散射分量组成一组归一化的特征矢量,进行FCM聚类分析。并且用日本机载L波段PiSAR数据验证了该算法具有较高的分类精度和较好的视觉效果。     

Change detection using distance-based algorithms between synthetic aperture radar polarimetric decompositions

The present study introduces distance based change detection (CD) algorithms in polarimetric synthetic aperture radar (PolSAR) data. PolSAR images, due to interactions between electromagnetic waves and target and because of the high spatial resolution, can be used to study changes in the Earth’s surface. The purpose of this paper is to use features extracted from the fully-polarimetry imaging radar that involved Yamaguchi four-component and H/A/α decomposition based on the distance between the vectors of features for CD. We first extract features from polarimetric decompositions of multi-looked covariance (or coherency) matrix data. We then use two well-known distance measures namely Canberra and Euclidean distances for measuring the similarity between the vectors of polarimetric decompositions at different times. Assessment of incorporated methods is performed using different criteria, such as overall accuracy, area under the receiver operating characteristic curve, and false alarms rate. The results of the experiments show that Canberra distance has better performance with high overall accuracy and low false alarm rate than Euclidean distance and other compared algorithms to detect changes.     

Radar response of periodic vegetation canopies

Vector radiative transfer theory is used to model the scattered intensity from a layer of randomly oriented particles over a periodic rough surface. To account for the periodic nature of row-structured vegetation, the number density of particles within the layer is assumed to be varying periodically in the horizontal direction. Using Fourier series expansions and orthogonality properties, the radiative transfer equation is solved for the transformation matrix relating the incident and scattered intensities, from which the backscattering coefficient of the layer can be computed for any incidence direction and polarization configuration. The experimental component of this investigation consisted of radar observations at 1–5,4–75, and 9–5 GHz made by a truck-mounted system for a field of corn under three conditions: (a) full, which means that the corn plants were in their natural state, (b) defoliated, which was accomplished by stripping off the leaves and removing them, thereby leaving behind only bare vertical stalks, and (c) bare soil, which corresponds to the soil surface after having removed the stalks. The soil surface is modelled as a composite consisting of a deterministic periodic component and a random roughness component. A two-scale polarimetric scattering model is formulated and used to compare with the experimental observations. Excellent agreement between theory and measurements is realized as a function of both incidence and azimuth angles at all three microwave frequencies. The canopy model was then applied to the corn canopy under the two other conditions: stalks alone and full canopy. The model results were compared with radar backscatter measurements made for each of three look directions, including perpendicular and parallel to the row direction and at 45° relative to the row direction. For the stalk canopy, it was observed that the quasi-periodic arrangement of the stalks within the row enhances the backscatter at L-band when looking perpendicular to the row direction, which is attributed to a coherent-scattering effect associated with the stalks. A heuristic approach is used to model the quasi-periodic structure of the stalks by deriving a coherency factor which multiplies the first-order radiative solution for randomly located stalks. A similar coherency factor was also introduced for the leaves of the full canopy. The modified model was found to provide good agreement with experimental observations at L-band for all polarizations and at all look directions.     

基于Krogager分解和SVM的极化SAR图像分类

目标分解包括基于Sinclair矩阵的相干目标分解和基于Mueller矩阵的部分相干目标分解,Krogager分解即属于相干目标分解,它可以将任一对称Sinclair矩阵分解为球散射体、二面角散射体和螺旋体3个分量,这是极化合成孔径雷达(Synthetic Aperture Radar,SAR)图像特征提取的有效途径。把3个分量的分解系数作为极化散射特征,由其组成样本向量,运用基于统计学习理论的支持向量机(Support Vector Machines,SVM)设计多类分类器,提出了一种极化SAR图像分类算法,并对实测极化SAR数据进行分类实验。结果表明,将Krogager分解和SVM分类器结合起来,对极化SAR图像进行分类是可行和有效的,并且选择不同的参数得到的分类结果差别很大,验证了参数选择在SVM分类器中的重要作用。     

The application of compact polarimetric decomposition algorithms to L-band PolSAR data in agricultural areas

ABSTRACT

In this paper, the applicability of the recently developed compact polarimetric decomposition and inversion algorithm to estimate soil moisture under low agricultural vegetation cover is investigated using simulated L-band compact polarimetric synthetic aperture radar (PolSAR) data. The surface scattering component is separated from the volume component of the vegetation through a model-based compact polarimetric decomposition (m-α) under the assumption of randomly orientated vegetation volume and reflection symmetry. The extracted surface scattering component is compared with two physics-based, low frequency surface scattering models such as extended Bragg (X-Bragg) and polarimetric two scale model (PTSM) in order to invert soil moisture for corresponding model- and data-derived surface scattering mechanism parameter α_s. In addition to the parameter α_s from m-α decomposition, the applicability of other scattering mechanism parameters, such as δ (relative phase) and χ (degree of circularity) from m-δ and m-χ decompositions are also investigated for their suitability to invert soil moisture. The algorithm is applied on a time series of simulated L-band compact polarimetric E-SAR data from the AgriSAR’2006 campaign over the Görmin test site in Northern Germany. The compact PolSAR-derived soil moisture is validated against in situ time-domain reflectometry (TDR) measurements. Including various growth stages of three different crop types, the estimated soil moisture values indicate an overall root mean square error (RMSE) of 9–12 and 9–15 vol.% using the X-Bragg model and the PTSM, respectively. The inversion rate for vegetation covered soils ranges from 5% to 40% including all phenological stages of the crops and different soil moisture conditions (range from 4 to 34 vol.%). The time series of soil moisture inversion results using compact polarimetry reveal that the developed algorithm is less sensitive to wet soils under growing agriculture crops due to less sensitivity of scattering mechanism parameters α_s and χ for ε_s > 20. Thus, further developments and investigations are needed to invert soil moisture for compact PolSAR data with high inversion rates and consistently less RMSE (<5 vol.%) over the various crop growing season.     

Radar modelling of coniferous forest using a tree growth model

The use of a tree growth model to provide statistical information about the microwave scattering components of boreal-type forests (in this case, Scots pine and Norwegian spruce), as an alternative to data obtained through intensive fieldwork, is described. The total backscatter from six test stands at C- and L-band frequency for three polarization combinations (HH, VV and HV) was predicted. Differences between measured C- and L-band data from a polarimetric airborne Synthetic Aperture Radar (EMISAR) and simulated backscatter values compare favourably with previous studies, with like- and cross-polarization differences generally less than 2.5 dB. Modelled backscatter values were consistently less than those observed. A likely explanation for such a discrepancy is the unrealistic manner in which the model incorporates the spatial distribution of tree needles.     

水面船只SAR探测的极化方式研究

利用SIR-C SAR的C和L波段全极化数据,分析水面船只的极化散射特性和船只与背景海面雷达后向散射的信噪比特性,研究水面船只SAR探测的最优极化方式。结果显示,二面角散射是水面船只SAR成像的主要机理。线性极化中,HV极化具有最大的船只与背景海面雷达后向散射信噪比。与线性极化相比,圆极化的雷达后向散射信噪比更优。C波段和L波段的水面船只的极化散射特性存在较大的差异,L波段的信噪比大于C波段的信噪比。水面船只的雷达后向散射特性表明,L波段的圆极化是水面船只探测的最优极化方式。     

Multi-chromatic analysis polarimetric interferometric synthetic aperture radar (MCA-PolInSAR) for urban classification

One of the problems of Synthetic Aperture Radar (SAR) polarimetric decomposition, is that oriented urban areas and vegetation signatures are decomposed into the same volume scattering mechanism. Such indetermination makes it difficult to distinguish vegetation from the oblique urban areas with respect to the radar illumination direction within the volume scattering mechanism. This event occurs because oriented targets exhibit similar polarimetric responses. This paper presents an improvement of the PolSAR decomposition scheme which permits the performing of more accurate classification. The method uses the information existing form the interference generated between two Doppler sub-aperture SAR images. This interferometric polarimetric SAR (PolInSAR) multi-chromatic analysis (MCA-PolInSAR) signal processing method permits the efficient separation of oriented buildings from vegetation yielding considerably improved results in which oriented urban areas are recognized, from volume scattering, as double-bounce objects. Results also show a considerable improvement in the robustness of classification and also in terms of definition and precision.     

On the relationship between the circumference of rubber trees and L-band waves

Despite substantial research conducted within the forestry domain, detailed assessments to monitor plantations and support their sustainable management have been understudied. This article attempts to fill this gap through coupling fully polarimetric L-band data and contemporary data mining methods for the estimation of tree circumference as: (1) a primary dataset for biomass accumulation studies; and, (2) critical information for operational management in rubber plantations. We used two rubber plantation sites in Subang (West Java) and Jember (East Java), Indonesia, to evaluate the capability of L-band radar data. Although polarimetric features derived from polarimetric decomposition theorems have been advocated by others, we show that backscatter coefficients, especially HV polarization, remain an important dataset for this research domain. Using Subang data to build the model, we found that modern machine learning methods do not always deliver the best performance. It appears that the data being ingested plays a significant role in obtaining a good model, hence careful selection of datasets from multiple forms of polarimetric SAR data needs to be further considered. The highest coefficient of determination (R² = 0.79) was achieved by Yamaguchi decomposition features with the aid of partial least squares regression. Nonetheless, we note that the R² gap was insignificant to the backscatter coefficient when random forests regression was used (R² = 0.78). Overall, only the backscatter coefficient dataset delivered fairly consistent results with any regression model, with the average R² being about 0.67. When tuning parameters were not assessed, random forests consistently outweighed support vector regressions in all forms of datasets. The latter generated a substantial increase in R² when a linear kernel was used instead of the popular radial basis function. The issue of transferability of the model is also addressed in this article. It appears that similarity of terrain characteristics substantially influences the model’s performance. Models developed in Subang, which has gentle slopes, seem valid only in plantations with similar terrain. Validation attempts in very flat terrain within two plantation sectors in Jember delivered a poor result, although they have similar elevations to the Subang site. In contrast, validation in a plantation sector with similar, gently sloping terrain achieved an R² of about 0.6 using some datasets.     

Investigation of the effect of the incidence angle on land cover classification using fully polarimetric SAR images

Incidence angle is one of the most important imaging parameters that affect polarimetric SAR (PolSAR) image classification. Several studies have examined the land cover classification capability of PolSAR images with different incidence angles. However, most of these studies provide limited physical insights into the mechanism how the variation of incidence angle affects PolSAR image classification. In the present study, land cover classification was conducted by using RADARSAT-2 Wide Fine Quad-Pol (FQ) images acquired at different incidence angles, namely, FQ8 (27.75°), FQ14 (34.20°), and FQ20 (39.95°). Land cover classification capability was examined for each single-incidence angle image and a multi-incidence angle image (i.e., the combination of single-incidence angle images). The multi-incidence angle image produced better classification results than any of the single-incidence angle images, and the different incidence angles exhibited different superiorities in land cover classification. The effect mechanisms of incidence angle variation on land cover classification were investigated by using the polarimetric decomposition theorem that decomposes radar backscatter into single-bounce scattering, double-bounce scattering and volume scattering. Impinging SAR easily penetrated crops to interact with the soil at a small incidence angle. Therefore, the difference in single-bounce scattering between trees and crops was evident in the FQ8 image, which was determined to be suitable for distinguishing between croplands and forests. The single-bounce scattering from bare lands increased with the decrease in incidence angles, whereas that from water changed slightly with the incidence angle variation. Consequently, the FQ8 image exhibited the largest difference in single-bounce scattering between bare lands and water and produced the fewest confusion between them among all the images. The single- and double-bounce scattering from urban areas and forests increased with the decrease in incidence angles. The increase in single- and double-bounce scattering from urban areas was more significant than that from forests because C-band SAR could not easily penetrate the crown layer of forests to interact with the trunks and ground. Therefore, the FQ8 image showed a slightly better performance than the other images in discriminating between urban areas and forests. Compared with other crops and trees, banana trees caused stronger single- and double-bounce scattering because of their large leaves. As a large incidence angle resulted in a long penetration path of radar waves in the crown layer of vegetation, the FQ20 image enhanced the single- and double-bounce scattering differences between banana trees and other vegetation. Thus, the FQ20 image outperformed the other images in identifying banana trees.     

Feature selection and weighted SVM classifier-based ship detection in PolSAR imagery

Target decomposition is an important method for ship detection in polarimetric synthetic aperture radar (SAR) imagery. Parameters such as the polarization entropy and alpha angle deduced from the coherency matrix eigenvalue decomposition capture the differences between the target and background from different views separately. However, under the conditions of a relatively high resolution and a rough sea, the contrast between ship and sea reduces in the aforementioned space. Based on the analyses of target decomposition theory and the target’s scattering mechanism, multi-polarization parameters can be used to characterize different scattering behaviours of the ship target and sea clutter. Moreover, each parameter has its own diverse significance in the practical detection problem. This article proposes a feature selection and weighted support vector machine (FSWSVM) classifier-based algorithm to detect ships in polarimetric SAR (PolSAR) imagery. First, the method constructs a feature vector that consists of multi-polarization parameters. Then, different polarization parameters are refined and weighted according to their significance in the support vector machine (SVM) classifier. Finally, ships are classified from the sea background and other false alarms by the classifier. The validation results on National Aeronautics and Space Administration/Jet Propulsion Laboratory (NASA/JPL) airborne synthetic aperture radar (AIRSAR) and Radarsat-2 quad polarimetric data illustrate that the method detects ship targets more precisely and reduces false alarms effectively.     

全极化SAR影像城区建筑密度分析

针对全极化SAR影像的建筑区特性,提出了一种基于极化特征共生矩阵的城区建筑密度分析方法。首先将极化特征与共生矩阵结合,在考虑建筑区极化散射机理和建筑朝向作用的同时,兼顾了建筑区的空间排列信息,在此基础上为了增强建筑密度的局部区域特性,将共生矩阵特征进行K-means聚类,结合图像分块形成标号直方图统计矢量,进而对该直方图统计矢量进行矢量量化实现SAR影像城区的建筑密度分级。RadarSat-2全极化SAR影像城区建筑密度分析的实验表明,该方法既适用于建筑朝向复杂城区也适用于建筑排列整齐城区的密度信息提取。