Home > GRASS, ushcn > Three's Company

Three's Company

2010 October 27

Yesterday I posted the simple unsupervised cluster analysis of USHCN station 051294, Canon City, CO. It neatly divided into two classes which seemed to make a good match with the visible appearance of man-made features -v- vegetated landscape with a few bare ground spots thrown in to the man-made side of the classification.

But it doesn’t always work so easily. Before I tried USHCN 051294, I did a test run on station picked at random, USHCN 199316, which is West Medway, Massachusetts.

USHCN 199316

Dividing this station into two classifications, similar to the exercise yesterday, and we see that many areas within the natural vegetation are marked the same as obviously man-made areas. Indeed, this image is divided into almost equally into two classes: 51.37% and 48.63%

USHCN 199316 class 2

So we take the clustering code we looked at yesterday, which divided the image into two classes, and divide it into 3 classes instead.

i.cluster group=51294_18 subgroup=51294_18 classes=3 sigfile=51294_18_sig.txt reportfile=51294_18_rpt.txt
i.maxlik group=51294_18 subgroup=51294_18 sigfile=51294_18_sig.txt class=51294_18c_class reject=51294_18c_reject

USHCN 199316 class 3

Taking a quick look and we see that it looks much more like we expect for a natural -v- man-made division. Most of the man-made stuff is in class 3. But it’s not perfect. The house in the lower left corner is mostly classed “natural” while the surrounding lawn is classed “man-made.” Note: the unsupervised classification is NOT dividing things by any predefined ‘natural’ or ‘man-made’ clusters – its just that the cluster analysis of the location of certain colors tends to lump MOST man-made landscapes into a single class – given the right number of classes.

USHCN 199316 class 3

Running the stats on the 3×3 smoothed version, we get these three classification coverages:
1 32.92%
2 52.70%
3 14.38%

Which seems about right.

USHCN 199316 class mask

A confounding factor in the above is the presence of dark shingles and shaded trees which works against a clean separation of classes.

  1. 2010 November 9 at 1:50 am

    None of the pictures seem to work for me (I just get the “could not load picture” box).

  2. Steven Mosher
    2010 November 9 at 11:35 am

    me either Ron.

    I was also thinking of looking at doing some edge detection to see if we count the number of buildings in an area.

    I’m busy tidying up some of the other work and will be posting on it when time permits.

    Thanks for all your help

  3. 2010 November 27 at 11:34 am

    Where have you been Ron? I was particularly enjoying the Climate Blog and News recaps you were doing regularly for a while there…can we expect them back sometime in the near future?

  1. No trackbacks yet.
Comments are closed.