Page updated: April 12, 2020
Author: Curtis Mobley

Thematic Mapping

Thematic mapping refers to the determination and display of a particular type of information (the theme). In terrestrial and oceanic remote sensing, a common theme is the type of surface material. On land, a thematic map might display the land areas covered by forest, grassland, water, crops, bare soil, pavement, etc. In shallow waters, the thematic map might distinguish bottom areas covered by mud, sand, rock, sea grass, coral, etc. Much work has been done recently on mapping bathymetry, bottom type, and water IOPs as extracted from hyperspectral imagery. This page compares the supervised classification technique used for terrestrial thematic mapping with spectrum matching techniques (e.g., Mobley et al. (2005); Dekker et al. (2011)) for shallow-water mapping of bottom type.

The simultaneous retrieval of bathymetry, bottom classification, and water IOPs is a much more difficult task than traditional thematic mapping to determine land surface type, as used in terrestrial remote sensing. In terrestrial thematic mapping, only the type of land surface must be deduced from an atmosphericly corrected image spectrum; there are no confounding influences by water IOPs and depth. We will see that terrestrial techniques for supervised classification are not well suited to the oceanic problem because of the additional complications of bottom depth and water optical properties, neither of which are present in terrestrial remote sensing.

Supervised classification

In supervised classification the object is to associate a given image spectrum with one of several pre-determined classes of spectra. In terrestrial remote sensing these classes are typically defined as soil, grass, trees, water, pavement, etc. A thematic map of earth surface features is then generated by classifying the spectrum from each image pixel into one of the pre-determined classes.

One approach to supervised classification is to compute the mean spectrum for each class and a corresponding covariance matrix that defines the ”size” of each class of spectra about its mean. The image spectrum is then compared only with the mean spectrum and size for each class, and the image spectrum is statistically associated with the class it is most likely to belong to according to some metric for distance between the image and mean spectra and user-specified assumptions about the statistical properties of the class members.

This page considers the terrestrial and oceanic problems in more detail and shows that the standard terrestrial thematic mapping methodology based on supervised classification is not easily applied to the ocean remote sensing problem.

Covariance and correlation matrices

Consider a collection of N remote sensing reflectance spectra Rrs, each with K wavelengths, which we denote by Rn(λk),n = 1,...,N and k = 1,...,K (dropping the rs subscript on Rrs for convenience). The spectra can be regarded as column vectors:

Rn = Rn(λ1) Rn(λ2) Rn(λK) = [Rn(λ1),Rn(λ2),,Rn(λK)]T , (1)

where bold type indicates a vector or matrix, and superscript T indicates transpose. In the spectrum matching technique described previously, these spectra are the database spectra, N is usually 105 or more, and K would be 75 for spectra from 380 to 750 nm with 5 nm resolution. Let

I = [R(λ1),R(λ2),,R(λK)]T

be the image spectrum that is to be classified.

Now consider subsets of the entire database that define various classes of spectra. To be specific in the illustrative computations below, we chose four classes of spectra: Rrs for 10 sand and sediment spectra seen through 0.01 m of water, 10 coral spectra seen through 0.01 m of water, and the same sand and coral spectra seen through 10 m of the same water. The water IOPs were based on measurements of the very clear water in the Bahamas. The sand and sediment spectra range from clean ooid sand to heavily biofilmed, darker sand. The coral spectra are different species of corals. Figure 1 shows the individual spectra in these four classes. To minimize the array sizes for the printout of Table 1 below, we subsampled the spectra to wavelengths of 400, 450, ..., 650, 700 nm, so that K = 7. The subsampled spectra are shown in Fig. 2.


Figure 1: Fig. 1. Example database Rrs spectra defining the four classes; each class has 10 spectra.


Figure 2: Fig. 2. The spectra of Fig. 1 as re-sampled at 50 nm intervals for use in the illustrative computations (light lines). The heavy lines are the class mean spectra.

These spectra are obviously correlated in wavelength. The amount of correlation between one wavelength and another is quantified by the covariance and correlation matrices, which are computed as follows. Let m = 1,...,M label the class, with M being the total number of classes (here 4). Class m contains Nm spectra (here, Nm = 10 for each class). Then the mean or average spectrum for each class is defined by

R¯m(λi) = 1 Nm n=1Nm Rn(λi), (2)

