Next Article in Journal
Taxonomic Review of Fossil Coleopterous Families (Insecta, Coleoptera). Suborder Archostemata: Superfamilies Coleopseoidea and Cupedoidea
Previous Article in Journal
Cascading Dynamics of the Hydrologic Cycle in California Explored through Observations and Model Simulations
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Reef Mapping Using Different Seabed Automatic Classification Tools

1
Laboratório de Geociências Marinhas, Dept Oceanografia, Universidade Federal do Espírito Santo, Vitoria-ES 29075-910, Brazil
2
Bolsista GEF MAR ICMBio, Parque Nacional Marinho dos Abrolhos, Caravelas–BA 45900-000, Brazil
3
SAGE/COPPE-Instituto de Biologia, Universidade Federal do Rio de Janeiro, Rio de Janeiro 21941-902, Brazil
4
Geological Survey of Norway, Postal Box 6315 Torgarden, NO-7491 Trondheim, Norway
*
Author to whom correspondence should be addressed.
Geosciences 2020, 10(2), 72; https://doi.org/10.3390/geosciences10020072
Submission received: 18 January 2020 / Revised: 11 February 2020 / Accepted: 12 February 2020 / Published: 15 February 2020

Abstract

:
There is a great demand to develop new acoustic techniques to efficiently map the seabed and automate the interpretation of acoustic, sedimentological, and imaging data sets, eliminating subjectivity. Here, we evaluate the potential, limitations and complementariety of distinct supervised and automatic classification techniques in the mapping of reefs by comparing these results with a reference map. The study was carried out in the Abrolhos Continental Shelf (Eastern Brazilian Continental Margin) using a multibeam echosounder and side scan sonar (SSS) dataset. Two automatic supervised techniques were applied. A reference map was derived by detailed manual interpretation carried out by three experts. The two supervised classification techniques were: benthic terrain modeler (BTM), a morphometric classification with focus on spatial analyses of the bathymetric grid derivatives, and object-based image analysis (OBIA), a segmentation applied to the backscatter data from the SSS mosaic. Both automatic techniques obtained similar values of reef coverage area, but overestimated the reef area when compared with the reference map. The agreement between BTM and OBIA results and the reference map was 69% and 67%, respectively. Disagreement was mainly due to quantity of reef (both methods over-estimated reef), while the disagreement in spatial allocation was relatively low, it indicates that both methods are reasonable representation of the spatial patterns of reef. Efficient mapping of reef in the wider area of the Abrolhos Continental Shelf will be best achieved by a further development of automatic methods tested against reference maps obained from representative areas of the seabed. By combining the results of the two automatic methods, it was possible to create an ensemble map, which achieved better agreement with the reference dataset.

1. Introduction

Terrestrial remote sensing has a comparatively long history since the 1970s and benefits from the availability of land surface imagery at various spatial and spectral resolutions, supplemented by ancillary data such as environmental data and digital terrain models. This enabled the development of sophisticated classification methods and algorithms that extract critical information for environmental management from remote sensing imagery [1,2]. Conversely, seabed automatic classification efforts are more recent [3,4,5] and mostly based on acoustic mapping technologies, since the application of electromagnetic methods has severe limitations due to the rapid attenuation of light in water.
The development of new technologies and classification tools is in great demand from the industry, government agencies and scientific institutions. The increasing interest and necessity for seabed habitat mapping in marine spatial planning initiatives has forced the industry of sonar makers, software developers and scientific research institutions to invest in new technology capable of producing seabed classification over regional areas [6]. Fisheries, oil and gas exploration, offshore engineering, as well as shallow and deep sea mining are examples of industries that are requiring seabed classification to support environmental risk assessments and spatial planning [6,7,8].
Side scan sonars (SSS) and multibeam echosounders (MBES) are among the most widely applied acoustic techniques for seabed mapping. Both systems have high sampling rates, which enables seabed mapping at a broad spatial scale, producing high resolution data sets in a relatively short period of time [9,10,11,12].
A series of classification models have been developed to reduce the subjectivity and automate the interpretation of geophysical, sedimentological and imaging data sets obtained from the seabed. Supervised and unsupervised classifications have improved the seabed classification [13,14] with the potential to also compile statistical analyses of those datasets.
The study presented here attempts to automatically classify and map submerged reefs. The identification and mapping of reefs through automatic classification from swath-mapping systems has great potential and applicability due to the global importance of these complex habitats, regarding planning and management of marine protected areas (MPA) and fisheries resources [15,16,17].
The analyses were carried out based on a dataset collected in the Abrolhos Continental Shelf (Eastern Brazilian Continental Margin), where diverse reef morphologies, including pinnacles, reef banks and rhodoliths, have been observed [18,19]. In order to automatically classify the distinct reef types, two supervised classification techniques were applied, taking into account that each of them has a different type of data input and distinct architecture. One technique is a morphometric classification tool (Benthic Terrain Modeler, BTM), with focus on the spatial analyses of the bathymetry derivative properties. The second technique is an object-based image analysis (OBIA, using eCognition software) applied only to the backscatter data from the SSS.
We have also compared the different methods using an assessment of map accuracy. To achieve this, the classification maps obtained from BTM and eCognition were compared against a reference dataset. In this context, the aim of this contribution is to analyze the potential benefits and limitations of the distinct techniques in the recognition and automatic classification of reefs through the complementary use of MBES or SSS or both data sets.

