Unsupervised Classification

Automatically groups pixels with similar multiband values without labelled training data. The result is a cluster map that must be interpreted before it becomes a semantic land-cover map.

Remote Sensing CoreEnd-User DocumentationTheory + Formula + Parameters

1. What This Feature Does

Automatically groups pixels with similar multiband values without labelled training data. The result is a cluster map that must be interpreted before it becomes a semantic land-cover map.

Cloud workflow: choose data → configure the scientific method → run → review the resulting layer/report. No programming is required.

2. Recommended Workflow

Select informative predictor bands. ↓ Keep standardization on when numeric scales differ. ↓ Choose a clustering algorithm. ↓ Set the class count or adaptive bounds. ↓ Run the clustering. ↓ Interpret and, if necessary, merge/rename spectral clusters.

3. Theory, Methods & Equations

K-Means

Minimizes within-cluster squared distance to K centroids.

J=Σ||xi−μc(i)||²
μk=(1/|Ck|)Σxi
Best used when: a fast general baseline with a fixed class count is desired.

ISODATA

Adaptive centroid clustering that can remove small clusters, merge nearby centers, and split high-variance clusters.

Merge if ||μi−μj|| < Tmerge
Split when max(σk,b) ≥ Tsplit
Best used when: the natural number of clusters is uncertain.

Self-Organizing Map (SOM)

Learns a 2-D topology of representative feature vectors, then groups map units into output classes.

BMU(x)=argminj||x−wj||²
hbj=exp(−dgrid²/(2σ²))
Best used when: complex spectral topology is worth exploring.

Gaussian Mixture Model (GMM)

Models data as a mixture of Gaussian components with configurable covariance.

p(x)=Σπk𝒩(x|μkk)
Best used when: clusters overlap or have elliptical covariance structure.

Fuzzy C-Means

Assigns soft membership to all clusters during fitting, then produces a class map from the strongest membership.

uik=1/Σj(dik/dij)2/(m−1)
ck=Σuikmxi/Σuikm
Best used when: mixed pixels and gradual spectral transitions are common.

BIRCH

Builds compact Clustering Feature summaries in a hierarchy, then performs global grouping.

CF=(N,LS,SS)
μ=LS/N
Best used when: efficient clustering of large feature samples is needed.

4. Input Data

InputTypeRequirementDescription
Input RasterrasterRequired

5. Parameters Available in the Application

ParameterDefaultChoices / RangeHow to Use It
Classification AlgorithmK-Means, ISODATA, Self-Organizing Map (SOM), Gaussian Mixture Model (GMM), Fuzzy C-Means, BIRCHChoose the clustering method.
Predictor BandsAll bandsSelect raster bands used to form each pixel's feature vector.
Number of Classes52 to 255Requested classes; for ISODATA this is the initial class count.
Training Pixel Sample100000100 to —Maximum valid pixels sampled to learn the cluster structure.
Standardize Predictor BandsOnRecommended when bands or stacked variables have different numeric ranges.
Prediction Batch Size2500001000 to 2000000Controls memory used while assigning the final raster classes.
Random State42Keep fixed to make comparative runs more reproducible.
Output Nodata00 to 32767Value reserved for invalid / nodata output pixels.
K-Means Mini-Batch Size409664 to 262144
Shown for: kmeans
K-Means Maximum Iterations1001 to 2000
Shown for: kmeans
K-Means Initializations31 to 50
Shown for: kmeans
ISODATA Minimum Classes22 to 255
Shown for: isodata
ISODATA Maximum Classes122 to 255
Shown for: isodata
ISODATA Maximum Iterations301 to 500
Shown for: isodata
ISODATA Minimum Pixels per Cluster501 to 1000000
Shown for: isodata
ISODATA Merge Distance Threshold0.50 to —Small distances encourage merging of similar clusters.
Shown for: isodata
ISODATA Split Standard Deviation Threshold1.00 to —High within-cluster variation can trigger splitting.
Shown for: isodata
ISODATA Convergence Tolerance0.011e-07 to —
Shown for: isodata
SOM Grid Rows51 to 64
Shown for: som
SOM Grid Columns51 to 64
Shown for: som
SOM Epochs301 to 1000
Shown for: som
SOM Initial Learning Rate0.51e-06 to 2.0
Shown for: som
SOM Initial Neighborhood Sigma2.00.1 to 100
Shown for: som
SOM Mini-Batch Size102416 to 65536
Shown for: som
GMM Covariance TypeFullFull, Tied, Diagonal, Spherical
Shown for: gmm
GMM Maximum Iterations2001 to 5000
Shown for: gmm
GMM Convergence Tolerance0.0011e-08 to —
Shown for: gmm
GMM Covariance Regularization1e-060 to —
Shown for: gmm
GMM Initializations11 to 50
Shown for: gmm
Fuzzy C-Means Fuzziness (m)2.01.0001 to 10
Shown for: fuzzy_cmeans
Fuzzy C-Means Maximum Iterations1501 to 5000
Shown for: fuzzy_cmeans
Fuzzy C-Means Convergence Tolerance0.00011e-08 to —
Shown for: fuzzy_cmeans
Fuzzy C-Means Initializations11 to 20
Shown for: fuzzy_cmeans
BIRCH Threshold0.51e-06 to —
Shown for: birch
BIRCH Branching Factor502 to 10000
Shown for: birch

6. Output & Interpretation

The result should be interpreted according to the selected method and the physical meaning of the input data. Preserve source units, coordinate reference information, NoData meaning, acquisition date, and preprocessing level when comparing results.

Scientific interpretation: an algorithm can produce a numerically valid result even when the input data are unsuitable. Always check masks, units, sensor characteristics, spatial resolution and reference data.
Important: Cluster IDs are arbitrary. Class 1 in one run does not necessarily represent the same surface type in another run.

7. Best Practices

  • Use analysis-ready inputs and remove invalid/cloud/noise artifacts that are not part of the target phenomenon.
  • Choose parameters from the physical scale of the data, not only from visual appearance.
  • Keep categorical and continuous rasters conceptually separate when selecting interpolation/resampling methods.
  • Compare the result with the source image and independent reference information.
  • Document the settings used when results will be compared across dates, sensors or study areas.