Startseite Consideration of Landsat-8 Spectral Band Combination in Typical Mediterranean Forest Classification in Halkidiki, Greece
Artikel Open Access

Consideration of Landsat-8 Spectral Band Combination in Typical Mediterranean Forest Classification in Halkidiki, Greece

  • Mohamed Elhag EMAIL logo
Veröffentlicht/Copyright: 13. Oktober 2017
Veröffentlichen auch Sie bei De Gruyter Brill

Abstract

Utilization of satellite data by means of remote sensing practices creates a wider window of opportunities to conduct robust outcomes Landsat OLI-8 data acquired over the peninsula of Halkidiki in Greece for the summer period of 2013, were analyzed to determine their utility to classify natural resources categories based on the determination of spectral bands combination. Consequently, spectral bands combination can be used to classify various categories based on their higher overall accuracy assessment. Spectral information contained in each of the utilized channels of Landsat OLI-8, statistics defining 12 cover type classes of interest were calculated and used as a basis for classification of the designated study area. Spectral bands combinations (4-5-6/ 3-4-5/ 3-4-5-6/ 4-5-6-7) fulfill the required accuracy for forest purpose as they conducted overall accuracies ranging between 7.85% and 96.29%. The spectral band combination 3-4-5-6-7 conducted the highest overall accuracy (98.15%). The poorest results were obtained by the analysis of spectral bands combination 4-6 (68.52%). Single spectral band 6 conducted the best overall results; however, spectral bands 5 and 6 seemed to be the most useful spectral bands combination for each category.

1 Introduction

The general term of image processing is to extract information either with no interference with the observed data (image preprocessing) and/or with adequate interference to improve the results concerning the observed phenomena [1, 2]. The most common practice in image processing is to recognize a pattern. Tou and Gonzalez [3] has defined pattern recognition as the categorization of input data into identifiable classes, via the extraction of significant features or attributes of the data from background irrelevant details [4, 5].

The first step in the analysis sequence involves the selection of data samples from which the computer may drive the necessary statistical parameters for the “training” of the classification algorithm [6, 7]. The training samples are based on “ground truth” information. This means that the used categories are defined based on the separability index [8, 9].

There are two major approaches, the unsupervised and the supervised model. The unsupervised model is useful when the spectral info is not identified [10]. An important criterion for clustering is the minimization of the sum of square error for determining the spectral class composition of the data, prior to detailed analysis by the method of supervised classification [2, 11]. Supervised classification assumes that each spectral class can be described by the probability distribution in the multispectral space; this will be a multivariate distribution with many variables as dimensions of space [12, 13].

Classification costs increase by the number of features which are used to describe pixel vectors in a multispectral space e.g. “with the number of spectral bands associated with a pixel” [14, 15]. Each class has a statistical spread associated with it, through the data point in each spectral band. As such probability distributions, can overlap, to greater or lesser extent; the class separability becomes a function of both the separation of the means and the statistical distribution of data points, within each class, for each dimension [8, 16].

Each of OLI-8 spectral bands can potentially contribute to the multispectral land cover classification. Although, there is a redundancy of information due to high inter-spectral band correlation [17]. Two basic approaches are used to identify subsets of spectral bands, which might be suitable classification of the multispectral data. These are separability analysis and evaluation of eigenvector and eigenvalue data derived from class statistics [18, 19]. There are several measures of separability available to predict best spectral bands combination for classification which are based on measurements of the statistical distance between spectral classes of interest [2022].

According to Gausman et al. [23], the wavelength of 0.68, 0.85, 1.65 and 2.20 µm is useful for monitoring vegetation. Lathrop et al. [24] reported that the most informative spectral intervals for monitoring for natural materials were 0.54-0.56, 0.66-0.68 and 0.78-0.82 µm. Many researchers are trying to evaluate the best single or combination of spectral bands for displaying, mapping and classifying forest lands [25].

Murtha and Watson [26] stated that a combination of spectral bands, primarily TM 5 and 6 are required for the interpretation and mapping forest clear cuts. Coggeshall and Hoffer [27] verified that at least one spectral band in the near infrared or middle infrared is necessary to accurately discriminate broadleaf and conifers forests. Spanner et al. [4] using principal component analysis found that the four optimal spectral bands of the Thematic Mapper Simulator (TMS) data for forest vegetation analysis, over a Northern Idaho study site, were the near infrared (TMS4), thermal infrared (TMS6), mid-infrared (TMS5) and red (TMS3). Latty and Hoffer [1], using average transformed divergence, found the highest separability between South Carolina forest classes in four spectral bands subset consisting of TM 1, –3, –4 and 5. They noted that such waveband selection results were “highly data and application dependent”. Duda et al. [28] pointed out that beyond a certain point the inclusion of additional features leads to worse rather than better performance. Nelson et al. [1] using stepwise discrimination analysis of TMS data, appeared to provide forest classification accuracies that are not significantly different from similar scholarly work [29, 30].

The main purpose of the study is to increase the knowledge about the use of satellite data in vegetation ecosystems. More specifically the objective is to determine the spectral bands combination that gives the highest overall accuracy that can be used for better classification of various categories of typical Mediterranean land cover.

2 Materials and Methods

2.1 Study Area

The study area is the peninsula Sithonia of the prefecture of Halkidiki (Figure 1). Study area selection was based on the fact that the designated study area belongs to a typical Mediterranean ecosystem. Therefore, most of the major forest cover types found in Greece are presented in the area [31]. Moreover, there are available ancillary data giving information about the vegetation, soil, and topography. Sithonia is situated on the middle of the three peninsulas of Halkidiki and occupies a place with longitude between 23° 36’ up to 24° 00’ E and latitude between 39° 56’ up to 40° 14’ N. The peninsula, with 43 km in length and 6 km in width in the northern part and 18 km in the center, is considered as a continuation of the mountain Holomonta. Its acreage is about 450 km2, half of which are covered by forests of Aleppo pine (Pinus halepensis), which settles its optimum development in Sithonia. The elevation is up to 823 m (hill Polielaios). The relief is gently looping, but there are places with an inclination of 50-60%. The peninsula has only small numerous water streams with seasonal activity, which create gully erosion phenomena. The pressure of the biotic factors of the area has affected the soil of Sithonia. The main rocks found are granites, metamorphic rocks of gneiss, sandstones, and alluvial rocks, semi-metamorphic rocks of phyllites and calc-schists and finally deposits of alluvial and sand-alluvial rocks [31, 32]. On the east side of Sithonia, where silicate rocks predominate, and the soils are acidic, shallow, with numerous rocks and littler fertility. On these soils, the stand of Aleppo pine of Stone pine (Pinus pinea) and Black pine (Pinus nigra) have been developed [33]. The understory is composed of acid-friendly bushes of Ericaceae (tree-heath and strawberry tree) and Cistaceae family. On non-silica soils, the maquis has been developed. Moreover, the development of Aleppo pine on such soils is quite better in contrast to the Stone pine [2, 34].