2. Materials and Methods

2.1. Acquisition, Processing and Analysis of Geophysical Data

The acoustic survey was carried out on the Abrolhos Continental Shelf (east Brazil), along the Abrolhos outer reef arc, known as California Reefs (Figure 1). The California Reefs are within the Abrolhos Marine National Park and the water depth ranges from 12 m (top of pinnacles) to 35 m deep (seabed).
MBES (Reson 7101-240 kHz) and SSS data (EdgeTech 4100-500 kHz) were acquired covering an area of approximately 20 km2. Both systems acquired data simultaneously. Bathymetric data were processed using Caris HIPS and SIPS 9.1. The applied workflow included patch test correction, correction for tidal elevation and measured sound velocity and spikes removal. A 5 × 5 m resolution grid was exported as an XYZ file and a digital bathymetric model (DBM) was produced in a GIS platform (ArcMAP 10.6) and AUTOCAD 3D. The SSS data were processed with SonarWiz 7 software (Chesapeake Technology), applying bottom tracking and the EGN (empirical gain normalization) filter workflow. A 1 × 1 m resolution georeferenced sonographic mosaic was exported as a geotiff image.

2.2. Seabed Classification Techniques

Two automatic supervised techniques were applied to classify reef occurrence and distribution. The results were compared against a reference map to assess agreement with reference data and estimate map error. Moreover, an analysis of overlap among the results from the different classifications was carried out. The following flow diagram (Figure 2) summarizes the analysis steps. The analyses were executed by different experts under the supervision of the same third expert. Thus, the final analysis attempted to standardize all reef classes in a single class (reef structures).
Despite the difference of classes names for each method, all of them were subcategorized in “Reef presence” or “Reef absence”. In terms of nomenclature, all classes categorized as “Reef absence” are defined as inter-reef areas.

2.2.1. Reference Map

The reference map was derived by manual interpretation of both the SSS backscatter and bathymetric data. Manual digitalisation of reef tops was guided by high backscatter, typical geomorphological shape and roughness produced by these conspicuous seabed features. To limit uncertainty in the reference map, the features were interpreted and reviewed by three experts to map the spatial occurrence of submerged reefs in the survey area. Extensive background knowledge [18,19] and the conspicuous nature of the reefs allowed us to confidently map reefs in the Abrolhos Continental Shelf even in the absence of ground-truth data.

2.2.2. Morphometry-Based Classification

The BTM is an ArcGIS® toolbox capable of analyzing digital terrain models derived from the different bathymetric derivative parameters, such as the Bathymetric Positioning Index (BPI) and the terrain slope [20]. The BPI is a bathymetric derivative that evaluates the elevation differences between a focal point and the mean elevation of its surrounding cells under broad and fine scales [21]. BPI and slope were calculated from the 5 × 5 m resolution grid. The inner and outer radius used for BPI fine was 5–50, and for BPI broad was 300–700. A final BTM map is generated using a rule-based classification of the morphometry based on a dictionary, where BPI, slope and depths are categorized into classes according to values designated by the user (Table 1). The dictionary was derived by an exhaustive testing of different parameter combinations and a visual assessment of the results. The categories within the dictionary work within a lower and upper bound. For BPI, the number 100 means the number of grid units that was used. Negative values mean below the standard deviation while positive values mean above it. For slope, we set the angle threshold as between 2 and 13°. Finally, in terms of depth, low relief reef banks were within 26 and 32 m water depth. All the 9 designated classes on the final result were elected due to the geomorphological heterogeneity of the seabed.

2.2.3. Object-Based Image Analysis Classification

The eCognition 9.1 software was used to provide tools for image segmentation and classification. The segmentation methodology from eCognition is based on image characteristics such as color, texture and shape of the feature. This technique can be considered as Object-Based Image Analysis (OBIA) and has been widely used in terrestrial remote sensing applications but also successfully applied in marine habitat mapping research [17,22]. The 1 × 1 m resolution SSS mosaic was used for the image segmentation. The supervised analysis of the sonographic mosaic involved three main steps: (i) segmentation according to defined parameters (multiresolution segmentation algorithm); (ii) determination of the number of classes based on the samples trained by user. The definition of the classes is a procedure performed by the analyst. In this case, they were determined according to the segmentation result. The sample training process is also a supervised process, in which the analyst chooses one or more segments to represent each defined class; (iii) automated classification. The eCognition scale parameter applied was 150 and a number of classes was determined as 5, according to a visual analysis of the mosaic.

2.2.4. Map Agreement and Error

The derived maps were simplified into reef presence/absence maps prior to further analysis. We carried out a direct map-to-map comparison in R [23] utilising the packages raster [24] and diffeR [25]. To that end, a contingency table between a comparison map and a reference map is created (Table 2). Disagreement between the comparison map and the reference map can be expressed with two simple metrics, quantity disagreement and allocation disagreement [26]. Quantity disagreement is defined as the amount of difference between the reference map and a comparison map that is due to the less than maximum match in the proportions of the categories (here: reef and non-reef). Allocation disagreement is defined as the amount of difference between a reference map and a comparison map that is due to the less than maximum match in the spatial allocation of the categories, given the proportions of the categories in the reference and comparison maps. The map agreement can be calculated by substracting the sum of quantity disagreement and allocation disagreement (both as fractions) from 1.
We also estimated the class-specific error of omission and error of commision based on the contingency table. Omission error refers to observations in the reference map that were classified in a category other than their true or known category. Commission error refers to observations that were incorrectly classified and do not belong in the category in which they were assigned according to the classification. Based on a 2 × 2 contingency table, errors of omission (omission absence— O A , and omission presence— O P ) and commision (comission absence— C A , and comission presence— C P ) can be calculated in the following way:
O A = F P T N + F P
O P = F N F N + T P
C A = F N T N + F N
C P = F P F P + T P

