Environmental and sensor challenges pose difficulties for the development of computer-assisted algorithms to segment synthetic aperture radar (SAR) sea ice imagery. In this research, in support of operational activities at the Canadian Ice Service, images containing visually separable classes of either ice and water or multiple ice classes are segmented. This work uses image intensity to discriminate ice from water and uses texture features to identify distinct ice types. In order to seamlessly combine image spatial relationships with various image features, a novel Bayesian segmentation approach is developed and applied. This new approach uses a function-based parameter to weight the two components in a Markov random field (MRF) model. The devised model allows for automatic estimation of MRF model parameters to produce accurate unsupervised segmentation results. Experiments demonstrate that the proposed algorithm is able to successfully segment various SAR sea ice images and achieve improvement over existing published methods including the standard MRF-based method, finite Gamma mixture model, and K-means clustering.

This paper presents a map-guided sea ice classification system built to work in parallel with the Canadian Ice Service (CIS) operations to produce pixel-based ice maps that complement actual "egg code" maps produced by CIS. The system uses the CIS maps as input to guide classification by providing information on the number of ice types and their final label for specific regions. Segmentation is based on a modified adaptive Markov random field (MRF) model that uses synthetic aperture radar (SAR) intensities and texture features as input. The ice type labeling is performed automatically by gathering evidences based on a priori information on one or two classes and deducing the other labels iteratively by comparing distributions of segments. Three methods for comparing the segment distributions (Fisher criterion, Mahalanobis distance, and Kolmogorov-Smirnov test) were implemented and compared. The system is fully described with special attention to the labeling procedure. Examples are presented in the form of two CIS SAR-based ice maps from the Gulf of Saint Lawrence region and one example from the Beaufort Sea. The results indicate that when the segmentation is good, the labeling attains best results (between 71% and 89%) based on evaluation by a sea ice analyst. Some problems remain to be assessed which are primarily attributable to discrepancies in the information provided by the egg code and what is actually visible in the SAR image. Subscale information on floe size and shape available to human analysts, but not in this classification system, also appear to be a critical information for separating some ice types.

Features based on Markov random field (MRF) models are sensitive to texture rotation. This paper develops an anisotropic circular Gaussian MRF (ACGMRF) model for retrieving rotation-invariant texture features. To overcome the singularity problem of the least squares estimate method, an approximate least squares estimate method is designed and implemented. Rotation-invariant features are obtained from the ACGMRF model parameters using the discrete Fourier transform. The ACGMRF model is demonstrated to be a statistical improvement over three published methods. The three methods include a Laplacian pyramid, an isotropic circular GMRF (ICGMRF), and gray level cooccurrence probability features.

The features based on Markov random field (MRF) models are usually sensitive to the rotation of image textures. The paper develops an anisotropic circular Gaussian MRF (ACGMRF) model for modeling rotated image textures and retrieving rotation-invariant texture features. To overcome the singularity problem of the least squares estimate (LSE) method, an approximate least squares estimate (ALSE) method is proposed to estimate the parameters of ACGMRF model. The rotation-invariant features can be obtained from the parameters of the ACGMRF model by the one-dimensional (1D) discrete Fourier transform (DFT). Significantly improved accuracy can be achieved by applying the rotation-invariant features to classify SAR (synthetic aperture radar) sea ice and Brodatz imagery.