The standard isodata clustering algorithm [Richards] was chosen as the unsupervised classification method for this dataset. The algorithm starts by randomly selecting cluster centers in the multidimensional input data space. Each pixel is then grouped into a candidate cluster based on the minimization of a distance function between that pixel and the cluster centers. After each iteration, the cluster means are updated, and clusters are possibly spilt or merged depending on the size and spread of the data points in the clusters.