3. Results

3.1. Seabed Morphology and Geomorphometric Analysis with BTM

The DBM of the study site is shown in Figure 3. Table 3 presents the reef classes defined in each technique used. The seabed map presents a great variability of reef morphologies, showing pinnacles (Figure 3c,d) and banks (Figure 3e,f). Banks are described here as Low Relief Reef Banks (LRRB), occurring in areas deeper than 30 m and generally flat seabed. The California Reef area presents a morphological variation from west to east, as shown in Figure 3b. The reef tops range in depth from 9 to 34 m.
Results of the morphometric (BTM) analysis are shown in Figure 4. Nine classes were defined based on the dictionary: “Deepest inter-reef shelf”, “East deep inter-reef shelf”, “West deep inter-reef shelf” and “Inter-reef Shelf” are associated with the deepest and the lowest roughness regions, which can also be interpreted by the slope and BPI parameters, respectively; “Inter-Reef Structures”, “Pinnacles/Irregular Relief Reef Banks (Irregular Relief Reef Banks-IRRB)”, ”Pinnacles/IRRB (artifacts)”, “Low Relief Reef Banks-LRRB” and “Edges of LRRB” are associated with reefs, having higher slope values and are also depicted in the BPI maps.
In general, the classes described as Inter-reef Shelf and Deep inter-reef shelf (west and east) predominate over the area. Both classes are characterized by relatively low BPI and bathymetric slope (Figure 4a–c). The BPI classifies them as depressed regions. The Deep Inter-Reef Shelf was separated in East and West, and although these classes are morphologically similar, they are placed in regions with distinct shelf morphologies (higher BPI/more rugged relief in the West than the East).
Reefs are well defined from the morphometric analysis. BPI fine and broad scale can easily depict reef distribution considering their elevation in relation to the adjacent seafloor. The LRRB class was also defined by the applied depth boundary (between 26 and 32 m). The LRRB edge was recognized as a class due to positive values of BPI and the limitation of the bathymetric slope values. The Edges of the LRRB class has a steeper slope value than the adjacent shelf and the tops are flatter than the pinnacles. In class LRRB, it is possible to note a limitation in the BTM application in the west-northwest sector of the study area, where there is a greater concentration of Pinnacles/IRRB and Pinnacles/IRRB (artifacts). LRRB were recognized due to both the relatively high values of BPI from the adjacent seabed and the bathymetric range. This may represent a transition from inter-reef shelf to reef dominated shelf.
The “Pinnacles/IRRB” class (dark red in Figure 4) corresponds to the structures with the highest absolute values of BPI (broad and fine scale) and is also highlighted by high bathymetric slope values. Therefore, this class represents either isolated reef structures or merged pinnacles. The other classes (pink and yellow) are in the northwest sector, a region where there is a high concentration of pinnacles and irregular relief reef banks. It is inferred that these regions were distinguished due to the interpolation of the bathymetric data in a fairly irregular region. This can’t be considered as a limitation of the classification technique, but rather a data input gap (from acquisition survey) to the construction of the model.
In total, approximately 6.85 km2, equalling 34.5% of the total area, were mapped as reef.

3.2. Seabed Backscatter and OBIA

The sonographic mosaic (Figure 5) presents three major backscatter features that were interpreted as: inter-reef shelf (more homogeneous)—framed in yellow in Figure 5a, pinnacles (isolated or merged)—framed in blue in Figure 5a, and low relief reef banks—framed in red in Figure 5a (Table 3).
The inter-reef shelf presents a more homogeneous backscatter pattern, with low roughness, representing a flat area. Although the backscatter seems to be more homogeneous and lower values than the reef areas, slight variations in signal intensity can be observed. Reefs are well characterized by a high backscatter and roughness, creating clear acoustic shadows that point to distinct morphologies (Figure 5b–e).
Aiming to adjust the classification to the spatial scale of this study, the spatial parameters of the segmentation needed to be adequate in a way that would make it possible to detail and recognize isolated reef structures such as pinnacles. In other words, it was necessary to set up a segment small enough for identification of pinnacular structures. Seeking to make the comparison easier, we choose to maintain the greatest similarity among the classes names used in each classification technique. After the segmentation, the sonographic mosaic was divided in the following classes: Inter-reef Shelf, Pinnacles, LRRB and IRRB (Figure 6).
In total, approximately 7.61 km2, equalling 38.3% of the total area, was mapped as reef.

3.3. Reference Map

The results of the manual interpretation using the backscatter data from the SSS mosaic (1m resolution) and the digital bathymetric model are shown in Figure 7. The analysis was focused on reef recognition. The result shows three main classes distinguished: pinnacles, reef bank and mixed reef structures (pinnacle and bank). Pinnacles appear isolated in central region, as well as in great concentration as seen in northwest area. Reef Banks (similar to LRRB) are mainly distributed on the east region of the study site.
The merged reef area is 2.96 km2, which represents approximately 15% of the total area.