where the sum is over the spectra belonging to class m. In vector notation this is

R¯m = 1 Nm n=1Nm Rn. (3)

The mean spectra for the example four classes are shown by the heavy lines in Fig. 2.

The elements of the K × K class covariance matrices Σm are defined by

Σm(i,j) = 1 Nm 1 n=1Nm [Rn(λi) R¯m(λi)][Rn(λj) R¯m(λj)]. (4)

Σm(i,j) expresses the covariance of the class spectra at wavelength λi with λj; Σm(i,i) is the variance of the class spectra at λi. For remote-sensing reflectance spectra Rrs with units of sr1, the units of Σm(i,j) are sr2. If we arrange the spectrum column vectors for class m in a K × Nm matrix with the class mean removed,

R(m) = R1(λ1) R¯m(λ1) RNm(λ1) R¯m(λ1) R1(λK) R¯m(λK)RNm(λK) R¯m(λK) , (5)

then the covariance matrix for class m can be compactly written as

Σm = 1 Nm 1R(m)R(m)T . (6)

The elements of the K × K correlation matrix ρm for class m are defined from the class covariance matrix Σm by

ρm(i,j) = Σm(i,j) Σm (i, i)Σm (j, j). (7)

Table 1 shows the class covariance and correlation matrices computed by these equations for the four classes of spectra shown in Fig. 2.

Table 1. Covariance and correlation matrices for the four classes of spectra seen in Fig. 2. Wavelength 1 (400 nm) is at the upper left and wavelength 7 (700 nm) is at the lower right of each array. Units for Σ are sr2; ρ is non-dimensional.

Σ(sandat0.1m) = 2.462e43.049e43.612e43.797e44.141e43.994e43.625e43.049e4 4.547e4 5.361e4 5.447e4 5.712e4 5.639e4 4.569e4 3.612e4 5.361e46.338e46.462e46.787e46.692e45.449e43.797e4 5.447e4 6.462e4 6.658e4 7.046e4 6.912e4 5.780e4 4.141e4 5.712e46.787e47.046e47.546e47.364e46.340e43.994e4 5.639e4 6.692e4 6.912e4 7.364e4 7.230e4 6.179e4 3.625e4 4.569e45.449e45.780e46.340e46.179e45.856e4

ρ(sandat0.1m) = 1.0000.9110.9140.9380.9610.9470.9550.911 1.000 0.999 0.990 0.975 0.983 0.885 0.914 0.9991.0000.9950.9810.9890.8940.938 0.990 0.995 1.000 0.994 0.996 0.926 0.961 0.9750.9810.9941.0000.9970.9540.947 0.983 0.989 0.996 0.997 1.000 0.950 0.955 0.8850.8940.9260.9540.9501.000

Σ(coralat0.1m) = 5.793e54.952e56.378e51.143e41.389e41.129e41.552e44.952e5 4.442e5 6.222e5 1.080e4 1.299e4 1.059e4 1.417e4 6.378e5 6.222e51.070e41.712e42.007e41.655e42.057e41.143e4 1.080e4 1.712e4 3.010e4 3.578e4 2.961e4 3.775e4 1.389e4 1.299e42.007e43.578e44.432e43.774e44.950e41.129e4 1.059e4 1.655e4 2.961e4 3.774e4 3.300e4 4.306e4 1.552e4 1.417e42.057e43.775e44.950e44.306e46.132e4

ρ(coralat0.1m) = 1.0000.9760.8100.8650.8670.8170.8230.976 1.000 0.903 0.934 0.926 0.875 0.859 0.810 0.9031.0000.9540.9220.8810.8030.865 0.934 0.954 1.000 0.980 0.940 0.879 0.867 0.9260.9220.9801.0000.9870.9490.817 0.875 0.881 0.940 0.987 1.000 0.957 0.823 0.8590.8030.8790.9490.9571.000

Σ(sandat10m) = 5.792e7 2.174e6 3.341e6 1.880e6 4.235e8 2.783e9 2.886e122.174e6 1.006e5 1.542e5 8.357e6 1.805e7 1.214e8 1.047e11 3.341e6 1.542e5 2.369e5 1.288e5 2.785e7 1.872e8 1.645e111.880e6 8.357e6 1.288e5 7.086e6 1.545e7 1.033e8 9.641e12 4.235e8 1.805e7 2.785e7 1.545e7 3.409e9 2.267e102.151e132.783e9 1.214e8 1.872e8 1.033e8 2.267e10 1.517e11 1.457e14 2.886e121.047e111.645e119.641e122.151e131.457e142.776e17

