Mapping and assessment of mangrove environment are crucial since the mangrove has an important role in the process of human-environment interaction. In Indonesia alone, 25% of South East Asia's mangroves available are under threat. Recognizing the availability and the ability of new sensor of Landsat data, this study investigates the use of Landsat ETM + 7 and Landsat 8, acquired in 2002 and 2013 respectively, for assessing the extent of mangroves along the South Sulawesi’s coastline. For each year, a supervised classification of the mangrove was performed using open source GRASS GIS software. The resulting maps were then compared to quantify the change. Field work activities were conducted and confirmed with the changes that occurred in the study area. Considering the accuracy assessment, our study shows that the RGB composite color-supervised classification is better than band ratio-supervised classification methods. By linking the open source software with the Landsat data and Google Earth satellite imagery that is available in public domain, mangroves forest conversion and changes can be observed remotely. Ground truth surveys confirmed that, decreases of mangroves forest is due to the expansion of fishpond activity. This technique could potentially allow rapid, inexpensive remote monitoring of cascading, indirect effects of human activities to mangroves forest.
Mangrove forests are among the most productive and biologically important ecosystems of the world because they provide important and unique ecosystem goods and services to human society, coastal and marine systems. The forests help stabilize shorelines and reduce the devastating impact of natural disasters such as tsunamis [
The mangrove biome dominates tropical and sub-tropical coastlines between latitudes 32˚N and 38˚S and covers approximately 22 million hectares. Around 28% of global mangroves are located in Southeast Asia, with Indonesia alone accounting for 25% [
Compared to other countries in Southeast Asia, Indonesia has the highest diversity with 43 true mangrove species [
The status of the mangrove environment in Indonesia has been assessed in some studies ( [
Unfortunately, it has been reported that worldwide, including in Indonesia, the area of mangrove forests have declined at a faster rate than inland tropical forests and coral reefs [
Other major threats to mangrove forest in Indonesia are oil spill from oil refinery industries and the relative sea-level rise [
For the protection of the mangrove forest the Indonesian government passed the regulation that a 50 - 200 m wide belt of mangroves has to be retained along the coast to preserve the ecological functions and ensure the natural regeneration process [
To preserve and protect the mangrove forest area from further illegal conversion, intensive mapping and monitoring is required. Yet considering the vast area, direct monitoring will require high cost and long period of time. Remote sensing provides lots of options for continuous monitoring of mangrove over different periods. On this note, mapping and monitoring techniques for the mangrove environment has been widely studied.
A work by Souza-Filho et al. [
Another study by Li [
Nascimento et al. [
A monitoring and mapping method for mangrove area which is cost-effective, e.g. applying free software and public data, will then be of great advantage. The present study aimed to develop an accurate and inexpensive classification method that will enable a more sustainable management of mangrove conservation.
Accessibility of Landsat data is freely available, and the recent launch of Landsat 8 will ensure the continued acquisition and availability of Landsat data utilizing a two-sensor payload, the Operational Land Imager (OLI) and the Thermal InfraRed Sensor (TIRS). Respectively, these two instruments will collect image data for nine shortwave bands and two longwave thermal bands.
Providing moderate-resolution imagery, from 15 meters to 100 meters, of Earth’s land surface and polar regions, Landsat 8 operates in the visible, near-infrared, short wave infrared, and thermal infrared spectrum. It also captures approximately 400 scenes a day, an increase from the 250 scenes a day on Landsat 7 [
OLI sensor collects data from nine spectral bands. Seven of the nine bands were consistent with the Thematic Mapper (TM) and Enhanced Thematic Mapper Plus (ETM+) sensors found on earlier Landsat satellites, which provides compatibility with the historical Landsat data, while also improving measurement capabilities. Two new spectral bands, a deep blue coastal/aerosol band and a shortwave-infrared cirrus band will be collected, allowing scientists to measure water quality and improve detection of high, thin clouds [
This study utilizes the benefit of new Landsat 8 data to develop a method of mangrove mapping and assessment. A supervised classification technique was developed for mapping and assessing the mangrove environment along South Sulawesi coastline, with specific focus along coastal environment in Maros Regency, the South Sulawesi Province, Indonesia. The study area was chosen, because the coastal environment in Maros Regency is one of the places that are under threat from anthropogenic activity.
The method was performed using free open source GRASS GIS and Quantum GIS software (available at http://www.qgis.org/en/site/forusers/download.html). Both are distributed as free software under the GNU Public License. GRASS GIS and Quantum GIS software presently are integrated in one graphical user interface, enabling the seamless combination of both software.
Multispectral images of Landsat 8 OLI (for year 2013) and Landsat 7 ETM+(for year 2002) acquired from http://earthexplorer.usgs.gov/, both at the spatial resolution of 15 m after fusion technique employed, using panchromatic band, were used in this study. The multispectral imagery consists of 11 spectral bands and 8 spectral bands respectively, for Landsat 8 OLI and Landsat 7 ETM+. The imageries were acquired over the Makassar Strait region (
Pre-processing, such as radiometric and atmospheric corrections, which is necessary for analysis of vegetation indices and land use/land cover parameters, were conducted. The Landsat 8 and 7 ETM+ level 1G images were geometry-corrected products.
The Landsat 8 and 7 ETM+ imagery were converted from digital numbers (DN) to Top-of-atmospheric (TOA) reflectance using a two-step process. The initial step is converting the DN values into radiance values using the Lmin and Lmax spectral radiance scaling factors. The values are specific to the individual scene and are produced in the image header file. However, the new .MTL file of Landsat 7 ETM+ data has “mult” and “add” values that are similar to Landsat 8. These values can be read as gain and bias for TOA radiance calibration. “Mult” is a multiplicative rescaling factor, while “add” is an additive rescaling factor for each band.
The second step is converting the radiance data into reflectance using the equations from the Landsat 7 Science Data Users Handbook to the ETM+ image [
The calculation method for image processing in GRASS GIS employs mathematical operators and functions in expressions involving maps or images. The maps and images which are stored in raster grid-cell format are two-dimensional matrices of integer values. The syntax for the algebra is result = expression where expression is built using maps and images, mathematical operators, functions, and temporary variables [
would produce a map Result where each cell is the multiplication of the values of the corresponding cells in map 1 and map 2.
To convert the DN values into reflectance in GRASS GIS software, we employed the “r.mapcalc” expression in the command console. In the command console of GRASS GIS, the formula to convert the DN of the Landsat 8 is written as follows (Equation (1)):
r.mapcalc
where,
B(i)TOA = Name for output of TOA planetary reflectance of Band (i).
REFLECTANCE_MULT_BAND(i) = Band (i) multiplicative rescaling factor from the metadata (.MTL).
B(i) = Quantized and calibrated standard product pixel values (DN) of Band (i).
REFLECTANCE_ADD_BAND(i) = Band (i) additive rescaling factor from the metadata.
The NDVI (Normal Difference Vegetation Index) [
NDVI has been widely monitored the quality and distribution of vegetation. This index can be computed using a simple formula (Equation (2)),
where Red is the reflectance or radiance in a visible wavelength channel (0.630 - 0.690 μm), and corresponds to band 3 and band 4 (0.630 - 0.680 μm) for ETM+ and OLI images, respectively. In the command console of GRASS GIS, the formula to calculate NDVI is written as follows (Equation (3)):
NDWI is sensitive to changes in the water content of vegetation canopies. It is considered as an independent vegetation index that was developed to delineate vegetation water content features, and to enhance their presence in remotely-sensed digital imagery. NDWI is expressed using the following (Equation (4)):
where,
NIR is the reflectance or radiance in a near infra-red wavelength channel (0.760 - 0.900 μm), and SWIR is the reflectance or radiance in a short wave infrared wave-length channel (1.550 - 1.750 μm).
NIR and SWIR correspond to bands 4 and 5 for ETM+ images, respectively. While for OLI images, NIR and SWIR correspond to band 5 (NIR, 0.845 - 0.885 μm) and band 6 (SWIR, 1.560 - 1.660 μm). In the command console of GRASS GIS, the formula to calculate NDWI of the Landsat 8 OLI is written as follows (Equation (5)):
An RGB (Red-Green-Blue) colour composite was produced using a blue band of Landsat 7 ETM+ and a deep blue coastal band of Landsat 8 OLI respectively in the Red channel, NDVI, and NDWI images were loaded into the green and blue channels, respectively. To produce RGB colour composite images, we used the “r.composite” expression in GRASS GIS software (Equation (6)).
where,
r = Red channel.
g = Green channel.
b = Blue channel.
output = The resulting raster map name of colour composite.
A data fusion technique was performed using the RGB image and the panchromatic band of Landsat 8 and 7 ETM+. Image fusion was performed to prepare data for visual interpretation and enhance the classification result. The high spatial resolution of the panchromatic image (Band 8) with the multispectral content of multiband images was fused into a single band. Since the Landsat 8 OLI store values as “DN” from zero to 65,535, we need to rescale it into zero to 255 to get an acceptable color-balanced composite image after pan-sharpening. Rescale method can be done using the “r.rescale” expression on the command console of GRASS GIS (Equation (7)). This study employed the Brovey spectral sharpening technique [
where,
ms1 = Name of input raster map (Green: B2).
ms2 = Name of input raster map (NIR: B4).
ms3 = Name of input raster map (MIR: B5).
Pan = Name of input raster map (Pancromatic: B8).
output prefix = Name for output raster map prefix (e.g. “fused7”).
in = The name of the raster map to be rescaled.
out = The resulting raster map name.
from = min,max (the input data range to be rescaled).
to = min,max (the output data range).
Image classification was performed on the fused image result. A maximum likelihood decision rule method was used to extract a mangrove map. Two classes were defined for the mangrove map: mangrove forest and non- mangrove forest (water, fishpond, built-up area). After the classification, the raster classification result was then converted into vector format.
Prior to classification, satellite images were subsetted to include only areas where mangrove forest is likely to occur. Based on the field work survey, mangroves 10 meters above sea level (masl) and 2 km of coastline were observed. The image was a subset by elevation using a digital elevation model (DEM) extracted from ASTER GDEM data to confine the mapping area below 10 masl. We performed a neighborhood filter on the ASTER GDEM data using the “r.neighbors” expression (Equation (9)), and extracted contour line using the “r.contour” expression (Equation (10)) in GRASS GIS software.
where,
input = Name of input raster map.
output = Name for output raster/vector map.
minlevel = Minimum contour level.
maxlevel = Maximum contour level.
step = Interval contour.
To generate a supervised classification map using specifically digitised training areas; we need to create training areas, but before training areas can be digitised, a vector map layer has to be created and added to the layer manager. The process of creating a vector map layer ready for training area digitisation is available from the graphical user interface of GRASS GIS; vector-develop vector map-create a new vector map. The vector training areas then can be converted into raster format using the “v.to.rast” expression (Equation (11)).
The “i.group” expression (Equation (12)) was used to select all relevant bands (blue band of Landsat 7 ETM+, a deep blue coastal band of Landsat 8 OLI, NDVI images, and NDWI images) of the raster images.
To generate statistics from the training area, the “i.gensig” expression (Equation (13)) was used. The input into the “i.gensig” for creation of signature statistics needs a ground truth training area image that is in raster format, which the conversion format has done using Equation (11).
The classification algorithm was employed using the “i.maxlik” expression (Equation (14)) in GRASS GIS. The classification result was converted into vector format using the “r.to.vect” expression (Equation (15)). Quantum GIS software was then used for the final analysis of mangrove forest changes and cartography design. The overall algorithm that was used in this study was shown in
where,
input = Name of input vector/raster map.
output = Name for output vector/raster map.
group = Name of imagery group.
subgroup = Name of imagery subgroup.
trainingmap = Name of training area in raster format.
signaturefile = Name of signature file.
class = Name for output classification result.
-s: smoothing the vector data.
feature = Feature type (area).
The classification was based on the band ratio and so the subsequent supervised classification had the possibility of producing the best result if the canopy spectra of mangrove vegetation. In the study of mangroves, Rahman [
However, in the present study, the same band ratio pair cannot be used, since the PCA analysis is dependent to the specific imagery, and that each new image requires a new recalculation. Therefore, we recalculated and combined the band ratio analysis with the principal component analysis (PCA), a common technique for extracting essential features of the data with high dimensionality ( [
PCA is a more powerful approach for feature extraction. PCA identifies the optimum linear combinations of the original channels that can account for the variations of pixel values within an image. Optimum values for coefficients are calculated by a procedure that ensures that the values they produce account for maximum variation within the entire dataset. This set of coefficients provides the maximum information that can be conveyed by any single channel formed by a linear combination of the original channels. If we make an image from all the values formed by applying this procedure to an entire image, we generate a single band of data that provides an optimum depiction of the information. However, the effectiveness of this procedure depends on the calculation of the optimum coefficients. More complete explanation requires the level of detail provided by Davis [
To perform PCA, this study calculated a variety of variables from all data sets: including eigenvalues, eigenvectors, and percent of importance. These values were calculated for each principal component (PC), and the results were shown in
where,
input = All band that used in this study.
output_prefix = Name for output raster map.
From the eigenvectors (
where,
B1-5 & B7 = Digital Numbers for the Landsat 7 ETM+ and Landsat 8 OLI, bands 1 - 5 & 7.
Since PC1 has the most positive (negative) contribution from B7 (B1), the band ratio derived from PC1 to Landsat 7 ETM+ is B1/B7. Similarly, the results for PC2 - PC7 lead to the choice of the ratios B1/B5, B7/B4, B3/B1, B5/B7, and B3/B2, respectively.
The band ratio derived from PC1 to Landsat 8 OLI is B5/B1. Similarly, the results for PC2 - PC7 lead to the choice of the ratios B5/B2, B5/B7, B7/B4, B7/B2, and B2/B1, respectively.
In accordance with the result of PCA analysis, we then employed GRASS GIS software to make a band ratio image using the “r.mapcalc” expression (Equation (19)).
PC1 | PC2 | PC3 | PC4 | PC5 | PC7 | |
---|---|---|---|---|---|---|
Band 1 | −0.455 | −0.596 | 0.405 | 0.458 | 0.053 | −0.244 |
Band 2 | −0.437 | −0.255 | −0.036 | −0.297 | −0.063 | 0.807 |
Band 3 | −0.445 | −0.168 | −0.446 | −0.524 | −0.173 | −0.520 |
Band 4 | −0.346 | 0.488 | 0.647 | −0.376 | 0.258 | −0.123 |
Band 5 | −0.433 | 0.521 | −0.128 | 0.450 | −0.566 | 0.046 |
Band 7 | −0.311 | 0.204 | −0.448 | 0.288 | 0.759 | 0.041 |
Note: The highest positive and negative loading in every PC column is marked as red and yellow shadings, respectively.
PC1 | PC2 | PC3 | PC4 | PC5 | PC7 | |
---|---|---|---|---|---|---|
Band 1 | −0.203 | 0.404 | −0.232 | −0.437 | 0.281 | 0.654 |
Band 2 | −0.221 | 0.422 | −0.255 | −0.304 | 0.284 | −0.729 |
Band 3 | −0.268 | 0.343 | −0.196 | 0.372 | −0.143 | 0.187 |
Band 4 | −0.287 | 0.390 | −0.077 | 0.516 | −0.354 | −0.052 |
Band 5 | −0.703 | −0.596 | −0.349 | −0.112 | −0.116 | −0.010 |
Band 7 | −0.279 | 0.179 | 0.583 | −0.472 | −0.562 | −0.031 |
Note: The highest positive and negative loading in every PC column is marked as red and yellow shadings, respectively.
Extensive field work was carried out September 1st, 2013. Location data are collected using Trimble Juno GPS for accuracy assessments. Local photographs and Google Earth imagery were used for visual inspection. Vector points data and fused images were transformed into a KML/KMZ format to enable land cover validation in Google Earth imagery. One hundred points were selected randomly from each classification result, and were then superimposed into Google Earth and a visual assessment was then completed (
To generate a kappa value, we needed to import shapefile data of ground truth into GRASS GIS with “v.in.ogr” expression (Equation (20)), and thereafter convert again the vector to a raster with “v.to.rast” expression (Equation (21)). Then we employed “r.kappa” expression (Equation (22)) in GRASS GIS command console.
where,
layer = Name for input shapefile vector data.
input = Name for input vector map.
output = Name for output vector/raster map.
classification = Name of raster map containing classification result.
reference = Name of raster map containing reference classes.
output = Name for output file containing error matrix and kappa.
From the field work activity, we found that there are two types of mangrove trees that exist in the study area, Rhizophora and Bruguiera as shown in the following
Furthermore, to enhance the analysis, we use the Global Distribution of Mangroves USGS [
This study utilized ground truth data to evaluate the accuracy of the classified mangrove forest obtained from Landsat 7 ETM+ and Landsat 8 OLI. The field validation and visual assessment through Google Earth imagery was used as reference to check the accuracy of the classified images.
The maximum likelihood classification decision rule was implemented quantitatively to consider several classes and several spectral channels simultaneously as well as form a powerful classification technique. However, the training area needs to be carefully selected; otherwise the classification result may introduce error. The maximum likelihood classification uses the training area as a means of estimating means and variance of the classes, which are then used to estimate the probabilities. Maximum likelihood classification considers not only the mean, or average, values in assigning classification, but also the variability of the brightness value of each class [
It can be seen from the RGB composite colors (
From the analysis using the Global Distribution of Mangroves as the reference data, the study provided the perfect strength of agreement, with the kappa statistic of 0.98 as the classification image of year 2013 using the RGB composite color-supervised classification (maximum likelihood decision rule).
Class/Region | Water | Mangrove | Fishpond | Built-up |
---|---|---|---|---|
Water | 98.00% | 0.000% | 7.000% | 0.000% |
Mangrove | 0.000% | 98.00% | 0.000% | 4.000% |
Fishpond | 1.00% | 0.000% | 92.00% | 0.000% |
Built-up | 0.000% | 2.000% | 0.000% | 95.00 |
Overall kappa = 0.95 | ||||
Overall classification accuracy = 96.47% |
To select the bands for band ratio classification, the eigenvalues and percentage of importance that was produced by PCA method showed that the first three eigenvalues contain approximate 99.10% and 99.68% of the total bands of Landsat 7 ETM+ and Landsat 8 OLI respectively. The band ratios of B1/B7, B1/B5, and B7/B4 were found to yield the major features of the classification from Landsat 7 ETM+ images. B5/B1, B5/B2, and B5/B7 yielded the major features of the classification from Landsat 8 OLI images. The kappa analysis of the Landsat 8 OLI classified images using band ratios yielded a value = 0.83
Furthermore, from the field work activities, mangrove area was observed to decrease due to the expansion of fishpond activity. However, the decrease rate is low (11 ha in period 2002-2013). To promote sustainable development in mangrove forest in the region, community participation in rehabilitation, conservation and management of mangroves had been carried out through mangrove rehabilitation programs supported by the Department of Forestry. Mangrove conservation and rehabilitation were initiated and promoted collaboratively by both local people and governmental institutions [
The results of the supervised classification methods are provided in
The spectral characteristics that were obtained from the ETM+ and OLI datasets, based on the physical properties of mangrove forests, in combination with ASTER GDEM data, increased the accuracy of mapping the mangrove forests. In addition, the RGB combination of the NDWI, NDVI, and the band 1 of Landsat 8 OLI, increased the accuracy of the classification. Using this method, mangrove forests have been successfully separated from other coastal land uses, located in the region.
In many cases, human eye interpretation of images was found as one of the best ways of obtaining information. The findings of the present study highlight the importance of integrating types of methods and visual and automated method.
Many studies have been conducted for mangroves forest mapping and assessment, for example Kiruiet [
Rahman [
Furthermore, Santos [
This study highlights the importance of freely available and less economically geospatial software and data available in the public domain, such as GRASS GIS software and Landsat data, which can be used for such purposes and, consequently, producing information to aid in decision making and elaboration of management plans for coastal conservation.
Freely-available satellite imagery of the entire Earth’s surface via the Landsat data and Google Earth allows examination of landscape features in even the most remote areas, including difficult-to-access habitats within them. Here we demonstrated that by linking the open source software with this data available in the public domain. Therefore, it is possible to remotely observe the mangroves forest changes and conversion for less economically developing nations/organization to employ remote sensing methods.
A limitation of this method is the inability to discriminate between different species of mangroves. Furthermore, the frequent existence of clouds in tropical and subtropical areas is considered as the main limitation in using optical remote sensing data, which must be taken into consideration when producing regional maps of mangrove forests.
This study has shown that Landsat 7 ETM+ and Landsat 8 OLI are suitable for the regional mapping and assessment of mangrove forests. The algorithm proposed in this study provides an inexpensive framework for the application of mapping and assessing mangroves on a large scale.
The utilisation of the deep blue coastal band in Landsat 8 OLI was useful for highlighting mangroves and separating them from water bodies when combined with the RGB color composite. In the present study, the choice of three bands (deep blue coastal band, NDVI, and NDWI) was suitable, resulting in an accuracy of 96.47%. The band ratios approach resulted in an accuracy of 86.73%.
GRASS GIS opens source software is freely available for use by universities, government agencies, and other organizations for any GIS applications. The capability of GRASS GIS is comparable with commercial software which has a high cost and is therefore not affordable for developing countries. In this study GRASS GIS has been applied to the coastline environment to monitor and assess the changes of mangrove’s foliage.
By linking the open source software with freely-available, the Landsat data and Google Earth satellite imagery, mangroves forest conversion and changes can be observed remotely. By using sequential, the Landsat data and Google Earth images of specific locations over time, this technique could potentially allow rapid, inexpensive remote monitoring of cascading, indirect effects of human activities in mangroves forest nearly anywhere on earth.
Ground truth surveys confirmed that, decreases in the mangrove forest are due to the expansion of fishpond activity. Regular monitoring and assessment is needed in order to monitor the changes in the mangrove environment. Furthermore, understanding the processes of change is also needed in order to minimize the negative impact from the human-environment linkage. This approach has the potential to lead to the sustainable development of the coastal environment.
The author would like to thank the team at Hasanuddin University, Makassar, South Sulawesi, Indonesia for conducting field work to support this study. We would also like to thank the reviewers of this paper and their valuable contribution towards refining the manuscript.