3.4. Map Agreement and Error

Full details of the analysis can be found in the Supplementary information (S1 and S2). A summary of the analysis of agreement is provided in Table 4. Both methods show similar trends: Agreement with the reference map is 69.1% for BTM and 67.1% for OBIA. The quantity component of disagreement is more pronounced than the allocation component. The disparity is more pronounced in the case of the OBIA result. This is mostly because both methods over-estimate the reef area (BTM: 6.85 km2 and OBIA: 7.61 km2) by a factor of >2. Low allocation disagreement does, however, indicate that the spatial configurations approximate the reference map.
The class-specific error analysis for the BTM map (Table 5) reveals that approximately 30% of reef absences and 38% of reef presences are omitted from the map. The error of commission is low for reef absences (9%), but high for reef presences (73%). Again, this can be linked to the over-estimation of reef. The above described patterns are similar for the analysis of the OBIA map (Table 6), although the error of omission for reef presence is slightly lower (32%).

4. Discussion

We have presented the results of two automated methods (BTM and OBIA) for mapping reefs in a test site situated in the Abrolhos outer reef arc. The accuracy of these methods in detecting reef against a background of non-reef was estimated by comparing them with a reference map that was obtained by detailed manual interpretation of SSS backscatter and MBES bathymetry data. In the following, we will discuss the obtained results qualitatively and quantitatively.
The estimation of map accuracy requires a reference dataset that can be acquired in different ways, but needs to fulfill at least one of the following two conditions: (1) The reference source data has to be of higher quality than the data used for map classification, or (2) if using the same source data for both map and reference classification, the process to create the reference data needs to be more accurate than the process used to generate the map classification to be evaluated [27]. Here, we use the same datasets, but both SSS and MBES data were used to create the reference map, while the two methods utilised either one or the other. Also, we are satisfied that the process of creating the reference map is more accurate due to the following reasons: (1) There exists substantial expertise from previous studies [18,19] on which the interpretation is based, (2) the reference map is the consensus of the interpretation of three skilled experts, and (3) reefs are conspicuous features that are easily recognizable in the acoustic data. We acknowledge that the reference map is not “the truth”, but under the given circumstances we are satisfied that it is close to reality and can be used as a reference against which the other methods are assessed.
The quantitative analysis of agreement with the baseline map showed that both methods, despite being very different in terms of input data and methodology, gave comparable results regarding map agreement and error. The map agreement of 67% (OBIA) and 69% (BTM) is acceptable, although lower than what has been reported in some recent studies that mapped reef [28,29,30]. These differences might be related to various factors, e.g., shallow warm-water coral reefs can be mapped with satellite or drone-based optical sensors, which give more spectral information than a one-band SSS backscatter dataset. However, this also indicates that there is room for future improvements and the potential to transfer mapping approaches from the optical remote sensing community [4]. The use of the accuracy statistics proposed by [26] gives some interesting insights: for both methods, disagreement is predominantly associated with quantity, i.e., the proportions of predicted reef and non-reef differ from the reference. This is also reflected in the high errors of commission for reef presence and is due to the over-estimation of reef area by both approaches, mainly attributable to reef pinnacles (see Figure 8a,b). Conversely, allocation disagreement is relatively low (around 10%) for both methods, indicating that the overall patterns and trends are mapped correctly.
Automated OBIA classification based on SSS backscatter intensity did not show seabed differences among inter-reef areas. Conversely, this was possible with the morphometric BTM classification. In these regions, the image characteristics (such as texture) did not respond to the differentiation presented by the BTM. Knowing that distinct morphometric characteristics are not necessarily reflected through SSS backscatter intensity, this technique’s complementarity was essential to achieve a better result, even with this limitation in the image classification.
Regarding reef structure classification (Figure 8d), the results by both automatic techniques were quantitatively similar considering the area covered by reef and disregarding reef sub-type (for example pinnacle or low relief bank). The BTM classification resulted in a coverage area of reef structures (classes: “Inter Reef Structures”, “Pinnacles/IRRB”, “Pinnacles/IRRB-artifacts”, “Edges of LRRB”, “LRRB”) of 6.85 km2, while the OBIA classification gave an area value of 7.61 km2 (classes: “Pinnacles”, “LRRB” and “IRRB”).
Reef classes, such as pinnacles, IRRB and LRRB, were not completly differentiated by both techniques. The BTM approach did not differentiate the isolated Pinnacles from IRRB (the value for the sum of the classes “Inter Reef Structures”, “Pinnacles/IRRB-Artifacts” and “Pinnacles/IRRB“ was 3.04 km2), while the OBIA method gave values of 1.57 km2 and 1.75 km2 for Pinnacles” and IRRB, respectively. The OBIA classification provides a better differentiation of these reef features due to their textural characteristics and form being more easily distinguished.
The LRRBs were spatially overestimated by both automatic techniques. The key factor for the overestimation with BTM was the BPI, since some regions presented this index in a similar way to the morphometric characteristics of those banks. When considering the image segmentation in the OBIA process, the selected scale parameter probably contributed to the overestimation of LRBB. In other words, since one of the objectives was to recognize isolated pinnacles, a large number of segments were generated in the image and this made some of these segments more susceptible to small changes in image characteristics (shape, color, and texture), and therefore the classification became particularly vulnerable to these slight differences.
While under the given circumstances, i.e., mapping of reefs, which exhibit specific morphometric and textural characteristics against a relatively flat and featureless background (non-reef areas) manual mapping can be an appropriate approach, it must be recognized that this is nevertheless a slow and painstaking process, especially when carried out by three analysts as it was done here. It is therefore unlikely that such a mapping approach is a viable option for larger seafloor areas. In this study, we mapped approximately 3 km2 of reef, while the reef area in the wider Abrolhos Shelf was estimated to 8844 km2 [18]. By contrast, automatic techniques can deal more efficiently with larger datasets. It might still be necessary to adapt classification dictionaries and rule-set parameters, but the general workflows could still be applied to other (larger) datasets. A viable process might be to develop reference maps for selected test areas, which are characteristic for the wider area. Then, automatic methods can be trialled and improved, with quantitative assessments of agreement and error as tools to decide on the most promising method.
Alternatively, the results of different automatic approaches could be combined to yield an ensemble map [31]. Such an approach is most efficient with complementary methods, which was the case in this study. By combining the results of the two methods, it was possible to create a new map, which does only include areas classified as reef by both methods. The new map was derived by multiplying BTM with OBIA. For this purpose, we consider reef presence as a value of 1, and absence as a value of 0. This final analysis (full details in Supplementary information S3) shows the over-estimation of reef extent is reduced (Agreement = 75.3%), and the argument of the complementarity of the methods is strengthened.
Future directions for comparing automatic mapping with a manual reference map should consider the combined use of bathymetric and backscatter dataset. Here, we used two different approaches considering the data availability and computer processing limitations, but OBIA could also be tested using both datasets and a higher resolution grid, e.g., 1-m grid cell. Further investigations in this reseacrh should consider the biological description in order to test if the reef morphology could be used as a proxy for reef benthic community changes along water depth. Also, the use of a higher resolution multibeam sonar could provide a better detail mapping of the Abrolhos reefs. The pinnacles in Abrolhos present a specific morphology that changes with depth [32]. Shallower reefs (<10 m deep) have wide tops, forming a mushroon-like morphology. Deepward, reef tops decrease in width and can form a single column at depths around 25–30 m. This was observed using scuba diving, but the resolution of the multibeam sonar used in this study was not able to map these detailed changes in pinnacle morphology.