ρ(sandat10m) = 1.0000.9000.9020.9280.9530.9390.7200.900 1.000 0.999 0.990 0.975 0.983 0.626 0.902 0.9991.0000.9940.9800.9870.6410.928 0.990 0.994 1.000 0.994 0.996 0.687 0.953 0.9750.9800.9941.0000.9970.6990.939 0.983 0.987 0.996 0.997 1.000 0.710 0.720 0.6260.6410.6870.6990.7101.000

Σ(coralat10m) = 1.956e7 5.520e7 9.484e7 9.105e7 2.255e8 1.305e9 1.714e125.520e7 1.638e6 3.068e6 2.868e6 7.055e8 4.100e9 4.858e12 9.484e7 3.068e6 7.054e6 6.135e6 1.484e7 8.731e9 8.282e129.105e7 2.868e6 6.135e6 5.845e6 1.424e7 8.393e9 9.001e12 2.255e8 7.055e8 1.484e7 1.424e7 3.620e9 2.197e102.322e131.305e9 4.100e9 8.731e9 8.393e9 2.197e10 1.372e11 1.359e14 1.714e124.858e128.282e129.001e122.322e131.359e142.221e17

ρ(coralat10m) = 1.0000.9750.8070.8520.8470.7970.8230.975 1.000 0.903 0.927 0.916 0.865 0.805 0.807 0.9031.0000.9550.9280.8880.6620.852 0.927 0.955 1.000 0.979 0.937 0.790 0.847 0.9160.9280.9791.0000.9860.8190.797 0.865 0.888 0.937 0.986 1.000 0.778 0.823 0.8050.6620.7900.8190.7781.000

These specific examples make it clear that

  • For a given class, Rrs at one wavelength is highly correlated with Rrs at another wavelength, as expected.
  • The covariance and correlation matrices are different for each class. These matrices depend not only on bottom type (sand vs coral) but also on bottom depth (and water IOPs, not explicitly shown here). In other words, the wavelength covariances carry information about both bottom type and water depth and IOPs.

Spectrum matching vs. statistical classification

One metric for comparing two spectra is the simple Euclidean metric, which measures the squared distance (in units of sr2) between an image spectrum I and each Rn in the database:

DE2(n) = i=1K[I(λ i) Rn(λi)]2 = [I R n]T [I R n] (8)

The spectrum Rn giving the minimum distance DE2(n) of all N database spectra determines the closest match to the image spectrum I. Note that this is not a statistical estimate in the sense that no probability model is involved. Note also that the image spectrum is being compared with every spectrum in the database, not just with pre-defined class mean spectra.

In traditional thematic classification, an image spectrum I is compared only with the mean spectrum and ”size” for each class, as expressed by the class mean R¯m and covariance Σm. Here ”size” is used in the sense that the variances and covariances in Σm are larger when the spread of Rrs spectra is greater. Inspect, for example, the elements of Σm for the class of sand at 0.01 m compared to sand at 10 m, for which the spectra are all much closer together (especially at blue and red wavelengths) and thus have smaller covariances. The class covariance matrix defines the size of the ”swarm of points” surrounding the centroid (mean class spectrum) representing the class in K-dimensional Rrs space. The image spectrum is assigned to a particular class according to a statistical model (often based on the assumption of a multivariate normal distribution of the swarm of points) that determines the probability that the image spectrum belongs to a particular the swarm of points defining a given class. The class spectra (K-dimensional swarms of points) generally overlap, so that an unambiguous, non-probabilistic association of I with a given class is not possible.

In maximum likelihood estimation (MLE; see Richards and Jia (1996) for an excellent discussion of this whole business), the distance metric is

DMLE2(m) = ln|Σ m| + [I R¯m]T Σ m1[I R¯ m], (9)

where |Σm| denotes the determinant of Σm and Σm1 denotes the inverse. |Σm| and Σm1 are of course pre-computed for each class before doing the spectrum matching. The image spectrum I is assigned to the class m having the smallest value of DMLE2(m). Note that now the image spectrum is compared only with the class mean spectra R¯m. The assignment of the image spectrum to a particular class is based on its closeness to the class mean and the spread of the ”swarm of points” surrounding the mean, as described by the covariance matrix. This metric involves matrix multiplication, which is computationally expensive, but the number of classes is generally small, so in practice this may not be a problem.