Figure 1 Location of the study area
Figure 1

Location of the study area

2.2 Data Acquisition

Landsat OLI-8 image was acquired on the 30th of April 2013. It has a nominal center path of 183 and row of 32 according to the Worldwide Reference System. During the processing of the image, Landsat OLI-8 band 1 and the thermal spectral bands were excluded. Spectral band 1 has been omitted because it is highly correlated with spectral band 2, and it’s affected by the atmospheric haze [16, 35].

Thermal spectral bands were excluded from further analysis because they give poor classification results due to their coarse spatial resolution and lack contrast [3, 36, 37]. Ancillary data in the form of thematic maps, orthophotographs, orthophoto maps and aerial photographs were obtained by the Hellenic Army Geographical Service, The Ministry for the Environment, Physical Planning and Public Works and from the Hellenic Forest Service. The vegetation and geological thematic maps (scale 1:500.000) used in this study have been developed by the Directorate of Forest Resources Development, Hellenic Forest Service and by the Institute of Geology and Mineral Exploration.

Field observations have been done during the spring of 2013. The most serious part of the field observations was the accurate recognition of the training sites in the field, in the analog photographs and afterward in the digital images.

2.3 Development of the classification scheme

The classification system proposed by Bonazountas et al. [38] has been used for this study. it has been properly adjusted to the Coordination of Information on the Environment “CORINE” Land Cover System (2006). The following classification system was required to establish a list of 12 categories suitable for the study as it illustrated in Figure 2.

Figure 2 The classification scheme
Figure 2

The classification scheme

2.4 Compilation of classification procedures

2.4.1 Unsupervised classification

An unsupervised classification is an analytical procedure based upon clustering, using different algorithms. The whole image of the study area has been segmented into 12 spectral categories. Various techniques have been used to get more detailed information. in this study we applied the Optimum Index Factor (OIF). The algorithm used to compute OIF for any subset of spectral bands is following Chavez et al. [39]:

OIF=k=1skj=1Abs(rj)(1)

Where:

sk = is the standard deviation from spectral band k.

rj = is the absolute value of the correlation coefficient between any of the two spectral bands being evaluated.

The spectral bands combination with the largest OIF has the most information, as measured by variance, with the least amount of duplication, as measured by correlation [40, 41].

2.4.2 Supervised classification

A supervised classification is described by the potential distribution of each class in the multispectral dimension. Therefore, ti’s will be a multivariable distribution with as many variables as the dimensions of the data [42]. A multidimensional normal distribution is described as a function of vector location in multispectral space by the following algorithm of Richards and Richards [12]:

p(x)=1(2π)N2||12exp(12(xm)t.1(xy)(2)

Where:

x = is a vector location is N-dimensional pixel space.

m = is the mean position of the spectral class.

∑= is the covariance matrix of the distribution, which describes its spread directionally in the pixel space.

If m and ∑ are known for every spectral class in the image, every pixel can be examined and categorized in the most likely class, based on the probabilities computed for the location. In this study, the maximum likelihood classifier has been used, because of its higher accuracy [43, 44].

To do this, an assumption is made that the cloud points forming the category training data is normally distributed [45, 46]. The position of the pixel points in a multispectral space can be described by vectors, whose components are the individual spectral responses in each spectral band. The mean position of the pixels in space is defined by the expected value of the pixel vector “x” according to Richards and Richards [12]:

m=E(x)=1kj=1kxj(3)

Where:

m = mean pixel vector.

Xj = individual pixel vectors of total number K.

E= expectation operator.

It is of value to have available a means by which the mean vector defines the average position of the pixels in multispectral vector space. This is the role of the covariance matrix which is defined as:

x=1k1j=1k(xjm)(xjm)t(4)

Where:

t = denote vector transpose.

The covariance matrix is one of the most important mathematical concepts in the analysis of multispectral remote sensing data [47]. So, if there is a correlation between the responses in a pair of spectral spectral bands, the corresponding off-diagonal element in the covariance matrix will be large by comparison to the diagonal terms. On the other hand, if there is a little correlation, the off-diagonal terms will be close to zero [48].

2.5 Accuracy assessment

The final step in the digital image analysis is the evaluation of the accuracy of the classification results. These results are expressed in tabular form are known as a confusion matrix. The statistical analysis presents the accuracy and the misfire of the classification based on the ground truth observations. A distinction is made between errors of omission and errors of commission. Errors of omission corresponded to those pixels belonging to the class of interest that the classifier has failed to recognize, whereas errors of commission are those that corresponded to pixels from other classes that the classifier has labeled as belonging to the class of interest [49]. The omission refers to columns of the confusion matrix, whereas the commission refers to rows.

The following diagram simplifies the methodological framework adopted in the current research study as it’s illustrated in Figure 3.

Figure 3 Schematic description of the followed procedures
Figure 3

Schematic description of the followed procedures

3 Results and discussion

Unsupervised and supervised procedures based on the spectral bands combinations have been used for the classification of twelve categories in Sithonia peninsula. The unsupervised classification showed low overall accuracy while the supervised showed acceptable overall accuracy for forest purposes. Spectral bands combinations were used for the best discrimination amongst the individual categories.

Once the training fields had been identified, there were grouped according to cover classes. The cover class groups of training fields were then individually clustered to resolve the cover classes into a set of spectral classes. The identified percentages of the quotient classified classes were demonstrated in Figure 4.

Figure 4 The contribution of the training fields in the classification
Figure 4

The contribution of the training fields in the classification

The unsupervised classification based on the Iterative Self-Organizing Data Analysis Technique (ISODATA) algorithm was performed using different spectral bands combinations. Table 1 shows the analysis of the classification results based on OIF that were produced. The evaluation of different spectral bands combinations is shown in Figures 5, 6, 7 and 8. In Figure 9, the best spectral bands combinations, based on OIF for the designated study area was reported. Maximum likelihood supervised classification was performed using the best spectral bands combinations based on OIF is shown in Figure 10.

Figure 5 Two spectral bands combination
Figure 5

Two spectral bands combination

Figure 6 Three spectral bands combination
Figure 6

Three spectral bands combination

Figure 7 Four spectral bands combination
Figure 7

Four spectral bands combination

Figure 8 Five spectral bands combination
Figure 8

Five spectral bands combination

Figure 9 Best spectral bands combination, based on OIF
Figure 9

Best spectral bands combination, based on OIF

Figure 10 Maximum likelihood classification
Figure 10

Maximum likelihood classification

Table 1

Matrix statistical separability for 3-5 spectral bands combination

CategoryPHPNMCMAM+RBRRAUAOLAGSAWA
PH0.0098.899.267.688.595.792.593.980.878.680.899.9
PN98.80.0075.397.598.699.799.710099.498.099.999.9
MC99.275.30.0097.398.698.699.399.999.298.799.9100.0
MA67.697.597.30.000.0094.289.495.467.667.497.8100.0
M+R88.598.698.680.796.396.374.383.475.776.995.6100.0
BR95.799.798.694.274.30.0095.798.891.996.499.4100.0
RA92.599.799.389.483.495.70.0085.368.890.296.6100.0
UA93.9100.099.995.483.498.885.30.0086.493.988.0100.0
OL80.899.498.677.375.791.968.886.40.0082.496.4100.0
AG78.698.098.767.476.996.490.293.982.40.0098.0100.0
SA97.099.999.097.895.699.496.688.096.498.00.00100.0
WA99.0100.0100.0100.0100.0100.0100.0100.0100.0100.0100.00.00
Sum of deviation0.20135.5107.5142.231.575.0108.333.434.5131.3119.533.2

The overall accuracy and the one for each individual category were calculated for each set of analyzed data. This was achieved by using test fields in the study area and random sampling procedure. In this phase, 270 test sites were visited and the data were entered in confusion matrices. This type of representation enables to estimate the omission and commission errors and the overall accuracy of each category of the study area. The results of the overall accuracy are shown in Figure 11.

Figure 11 Overall accuracy assessment of various spectral bands combinations
Figure 11

Overall accuracy assessment of various spectral bands combinations

The overall accuracy assessments ranged from 68.52% to 98.15% and were considered reasonable for forestry purposes by Congalton [49] and Congalton and Green [10]. When spectral bands combination 3-5 (Table 1) was used, the overall accuracy was the lowest (68.52%). The category of conifers was classified with an accuracy ranging from 62.86% (Pinus nigra) to 84.44% (Pinus halepensis). Also, there were misclassification errors between the categories of Pinus halepensis and Pinus nigra and mixed conifers, maquis and maquis + rocks, range and agriculture, sand and urban areas, range and maquis according to the statistical separability index. This combination gave the lowest producer’s accuracy for the categories of Pinus nigra, maquis + rocks and range (62.86, 46.67 and 54.84 respectively) amongst various spectral bands combinations. Most of the forest categories were classified with user’s accuracy of less than 80% as the category of maquis showed the lowest accuracy (46.67%).

The spectral bands combination 3-7 (Table 2) gave an overall accuracy of 71.11%. Only a little lower than the spectral bands combinations 3-4-5-6 (71.85%). In this combination, the categories of Pinus halepensis and maquis, Pinus nigra and mixed conifers, olive trees, and agriculture gave the highest misclassification errors with the lowest producer’s accuracy amongst the other spectral bands combinations. The user’s accuracy was not satisfactory for 50% of the categories, except the categories of Pinus halepensis (90.48%), maquis + rocks (85.00%), broadleaf’s (93.33%), range (80.00%), urban areas (94.44%) and water (100.00%).

Table 2

Matrix statistical separability for 3-7 spectral bands combination

CategoryPHPNMCMAM+RBRRAUAOLAGSAWA
PH0.0096.992.254.070.271.176.993.267.791.297.498.1
PN96.90.0077.596.298.499.399.8100.099.695.699.9100.0
MC92.277.50.0089.695.896.299.1100.098.493.199.898.9
MA54.096.296.20.0071.869.582.495.473.292.598.098.9
M+R70.298.498.471.80.0084.867.081.163.897.394.099.8
BR71.199.396.269.584.80.0090.798.483.893.598.899.2
RA76.999.899.182.481.190.70.0081.161.399.195.4100.0
UA93.2100.0100.095.467.098.481.10.0086.398.385.7100.0
OL67.799.698.473.281.183.861.386.30.0096.496.299.9
AG91.295.693.192.563.893.599.198.396.40.0098.597.2
SA97.499.999.898.094.098.995.485.796.298.50.00100.0
WA98.1100.098.998.899.899.2100.0100.099.997.2100.00.00
Sum of deviation191.937.060.4178.9176.5114.3147.380.6173.747.436.59.90

The spectral bands combination 4-6 (Table 3) gave the second lowest accuracy assessment. Misclassification errors appeared between the categories Pinus Halepensis and maquis, olive trees and agriculture, sand, and urban areas. The categories of the maquis, broadleaf’s and sand were classified with the lowest accuracies of the various spectral bands combinations. The user’s accuracy was satisfactory only for 42% of the categories (Pinus nigra, broadleaf’s, range, urban areas, and water).

Table 3

Matrix statistical separability for 4-6 spectral bands combination

CategoryPHPNMCMAM+RBRRAUAOLAGSAWA
PH0.0095.589.856.570.0100.086.693.177.492.495.799.7
PN96.00.0077.295.597.772.599.9100.099.596.499.8100.0
MC92.477.20.0089.895.282.399.5100.098.596.899.6100.0
MA56.595.589.80.0070.683.188.295.977.494.599.699.8
M+R70.097.799.570.60.0097.968.483.967.397.796.6100.0
BR82.398.093.572.580.80.0092.498.583.198.990.3100.0
RA86.699.999.588.268.492.40.0078.067.599.497.9100.0
UA93.1100.0100.095.983.998.078.00.0087.297.788.7100.0
OL77.499.598.577.467.393.567.587.20.0097.393.0100.0
AG92.496.496.894.597.780.899.497.797.30.0098.8100.0
SA95.799.899.696.690.398.988.781.893.098.90.00100.0
WA99.7100.0100.099.8100.00.0100.0100.0100.0100.0100.00.00
Sum of deviation15840.057.3162.6178.1102.0131.483.9152.029.957.00.60

The spectral bands combination 3-4-5 (Table 4) gave overall accuracy (88.52%) while the classification errors appeared between the categories of sand and urban areas, olive trees and agriculture. The category of broadleaf gave the highest producers accuracy (100.00%) as the same percentage appeared with spectral bands combination 4-5-6 and 3-4-5-6-7. The user’s accuracy was satisfactory for 42.00% pf the categories but it is noticeable that the categories of conifers gave user’s accuracy between 93.48% (Pinus halepensis) and 97.96% (mixed conifers), which is satisfactory for forest purposes [17].

Table 4

Matrix statistical separability for 3-4-5 spectral bands combination

CategoryPHPNMCMAM+RBRRAUAOLAGSAWA
PH0.0099.099.469.990.497.796.195.985.495.797.5100.0
PN99.00.0078.798.099.399.899.9100.099.798.2100.0100.0
MC99.478.70.0097.598.999.099.7100.099.299.5100.0100.0
MA69.998.097.50.0081.995.892.896.680.395.598.1100.0
M+R90.499.398.981.90.0096.880.786.177.998.096.3100.0
BR97.799.899.095.896.70.0096.199.293.3100.099.5100.0
RA96.199.999.792.880.796.10.0089.571.599.897.4100.0
UA95.9100.0100.096.686.199.289.50.0089.398.089.2100.0
OL85.499.799.280.377.993.371.589.30.0097.996.7100.0
AG95.798.299.595.598.0100.099.898.097.90.0099.4100.0
SA97.5100.0100.098.196.399.597.489.296.799.40.00100.0
WA100.0100.0100.0100.0100.0100.0100.0100.0100.0100.0100.00.00
Sum of deviation73.127.528.393.593.822.976.556.210818.125.90.00

The spectral bands combinations 4-5-6 (Table 5) gave an overall accuracy of 96.29%, a little lower than the spectral bands combinations 3-4-5-6-7 (98.15%). The categories of mixed conifers, maquis, maquis + rocks, broadleaf, urban areas, and sand gave the highest procedure’s accuracy among the various spectral bands combinations. Most (83.33%) of the categories gave user’s accuracy more than 89% which is satisfactory [6, 19].

Table 5

Matrix statistical separability for 4-5-6 spectral bands combination

CategoryPHPNMCMAM+RBRRAUAOLAGSAWA
PH0.0098.899.169.988.696.993.295.484.895.895.8100.0
PN98.80.0080.598.099.399.6100.0100.099.798.0100.0100.0
MC99.180.50.0097.298.998.999.8100.099.199.399.9100.0
MA69.998.097.20.0079.895.291.296.880.795.897.0100.0
M+R88.699.398.979.80.0096.270.985.776.897.994.1100.0
BR96.699.698.995.296.20.0096.699.292.499.999.1100.0
RA93.2100.099.891.270.996.60.0079.871.899.592.4100.0
UA95.4100.0100.096.885.799.279.80.0088.798.286.4100.0
OL84.899.799.180.776.892.471.888.70.0098.094.3100.0
AG95.898.099.895.897.999.999.598.298.00.0099.3100.0
SA95.8100.099.997.094.199.192.486.494.399.30.00100.0
WA100.0100.0100.0100.0100.0100.0100.0100.0100.0100.0100.00.00
Sum of deviation81.826.127.398.2111.625.9104.969.9113.718.341.70.00

The spectral bands combinations 3-4-5-6 (Table 6) gave almost the same overall accuracy as spectral bands combination 3-7 (71.85 and 71.11% respectively). Classification errors appeared between the categories of maquis and maquis + rocks, urban areas, and agriculture. The category urban areas gave the lowest producer’s accuracy amongst the various spectral bands combinations.

Table 6

Matrix statistical separability for 3-4-5-6 spectral bands combination

CategoryPHPNMCMAM+RBRRAUAOLAGSAWA
PH0.0099.399.270.990.897.997.696.286.096.098.6100.0
PN99.30.0081.298.499.499.7100.0100.099.898.6100.0100.0
MC99.481.20.0097.799.099.099.8100.099.299.5100.0100.0
MA70.998.497.70.0082.796.296.597.181.896.098.6100.0
M+R90.899.499.082.70.0096.891.187.678.798.296.8100.0
BR97.999.899.096.296.80.0097.299.493.5100.099.7100.0
RA97.6100.099.896.591.197.20.0093.185.699.998.7100.0
UA96.2100.0100.097.187.699.493.10.0091.098.392.6100.0
OL86.099.899.281.878.793.585.691.00.0098.197.9100.0
AG96.098.699.596.098.2100.099.998.398.10.0099.6100.0
SA98.6100.0100.098.696.899.798.792.697.799.600.0100.0
WA100.0100.0100.0100.0100.0100.0100.0100.0100.0100.0100.00.00
Sum of deviation67.423.525.184.279.020.540.544.788.515.817.60.00

The spectral bands combination 4-5-6-7 (Table 7) gave higher accuracy than the spectral bands combination 3-4-5-6 (75.18 and 71.85% respectively). Classification errors appeared between the categories of maquis and maquis + rocks, agriculture and olive trees.

Table 7

Matrix statistical separability for 4-5-6-7 spectral bands combination

CategoryPHPNMCMAM+RBRRAUAOLAGSAWA
PH0.0099.499.577.090.697.895.796.589.097.097.3100.0
PN99.40.0083.998.599.599.7100.0100.099.898.2100.0100.0
MC99.583.90.0097.598.999.199.9100.099.299.3100.0100.0
MA77.098.597.50.0080.295.893.797.983.196.197.9100.0
M+R90.699.598.980.20.0096.679.990.179.698.195.6100.0
BR97.899.799.195.896.60.0097.299.693.4100.099.3100.0
RA95.7100.099.993.779.997.20.0089.973.599.694.4100.0
UA96.5100.0100.09.890.199.689.90.0093.299.489.0100.0
OL89.099.899.283.179.693.473.593.20.0098.395.5100.0
AG97.098.299.396.198.1100.099.699.498.30.0099.6100.0
SA97.3100.0100.097.995.699.394.489.095.599.60.00100.0
WA100.0100.0100.0100.0100.0100.0100.0100.0100.0100.0100.00.00
Sum of deviation60.321.128.882.390.821.676.244.495.314.331.50.00

The spectral bands combination 3-4-5-6-7 (Table 8) gave the highest overall accuracy (98.15%). The categories of Pinus halepensis, Pinus nigra, mixed conifers, maquis, broadleaf’s, range, olive trees, agriculture, and sand gave the highest producer’s accuracy amongst the various spectral bands combinations.

Table 8

Matrix statistical separability for 3-4-5-6-7 spectral bands combination

CategoryPHPNMCMAM+RBRRAUAOLAGSAWA
PH0.0099.699.777.692.598.498.197.189.797.298.8100.0
PN99.60.0084.498.899.599.9100.0100.099.898.8100.0100.0
MC99.784.40.0097.999.099.299.9100.099.399.6100.0100.0
MA77.699.897.70.0083.296.697.098.183.896.299.0100.0
M+R92.599.599.083.20.0097.192.291.581.198.597.6100.0
BR98.499.999.296.697.10.0097.699.794.4100.099.8100.0
RA98.1100.099.997.092.297.60.0096.585.999.999.3100.0
UA97.1100.0100.098.191.599.796.50.0094.999.593.6100.0
OL89.799.899.383.881.194.485.994.90.0098.498.5100.0
AG97.298.899.696.298.5100.099.999.598.40.0099.7100.0
SA98.8100.0100.099.097.699.899.393.698.599.70.00100.0
WA100.0100.0100.0100.0100.0100.0100.0100.0100.0100.0100.00.00
Sum of deviation58.619.726.397.985.219.575.142.591.613.127.40.00

Results from the analysis indicated that all the three spectral regions are suitable for forest cover classification and mapping [13, 34]. The combination of the visible spectral bands with the near and middle infrared bands give more accurate mapping results [8, 24].

The visible range of the spectrum seemed separating forests from other categories, but differentiating results broadleaf and coniferous forests [1, 9]. The near and middle infrared give quite satisfactory results in separating the two forest categories from each other (Pinus halepensis and Pinus nigra), with some failures between the broadleaf’s stands and the range [16, 18] (Figure 12).

Figure 12 Classification accuracy assessment of different categories of various spectral bands combinations
Figure 12

Classification accuracy assessment of different categories of various spectral bands combinations

4 Conclusions

Some authors claim to achieve better results with supervised than an unsupervised classification [7, 50]. It is easier to delineate well-defined classes by appropriately supervising techniques, but this requires a certain familiarity with the region. For areas, as large and complex as this study area, an unsupervised method must be used initially to obtain an indication of the spectral variability of the various categories. Once the unsupervised classification gives a general knowledge of the study area, then the supervised classification will yield greater accuracy in separating forest types, which most of the time have reflectances values close to each other. By utilizing the visible, near and mid-infrared parts of the electromagnetic spectrum gave an overall accuracy more than 98.10% can be achieved. High classification accuracies ranging from 82.61% to 98.15% for the Pinus halepensis category resulted from the analysis of the various spectral bands combination. For the category of Pinus nigra, the classification overall accuracy ranged between 62.86% and 96.77% while the mixed conifers between 57.14% and 95.45%. The category of maquis can be classified and mapped with accuracies ranging from 72.22% to 100%; the category of maquis + rocks lead to misclassification results, due to the presence of different amount of rocks which affected the spectral behavior of this category. The broadleaf’s classified, with overall accuracy ranging between 52.63% and 100%. The range category can be classified with accuracies ranging from 54.84% to 96.15%; the lower limit was due to classification errors between this category and maquis category, as in Sithonia there are not improved pasture, but only degraded forested lands where maquis predominate. The olive trees and agriculture categories lead several times to classification errors, due to the similar spectral responce. The agricultural lands showed the lowest overall accuracy (40% in spectral bands combination 3-7). The category sand showed the lowest accuracy (40%) in spectral bands combination 3-5 and 4-6.This category was confused with the category of urban areas due to the course resolution of the image (mixed pixel)‥ Water, on the contrary, was classified with the highest accuracy (100%). Finally, it seems that, spectral bands 3-4-5-6-7 are the most suitable for forest inventories. It should be noticed that more detailed inventory during the development of the training sites should have been done because the classification proved the presence of broadleaf’s in ravines, where no training data were collected.

Acknowledgement

This article was funded by the Deanship of Scientific Research (DSR) at King Abdulaziz University, Jeddah. The authors, therefore, acknowledged with thanks DSR for technical and financial support.

The author declares no conflict of interest regarding the publication of this manuscript.

References

[1] Nelson, R.F., R. Latty, and G. Mott, Classifying northern forests using Thematic Mapper simulator data. 1984.Suche in Google Scholar

[2] Yuan, F., et al., Multi-level land cover mapping of the Twin Cities (Minnesota) metropolitan area with multi-seasonal Landsat TM/ETM+ data. Geocarto International, 2005. 20(2): p. 5-13.10.1080/10106040508542340Suche in Google Scholar

[3] Tou, J. and R. Gonzales, Pattern Recognition Principles. Addison-Weley. Reading, MA, 1974. 377.Suche in Google Scholar

[4] Spanner, M.A., J.A. Brass, and D.L. Peterson, Feature selection and information content of thematic mapper simulator data for a forested environment. 1983.Suche in Google Scholar

[5] Clinton, N., et al., Accuracy assessment measures for objectbased image segmentation goodness. Photogrammetric Engineering and remote sensing, 2010. 76(3): p. 289-299.10.14358/PERS.76.3.289Suche in Google Scholar

[6] Rogan, J., et al., Land-cover change monitoring with classification trees using Landsat TM and ancillary data. Photogrammetric Engineering & Remote Sensing, 2003. 69(7): p. 793-804.10.14358/PERS.69.7.793Suche in Google Scholar

[7] Elhag, M., A. Psilovikos, and M. Sakellariou, Detection of land cover changes for water recourses management using remote sensing data over the Nile Delta Region. Environment, Development and Sustainability, 2013. 15(5): p. 1189-1204.10.1007/s10668-013-9433-5Suche in Google Scholar

[8] Friedl, M.A., et al., MODIS Collection 5 global land cover: Algorithm refinements and characterization of new datasets. Remote sensing of Environment, 2010. 114(1): p. 168-182.10.1016/j.rse.2009.08.016Suche in Google Scholar

[9] Ghimire, B., J. Rogan, and J. Miller, Contextual land-cover classification: incorporating spatial dependence in land-cover classification models using random forests and the Getis statistic. Remote Sensing Letters, 2010. 1(1): p. 45-54.10.1080/01431160903252327Suche in Google Scholar

[10] Congalton, R.G. and K. Green, Assessing the accuracy of remotely sensed data: principles and practices. 2008: CRC press.10.1201/9781420055139Suche in Google Scholar

[11] Townshend, J.R., et al., Global characterization and monitoring of forest cover using Landsat data: opportunities and challenges. International Journal of Digital Earth, 2012. 5(5): p. 373-397.10.1080/17538947.2012.713190Suche in Google Scholar

[12] Richards, J.A. and J. Richards, Remote sensing digital image analysis. Vol. 3. 1999: Springer.10.1007/978-3-662-03978-6Suche in Google Scholar

[13] Gislason, P.O., J.A. Benediktsson, and J.R. Sveinsson, Random forests for land cover classification. Pattern Recognition Letters, 2006. 27(4): p. 294-300.10.1016/j.patrec.2005.08.011Suche in Google Scholar

[14] Betts, R.A., et al., Biogeophysical effects of land use on climate: Model simulations of radiative forcing and large-scale temperature change. Agricultural and forest meteorology, 2007. 142(2): p. 216-233.10.1016/j.agrformet.2006.08.021Suche in Google Scholar

[15] Chapman, D.S., et al., Random Forest characterization of upland vegetation and management burning from aerial imagery. Journal of Biogeography, 2010. 37(1): p. 37-46.10.1111/j.1365-2699.2009.02186.xSuche in Google Scholar

[16] Hansen, M.C., et al., High-resolution global maps of 21st-century forest cover change. science, 2013. 342(6160): p. 850-853.10.1126/science.1244693Suche in Google Scholar PubMed

[17] Gong, P., et al., Finer resolution observation and monitoring of global land cover: First mapping results with Landsat TM and ETM+ data. International Journal of Remote Sensing, 2013. 34(7): p. 2607-2654.10.1080/01431161.2012.748992Suche in Google Scholar

[18] Lunetta, R.S. and M.E. Balogh, Application of multi-temporal Landsat 5 TM imagery for wetland identification. Photogrammetric Engineering and Remote Sensing, 1999. 65(11): p. 1303-1310.Suche in Google Scholar

[19] Li, C., et al., Comparison of classification algorithms and training sample sizes in urban land classification with Landsat thematic mapper imagery. Remote Sensing, 2014. 6(2): p. 964-983.10.3390/rs6020964Suche in Google Scholar

[20] Davis, S.M., et al., Remote sensing: the quantitative approach. New York, McGraw-Hill International Book Co., 1978. 405 p., 1978. 1.Suche in Google Scholar

[21] Manolakis, D. and G. Shaw, Detection algorithms for hyperspectral imaging applications. IEEE signal processing magazine, 2002. 19(1): p. 29-43.10.21236/ADA399744Suche in Google Scholar

[22] Mallinis, G., et al., Forest parameters estimation in a European Mediterranean landscape using remotely sensed data. Forest Science, 2004. 50(4): p. 450-460.Suche in Google Scholar

[23] Gausman, H., et al., Reflectance discrimination of cotton and corn at four growth stages. Agronomy Journal, 1973. 65(2): p. 194-198.10.2134/agronj1973.00021962006500020003xSuche in Google Scholar

[24] Lathrop, R., T.M. Lillesand, and B.S. Yandell, An evaluation of Thematic Mapper data for forest cover mapping in northern Wisconsin. Proc. 11th Pecora Sym., Sioux Falls, SD, 1987: p. 386-393.Suche in Google Scholar

[25] Elhag, M., Evaluation of different soil salinity mapping using remote sensing techniques in arid ecosystems, Saudi Arabia. Journal of Sensors, 2016. 2016.10.1155/2016/7596175Suche in Google Scholar

[26] Murtha, P. and E. Watson. Mapping of forest clear-cutting, south Vancouver Island, from Landsat imagery. in Canadian Symposium on Remote Sensing, 3 rd, Edmonton, Alberta, Canada. 1976.Suche in Google Scholar

[27] Coggeshall, M.E. and R.M. Hoffer, Basic forest cover mapping using digitized remote sensor data and automated data processing techniques. 1973.Suche in Google Scholar

[28] Duda, R.O., P.E. Hart, and D.G. Stork, Pattern classification. Vol. 2. 1973: Wiley New York.Suche in Google Scholar

[29] Briem, G.J., J.A. Benediktsson, and J.R. Sveinsson, Multiple classifiers applied to multisource remote sensing data. IEEE transactions on geoscience and remote sensing, 2002. 40(10): p. 2291-2299.10.1109/TGRS.2002.802476Suche in Google Scholar

[30] Lippitt, C.D., et al., Mapping Selective Logging in Mixed Deciduous Forest. Photogrammetric Engineering & Remote Sensing, 2008. 74(10): p. 1201-1211.10.14358/PERS.74.10.1201Suche in Google Scholar

[31] Franklin, S. and M. Wulder, Remote sensing methods in medium spatial resolution satellite data land cover classification of large areas. Progress in Physical Geography, 2002. 26(2): p. 173-205.10.1191/0309133302pp332raSuche in Google Scholar

[32] Bonan, G.B., Forests and climate change: forcings, feedbacks, and the climate benefits of forests. science, 2008. 320(5882): p. 1444-1449.10.1126/science.1155121Suche in Google Scholar

[33] Steele, B.M., Combining multiple classifiers: An application using spatial and remotely sensed information for land cover type mapping. Remote sensing of environment, 2000. 74(3): p. 545-556.10.1016/S0034-4257(00)00145-0Suche in Google Scholar

[34] Rogan, J., et al., Mapping land-cover modifications over large areas: A comparison of machine learning algorithms. Remote Sensing of Environment, 2008. 112(5): p. 2272-2283.10.1016/j.rse.2007.10.004Suche in Google Scholar

[35] Huang, H., et al., Reduction of atmospheric and topographic effect on Landsat TM data for forest classification. International Journal of Remote Sensing, 2008. 29(19): p. 5623-5642.10.1080/01431160802082148Suche in Google Scholar

[36] Coppin, P.R. and M.E. Bauer, Digital change detection in forest ecosystems with remote sensing imagery. Remote sensing reviews, 1996. 13(3-4): p. 207-234.10.1080/02757259609532305Suche in Google Scholar

[37] Hopkins, P.F., A.L. Maclean, and T.M. Lillesand, Assessment of Thematic Mapper imagery for forestry applications under Lake States conditions. Photogrammetric Engineering and Remote Sensing, 1988. 54(1): p. 61-68.Suche in Google Scholar

[38] Bonazountas, M., et al., A decision support system for managing forest fire casualties. Journal of Environmental Management, 2007. 84(4): p. 412-418.10.1016/j.jenvman.2006.06.016Suche in Google Scholar

[39] Chavez Jr, P.S., G.L. Berlin, and L.B. Sowers, Statistical Method for Selecting Landsat MSS. Journal of applied photographic engineering, 1982. 8(1): p. 23-30.Suche in Google Scholar

[40] DeFries, R. and J.C.-W. Chan, Multiple criteria for evaluating machine learning algorithms for land cover classification from satellite data. Remote Sensing of Environment, 2000. 74(3): p. 503-515.10.1016/S0034-4257(00)00142-5Suche in Google Scholar

[41] Cutler, D.R., et al., Random forests for classification in ecology. Ecology, 2007. 88(11): p. 2783-2792.10.1890/07-0539.1Suche in Google Scholar PubMed

[42] Elhag, M. and J.A. Bahrawi, Conservational use of remote sensing techniques for a novel rainwater harvesting in arid environment. Environmental earth sciences, 2014. 72(12): p. 4995-5005.10.1007/s12665-014-3367-6Suche in Google Scholar

[43] Sesnie, S.E., et al., Integrating Landsat TM and SRTM-DEM derived variables with decision trees for habitat classification and change detection in complex neotropical environments. Remote Sensing of Environment, 2008. 112(5): p. 2145-2159.10.1016/j.rse.2007.08.025Suche in Google Scholar

[44] Foody, G.M., Sample size determination for image classification accuracy assessment and comparison. International Journal of Remote Sensing, 2009. 30(20): p. 5273-5291.10.1080/01431160903130937Suche in Google Scholar

[45] Petropoulos, G., et al., A global Bayesian sensitivity analysis of the 1d SimSphere soil–vegetation–atmospheric transfer (SVAT) model using Gaussian model emulation. Ecological Modelling, 2009. 220(19): p. 2427-2440.10.1016/j.ecolmodel.2009.06.006Suche in Google Scholar

[46] Elhag, M., Inconsistencies of SEBS Model Output Based on the Model Inputs: Global Sensitivity Contemplations. Journal of the Indian Society of Remote Sensing, 2016. 44(3): p. 435-442.10.1007/s12524-015-0502-0Suche in Google Scholar

[47] Psilovikos, A. and M. Elhag, Forecasting of remotely sensed daily evapotranspiration data over Nile Delta region, Egypt. Water Resources Management, 2013. 27(12): p. 4115-4130.10.1007/s11269-013-0368-2Suche in Google Scholar

[48] Elhag, M. and J.A. Bahrawi, Soil salinity mapping and hydrological drought indices assessment in arid environments based on remote sensing techniques. Geoscientific Instrumentation, Methods and Data Systems, 2017. 6(1): p. 149.10.5194/gi-6-149-2017Suche in Google Scholar

[49] Congalton, R.G., A review of assessing the accuracy of classifications of remotely sensed data. Remote sensing of environment, 1991. 37(1): p. 35-46.10.1016/0034-4257(91)90048-BSuche in Google Scholar

[50] Elhag, M. and J.A. Bahrawi, Realization of daily evapotranspiration in arid ecosystems based on remote sensing techniques. Geoscientific Instrumentation, Methods and Data Systems, 2017. 6(1): p. 141.10.5194/gi-2016-40-AC1Suche in Google Scholar

Received: 2017-3-16
Accepted: 2017-7-10
Published Online: 2017-10-13

© 2017 M. Elhag

This work is licensed under the Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 License.

Artikel in diesem Heft

  1. Regular Articles
  2. Two types of gabbroic xenoliths from rhyolite dominated Niijima volcano, northern part of Izu-Bonin arc: petrological and geochemical constraints
  3. Regular Articles
  4. CRSP, numerical results for an electrical resistivity array to detect underground cavities
  5. Regular Articles
  6. Magma evolution inside the 1631 Vesuvius magma chamber and eruption triggering
  7. Regular Articles
  8. Probabilities of Earthquake Occurrences along the Sumatra-Andaman Subduction Zone
  9. Regular Articles
  10. Modelling of carrying capacity in National Park - Fruška Gora (Serbia) case study
  11. Regular Articles
  12. Knickzone Extraction Tool (KET) – A new ArcGIS toolset for automatic extraction of knickzones from a DEM based on multi-scale stream gradients
  13. Regular Articles
  14. When the display matters: A multifaceted perspective on 3D geovisualizations
  15. Regular Articles
  16. Dependence of Gully Networks on Faults and Lineaments Networks, Case Study from Hronska Pahorkatina Hill Land
  17. Regular Articles
  18. Environmental Geochemistry of Geophagic Materials from Free State Province in South Africa
  19. Regular Articles
  20. Neotectonic interpretations and PS-InSAR monitoring of crustal deformations in the Fujian area of China
  21. Regular Articles
  22. Dual-shale-content method for total organic carbon content evaluation from wireline logs in organic shale
  23. Regular Articles
  24. The dolerite dyke swarm of Mongo, Guéra Massif (Chad, Central Africa): Geological setting, petrography and geochemistry
  25. Regular Articles
  26. Seismic data filtering using non-local means algorithm based on structure tensor
  27. Regular Articles
  28. Pore Distribution Characteristics of the Igneous Reservoirs in the Eastern Sag of the Liaohe Depression
  29. Regular Articles
  30. Three-dimensional structural model of the Qaidam basin: Implications for crustal shortening and growth of the northeast Tibet
  31. Regular Articles
  32. Failure Mode of the Water-filled Fractures under Hydraulic Pressure in Karst Tunnels
  33. Regular Articles
  34. Creating a low carbon tourism community by public cognition, intention and behaviour change analysisa case study of a heritage site (Tianshan Tianchi, China)
  35. Regular Articles
  36. Mapping Mangrove Density from Rapideye Data in Central America
  37. Regular Articles
  38. Marine sediment cores database for the Mediterranean Basin: a tool for past climatic and environmental studies
  39. Regular Articles
  40. Retrofitting the Low Impact Development Practices into Developed Urban areas Including Barriers and Potential Solution
  41. Regular Articles
  42. Spatial uncertainty of a geoid undulation model in Guayaquil, Ecuador
  43. Regular Articles
  44. Structure and Filling Characteristics of Paleokarst Reservoirs in the Northern Tarim Basin, Revealed by Outcrop, Core and Borehole Images
  45. Regular Articles
  46. Ground volume assessment using ’Structure from Motion’ photogrammetry with a smartphone and a compact camera
  47. Regular Articles
  48. Classification of coal seam outburst hazards and evaluation of the importance of influencing factors
  49. Regular Articles
  50. Geochemical characterization of Neogene sediments from onshore West Baram Delta Province, Sarawak: paleoenvironment, source input and thermal maturity
  51. Regular Articles
  52. Influence of Social-economic Activities on Air Pollutants in Beijing, China
  53. Regular Articles
  54. Spectral properties of weathered and fresh rock surfaces in the Xiemisitai metallogenic belt, NW Xinjiang, China
  55. Regular Articles
  56. Geochemistry of sandstones and shales from the Ecca Group, Karoo Supergroup, in the Eastern Cape Province of South Africa: Implications for provenance, weathering and tectonic setting
  57. Regular Articles
  58. Petrology and geochemistry of meta-ultramafic rocks in the Paleozoic Granjeno Schist, northeastern Mexico: Remnants of Pangaea ocean floor
  59. Regular Articles
  60. Distal turbidite fan/lobe succession of the Late Oligocene Zuberec Fm. – architecture and hierarchy (Central Western Carpathians, Orava–Podhale basin)
  61. Regular Articles
  62. Fourier Transform Infrared Spectroscopy of Clay Size Fraction of Cretaceous-Tertiary Kaolins in the Douala Sub-Basin, Cameroon
  63. Regular Articles
  64. Optimized AVHRR land surface temperature downscaling method for local scale observations: case study for the coastal area of the Gulf of Gdańsk
  65. Regular Articles
  66. New non-linear model of groundwater recharge: Inclusion of memory, heterogeneity and visco-elasticity
  67. Regular Articles
  68. “Urban geosites” as an alternative geotourism destination - evidence from Belgrade
  69. Regular Articles
  70. A customized resistivity system for monitoring saturation and seepage in earthen levees: installation and validation
  71. Regular Articles
  72. Consideration of Landsat-8 Spectral Band Combination in Typical Mediterranean Forest Classification in Halkidiki, Greece
  73. Regular Articles
  74. Coda Wave Attenuation Characteristics for North Anatolian Fault Zone, Turkey
  75. Regular Articles
  76. Modal composition and tectonic provenance of the sandstones of Ecca Group, Karoo Supergroup in the Eastern Cape Province, South Africa
  77. Regular Articles
  78. Quantitative studies of the morphology of the south Poland using Relief Index (RI)
  79. Regular Articles
  80. Interpretation of sedimentological processes of coarse-grained deposits applying a novel combined cluster and discriminant analysis
  81. Regular Articles
  82. Utilizing borehole electrical images to interpret lithofacies of fan-delta: A case study of Lower Triassic Baikouquan Formation in Mahu Depression, Junggar Basin, China
  83. Regular Articles
  84. Grain size statistics and depositional pattern of the Ecca Group sandstones, Karoo Supergroup in the Eastern Cape Province, South Africa
  85. Regular Articles
  86. Carbonate stable isotope constraints on sources of arsenic contamination in Neogene tufas and travertines of Attica, Greece
  87. Regular Articles
  88. Appreciation of landscape aesthetic values in Slovakia assessed by social media photographs
  89. Regular Articles
  90. Geochemistry of Selected Kaolins from Cameroon and Nigeria
  91. Regular Articles
  92. Spatial pattern of ASG-EUPOS sites
  93. Regular Articles
  94. A Stream Tilling Approach to Surface Area Estimation for Large Scale Spatial Data in a Shared Memory System
  95. Regular Articles
  96. A location-based multiple point statistics method: modelling the reservoir with non-stationary characteristics
  97. Regular Articles
  98. Water Inrush Analysis of the Longmen Mountain Tunnel Based on a 3D Simulation of the Discrete Fracture Network
  99. Regular Articles
  100. A Computer Program for Practical Semivariogram Modeling and Ordinary Kriging: A Case Study of Porosity Distribution in an Oil Field
  101. Regular Articles
  102. Imaging and locating paleo-channels using geophysical data from meandering system of the Mun River, Khorat Plateau, Northeastern Thailand
  103. Regular Articles
  104. Rare earth element contents of the Lusi mud: An attempt to identify the environmental origin of the hot mudflow in East Java – Indonesia
  105. Regular Articles
  106. Is Nigeria losing its natural vegetation and landscape? Assessing the landuse-landcover change trajectories and effects in Onitsha using remote sensing and GIS
  107. Regular Articles
  108. Methodological approach for the estimation of a new velocity model for continental Ecuador
Heruntergeladen am 20.9.2025 von https://www.degruyterbrill.com/document/doi/10.1515/geo-2017-0036/html
Button zum nach oben scrollen