5. Conclusions

The combined use of the bathymetric and sonographic mapping was fundamental for the characterization and comprehension of seabed geomorphology of the study area. The results provided important information about reef morphologies and have great potential to be integrated with studies aiming to examine the associated biodiversity, mainly in terms of mesophotic reefs.
The presented analysis employing a reference map and two automatic techniques for seabed classification (morphometry and object-based image segmentation) have shown limitations and potential for each one. Morphological variations in inter-reef areas were better captured by the morphometric analysis. Conversely, the object-based image classification achieved a greater differentiation between the types of reef morphologies. The assessment of agreement and error also highlighted that both automatic approaches over-estimated the reef area, while at the same time giving a reasonable representation of the spatial patterns of reef distribution. Both automatic techniques obtained similar values of agreement with the reference map and reef coverage area. The results confirm the usefulness of the two methods for reef mapping purposes in potential future studies. The purpose in using these three distinct techniques was not only to compare their effectiveness, but rather to achieve a better recognition of the reefs and inter-reef areas through their potential complementarity. The final analysis achieved a better agreement with the reference dataset through combining the results of the two methods.

Supplementary Materials

The following are available online at https://www.mdpi.com/2076-3263/10/2/72/s1, S1: Map Comparison BTM, S2: Map Comparison OBIA, S3: Map Comparison.

Author Contributions

Conceptualization, P.S.M., G.B. and A.C.B.; methodology, P.S.M., G.B., F.V.V., A.C.L., L.C.F., M.D.; software, P.S.M., G.B., F.V.V., A.C.L., L.C.F. and M.D.; formal analysis, P.S.M., G.B., F.V.V., A.C.L., L.C.F., A.C.B., R.M., M.D.; investigation, G.B., L.C.F.; resources, A.C.B., R.L.M.; data curation, A.C.B.; writing—original draft preparation, P.S.M., G.B.; writing—review and editing, A.C.B., R.L.M.; visualization, P.S.M., G.B., F.V.V., A.C.L., L.C.F., R.L.M.; supervision, A.C.B.; project administration, A.C.B.; funding acquisition, A.C.B., R.L.M. All the authors have contributed substantially to the work reported. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by CAPES grants Ciências do Mar II and IODP-Brasil; CAPES/CNPq/FAPES PELD Abrolhos grant. This is a contribution to Rede Abrolhos (www.abrolhos.org).

Acknowledgments