It is often said that the incorporation of Σm into the distance metric ”removes the effect of correlations between wavelengths.” This interpretation of the effect of Σm relates to the fact that covariance matrices are the foundation of principle component analysis (PCA; see Preisendorfer (1988)). In PCA the original independent, physical variables (here, the wavelengths) are transformed to obtain a new set of (generally unphysical) independent variables for which the data are uncorrelated. This transformation can be viewed as a rotation of the axes of the original, physical data space (here the wavelength axes used for plots in K-dimensional space) to generate new (generally unphysical) axes for which the data are uncorrelated.

If the class covariances are equal (or assumed to be equal), then ln|Σm| is the same for each class and can be ignored. The MLE metric then reduces to the Mahalanobis distance metric,

DM2(m) = [I R¯ m]T Σ1[I R¯ m], (10)

where Σ is the common value of Σm. The image spectrum I is then assigned to the class m having the smallest value of DM2(m).

We have seen by the specific examples of Fig. 2 and Table 1 that the covariance matrices are different for different classes of the sort that are relevant for ocean-bottom remote sensing. Indeed, Table 1 shows that the elements of the Σm can change by orders of magnitude as a function of water depth. This inequality of the Σm for different classes precludes use of the Mahalanobis metric for classes as defined here. For the retrievals needed for shallow-water mapping of bottom type, MLE (or something else) would have to be used with a different covariance matrix for each class.

However, it is not at all clear how meaningful classes should be defined for simultaneous retrievals of bottom type, water column IOPs, and bottom depth. Should one class be ”sand spectra at 5.25 m depth with a particular set of water absorption, scattering, and backscatter spectra,” and another class be ”sand spectra at 5.25 m depth with the same absorption and scattering spectra but a different backscatter fraction,” and another class be ”sand spectra at 5.50 m with the first set of IOPs,” and then another class be ”sea grass spectra at 7.50 m with yet another set of IOPs,” and so on? If so, then the number of classes quickly becomes as large as the number of depths, IOP sets, and pre-chosen classes of bottom type (sand, coral, sea grass, etc.). A database generated as previously described easily could have hundreds or thousands of classes (a database often has 50-100 bottom depths, several dozen to several hundred sets of IOPs, and more than 100 bottom reflectance spectra). With such a large number of classes, the validity of doing traditional thematic mapping becomes uncertain, not to mention the additional computational costs involved with the matrix multiplications.

Clearly spectrum matching for shallow-water applications addresses a much more complicated problem than classic terrestrial thematic mapping, which corresponds to retrieval of bottom type if there were no water present, i.e. no simultaneous retrieval of depth and IOPs. Because of the greater complexity of the oceanographic retrieval problem, and because of the difficulty in defining meaningful classes, shallow-water spectrum matching does not use statistical classification techniques such as MLE. The spectrum matching approach of Mobley et al. (2005) does not compare an image spectrum to a class mean spectrum. In that technique, an image spectrum is compared to every spectrum in a database to find the closest match by the chosen (Euclidean or some other) metric, which is appropriate in this case. In a manner of speaking, each database Rrs spectrum is a separate class corresponding to a particular depth, bottom reflectance spectrum, and set of IOPs. In such a situation (only one member in each class) the covariance matrix is undefined.

Moreover, for the present problem it is not even desirable to remove the effects of wavelength correlations, as can be done with the MLE or Mahalanobis metrics, because the wavelength correlations carry information that is critical to separating depth and IOPs effects from bottom type effects.

The spectrum-matching approach of Mobley et al. (2005) for shallow-water benthic mapping therefore avoids defining predetermined classes and finds the closest match from the entire database. This gives the highest possible resolution (in depth, bottom type, and water IOPs) of retrievals. This approach retrieves a particular bottom reflectance spectrum (which represents a particular bottom type), not just a generic bottom type such as sand or coral. If the user later wishes to group the particular spectra for the retrieved bottom types into broader classes such as corals vs. sediments, or to group the retrieved IOPs into low, medium, and high absorption bins, for example, then that is easily done from the full-resolution retrieval.

Comments for Thematic Mapping:

Loading Conversation