In this paper we investigate the performance of an algorithm for automatic segmentation of full polarimetric, synthetic aperture radar (SAR) sea ice scenes. The algorithm uses statistical and polarimetric properties of the backscattered radar signals to segment the SAR image into a specified number of classes. This number was determined in advance from visual inspection of the SAR image and by available in situ measurements. The segmentation result was then compared to ice charts drawn by ice service analysts. The comparison revealed big discrepancies between the charts of the analysts, and between the manual and the automatic segmentations. In the succeeding analysis, the automatic segmentation chart was labeled into ice types by sea ice experts, and the SAR features used in the segmentation were interpreted in terms of physical sea ice properties. Utilizing polarimetric information in sea ice charting will increase the efficiency and exactness of the maps. The number of classes used in the segmentation has shown to be of significant importance. Thus, studies of automatic and robust estimation of the number of ice classes in SAR sea ice scenes will be highly relevant for future work.

We provide sea ice classification maps of a sub-weekly time series of single (horizontal–horizontal, HH) polarization X-band TerraSAR-X scanning synthetic aperture radar (TSX SC) images from November 2019 to March 2020, covering the Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) expedition. This classified time series benefits from the wide spatial coverage and relatively high spatial resolution of TSX SC data and is a useful basic dataset for future MOSAiC studies on physical sea ice processes and ocean and climate modeling. Sea ice is classified into leads, young ice with different backscatter intensities, and first-year ice (FYI) or multiyear ice (MYI) with different degrees of deformation. We establish the per-class incidence angle (IA) dependencies of TSX SC intensities and gray-level co-occurrence matrix (GLCM) textures and use a classifier that corrects for the class-specific decreasing backscatter with increasing IAs, with both HH intensities and textures as input features. Optimal parameters for texture calculation are derived to achieve good class separation while maintaining maximum spatial detail and minimizing textural collinearity. Class probabilities yielded by the classifier are adjusted by Markov random field contextual smoothing to produce classification results. The texture-based classification process yields an average overall accuracy of 83.70 % and good correspondence to geometric ice surface roughness derived from in situ ice thickness measurements (correspondence consistently close to or higher than 80 %). A positive logarithmic relationship is found between geometric ice surface roughness and TSX SC HH backscatter intensity, similar to previous C- and L-band studies. Areal fractions of classes representing ice openings (leads and young ice) show prominent increases in middle to late November 2019 and March 2020, corresponding well to ice-opening time series derived from in situ data in this study and those derived from satellite synthetic aperture radar (SAR) and optical data in other MOSAiC studies.

We introduce the fully automatic design of a numerically optimized decision-tree algorithm and demonstrate its application to sea ice classification from SAR data. In the decision tree, an initial multi-class classification problem is split up into a sequence of binary problems. Each branch of the tree separates one single class from all other remaining classes, using a class-specific selected feature set. We optimize the order of classification steps and the feature sets by combining classification accuracy and sequential search algorithms, looping over all remaining features in each branch. The proposed strategy can be adapted to different types of classifiers and measures for the class separability. In this study, we use a Bayesian classifier with non-parametric kernel density estimation of the probability density functions. We test our algorithm on simulated data as well as airborne and spaceborne SAR data over sea ice. For the simulated cases, average per-class classification accuracy is improved between 0.5% and 4% compared to traditional all-at-once classification. Classification accuracy for the airborne and spaceborne SAR datasets was improved by 2.5% and 1%, respectively. In all cases, individual classes can show larger improvements up to 8%. Furthermore, the selection of individual feature sets for each single class can provide additional insights into physical interpretation of different features. The improvement in classification results comes at the cost of longer computation time, in particular during the design and training stage. The final choice of the optimal algorithm therefore depends on time constraints and application purpose.

This paper describes a flexible non-Gaussian statistical method used to model polarimetric synthetic aperture radar (POLSAR) data. We outline the theoretical basis of the well-know product model as described by the class of Scale Mixture models and discuss their appropriateness for modelling radar data. The statistical distributions of several Scale mixture models are then described, including the commonly used Gaussian model, and techniques for model parameter estimation are given. Real data evaluations are made using airborne fully polarimetric SAR studies for several distinct land cover types. Generic scale mixture of Gaussian features is extracted from the model parameters and a simple clustering example presented.

Synthetic Aperture Radar (SAR) systems perform signal acquisition under varying incidence angles and register an implicit intensity decay from near to far range. Due to the geometrical interaction between microwaves and the imaged targets, the rates at which intensities decay depend on the nature of the targets, thus rendering single-rate image correction approaches only partially successful. The decay, also known as the incidence angle effect, impacts the segmentation of wide-swath images performed on absolute intensity values. We propose to integrate the target-specific intensity decay rates into a non-stationary statistical model, for use in a fully automatic and unsupervised segmentation algorithm. We demonstrate this concept by assuming Gaussian distributed log-intensities and linear decay rates, a fitting approximation for the smooth, systematic decay observed for extended flat targets. The segmentation is performed on Sentinel-1, Radarsat-2 and UAVSAR wide-swath scenes containing open water, sea ice and oil slicks. As a result, we obtain segments connected throughout the entire incidence angle range, thus overcoming the limitations of modeling that does not account for different per-target decays. The model simplicity also allows for short execution times and presents the segmentation approach as a potential operational algorithm. In addition, we estimate the log-linear decay rates and examine their potential for a physical interpretation of the segments.