This study reviews the recent advances in data-driven polarimetric imaging technologies based on a wide range of practical applications.The widespread international research and activity in polarimetric imaging techni...This study reviews the recent advances in data-driven polarimetric imaging technologies based on a wide range of practical applications.The widespread international research and activity in polarimetric imaging techniques demonstrate their broad applications and interest.Polarization information is increasingly incorporated into convolutional neural networks(CNN)as a supplemental feature of objects to improve performance in computer vision task applications.Polarimetric imaging and deep learning can extract abundant information to address various challenges.Therefore,this article briefly reviews recent developments in data-driven polarimetric imaging,including polarimetric descattering,3D imaging,reflection removal,target detection,and biomedical imaging.Furthermore,we synthetically analyze the input,datasets,and loss functions and list the existing datasets and loss functions with an evaluation of their advantages and disadvantages.We also highlight the significance of data-driven polarimetric imaging in future research and development.展开更多
Accurate radar quantitative precipitation estimation(QPE)plays an essential role in disaster prevention and mitigation.In this paper,two deep learning-based QPE networks including a single-parameter network and a mult...Accurate radar quantitative precipitation estimation(QPE)plays an essential role in disaster prevention and mitigation.In this paper,two deep learning-based QPE networks including a single-parameter network and a multi-parameter network are designed.Meanwhile,a self-defined loss function(SLF)is proposed during modeling.The dataset includes Shijiazhuang S-band dual polarimetric radar(CINRAD/SAD)data and rain gauge data within the radar’s 100-km detection range during the flood season of 2021 in North China.Considering that the specific propagation phase shift(KDP)has a roughly linear relationship with the precipitation intensity,KDP is set to 0.5°km^(-1 )as a threshold value to divide all the rain data(AR)into a heavy rain(HR)and light rain(LR)dataset.Subsequently,12 deep learning-based QPE models are trained according to the input radar parameters,the precipitation datasets,and whether an SLF was adopted,respectively.The results suggest that the effects of QPE after distinguishing rainfall intensity are better than those without distinguishing,and the effects of using SLF are better than those that used MSE as a loss function.A Z-R relationship and a ZH-KDP-R synthesis method are compared with deep learning-based QPE.The mean relative errors(MRE)of AR models using SLF are improved by 61.90%,51.21%,and 56.34%compared with the Z-R relational method,and by 38.63%,42.55%,and 47.49%compared with the synthesis method.Finally,the models are further evaluated in three precipitation processes,which manifest that the deep learning-based models have significant advantages over the traditional empirical formula methods.展开更多
This paper focuses on the adaptive detection of range and Doppler dual-spread targets in non-homogeneous and nonGaussian sea clutter.The sea clutter from two polarimetric channels is modeled as a compound-Gaussian mod...This paper focuses on the adaptive detection of range and Doppler dual-spread targets in non-homogeneous and nonGaussian sea clutter.The sea clutter from two polarimetric channels is modeled as a compound-Gaussian model with different parameters,and the target is modeled as a subspace rangespread target model.The persymmetric structure is used to model the clutter covariance matrix,in order to reduce the reliance on secondary data of the designed detectors.Three adaptive polarimetric persymmetric detectors are designed based on the generalized likelihood ratio test(GLRT),Rao test,and Wald test.All the proposed detectors have constant falsealarm rate property with respect to the clutter texture,the speckle covariance matrix.Experimental results on simulated and measured data show that three adaptive detectors outperform the competitors in different clutter environments,and the proposed GLRT detector has the best detection performance under different parameters.展开更多
A medium wave(MW)640×512(25μm)Mercury Cadmium Telluride(HgCdTe)polarimetric focal plane array(FPA)was demonstrated.The micro-polarizer array(MPA)has been carefully designed in terms of line grating structure opt...A medium wave(MW)640×512(25μm)Mercury Cadmium Telluride(HgCdTe)polarimetric focal plane array(FPA)was demonstrated.The micro-polarizer array(MPA)has been carefully designed in terms of line grating structure optimization and crosstalk suppression.A monolithic fabrication process with low damage was explored,which was verified to be compatible well with HgCdTe devices.After monolithic integration of MPA,NETD<9.5 mK was still maintained.Furthermore,to figure out the underlying mechanism that dominat⁃ed the extinction ratio(ER),specialized MPA layouts were designed,and the crosstalk was experimentally vali⁃dated as the major source that impacted ER.By expanding opaque regions at pixel edges to 4μm,crosstalk rates from adjacent pixels could be effectively reduced to approximately 2%,and promising ERs ranging from 17.32 to 27.41 were implemented.展开更多
During the pre-summer rainy season,heavy rainfall occurs frequently in South China.Based on polarimetric radar observations,the microphysical characteristics and processes of convective features associated with extrem...During the pre-summer rainy season,heavy rainfall occurs frequently in South China.Based on polarimetric radar observations,the microphysical characteristics and processes of convective features associated with extreme rainfall rates(ERCFs)are examined.In the regions with high ERCF occurrence frequency,sub-regional differences are found in the lightning flash rate(LFR)distributions.In the region with higher LFRs,the ERCFs have larger volumes of high reflectivity factor above the freezing level,corresponding to more active riming processes.In addition,these ERCFs are more organized and display larger spatial coverage,which may be related to the stronger low-level wind shear and higher terrain in the region.In the region with lower LFRs,the ERCFs have lower echo tops and lower-echo centroids.However,no clear differences of the most unstable convective available potential energy(MUCAPE)exist in the ERCFs in the regions with different LFR characteristics.Regardless of the LFRs,raindrop collisional coalescence is the main process for the growth of raindrops in the ERCFs.In the ERCFs within the region with lower LFRs,the main mechanism for the rapid increase of liquid water content with decreasing altitude below 4 km is through the warm-rain processes converting cloud drops to raindrops.However,in those with higher LFRs,the liquid water content generally decreases with decreasing altitude.展开更多
The performance of different quantitative precipitation estimation(QPE) relationships is examined using the polarimetric variables from the X-band polarimetric phased-array radars in Guangzhou,China.Three QPE approach...The performance of different quantitative precipitation estimation(QPE) relationships is examined using the polarimetric variables from the X-band polarimetric phased-array radars in Guangzhou,China.Three QPE approaches,namely,R(ZH),R(ZH,ZDR) and R(KDP),are developed for horizontal reflectivity,differential reflectivity and specific phase shift rate,respectively.The estimation parameters are determined by fitting the relationships to the observed radar variables using the T-matrix method.The QPE relationships were examined using the data of four heavy precipitation events in southern China.The examination shows that the R(ZH) approach performs better for the precipitation rate less than 5 mm h-1, and R(KDP) is better for the rate higher than 5 mm h-1, while R(ZH,ZDR) has the worst performance.An adaptive approach is developed by taking the advantages of both R(ZH) and R(KDP) approaches to improve the QPE accuracy.展开更多
Compared to single-polarization synthetic aperture radar(SAR)data,fully polarimetric SAR data can provide more detailed information of the sea surface,which is important for applications such as shallow sea topography...Compared to single-polarization synthetic aperture radar(SAR)data,fully polarimetric SAR data can provide more detailed information of the sea surface,which is important for applications such as shallow sea topography detection.The Gaofen-3 satellite provides abundant polarimetric SAR data for ocean research.In this paper,a shallow sea topography detection method was proposed based on fully polarimetric Gaofen-3 SAR data.This method considers swell patterns and only requires SAR data and little prior knowledge of the water depth to detect shallow sea topography.Wave tracking was performed based on preprocessed fully polarimetric SAR data,and the water depth was then calculated considering the wave parameters and the linear dispersion relationships.In this paper,four study areas were selected for experiments,and the experimental results indicated that the polarimetric scattering parameterαhad higher detection accuracy than quad-polarization images.The mean relative errors were 14.52%,10.30%,12.56%,and 12.90%,respectively,in the four study areas.In addition,this paper also analyzed the detection ability of this model for different topographies,and the experiments revealed that the topography could be well recognized when the topography gradient is small,the topography gradient direction is close to the wave propagation direction,and the isobath line is regular.展开更多
A mobile C-band dual polarimetric weather radar J type (PCDJ), which adopts simultaneous transmission and simultaneous reception (STSR) of horizontally and vertically polarized signals, was first developed in Chin...A mobile C-band dual polarimetric weather radar J type (PCDJ), which adopts simultaneous transmission and simultaneous reception (STSR) of horizontally and vertically polarized signals, was first developed in China in 2008. It was deployed in the radar observation plan in the South China Heavy Rainfall Experiment (SCHeREX) in the summer of 2008 and 2009, as well as in Tropical Western Pacific Ocean Observation Experiments and Research on the Predictability of High Impact Weather Events from 2008 to 2010 in China (TWPOR). Using the observation data collected in these experiments, the radar systematic error and its sources were analyzed in depth. Meanwhile an algorithm that can smooth differential propagation phase (~Dp) for estimating the high-resolution specific differential phase (KDP) was developed. After attenuation correction of reflectivity in horizontal polarization (ZH) and differential reflectivity (ZDR) of PCDJ radar by means of KDP, the data quality was improved significantly. Using quality-controlled radar data, quantitative rainfall estimation was performed, and the resutls were compared with rain-gauge measurements. A synthetic ZH /KDp-based method was analyzed. The results the traditional ZH-based method when the rain suggest that the synthetic method has the advantage over rate is 〉5 mm h^-1. The more intensive the rain rates, the higher accuracy of the estimation.展开更多
The exact radar cross-section (RCS) measurement is difficult when the scattering of targets is low. Ful polarimetric cali-bration is one technique that offers the potential for improving the accuracy of RCS measurem...The exact radar cross-section (RCS) measurement is difficult when the scattering of targets is low. Ful polarimetric cali-bration is one technique that offers the potential for improving the accuracy of RCS measurements. There are numerous polarimetric calibration algorithms. Some complex expressions in these algo-rithms cannot be easily used in an engineering practice. A radar polarimetric coefficients matrix (RPCM) with a simpler expression is presented for the monostatic radar polarization scattering matrix (PSM) measurement. Using a rhombic dihedral corner reflector and a metal ic sphere, the RPCM can be obtained by solving a set of equations, which can be used to find the true PSM for any target. An example for the PSM of a metal ic dish shows that the proposed method obviously improves the accuracy of cross-polarized RCS measurements.展开更多
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 pro...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.展开更多
Multifrequency polarimetric SAR imagery provides a very convenient approach for signal processing and acquisition of radar image. However, the amount of information is scattered in several images, and redundancies exi...Multifrequency polarimetric SAR imagery provides a very convenient approach for signal processing and acquisition of radar image. However, the amount of information is scattered in several images, and redundancies exist between different bands and polarizations. Similar to signal-polarimetric SAR image, multifrequency polarimetric SAR image is corrupted with speckle noise at the same time. A method of information compression and speckle reduction for multifrequency polarimetric SAR imagery is presented based on kernel principal component analysis (KPCA). KPCA is a nonlinear generalization of the linear principal component analysis using the kernel trick. The NASA/JPL polarimetric SAR imagery of P, L, and C bands quadpolarizations is used for illustration. The experimental results show that KPCA has better capability in information compression and speckle reduction as compared with linear PCA.展开更多
Speckle filtering is an indispensable pre-processing step for applications of polarimetric synthetic aperture radar (POLSAR), such as terrain classification, target detection, etc. As one of the most typical methods...Speckle filtering is an indispensable pre-processing step for applications of polarimetric synthetic aperture radar (POLSAR), such as terrain classification, target detection, etc. As one of the most typical methods, the polarimetric whitening filter (PWF) can be used to produce a minimum-speckle image by combining the complex elements of the scattering matrix, but polarimetric information is lost after the filtering process. A polarimetric filter based on subspaze decomposition which was proposed by Cu et al specializes in retrieving principle scattering characteristics, but the corresponding mean value of an image after filtering is not kept well. A new filter is proposed for improving the disadvantage based on subspace decomposition. Under the constraint that a weighted combination of the polarimetric SAR images equals to the output of the PWF, the Euclidean distance between an unfiltered parameter vector and a signal space vector is minimized so that noises can be reduced. It is also shown that the proposed method is equivalent to the subspace filter in the case of no constraint. Experimental results with the NASA/JPL airborne polarimetric SAR data demonstrate the effectiveness of the proposed method.展开更多
Azimuth ambiguities (ghost targets) discrimination is of great interest with the development of a synthet- ic aperture radar (SAR). And the azimuth ambiguities are often mistaken as actual targets and cause false ...Azimuth ambiguities (ghost targets) discrimination is of great interest with the development of a synthet- ic aperture radar (SAR). And the azimuth ambiguities are often mistaken as actual targets and cause false alarms. For actual targets, HV channel signals acquired by a fully polarimetric SAR are approximately equal to a VH channel in magnitude and phase, i.e., the reciprocity theorem applies, but shifted in phase about ±π for the first-order azimuth ambiguities. Exploiting this physical behavior, the real part of the product of the two cross-polarized channels, i.e. (SHVSVH), hereafter called A12r, is employed as a new parameter for a target detection at sea. Compared with other parameters, the contrast of A12r image between a target and the surrounding sea surface will be obviously increased when A12r image is processed by mean filtering algo- rithm. Here, in order to detect target with constant false-alarm rates (CFARs), an analytical expression for the probability density function (pdf) ofA12r is derived based on the complexWishart-distribution. Because a value of A12r is greater/less than 0 for real target/its azimuth ambiguities, the first-order azimuth ambiguities can be completely removed by this A12r-based CFAR technology. Experiments accomplished over C-band RADARSAT-2 fully polarimetric imageries confirm the validity.展开更多
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 sca...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.展开更多
We propose a multi-sensor multi-spectral and bi-temporal dual-polarimetric Synthetic Aperture Radar(SAR) data integration scheme for dry/wet snow mapping using Sentinel-2 and Sentinel-1 data which are freely available...We propose a multi-sensor multi-spectral and bi-temporal dual-polarimetric Synthetic Aperture Radar(SAR) data integration scheme for dry/wet snow mapping using Sentinel-2 and Sentinel-1 data which are freely available to the research community. The integration is carried out by incorporating the information retrieved from ratio images of the conventional method for wet snow mapping and the multispectral data in two different frameworks. Firstly, a simple differencing scheme is employed for dry/wet snow mapping, where the snow cover area is derived using the Normalized Differenced Snow Index(NDSI). In the second framework, the ratio images are stacked with the multispectral bands and this stack is used for supervised and unsupervised classification using support vector machines for dry/wet snow mapping. We also investigate the potential of a state of the art backscatter model for the identification of dry/wet snow using Sentinel-1 data. The results are validated using a reference map derived from RADARSAT-2 full polarimetric SAR data. A good agreement was observed between the results and the reference data with an overall accuracy greater than 0.78 for the different blending techniques examined. For all the proposed frameworks, the wet snow was better identified. The coefficient of determination between the snow wetness derived from the backscatter model and the reference based on RADARSAT-2 data was observed to be 0.58 with a significantly higher root mean square error of 1.03 % by volume.展开更多
The wind retrieval performance of HY-2 A scanning scatterometer operating at Ku-band in HH and VV polarizations has been well evaluated in the wind speed range of 0–25 m s^-1.In order to obtain more accurate ocean wi...The wind retrieval performance of HY-2 A scanning scatterometer operating at Ku-band in HH and VV polarizations has been well evaluated in the wind speed range of 0–25 m s^-1.In order to obtain more accurate ocean wind field,a potential extension of dual-frequency(C-band and Ku-band)polarimetric measurements is investigated for both low and very high wind speeds,from 5 to 45 m s^-1.Based on the geophysical model functions of C-band and Ku-band,the simulation results show that the polarimetric measurements of Ku-band can improve the wind vector retrieval over the entire scatterometer swath,especially in nadir area,with the wind direction root-mean-square error(RMSE)less than 12?in the wind speed range of 5–25 m s^-1.Furthermore,the results also show that C-band cross-polarization plays a very important role in improving the wind speed retrieval,with the wind speed retrieval accuracy better than 2 m s^-1 for all wind conditions(0–45 m s^-1).For extreme winds,the C-band HH backscatter coefficients modeled by CMOD5.N(H)and the ocean co-polarization ratio model at large incidence are used to retrieve sea surface wind vector.This result reveals that there is a big decrease of wind direction retrieval RMSE for extreme wind fields,and the retrieved result of C-band HH polarization is nearly the same as that of C-band VV polarization for low-to-high wind speed(5–25 m s^-1).Thus,to improve the wind retrieval for all wind conditions,the dual-frequency polarimetric scatterometer with C-band and Ku-band horizontal polarization in inner beam,and C-band horizontal and Ku-band vertical polarization in outer beam,can be used to measure ocean winds.This study will contribute to the wind retrieval with merged satellites data and the future spaceborne scatterometer.展开更多
The UK Met Office Unified Model(UM) is employed by many weather forecasting agencies around the globe. This model is designed to run across spatial and time scales and known to produce skillful predictions for large...The UK Met Office Unified Model(UM) is employed by many weather forecasting agencies around the globe. This model is designed to run across spatial and time scales and known to produce skillful predictions for large-scale weather systems. However, the model has only recently begun running operationally at horizontal grid spacings of ~1.5 km [e.g.,at the UK Met Office and the Korea Meteorological Administration(KMA)]. As its microphysics scheme was originally designed and tuned for large-scale precipitation systems, we investigate the performance of UM microphysics to determine potential inherent biases or weaknesses. Two rainfall cases from the KMA forecasting system are considered in this study: a Changma(quasi-stationary) front, and Typhoon Sanba(2012). The UM output is compared to polarimetric radar observations in terms of simulated polarimetric radar variables. Results show that the UM generally underpredicts median reflectivity in stratiform rain, producing high reflectivity cores and precipitation gaps between them. This is partially due to the diagnostic rain intercept parameter formulation used in the one-moment microphysics scheme. Model drop size is generally both underand overpredicted compared to observations. UM frozen hydrometeors favor generic ice(crystals and snow) rather than graupel, which is reasonable for Changma and typhoon cases. The model performed best with the typhoon case in terms of simulated precipitation coverage.展开更多
An optimization of polarimetric contrast enhancement method is proposed to detect ships with lowship-to-clutter power ratio.The received power is calculated with Kennaugh matrix and an iterative algo-rithm is adopted ...An optimization of polarimetric contrast enhancement method is proposed to detect ships with lowship-to-clutter power ratio.The received power is calculated with Kennaugh matrix and an iterative algo-rithm is adopted to get the optimal polarimetric states.The optimization method depresses the power of o-cean clutter and increases the power of ship signal.With the double effects,the contrast of ship to oceanis dramatically increased.Thus small ship or weak signals of low ship-to-ocean power ratio can easily bedetected.Ship signals can be distinguished from speckle noise using the different variation trend after op-timization,and thus the threshold problem can be avoided.Moreover,the analyses of different ship'sKennaugh matrices give two implications.One is that the results are affected little by choosing differentKennaugh matrices of ships with strong intensity from Synthetic Aperture Radar(SAR)images.The otheris that ship's Kennaugh matrix chosen from real SAR images is more favorable than that of ideal scatter-ing.Finally,the optimization results are confirmed by polarimetric scattering angle and co-polarizationphase difference.展开更多
Polarimetric radar and 2D video disdrometer observations provide new insights into the precipitation microphysical processes and characteristics in the inner rainband of tropical cyclone(TC)Kajiki(2019)in the South Ch...Polarimetric radar and 2D video disdrometer observations provide new insights into the precipitation microphysical processes and characteristics in the inner rainband of tropical cyclone(TC)Kajiki(2019)in the South China Sea for the first time.The precipitation of Kajiki is dominated by high concentrations and small(<3 mm)raindrops,which contribute more than 98%to the total precipitation.The average mass-weighted mean diameter and logarithmic normalized intercept are 1.49 mm and 4.47,respectively,indicating a larger mean diameter and a lower concentration compared to the TCs making landfall in eastern China.The ice processes of the inner rainband are dramatically different among different stages.The riming process is dominant during the mature stage,while during the decay stage the aggregation process is dominant.The vertical profiles of the polarimetric radar variables together with ice and liquid water contents in the convective region indicate that the formation of precipitation is dominated by warm-rain processes.Large raindrops collect cloud droplets and other raindrops,causing reflectivity,differential reflectivity,and specific differential phase to increase with decreasing height.That is,accretion and coalescence play a critical role in the formation of heavy rainfall.The melting of different particles generated by the ice process has a great influence on the initial raindrop size distribution(DSD)to further affect the warm-rain processes.The DSD above heavy rain with the effect of graupel has a wider spectral width than the region without the effect of graupel.展开更多
For better interpretation of synthetic aperture radar(SAR) images,the speckle filtering is an important issue.In the area of speckle filtering,the proper averaging of samples with similar scattering characteristics ...For better interpretation of synthetic aperture radar(SAR) images,the speckle filtering is an important issue.In the area of speckle filtering,the proper averaging of samples with similar scattering characteristics is of great importance.However,existing filtering algorithms are either lack of a similarity judgment of scattering characteristics or using only intensity information for similarity judgment.A novel polarimetric SAR(PolSAR) speckle filtering algorithm based on the mean shift theory is proposed.As polarimetric covariance matrices or coherency matrices form Riemannian manifold,the pixels with similar scattering characteristics gather closely and those with different scattering characteristics separate in this hyperspace.By using the range-spatial joint mean shift theory in Riemannian manifold,the pixels chosen for averaging are ensured to be close not only in scattering characteristics but also in the spatial domain.German Aerospace Center(DLR) L-Band Experiment SAR(E-SAR) data and East China Research Institute of Electronic Engineering(ECRIEE) PolSAR data are used to demonstrate the efficiency of the proposed algorithm.The filtering results of two commonly used speckle filtering algorithms,refined Lee filtering algorithm and intensity driven adaptive neighborhood(IDAN) filtering algorithm,are also presented for the comparison purpose.Experiment results show that the proposed speckle filtering algorithm achieves a good performance in terms of speckle filtering,edge protection as well as polarimetric characteristics preservation.展开更多
基金support from the National Natural Science Foundation of China(Nos.62205259,62075175,61975254,62375212,62005203 and 62105254)the Open Research Fund of CAS Key Laboratory of Space Precision Measurement Technology(No.B022420004)the Fundamental Research Funds for the Central Universities(No.ZYTS23125).
文摘This study reviews the recent advances in data-driven polarimetric imaging technologies based on a wide range of practical applications.The widespread international research and activity in polarimetric imaging techniques demonstrate their broad applications and interest.Polarization information is increasingly incorporated into convolutional neural networks(CNN)as a supplemental feature of objects to improve performance in computer vision task applications.Polarimetric imaging and deep learning can extract abundant information to address various challenges.Therefore,this article briefly reviews recent developments in data-driven polarimetric imaging,including polarimetric descattering,3D imaging,reflection removal,target detection,and biomedical imaging.Furthermore,we synthetically analyze the input,datasets,and loss functions and list the existing datasets and loss functions with an evaluation of their advantages and disadvantages.We also highlight the significance of data-driven polarimetric imaging in future research and development.
基金supported by National Key R&D Program of China(Grant No.2022YFC3003903)the S&T Program of Hebei(Grant No.19275408D),the Key-Area Research and Development Program of Guangdong Province(Grant No.2020B1111200001)+1 种基金the Key Project of Monitoring,Early Warning and Prevention of Major Natural Disasters of China(Grant No.2019YFC1510304)the Joint Fund of Key Laboratory of Atmosphere Sounding,CMA,and the Research Centre on Meteorological Observation Engineering Technology,CMA(Grant No.U2021Z05).
文摘Accurate radar quantitative precipitation estimation(QPE)plays an essential role in disaster prevention and mitigation.In this paper,two deep learning-based QPE networks including a single-parameter network and a multi-parameter network are designed.Meanwhile,a self-defined loss function(SLF)is proposed during modeling.The dataset includes Shijiazhuang S-band dual polarimetric radar(CINRAD/SAD)data and rain gauge data within the radar’s 100-km detection range during the flood season of 2021 in North China.Considering that the specific propagation phase shift(KDP)has a roughly linear relationship with the precipitation intensity,KDP is set to 0.5°km^(-1 )as a threshold value to divide all the rain data(AR)into a heavy rain(HR)and light rain(LR)dataset.Subsequently,12 deep learning-based QPE models are trained according to the input radar parameters,the precipitation datasets,and whether an SLF was adopted,respectively.The results suggest that the effects of QPE after distinguishing rainfall intensity are better than those without distinguishing,and the effects of using SLF are better than those that used MSE as a loss function.A Z-R relationship and a ZH-KDP-R synthesis method are compared with deep learning-based QPE.The mean relative errors(MRE)of AR models using SLF are improved by 61.90%,51.21%,and 56.34%compared with the Z-R relational method,and by 38.63%,42.55%,and 47.49%compared with the synthesis method.Finally,the models are further evaluated in three precipitation processes,which manifest that the deep learning-based models have significant advantages over the traditional empirical formula methods.
基金supported by the National Natural Science Foundation of China(62371382,62071346)the Science,Technology&Innovation Project of Xiong’an New Area(2022XAGG0181)the Special Funds for Creative Research(2022C61540)。
文摘This paper focuses on the adaptive detection of range and Doppler dual-spread targets in non-homogeneous and nonGaussian sea clutter.The sea clutter from two polarimetric channels is modeled as a compound-Gaussian model with different parameters,and the target is modeled as a subspace rangespread target model.The persymmetric structure is used to model the clutter covariance matrix,in order to reduce the reliance on secondary data of the designed detectors.Three adaptive polarimetric persymmetric detectors are designed based on the generalized likelihood ratio test(GLRT),Rao test,and Wald test.All the proposed detectors have constant falsealarm rate property with respect to the clutter texture,the speckle covariance matrix.Experimental results on simulated and measured data show that three adaptive detectors outperform the competitors in different clutter environments,and the proposed GLRT detector has the best detection performance under different parameters.
基金Supported by the self-funded project of Kunming Institute of Physics。
文摘A medium wave(MW)640×512(25μm)Mercury Cadmium Telluride(HgCdTe)polarimetric focal plane array(FPA)was demonstrated.The micro-polarizer array(MPA)has been carefully designed in terms of line grating structure optimization and crosstalk suppression.A monolithic fabrication process with low damage was explored,which was verified to be compatible well with HgCdTe devices.After monolithic integration of MPA,NETD<9.5 mK was still maintained.Furthermore,to figure out the underlying mechanism that dominat⁃ed the extinction ratio(ER),specialized MPA layouts were designed,and the crosstalk was experimentally vali⁃dated as the major source that impacted ER.By expanding opaque regions at pixel edges to 4μm,crosstalk rates from adjacent pixels could be effectively reduced to approximately 2%,and promising ERs ranging from 17.32 to 27.41 were implemented.
基金primarily supported by the National Natural Science Foundation of China(Grant Nos.42025501,41905019,and 61827901)the National Key Research and Development Program of China(Grant 2018YFC1506404 and Grant 2017YFC1501703)。
文摘During the pre-summer rainy season,heavy rainfall occurs frequently in South China.Based on polarimetric radar observations,the microphysical characteristics and processes of convective features associated with extreme rainfall rates(ERCFs)are examined.In the regions with high ERCF occurrence frequency,sub-regional differences are found in the lightning flash rate(LFR)distributions.In the region with higher LFRs,the ERCFs have larger volumes of high reflectivity factor above the freezing level,corresponding to more active riming processes.In addition,these ERCFs are more organized and display larger spatial coverage,which may be related to the stronger low-level wind shear and higher terrain in the region.In the region with lower LFRs,the ERCFs have lower echo tops and lower-echo centroids.However,no clear differences of the most unstable convective available potential energy(MUCAPE)exist in the ERCFs in the regions with different LFR characteristics.Regardless of the LFRs,raindrop collisional coalescence is the main process for the growth of raindrops in the ERCFs.In the ERCFs within the region with lower LFRs,the main mechanism for the rapid increase of liquid water content with decreasing altitude below 4 km is through the warm-rain processes converting cloud drops to raindrops.However,in those with higher LFRs,the liquid water content generally decreases with decreasing altitude.
基金Guangzhou Science and Technology Plan Project(202103000030)Guangdong Meteorological Bureau Science and Technology Project(GRMC2020Z08)a project co-funded by the Development Team of Radar Application and Severe Convection Early Warning Technology(GRMCTD202002)。
文摘The performance of different quantitative precipitation estimation(QPE) relationships is examined using the polarimetric variables from the X-band polarimetric phased-array radars in Guangzhou,China.Three QPE approaches,namely,R(ZH),R(ZH,ZDR) and R(KDP),are developed for horizontal reflectivity,differential reflectivity and specific phase shift rate,respectively.The estimation parameters are determined by fitting the relationships to the observed radar variables using the T-matrix method.The QPE relationships were examined using the data of four heavy precipitation events in southern China.The examination shows that the R(ZH) approach performs better for the precipitation rate less than 5 mm h-1, and R(KDP) is better for the rate higher than 5 mm h-1, while R(ZH,ZDR) has the worst performance.An adaptive approach is developed by taking the advantages of both R(ZH) and R(KDP) approaches to improve the QPE accuracy.
基金The National Natural Science Foundation of China under contract Nos 51839002 and U2006207.
文摘Compared to single-polarization synthetic aperture radar(SAR)data,fully polarimetric SAR data can provide more detailed information of the sea surface,which is important for applications such as shallow sea topography detection.The Gaofen-3 satellite provides abundant polarimetric SAR data for ocean research.In this paper,a shallow sea topography detection method was proposed based on fully polarimetric Gaofen-3 SAR data.This method considers swell patterns and only requires SAR data and little prior knowledge of the water depth to detect shallow sea topography.Wave tracking was performed based on preprocessed fully polarimetric SAR data,and the water depth was then calculated considering the wave parameters and the linear dispersion relationships.In this paper,four study areas were selected for experiments,and the experimental results indicated that the polarimetric scattering parameterαhad higher detection accuracy than quad-polarization images.The mean relative errors were 14.52%,10.30%,12.56%,and 12.90%,respectively,in the four study areas.In addition,this paper also analyzed the detection ability of this model for different topographies,and the experiments revealed that the topography could be well recognized when the topography gradient is small,the topography gradient direction is close to the wave propagation direction,and the isobath line is regular.
基金funded by National Natural Science Foundation of China (Grant Nos. 40975013 and 40975014)Chinese Academy of Meteorological Sciences (CAMS) basic scientific and operational project:Observation and retrieval methods of microphysics and dynamic parameters of cloud and precipitation with multi-wavelength Remote Sensing,SCHeREX and TWPOR
文摘A mobile C-band dual polarimetric weather radar J type (PCDJ), which adopts simultaneous transmission and simultaneous reception (STSR) of horizontally and vertically polarized signals, was first developed in China in 2008. It was deployed in the radar observation plan in the South China Heavy Rainfall Experiment (SCHeREX) in the summer of 2008 and 2009, as well as in Tropical Western Pacific Ocean Observation Experiments and Research on the Predictability of High Impact Weather Events from 2008 to 2010 in China (TWPOR). Using the observation data collected in these experiments, the radar systematic error and its sources were analyzed in depth. Meanwhile an algorithm that can smooth differential propagation phase (~Dp) for estimating the high-resolution specific differential phase (KDP) was developed. After attenuation correction of reflectivity in horizontal polarization (ZH) and differential reflectivity (ZDR) of PCDJ radar by means of KDP, the data quality was improved significantly. Using quality-controlled radar data, quantitative rainfall estimation was performed, and the resutls were compared with rain-gauge measurements. A synthetic ZH /KDp-based method was analyzed. The results the traditional ZH-based method when the rain suggest that the synthetic method has the advantage over rate is 〉5 mm h^-1. The more intensive the rain rates, the higher accuracy of the estimation.
基金supported by the National Basic Research Program of China(973 Program)(2010CB731905)
文摘The exact radar cross-section (RCS) measurement is difficult when the scattering of targets is low. Ful polarimetric cali-bration is one technique that offers the potential for improving the accuracy of RCS measurements. There are numerous polarimetric calibration algorithms. Some complex expressions in these algo-rithms cannot be easily used in an engineering practice. A radar polarimetric coefficients matrix (RPCM) with a simpler expression is presented for the monostatic radar polarization scattering matrix (PSM) measurement. Using a rhombic dihedral corner reflector and a metal ic sphere, the RPCM can be obtained by solving a set of equations, which can be used to find the true PSM for any target. An example for the PSM of a metal ic dish shows that the proposed method obviously improves the accuracy of cross-polarized RCS measurements.
基金The National Natural Science Foundation of China under contract Nos 41306186,41076012 and 41276019the Youth Science Fund Project of State Oceanic Administration of China
文摘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 Specialized Research Found for the Doctoral Program of Higher Education (20070699013)the Natural Science Foundation of Shaanxi Province (2006F05)the Aeronautical Science Foundation (05I53076).
文摘Multifrequency polarimetric SAR imagery provides a very convenient approach for signal processing and acquisition of radar image. However, the amount of information is scattered in several images, and redundancies exist between different bands and polarizations. Similar to signal-polarimetric SAR image, multifrequency polarimetric SAR image is corrupted with speckle noise at the same time. A method of information compression and speckle reduction for multifrequency polarimetric SAR imagery is presented based on kernel principal component analysis (KPCA). KPCA is a nonlinear generalization of the linear principal component analysis using the kernel trick. The NASA/JPL polarimetric SAR imagery of P, L, and C bands quadpolarizations is used for illustration. The experimental results show that KPCA has better capability in information compression and speckle reduction as compared with linear PCA.
基金supported by the National Natural Science Foundation of China (40571099)the Research Fund for the Doctoral Program of Higher Education of China.
文摘Speckle filtering is an indispensable pre-processing step for applications of polarimetric synthetic aperture radar (POLSAR), such as terrain classification, target detection, etc. As one of the most typical methods, the polarimetric whitening filter (PWF) can be used to produce a minimum-speckle image by combining the complex elements of the scattering matrix, but polarimetric information is lost after the filtering process. A polarimetric filter based on subspaze decomposition which was proposed by Cu et al specializes in retrieving principle scattering characteristics, but the corresponding mean value of an image after filtering is not kept well. A new filter is proposed for improving the disadvantage based on subspace decomposition. Under the constraint that a weighted combination of the polarimetric SAR images equals to the output of the PWF, the Euclidean distance between an unfiltered parameter vector and a signal space vector is minimized so that noises can be reduced. It is also shown that the proposed method is equivalent to the subspace filter in the case of no constraint. Experimental results with the NASA/JPL airborne polarimetric SAR data demonstrate the effectiveness of the proposed method.
基金The National Natural Science Foundation of China under contract Nos 41376179 and 41106153
文摘Azimuth ambiguities (ghost targets) discrimination is of great interest with the development of a synthet- ic aperture radar (SAR). And the azimuth ambiguities are often mistaken as actual targets and cause false alarms. For actual targets, HV channel signals acquired by a fully polarimetric SAR are approximately equal to a VH channel in magnitude and phase, i.e., the reciprocity theorem applies, but shifted in phase about ±π for the first-order azimuth ambiguities. Exploiting this physical behavior, the real part of the product of the two cross-polarized channels, i.e. (SHVSVH), hereafter called A12r, is employed as a new parameter for a target detection at sea. Compared with other parameters, the contrast of A12r image between a target and the surrounding sea surface will be obviously increased when A12r image is processed by mean filtering algo- rithm. Here, in order to detect target with constant false-alarm rates (CFARs), an analytical expression for the probability density function (pdf) ofA12r is derived based on the complexWishart-distribution. Because a value of A12r is greater/less than 0 for real target/its azimuth ambiguities, the first-order azimuth ambiguities can be completely removed by this A12r-based CFAR technology. Experiments accomplished over C-band RADARSAT-2 fully polarimetric imageries confirm the validity.
基金The National Natural Science Foundation of China under contract Nos 41576170 and 41376179the Public Science and Technology Research Funds Projects of Ocean(Ocean University of China) under contract No.2013418025-2
文摘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.
基金partly supported by Project number DST-2016056, funded by the Department of Science and Technology, Government of India
文摘We propose a multi-sensor multi-spectral and bi-temporal dual-polarimetric Synthetic Aperture Radar(SAR) data integration scheme for dry/wet snow mapping using Sentinel-2 and Sentinel-1 data which are freely available to the research community. The integration is carried out by incorporating the information retrieved from ratio images of the conventional method for wet snow mapping and the multispectral data in two different frameworks. Firstly, a simple differencing scheme is employed for dry/wet snow mapping, where the snow cover area is derived using the Normalized Differenced Snow Index(NDSI). In the second framework, the ratio images are stacked with the multispectral bands and this stack is used for supervised and unsupervised classification using support vector machines for dry/wet snow mapping. We also investigate the potential of a state of the art backscatter model for the identification of dry/wet snow using Sentinel-1 data. The results are validated using a reference map derived from RADARSAT-2 full polarimetric SAR data. A good agreement was observed between the results and the reference data with an overall accuracy greater than 0.78 for the different blending techniques examined. For all the proposed frameworks, the wet snow was better identified. The coefficient of determination between the snow wetness derived from the backscatter model and the reference based on RADARSAT-2 data was observed to be 0.58 with a significantly higher root mean square error of 1.03 % by volume.
基金supported by the National Key R&D Program of China (No. 2016YFC1401006)the National Natural Science Foundation of China (Nos. 51279186, 51479183 and 41676169)+2 种基金the National Program on Key Basic Research Project (No. 2011CB013704)the 111 Project (No. B14028)the Marine and Fishery Information Center Project of Jiangsu Province (No. SJC2014 110338)
文摘The wind retrieval performance of HY-2 A scanning scatterometer operating at Ku-band in HH and VV polarizations has been well evaluated in the wind speed range of 0–25 m s^-1.In order to obtain more accurate ocean wind field,a potential extension of dual-frequency(C-band and Ku-band)polarimetric measurements is investigated for both low and very high wind speeds,from 5 to 45 m s^-1.Based on the geophysical model functions of C-band and Ku-band,the simulation results show that the polarimetric measurements of Ku-band can improve the wind vector retrieval over the entire scatterometer swath,especially in nadir area,with the wind direction root-mean-square error(RMSE)less than 12?in the wind speed range of 5–25 m s^-1.Furthermore,the results also show that C-band cross-polarization plays a very important role in improving the wind speed retrieval,with the wind speed retrieval accuracy better than 2 m s^-1 for all wind conditions(0–45 m s^-1).For extreme winds,the C-band HH backscatter coefficients modeled by CMOD5.N(H)and the ocean co-polarization ratio model at large incidence are used to retrieve sea surface wind vector.This result reveals that there is a big decrease of wind direction retrieval RMSE for extreme wind fields,and the retrieved result of C-band HH polarization is nearly the same as that of C-band VV polarization for low-to-high wind speed(5–25 m s^-1).Thus,to improve the wind retrieval for all wind conditions,the dual-frequency polarimetric scatterometer with C-band and Ku-band horizontal polarization in inner beam,and C-band horizontal and Ku-band vertical polarization in outer beam,can be used to measure ocean winds.This study will contribute to the wind retrieval with merged satellites data and the future spaceborne scatterometer.
基金supported by a research grant of “Development of a Polarimetric Radar Data Simulator for Local Forecasting Model (Ⅱ)” by the KMAsupport was provided by a NOAA Warn-on-Forecast grant (Grant No. NA16OAR4320115)a National Science Foundation grant (Grant No. AGS-1261776)
文摘The UK Met Office Unified Model(UM) is employed by many weather forecasting agencies around the globe. This model is designed to run across spatial and time scales and known to produce skillful predictions for large-scale weather systems. However, the model has only recently begun running operationally at horizontal grid spacings of ~1.5 km [e.g.,at the UK Met Office and the Korea Meteorological Administration(KMA)]. As its microphysics scheme was originally designed and tuned for large-scale precipitation systems, we investigate the performance of UM microphysics to determine potential inherent biases or weaknesses. Two rainfall cases from the KMA forecasting system are considered in this study: a Changma(quasi-stationary) front, and Typhoon Sanba(2012). The UM output is compared to polarimetric radar observations in terms of simulated polarimetric radar variables. Results show that the UM generally underpredicts median reflectivity in stratiform rain, producing high reflectivity cores and precipitation gaps between them. This is partially due to the diagnostic rain intercept parameter formulation used in the one-moment microphysics scheme. Model drop size is generally both underand overpredicted compared to observations. UM frozen hydrometeors favor generic ice(crystals and snow) rather than graupel, which is reasonable for Changma and typhoon cases. The model performed best with the typhoon case in terms of simulated precipitation coverage.
基金the High Technology Research and Development Programme of China(No.2002AA633120)Sharing and Opening Projects of ENVISAT ASAR Data
文摘An optimization of polarimetric contrast enhancement method is proposed to detect ships with lowship-to-clutter power ratio.The received power is calculated with Kennaugh matrix and an iterative algo-rithm is adopted to get the optimal polarimetric states.The optimization method depresses the power of o-cean clutter and increases the power of ship signal.With the double effects,the contrast of ship to oceanis dramatically increased.Thus small ship or weak signals of low ship-to-ocean power ratio can easily bedetected.Ship signals can be distinguished from speckle noise using the different variation trend after op-timization,and thus the threshold problem can be avoided.Moreover,the analyses of different ship'sKennaugh matrices give two implications.One is that the results are affected little by choosing differentKennaugh matrices of ships with strong intensity from Synthetic Aperture Radar(SAR)images.The otheris that ship's Kennaugh matrix chosen from real SAR images is more favorable than that of ideal scatter-ing.Finally,the optimization results are confirmed by polarimetric scattering angle and co-polarizationphase difference.
基金This work was primarily supported by the National Key Research and Development Program of China(Grant No.2018YFC1507304)the National Natural Science Foundation of China(Grant Nos.42075080,41975066 and 41865009).
文摘Polarimetric radar and 2D video disdrometer observations provide new insights into the precipitation microphysical processes and characteristics in the inner rainband of tropical cyclone(TC)Kajiki(2019)in the South China Sea for the first time.The precipitation of Kajiki is dominated by high concentrations and small(<3 mm)raindrops,which contribute more than 98%to the total precipitation.The average mass-weighted mean diameter and logarithmic normalized intercept are 1.49 mm and 4.47,respectively,indicating a larger mean diameter and a lower concentration compared to the TCs making landfall in eastern China.The ice processes of the inner rainband are dramatically different among different stages.The riming process is dominant during the mature stage,while during the decay stage the aggregation process is dominant.The vertical profiles of the polarimetric radar variables together with ice and liquid water contents in the convective region indicate that the formation of precipitation is dominated by warm-rain processes.Large raindrops collect cloud droplets and other raindrops,causing reflectivity,differential reflectivity,and specific differential phase to increase with decreasing height.That is,accretion and coalescence play a critical role in the formation of heavy rainfall.The melting of different particles generated by the ice process has a great influence on the initial raindrop size distribution(DSD)to further affect the warm-rain processes.The DSD above heavy rain with the effect of graupel has a wider spectral width than the region without the effect of graupel.
基金supported by the National Natural Science Foundation of China(61101180)the China Postdoctoral Science Foundation (20110490088)
文摘For better interpretation of synthetic aperture radar(SAR) images,the speckle filtering is an important issue.In the area of speckle filtering,the proper averaging of samples with similar scattering characteristics is of great importance.However,existing filtering algorithms are either lack of a similarity judgment of scattering characteristics or using only intensity information for similarity judgment.A novel polarimetric SAR(PolSAR) speckle filtering algorithm based on the mean shift theory is proposed.As polarimetric covariance matrices or coherency matrices form Riemannian manifold,the pixels with similar scattering characteristics gather closely and those with different scattering characteristics separate in this hyperspace.By using the range-spatial joint mean shift theory in Riemannian manifold,the pixels chosen for averaging are ensured to be close not only in scattering characteristics but also in the spatial domain.German Aerospace Center(DLR) L-Band Experiment SAR(E-SAR) data and East China Research Institute of Electronic Engineering(ECRIEE) PolSAR data are used to demonstrate the efficiency of the proposed algorithm.The filtering results of two commonly used speckle filtering algorithms,refined Lee filtering algorithm and intensity driven adaptive neighborhood(IDAN) filtering algorithm,are also presented for the comparison purpose.Experiment results show that the proposed speckle filtering algorithm achieves a good performance in terms of speckle filtering,edge protection as well as polarimetric characteristics preservation.