unsupervised classification with R

unsupervised classification with R

m

January 29, 2016

Here we see three simple ways to perform an unsupervised classification on a raster dataset in R. I will show these approaches, but first we need to load the relevant packages and the actual data. You could use the Landsat data used in the “Remote Sensing and GIS for Ecologists” book which can be downloaded here.

library("raster")  
library("cluster")
library("randomForest")

# loading the layerstack  
# here we use a subset of the Landsat dataset from "Remote Sensing and GIS for Ecologists" 
image <- stack("path/to/raster")
plotRGB(image, r=3,g=2,b=1,stretch="hist")

RGBimage

Now we will prepare the data for the classifications. First we convert the raster data in a matrix, then we remove the NA-values.

## returns the values of the raster dataset and write them in a matrix. 
v <- getValues(image)
i <- which(!is.na(v))
v <- na.omit(v)

The first classification method is the well-known k-means method. It separates n observations into  k clusters. Each observation belongs to the cluster with the nearest mean.

## kmeans classification 
E <- kmeans(v, 12, iter.max = 100, nstart = 10)
kmeans_raster <- raster(image)
kmeans_raster[i] <- E$cluster
plot(kmeans_raster)

Kmeans

The second classification method is called clara (Clustering for Large Applications). It work by clustering only a sample of the dataset and then assigns all object in the dataset to the clusters.

## clara classification 
clus <- clara(v,12,samples=500,metric="manhattan",pamLike=T)
clara_raster <- raster(image)
clara_raster[i] <- clus$clustering
plot(clara_raster)

clara

The third method uses a random Forest model to calculate proximity values. These values were clustered using k-means. The clusters are used to train another random Forest model for classification.

## unsupervised randomForest classification using kmeans
vx<-v[sample(nrow(v), 500),]
rf = randomForest(vx)
rf_prox <- randomForest(vx,ntree = 1000, proximity = TRUE)$proximity

E_rf <- kmeans(rf_prox, 12, iter.max = 100, nstart = 10)
rf <- randomForest(vx,as.factor(E_rf$cluster),ntree = 500)
rf_raster<- predict(image,rf)
plot(rf_raster)

randomForest

The three classifications are stacked into one layerstack and plotted for comparison.

class_stack <- stack(kmeans_raster,clara_raster,rf_raster)
names(class_stack) <- c("kmeans","clara","randomForest")

plot(class_stack)

Comparing the three classifications:

Looking at the different classifications we notice, that the kmeans and clara classifications have only minor differences.
The randomForest classification shows a different image.

 

want to read more about R and classifications? check out this book:

you may also like:

EO4Cam Meeting in Oberpfaffenhofen

EO4Cam Meeting in Oberpfaffenhofen

On 5 June 2025, representatives from the Bavarian State Ministry of the Environment and Consumer Protection (StMUV), the Bavarian Environment Agency (LfU), and the German Aerospace Center (DLR) met at DLR Oberpfaffenhofen as part of the EO4CAM (Earth Observation...

EORC at the International Africa Festival 2025 in Würzburg

EORC at the International Africa Festival 2025 in Würzburg

If you know Würzburg, you certainly know the International Africa Festival, Europe's largest and oldest festival for African music and culture. For 15 years in a row now, the university tent has been an integral part of the festival. This is where the...

EO4CAM meeting at LfU in Augburg on grassland in Bavaria

EO4CAM meeting at LfU in Augburg on grassland in Bavaria

As part of the EO4CAM project (Earth Observation Laboratory for Climate Adaption and Mitigation), representatives of the Bavarian Environment Agency (LfU), the German Aerospace Centre (DLR), and the Earth Observation Research Cluster (EORC) met at the LfU in...

New R Package Enhances UAS Research and Planning

New R Package Enhances UAS Research and Planning

We’re excited to share the development of a new R package created by our PhD student, Antonio Castaneda Gomez, whose contributions to Uncrewed Aerial Systems (UAS) research continue to impress. Known as the brain behind many of our UAS data collection...

Radio Bavaria BR2 covered our activities at the Africa-Festival

Radio Bavaria BR2 covered our activities at the Africa-Festival

Once again, our team proudly took part in the International Africa Festival in Würzburg, continuing our active participation within the University of Würzburg's exhibition—a tradition we’ve upheld for many years. This year’s event highlighted the ongoing commitment of...