## Servicios Personalizados

## Revista

## Articulo

## Indicadores

- Citado por SciELO
- Accesos

## Links relacionados

## Compartir

## Gayana (Concepción)

##
*versión impresa* ISSN 0717-652X*versión On-line* ISSN 0717-6538

### Gayana (Concepc.) v.68 n.2 supl.TIIProc Concepción 2004

#### http://dx.doi.org/10.4067/S0717-65382004000300021

Gayana 68(2): 432-438, 2004
Royal Military Academy, Brussels, Belgium
This paper presents a new method to perform sea/ice discrimination in single-pass ERS-2 scatterometer data. Existing methods are 1rst reviewed and compared in a consistent framework. Next, the ice probability according to the individual existing methods is learned through the use of a neural network. Finally, the individual criteria are combined together in order to increase the sea-ice discrimination accuracy. The proposed method is shown to provide an acceptable performance even on single-pass data, i.e., without requiring temporal averaging.
Besides the study of the ice coverage at the poles, discrimination between open sea and ice is relevant to sea-surface wind extraction using scatterometer data since no wind can be extracted over icecovered areas. Indeed, winds obtained from ERS-2 s , ^{0}_{f}s^{0}_{m},s) measured respectively by the fore, mid and aft antenna, are ambiguous and removing this ambiguity requires spatial coherence of the wind field. Hence, having erroneous wind vectors disturbs the wind ambiguity removal process. ^{0}_{a}Moreover, from an operational point of view and taking into account the near real-time dissemination Further author information: (Send correspondence to Xavier Neyt) E-mail: Xavier.Neyt@rma.ac.be, Telephone: +32 2 737 6473, Address: Signal and Image Center, Electrical Engineering Dept, Royal Military Academy, 30 av de la Renaissance, B-1000 Brussels, Belgium. Requirement, it is advantageous to be able to discriminate ice from open-sea using only the data acquired during the pass being processed. Several criteria for ice discrimination have been proposed in the literature. In Section 2, these different criteria are reviewed. Section 3 compares the performances of each of these criteria in order to assess the conditions under which they are discriminant. This is done using a neural network architecture. The network is trained to "learn" the ice-probability according to the considered training set. After training, the network is able to provide an ice probability based on one single measurement. In the same section, we propose to combine the output of individual criteria to a single result. Finally, the performance of the overall scheme is discussed.
The basis of the discrimination is that the statistics of the
Besides the raw
In (1), (2) a measure of the isotropy of the backscattering is proposed as discriminating criterion. The isotropy fator is definied as
where the s are provided in dB. Open sea typically results in a high isotropy factor, while sea-ice corresponds to rather low isotropy factor values. However, when the wind is blowing parallel or perpendicular to the ground track, a low isotropy value will be observed over sea. Indeed, in this case, the capillary waves generated by the wind are propagating perpendicular or parallel to the satellite ground track hence the fore and aft measurements will have similar values. Figure 3 (left) illustrates the repartition of isotropy factor between sea and ice for the various incidence angles. Clearly, the isotropy factor corresponding to sea-ice is typically small. However, as can be seen, small isotropy factor over open sea are possible too.
In,3 the derivative of the
where the so values are in dB and the incidence angle values in degrees. Notice that, due to the use of the zero-gyro mode following the gyroscopes anomaly, the incidence angles of the fore and of the aft beam cannot be assumed equal anymore [4] hence the modification of the original formula given in [3] Measurements over ice typically result in a lower D than measurements over sea, at least when the incidence angle is not too high. From figure 3 (right) illustrating the repartition of the derivative between sea and ice nodes for the various incidence angles, it is clear that the criterion is more discriminant for low incidence angles (low node number, near swath). 2.4. Distance to wind model The (Euclidean) distance to the wind model was also proposed [5] as a criterion. The reasoning behind this criterion is that measurements acquired over sea should lie close to the wind model, while measurements performed over ice might be located further away from the wind model. Figure 4. (left) shows that this criterion is hardly discriminating for medium incidence angle. This is related to the fact that for medium incidence angles, S
Several models of the so over ice where proposed [6], [7] The idea is to use the euclidean distance to the ice-model as discriminating factor. The models selected in this document is the ERS ice model described in.7 This model can however be generalized to other types of scatterometers. The model consists in an incidence angle-dependent line. For the details of the computation, the reader is referred to [7] As can be seen in figure 4 (right), this metric is essentially discriminant for low incidence angles. Indeed, for higher incidence angles, the distribution of ice and sea nodes tends to be closer to each other as can be seen in figures 1 and 2.
Typically, the criteria described in the previous section are thresholded to determine whether the measurement made corresponds to open sea or to ice [1, 7] That threshold is also typically dependent on the incidence angle. However, a binary threshold does not take into account the fact that, in ambiguous cases, it is not possible to make a decision based only on one single s-triplet measurement. In [7] a decision algorithm using 4 classes (sea, ice, mixed and not valid) is proposed. This algorithm is based on an incidence-angle dependent thresholding of the distances to the ice model and to the wind model. The existence of the mixed class acknowledges the fact that some nodes cannot clearly be classified as ice or sea. These decision methods provide a binary (quaternary) answer regarding the status of the considered measurement point (sea or ice). No information is provided on the (un)certainty of the classification and the classification accuracy trade off is difficult to master. In the following section, we define an ice probability based on the criteria defined in the previous section. This ice probability can then be thresholded to perform a classification. Furthermore, the criteria presented above can be combined together, as individual experts, to provide a combined ice probability.
Our goal is thus to compute the ice probability given by the measurements, i corresponds to ice m is the numerical value of the criterion _{c,i}c for measurement i and nis the across-track node number at which measurement i was made. _{i} It is well known [8] that this probability can be learned by a Multi-Layer Perceptron (MLP). The learning process consists in feeding the MLP with the measurements made ( H_{1} is true for measurement i and 0 else. In order to avoid biasing the MLP output due to a differing a priori probability we must ensure P(H_{1}) = P(H_{1}) over the training set. This imply that the number of measurements corresponding to ice in the training sets must be equal to the number of measurements corresponding to sea. We considered an MLP with two inputs, one for the value of the criterion n. It has one single hidden layer counting 5 neurons. The number of hidden layers and the number of neurones in these layers govern the complexity of the non-linear function that the MLP will be able to approximate. We are actually seeking to approximate reasonably "simple" functions, hence the single hidden layer and the small number of neurons in that hidden layer. _{i}The results of the learning of the probability by the MLP for the different criteria are shown in figure 5. For clarity, the 3D-surface was thresholded. As can be seen from this figure, the uncertain areas for the sigma anisotropy and for the distance to wind model criteria are quite large. In itself a large uncertain area is not bad, as long as no measurement samples fall inside that area. As can be seen by comparing figure 5 with the figures from Section 2, there are actually a lot of samples that fall inside that uncertain area for the Isotropy and for the Distance to the wind model criteria.
Consequently, these methods exhibit a high rate of "No Decision" answers.
The output of the MLP provides the probability that a given input corresponds to ice ( Figure 6 (left) shows the ROC for ice, where the independent variable is the threshold used in taking the decision. As can be seen, the "Distance to ice mode" criterion is the most discriminant, closely followed by the "Derivative of sigma" The two other criteria would imply a much lower True Ice rate if a low False sea rate was to be achieved. Figure 6 (right) shows the unclassified ice rate in function of the across track node-number (which corresponds approximatively to an incidence angle).
As can be seen, the "Distance to wind model" criterion fails to discriminate between ice and sea at mid swath. This is due to the fact that at those incidence angles, s corresponding to ice are also very close to the wind model. Also, the s "Isotropy" criterion does not prove to be very decisive at any incidence angle. This confirms the deductions made in Section 2.
When several sources of information are available, they can be combined to reduce imprecision and uncertainty and increase completeness [9] A good review of some of the existing combination methods is to be found in [9] It only makes sense to combine the "best" sources of information, hence we will combine the "Distance to ice model" and the "Derivative of sigma" criteria. We will compare two combination methods: the mean operator and the symmetrical associative sums. The mean operator simply outputs the mean of its inputs, i.e.
In this paper, we presented a Neural-Network based method to compare the various ice detection criteria found in the literature. We also showed that by combining the two best criteria, a yet better performance was achieved. This performance was obtained in a scenario where single observations are considered, i.e., excluding any kind of time-averaging. ACKNOWLEDGMENTS
This work was performed under European Space Agency (ESA) contracts. We would like to thank the ESA for the use of data and the provision of industry-confidential information. This work would not have been possible without the help of Pascal Lecomte and Raffaele Crapolicchio from the European Space Agency/ESRIN.
P. Lecomte,.Wind scatterometer processing requirements for ice detection,. Tech. Rep. ER-SA-ESA-SY- 1122, ESA ESRIN, Sept. 1993. ERS-2 A8. [ Links ] [1] A. Cavani'e, F. Gohin, Y. Quilfen, & P. Lecomte,. Identification of sea ice zones using the AMI-wind: Physical bases and applications to the FDP and CERSAT processing chains,. in Proc. Second ERS-1 Validation Workshop, pp. 1009.1012, ESTEC, (Hamburg, Germany), Oct. 1993. [ Links ] [2] F. Gohin & A. Cavani'e,.Some active and passive microwave signatures of antarctic sea ice from midwinter to spring 1991,. International Journal of Remote Sensing 16(11), pp. 2031.2054, 1995. [ [ Links ]3] X. Neyt, P. Pettiaux, & M. Acheroy, .Scatterometer ground processing review for gyro-less operations,. in Proceedings of SPIE: Remote Sensing of the Ocean and Sea Ice 2002, 4880, 2002. [ Links ] [4] F. Gohin & C. Maroni, .ERS scatterometer polar sea ice grids: User manual,. Tech. Rep. C2-MUT-W- 03-IF, IFREMER, Mar. 1998. [ Links ] [5] A. Cavani'e, .An empirical C-band backscatter model over arctic sea ice from ERS-1 AMI-wind data,. in Proc. of a Joint ESA-Eumetsat Workshop on Emerging Scatterometer Applications, pp. 99.106, ESTEC, ESTEC, (Noordwijk, The Netherlands), Nov. 1995. [ Links ] [6] S. de Haan & A. Stoffelen, .Ice discrimination using ERS scatterometer,. tech. rep., KNMI, Sept. 2001. [ Links ] [7] C. Bishop, Neural Networks for Pattern Recognition, Carendon Press, 1995. [ Links ] [8] I. Bloch, .Information combination operators for data fusion: A comparative review with classification,. IEEE Transactions on Systems, Man and Cybernetics 26, pp. 52.67, Jan. 1996. [ Links ] [9] |