We are thankful to the crew of the Armando vessel. ACB, RLM are CNPq fellows. We also take this moment to honor Gilberto Amado Filho, who passed away in the final stages of the manuscript preparation and was a great supporter and passionate of the use of seabed mapping for biodiversity studies.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Lu, D.; Weng, Q. A survey of image classification methods and techniques for improving classification performance. Int. J. Remote Sens. 2007, 28, 823–870. [Google Scholar] [CrossRef]
  2. Lecours, V.; Dolan, M.F.; Micallef, A.; Lucieer, V.L. A review of marine geomorphometry, the quantitative study of the seafloor. Hydrol. Earth Syst. Sci. 2016, 20, 3207–3244. [Google Scholar] [CrossRef] [Green Version]
  3. Anderson, J.T.; Holliday, D.V.; Kloser, R.; Reid, D.G.; Simard, Y. Acoustic seabed classification: Current practice and future directions. ICES J. Mar. Sci. 2008, 65, 1004–1011. [Google Scholar] [CrossRef]
  4. Diesing, M.; Mitchell, P.; Stephens, D. Image-based seabed classification: What can we learn from terrestrial remote sensing? ICES J. Mar. Sci. 2016, 73, 2425–2441. [Google Scholar] [CrossRef]
  5. Kendall, M.S.; Buja, K.; Menza, C.; Battista, T. Where, What, When, and Why Is Bottom Mapping Needed an On-Line Application to Set Priorities Using Expert Opinion. Geosciences 2018, 8, 379. [Google Scholar] [CrossRef] [Green Version]
  6. Lucieer, V.; Roche, M.; Degrendele, K.; Malik, M.; Dolan, M.; Lamarche, G. User expectations for multibeam backscatter data—Looking back into the future. In: Seafloor backscatter data from swath mapping echosounders: From technological development to novel applications. Mar. Geophys. Res. 2017, 39, 23–40. [Google Scholar] [CrossRef]
  7. Bourguignon, S.N.; Bastos, A.C.; Quaresma, V.S.; Vieira, F.V.; Pinheiro, H.; Amado-Filh, G.M.; De Moura, R.L.; Teixeira, J.B. Seabed Morphology and Sedimentary Regimes defining Fishing Grounds along the Eastern Brazilian Shelf. Geosciences 2018, 8, 91. [Google Scholar] [CrossRef] [Green Version]
  8. Mayer, L.; Jakobsson, M.; Allen, G.; Dorschel, B.; Falconer, R.; Ferrini, V.; Lamarche, G.; Snaith, H.; Weatherall, P. The Nippon Foundation—GEBCO Seabed 2030 Project: The Quest to See the World’s Oceans Completely Mapped by 2030. Geosciences 2018, 8, 63. [Google Scholar] [CrossRef] [Green Version]
  9. Walker, B.K.; Riegl, B.; Dodge, R.E. Mapping coral reef habitats in southeast Florida using a combined technique approach. J. Coast. Res. 2008, 24, 1138–1150. [Google Scholar] [CrossRef] [Green Version]
  10. Le Bas, T.P.; Huvenne, V.A.I. Acquisition and processing of backscatter data for habitat mapping—Comparison of multibeam and sidescan systems. Appl. Acoust. 2009, 70, 1248–1257. [Google Scholar] [CrossRef]
  11. Brown, C.J.; Smith, S.J.; Lawton, P.; Anderson, J.T. Benthic habitat mapping: A review of progress towards improved understanding of the spatial ecology of the seafloor using acoustic techniques. Estuar. Coast. Shelf Sci. 2011, 92, 502–520. [Google Scholar] [CrossRef]
  12. International Hydrographic Organization, Intergovernmental Oceanographic Commission. The IHO-IOC GEBCO Cook Book; IHO Publication B-11; IOC Manuals and Guides 63; IHO: Monaco, France, 2015; 429p. [Google Scholar]
  13. Stephens, D.; Diesing, M. A Comparison of Supervised Classification Methods for the Prediction of Substrate Type Using Multibeam Acoustic and Legacy Grain-Size Data. PLoS ONE 2014, 9, e93950. [Google Scholar] [CrossRef] [PubMed]
  14. Ierodiaconou, D.; Schimel, A.C.G.; Kennedy, D.; Monk, J.; Gaylard, G.; Young, M.; Diesing, M.; Rattray, A. Combining pixel and object based image analysis of ultra-high resolution multibeam bathymetry and backscatter for habitat mapping in shallow marine waters. Mar. Geophys. Res. 2018, 39, 271–288. [Google Scholar] [CrossRef]
  15. Locker, S.D.; Armstrong, R.A.; Battista, T.A.; Rooney, J.J.; Sherman, C.; Zawada, D.G. Geomorphology of mesophotic coral ecosystems: Current perspectives on morphology, distribution, and mapping strategies. Coral Reefs 2010, 29, 329. [Google Scholar] [CrossRef]
  16. Costa, B.M.; Battista, T.A. The semi-automated classification of acoustic imagery for characterizing coral reef ecosystems. Int. J. Remote Sens. 2013. [Google Scholar] [CrossRef]
  17. Diesing, M.; Thorsnes, T. Mapping of Cold-Water Coral Carbonate Mounds Based on geomorphometric Features: An Object-Based Approach. Geosciences 2018, 8, 34. [Google Scholar] [CrossRef] [Green Version]
  18. Moura, R.L.; Secchin, N.A.; Amado-Filho, G.M.; Francini-Filho, R.B.; Freitas, M.O.; Minte-Vera, C.V.; Teixeira, J.B.; Thompson, F.L.; Dutra, G.F.; Sumida, P.Y.G.; et al. Spatial patterns of benthic megahabitats and conservation planning in the Abrolhos Bank. Cont. Shelf Res. 2013, 70, 109–117. [Google Scholar] [CrossRef]
  19. Bastos, A.C.; Quaresma, V.S.; Marangoni, M.B.; D’Agostini, D.P.; Bourguignon, S.N.; Cetto, P.H.; Silva, A.E.; AmadoFilho, G.M.; Moura, R.L.; Collins, M. Shelf morphology as an indicator of sedimentary regimes: A synthesis from a mixed siliciclastic–carbonate shelf on the eastern Brazilian margin. J. S. Am. Earth Sci. 2015, 63, 125–136. [Google Scholar] [CrossRef] [Green Version]
  20. Walbridge, S.; Slocum, N.; Pobuda, M.; Wright, D.J. Unified Geomorphological Analysis Workflows with Benthic Terrain Modeler. Geosciences 2018, 8, 94. [Google Scholar] [CrossRef] [Green Version]
  21. Lundblad, E.R.; Wright, D.J.; Miller, J.; Larkin, E.M.; Rinehart, R.; Naar, D.F.; Donahue, B.T.; Anderson, S.M.; Battista, T. A Benthic Terrain Classification Scheme for American Samoa. Mar. Geodesy 2006, 29, 89–111. [Google Scholar] [CrossRef] [Green Version]
  22. Lucieer, V.; Lamarche, G. Unsupervised fuzzy classification and object-based image analysis of multibeam data to map deep water substrates, Cook Strait, New Zealand. Cont. Shelf Res. 2011, 31, 1236–1247. [Google Scholar] [CrossRef]
  23. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2018; Available online: https://www.R-project.org (accessed on 15 February 2020).
  24. Hijmans, R.J. Raster: Geographic Data Analysis and Modeling. R Package Version 3.0-7. 2019. Available online: https://CRAN.R-project.org/package=raster (accessed on 20 December 2019).
  25. Pontius, R.G., Jr.; Santacruz, A. diffeR: Metrics of Difference for Comparing Pairs of Maps or Pairs of Variables. R Package Version 0.0-6. 2019. Available online: https://CRAN.R-project.org/package=differ (accessed on 20 December 2019).
  26. Pontius, R.G.; Millones, M. Death to Kappa: Birth of quantity disagreement and allocation disagreement for accuracy assessment. Int. J. Remote Sens. 2011, 32, 4407–4429. [Google Scholar] [CrossRef]
  27. Olofsson, P.; Foody, G.M.; Herold, M.; Stehman, S.V.; Woodcock, C.E.; Wulder, M.A. Good practices for estimating area and assessing accuracy of land change. Remote Sens. Environ. 2014, 148, 42–57. [Google Scholar] [CrossRef]
  28. Collin, A.; Ramambason, C.; Pastol, Y. Very high resolution mapping of coral reef state using airborne bathymetric LiDAR surface-intensity and drone imagery. Int. J. Remote Sens. 2018, 39, 5676–5688. [Google Scholar] [CrossRef] [Green Version]
  29. Jarna, A.; Baeten, N.J.; Elvenes, S.; Bellec, V.K.; Thorsnes, T. Semi-automatic versus manual mapping of cold-water coral carbonate mounds located offshore Norway. Int. J. Geo-Inf. 2019, 8, 40. [Google Scholar] [CrossRef] [Green Version]
  30. Li, J.; Schill, S.R.; Knapp, D.E.; Asner, G.P. Object-Based Mapping of Coral Reef Habitats Using Planet Dove Satellites. Remote Sens. 2019, 11, 1445. [Google Scholar] [CrossRef] [Green Version]
  31. Diesing, M.; Stephens, D. A multi-model ensemble approach to seabed mapping. J. Sea Res. 2015, 100, 62–69. [Google Scholar] [CrossRef]
  32. Bastos, A.C.; Moura, R.L.; Moraes, F.C.; Vieira, L.S.; Braga, J.C.; Ramalho, L.V.; Amado-Filho, G.M.; Magdalena, U.R.; Webster, J.M. Bryozoans are major modern builders of South Atlantic oddly shaped reefs. Sci. Rep. 2018, 8, 9638. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Survey area location, along the eastern Brazilian coast. The detailed frame in the top right corner shows the survey plan (28.9 km of acquisition tracklines).
Figure 1. Survey area location, along the eastern Brazilian coast. The detailed frame in the top right corner shows the survey plan (28.9 km of acquisition tracklines).
Geosciences 10 00072 g001
Figure 2. Flow diagram of the analyses.
Figure 2. Flow diagram of the analyses.
Geosciences 10 00072 g002
Figure 3. (a) DBM of the study area; (b) bathymetric profile A-B; (c) 3D view of the isolated pinnacles; (d) 2D view of the isolated pinnacles; (e) 2D view of the LRRB; (f) 3D view of the LRRB.
Figure 3. (a) DBM of the study area; (b) bathymetric profile A-B; (c) 3D view of the isolated pinnacles; (d) 2D view of the isolated pinnacles; (e) 2D view of the LRRB; (f) 3D view of the LRRB.
Geosciences 10 00072 g003
Figure 4. Products generated with the BTM toolbox: (a) BPI broad scale; (b) BPI fine scale; (c) Slope; (d) BTM classification.
Figure 4. Products generated with the BTM toolbox: (a) BPI broad scale; (b) BPI fine scale; (c) Slope; (d) BTM classification.
Geosciences 10 00072 g004
Figure 5. (a) SSS mosaic over the survey area; (b,c) detailed backscatter showing reef pinnacles; and (d,e) backscatter pattern of low relief banks as viewed from above or in 3D.
Figure 5. (a) SSS mosaic over the survey area; (b,c) detailed backscatter showing reef pinnacles; and (d,e) backscatter pattern of low relief banks as viewed from above or in 3D.
Geosciences 10 00072 g005
Figure 6. Supervised SSS data segmentation.
Figure 6. Supervised SSS data segmentation.
Geosciences 10 00072 g006
Figure 7. Map showing baseline map of reef distribution derived by manual interpretation.
Figure 7. Map showing baseline map of reef distribution derived by manual interpretation.
Geosciences 10 00072 g007
Figure 8. (a) Intersection of reef structures mapped by BTM and OBIA; (b) Overlapping of reef structures mapped by “BTM and reference map” as well by OBIA and the reference map; (c) Intersection among the three techniques; (d) Spatial distribution of all mapped reefs indicated by the two automatic classification techniques and the reference map, as well as the intersection and spatial coverage indicated by both.
Figure 8. (a) Intersection of reef structures mapped by BTM and OBIA; (b) Overlapping of reef structures mapped by “BTM and reference map” as well by OBIA and the reference map; (c) Intersection among the three techniques; (d) Spatial distribution of all mapped reefs indicated by the two automatic classification techniques and the reference map, as well as the intersection and spatial coverage indicated by both.
Geosciences 10 00072 g008
Table 1. Dictionary classes.
Table 1. Dictionary classes.
ClassesZoneBroadBPI LowerBroadBPI UpperFineBPI LowerFineBPI UpperSlope LowerSlope UpperDepth LowerDepth Upper
1Deepest Inter-reef shelf −100 −100
2East deep Inter-reef shelf −100−100100
3West deep Inter-reef shelf−100100 −100
4Inter reef structures −100
5Pinnacles/Irregular relief reef Banks (IRRB)100 100
6Pinnacles/IRRB (artifacts)100 −100100
7Edges of LRRB−100100100 213
8Low relief reef banks (LRRB) 100−100100 −32−26
9Inter-reef shelf 100−100100
Table 2. Contingency table for reef presence-absence predictions.
Table 2. Contingency table for reef presence-absence predictions.
Reef AbsenceReef Presence
Predicted absenceTrue negative (TN)False negative (FN)
Predicted presenceFalse positive (FP)True positive (TP)
Table 3. Classes adopted for each classification method.
Table 3. Classes adopted for each classification method.
Classification MethodClasses
Reef PresenceReef Absence
Morphometry
  • “Inter-Reef Structures”
  • “Pinnacles/Irregular relief reef Banks (IRRB)”
  • “Pinnacles/IRRB (artifacts)”
  • “Low Relief Reef Banks-LRRB”
  • “Edges of LRRB”
  • Deepest Inter-reef shelf
  • East deep Inter-reef shelf
  • West deep Inter-reef shelf
  • Inter-reef shelf
OBIA
  • Pinnacles
  • LRBB (Low relief reef banks)
  • IRRB (Irregular Relief reef banks)
  • Inter-reef shelf
Manual interpretation
  • Pinnacle
  • Reef bank
  • Mixed (pinnacle and bank)
-
Table 4. Results of the agreement analysis.
Table 4. Results of the agreement analysis.
MethodAgreementQuantity DisagreementAllocation Disagreement
BTM0.6910.1960.113
OBIA0.6710.2350.094
Table 5. Class-specific errors of omission and commission for the BTM method.
Table 5. Class-specific errors of omission and commission for the BTM method.
ClassError of OmissionError of Commission
Reef absence0.2970.086
Reef presence0.3790.732
Table 6. Class-specific errors of omission and commission for the OBIA method.
Table 6. Class-specific errors of omission and commission for the OBIA method.
ClassError of OmissionError of Commission
Reef absence0.3310.076
Reef presence0.3160.735

Share and Cite

MDPI and ACS Style

Menandro, P.S.; Bastos, A.C.; Boni, G.; Ferreira, L.C.; Vieira, F.V.; Lavagnino, A.C.; Moura, R.L.; Diesing, M. Reef Mapping Using Different Seabed Automatic Classification Tools. Geosciences 2020, 10, 72. https://doi.org/10.3390/geosciences10020072

AMA Style

Menandro PS, Bastos AC, Boni G, Ferreira LC, Vieira FV, Lavagnino AC, Moura RL, Diesing M. Reef Mapping Using Different Seabed Automatic Classification Tools. Geosciences. 2020; 10(2):72. https://doi.org/10.3390/geosciences10020072

Chicago/Turabian Style

Menandro, Pedro S., Alex C. Bastos, Geandré Boni, Lucas C. Ferreira, Fernanda V. Vieira, Ana Carolina Lavagnino, Rodrigo L. Moura, and Markus Diesing. 2020. "Reef Mapping Using Different Seabed Automatic Classification Tools" Geosciences 10, no. 2: 72. https://doi.org/10.3390/geosciences10020072

APA Style

Menandro, P. S., Bastos, A. C., Boni, G., Ferreira, L. C., Vieira, F. V., Lavagnino, A. C., Moura, R. L., & Diesing, M. (2020). Reef Mapping Using Different Seabed Automatic Classification Tools. Geosciences, 10(2), 72. https://doi.org/10.3390/geosciences10020072

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop