# Photogrammetric Surveys and Geometric Processes to Analyse and Monitor Red Coral Colonies

Jean-Philip Royer, Mohamad Motasem Nawaf, Djamal Merad, Mauro

Saccone, Olivier Bianchimani, Joaquim Garrabou, Jean-Baptiste Ledoux,

Angel Lopez-Sanz, Pierre Drap

► To cite this version:

Jean-Philip Royer, Mohamad Motasem Nawaf, Djamal Merad, Mauro Saccone, Olivier Bianchimani, et al.. Photogrammetric Surveys and Geometric Processes to Analyse and Monitor Red Coral Colonies.

Journal of Marine Science and Engineering, 2018, 6 (2), pp.42. 10.3390/jmse6020042. hal-03509923

HAL Id: hal-03509923

https://amu.hal.science/hal-03509923

Submitted on 30 Aug 2023

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

Article

# Photogrammetric Surveys and Geometric Processes to Analyse and Monitor Red Coral Colonies

Jean-Philip Royer <fnref n="1" />,†, Mohamad Motasem Nawaf <fnref n="1" />,†, Djamal Merad <fnref n="1" />,†, Mauro Saccone <fnref n="1" />,†, Olivier Bianchimani <fnref n="2" />,†, Joaquim Garrabou <fnref n="3" />,†, Jean-Baptiste Ledoux <fnref n="3" />,<fnref n="4" />,†, Angel Lopez-Sanz <fnref n="3" />,† and Pierre Drap <fnref n="1" />,*,†

Received: 1 February 2018; Accepted: 26 March 2018; Published: 12 April 2018

<figure><img src="image_3.png" /></figure>

check for
updates

**Abstract:** This article describes the set of photogrammetric tools developed for the monitoring of Mediterranean red coral *Corallium rubrum* populations. The description encompasses the full processing chain: from the image acquisition to the information extraction and data interpretation. The methods applied take advantage of existing tools and new, innovative and specific developments in order to acquire data on relevant ecological information concerning the structure and functioning of a red coral population. The tools presented here are based on: (i) automatic orientation using coded quadrats; (ii) use of non-photorealistic rendering (NPR) and 3D skeletonization techniques; (iii) computation of distances between colonies from a same site; and (iv) the use of a plenoptic approach in an underwater environment.

**Keywords:** underwater photogrammetry; plenoptic; red coral; skeleton; NPR; orthophoto; coded targets; climate change; conservation biology

# 1. Introduction

## 1.1. *Ecological Context and Associated Project*

Red coral (*Corallium rubrum* [1]) is considered an emblematic species for the Mediterranean Sea [2]. Red coral is an invertebrate that grows forming arborescent colonies in rocky habitats from 10 to 800 m in depth (Costantini [3]). With a lifetime of up to 100 years, this species displays a slow growth rate (<1 mm in height per year) with a recorded maximum size up to 50 cm in height [4,5]. Red coral has been fished in the Mediterranean since antiquity for its use in jewelry. As a result of this secular fishing, the current population is characterized by colonies that rarely reach more than 10 cm in height, except in protected areas or in deep water where an unexploited population can still be found.

Moreover, red coral populations are affected by global warming, leading to an increase in its mortality rate. In this context, the red coral is included in the Annex III Protocol Concerning Specially Protected Areas and Biological Diversity in the Mediterranean within the Barcelona convention

<footnote n="1">¹ Aix Marseille Université, CNRS, ENSAM, Université de Toulon, LIS UMR 7020, Domaine Universitaire de Saint-Jérôme, Bâtiment Polytech, Avenue Escadrille Normandie Niemen, 13397 Marseille, France; jph.royer@gmail.com (J.-P.R.); mohamadmotasem.nawaf@lsis.org (M.M.N.); djamal.merad@univ-amu.fr (D.M.); mauro.saccone@gmail.com (M.S.)</footnote>

<footnote n="2">² Septentrion Environnement, 13007 Marseille, France; olivier.bianchimani@septentrion-env.com</footnote>

<footnote n="3">³ Institut de Ciencies del Mar-CSIC, 37-49, E-08003 Barcelona, Spain; garrabou@icm.csic.es (J.G.); jbaptiste.ledoux@gmail.com (J.-B.L.); alopez@icm.csic.es (A.L.-S.)</footnote>

<footnote n="4">⁴ CIIMAR/CIMAR, Centro Interdisciplinar de Investigação Marinha e Ambiental, Universidade do Porto, 4050-123 Porto, Portugal</footnote>

<footnote>* Correspondence: Pierre.Drap@univ-amu.fr; Tel.: +33-4-9182-8520</footnote>

<footnote>† These authors contributed equally to this work.</footnote>
adopted by the EU through Council Decision 77/585/CEE and classified as an endangered species in the International Union for Conservation of Nature Red List (IUCN 2015).

## 1.2. Contribution

In this article, we present a set of tools as well as a full methodology dedicated to the application of different technologies for in situ monitoring of Mediterranean red coral *Corallium rubrum* population. We will describe the pipeline that was set up to characterize the structure and functioning of red coral populations; data collection, either from a multi-view approach implying photogrammetry or from a plenoptic-based method, up to the information processing. In this study, the image acquisition is generally performed by scuba divers. Since diving time is limited and hard to carry out, we chose to use a lightweight setup: a single camera with two strobes for stereoscopic acquisition. We will discuss the different approaches used during in situ surveys that also determine the type of data acquisition. Finally, we will introduce the calibration and orientation methods adopted in each case. These tools provide an estimation of dimensional or geometric characteristics allowing us to explore 3D models of single colonies using sets of photographs (2 to 10) or a plenoptic based acquisition showing the same scene and orthophotographs covering several square meters.

The tools described in this paper are designed to contribute to define efficient management and conservation of red coral populations. Bearing in mind the species' patchy distribution and depth range as well as the colonies' complexity and fragility, it is extremely difficult to acquire sound information on key features of the population's structure and functioning.

To overcome these barriers, we adopted a multidisciplinary approach involving underwater photogrammetry (new study tools and image processing) and biologists (performing in situ acquisitions, statistical and genetic analysis) able to encompass relevant temporal and spatial scales and implementing an innovative analytic tools to obtain quantitative data on the structure and functioning of red coral populations for the very first time.

We contend that this approach will largely benefit the efficiency of management and conservation plans for red coral populations. Particularly in this paper, we will address three main goals:

* Characterization of population structure (scale 1 cm to 50 cm); see Section 3.

* Population dynamics and functioning (scale m²); see Section 4.

* Characterization of colony architecture/complexity (resolution scale 1 mm to 1 cm); see Section 5.

The first goal implies enhanced data acquisition on population size structure. The second goal allows to obtain accurate maps of the surveyed areas (up to 10 m²) [6]. The third goal addresses a dimensional study, as well as estimating statistical information on common colonies and nowadays rare large-sized colonies. Due to the fragility of corals, all contact should be excluded during the surveys. Thus, the methods developed to tackle this problem should meet several requirements: user-friendly, non-invasive, remote sensing and minimal dive time.

## 1.3. *Short Description of the Surveys Associated with the Presented Tools*

In the framework of this study, we have developed a number of tools for monitoring red coral populations and tackling the inherent difficulties through three kinds of surveys that address the previous three goals:

1. The first survey aims at performing measurements on randomly chosen populations to characterize the populations. In practice, the underwater site is randomly sampled. Defining an area of interest with a homemade printed quadrat, a plastic frame of size 20 × 20 cm, a scuba diver shoots several photos (distance 20 to 40 cm) with a high overlap (50% to 70%). We opted for a lightweight and practical approach, relying on a single device (a DSLR -Digital Single-Lens Reflex camera- camera) with a housing and two electronic strobes. This multi-view procedure minimizes the diver's tasks and enables them to easily adapt the measurement procedure to
the scene configuration. These photos are oriented in a global frame with an excellent accuracy, thanks to the coded targets drawn on the quadrat [7]. Measurements are directly performed on the photographs using a dedicated graphical interface and the corresponding 3D points are computed. They provide reliable results about the colony size (base diameter, height). Comments on the occurrence of necrosis, for example, can be added.

2. The second kind of survey enables the tracking of temporal changes in red coral populations. It gives a fine scale spatial distribution of red coral populations. For these surveys, permanent plots were setup, which allowed for sampling the same areas through time. The plots were marked with PVC (PolyVinyl Chloride) or metal screws embedded in holes of the rocky substratum with epoxy putty. The total covered area varies between 1 and 6 m². That area is large enough to monitor hundreds of colonies over time. We used different DSLR cameras embedded in a housing and including two electronic strobes to photograph the transects (from Nikon D70 to D700). We implemented two different sampling methods. The first one consisted of setting a cord attached to the screws to define the plots in each sampling period. Quadrats with coded targets were then positioned above and below the cord as reference. At each position, a few photos were shot and a local coordinate system was defined. The second method consisted of obtaining orthophotos from the permanent plots. Since, in order to build orthophotos, moving a quadrat including coded targets is not possible, we use a global orientation of the scene and for studies at different times, we rely on the whole transect, which always remains fixed, to keep the same coordinate system and compare orthophotos. If biologists really need to sample individual colonies and visualize at the same time the standard of size defined by the quadrat (see Section 3.1), it is still possible to move a virtual 20 × 20 cm quadrat in the scene, using our dedicated application.

3. The third type survey implemented was designed for the characterization of emblematic large colonies that have a high natural heritage value (last goal of the paper). The complex shape exhibited by large colonies as well as the entanglement of the branches lead to several problems. Likewise, the proximity of the substrate on which the colony is attached may disturb the analysis. The survey is designed to acquire and build a 3D model first in several ways and then perform automatic measurements on the colonies and determine relevant quantitative morphological information (e.g., like the number of branches). Non-photorealistic rendering based post-processing to smooth the shape, as well as skeletonization methods, were developed.

### 1.4. *Structure of the Document*

The paper is structured as follows: first, in Section 2, we will present recent works putting in relation photogrammetry and measurements on living organisms or on the seabed directly. In Sections 3–5, each survey will be more precisely described (the photogrammetry approach to measure dimensional data on colonies, the photogrammetric dimensional change monitoring to track the evolution of the populations and finally a study to classify large colonies using plenoptic and photogrammetric approaches). Finally, we will conclude with the contributions brought by the survey methodology developed, based on photogrammetric and plenoptic approaches.

## 2. Red Coral Surveys

### *Use of Photogrammetry for Marine Benthic Communities*

The surveys applied in this study on red coral are based on standard sampling methods used in marine ecology using quadrats and transects to quantify the abundance and diversity of coral reef communities. Transects are rectangular areas defined by a line deployed over the communities defining the length and quadrats or bars that help to set the width of the transect, the transects can be permanent when marks are setup in the substratum as reference or random if the transects are set in different areas for each sampling. Quadrats are plastic frames allowing for discretization of the observed surface in situ. For instance, Caroline S. Rogers in [8] cites this approach in 1983 and Arthur
Lyon Dahl refers to exploited transects since 1917. The use of permanent marks allows observations of change over time. This technique was used for example by Caroline S. Rogers to observe the impact of Hurricane David in 1979 [8].

Indeed, Arthur Lyon Dahl mentioned in [9]: "In 1917, Alfred G Mayor, made a detailed quantitative survey of the fringing reef in Aua in Pago Pago Harbour, American Samoa [10]. Mayor laid a transect across the reef, and counted all the corals and certain other animals within a series of 7.32 m quadrat along the line from the shore to the outer edge of the reef. The transect was repeated using an identical method in 1973 [11] and again in 1980 by the author. Two markers placed on the line in 1973 were still present in 1980, assuring the accurate position of the line". In this study, we apply the same methods whereas now photogrammetry allows us to be more accurate and exhaustive.

Biologists have integrated photogrammetric computer vision in their fields of research for some years to measure or count underwater living animals or work directly with the seabed. Sanchez et al. [12] used a light photogrammetric approach in shallow water, while Ref. [13] preferred a specifically designed autonomous underwater vehicle (AUV) and Refs. [14,15] relied on specific stereoscopic photography equipment in deep water habitats. Finally, many studies applied photogrammetry to measure benthic communities [16–19] and also [20]. Shortis [21] provided a clear and synthetic overview.

Underwater photogrammetry has a drawback since geolocation can be difficult. However, the need to use a geographic information system (GIS) is still a reality. Actually, a local coordinate system can be defined [22,23]. Photogrammetry and GIS approaches can be used to characterize reproductive units within red coral populations by studying family relationships by combining spatially explicit positions with genetic analysis [24]. Contrary to traditional methods used in marine biology, which can be destructive and difficult to be implemented, photogrammetry is non-invasive and can be performed in situ [25]. Moreover, the emergence of powerful photogrammetry processing software, like *PhotoScan* [26], has changed the way to handle these studies. This approach is now more user-friendly for non-specialist researchers [19]. Finally, photogrammetry is being perceived as a useful and affordable approach for an increasing number of research teams. See, for example, Refs. [27–29] or the study of structural complexity of underwater seafloor [30–32]. GIS software is also coupled with photogrammetry in [33] for a deep coral survey.

Some stereo tools, diver-operated or embedded in an AUV, have been developed since some years ago, in the field of underwater robotics in order to easily survey a site using photogrammetry with specific methods to manage the stereo systems and local bundle adjustment. These developments were intended to support inclusion of bathymetry in an underwater archeology context as well as in marine biology (see [32,34–36]). Similarly, we have also developed our own low-cost system, which is easier to use in complex environments [37]. The system is based on embedded and distributed low-power mini computers (Raspberry Pi). A set of software programs developed by our team ensures a general good precision with local bundle adjustment and real-time odometry. Moreover, a bridge was done with the commercial application Photoscan developed by Agisoft (Saint-Petersbourg, Russia) [26].

Today, many new and complex applications merge several media, lidar and photogrammetry as well as using mobile systems such as aerial drones, remotely operated underwater vehicles, etc. [38]. However, what we can see in this brief overview of methods used to make an underwater survey of coral reef colonies is that, even if the use of photogrammetry has dramatically increased over the last five years, the main goal of these surveys can be achieved only by human observations.

Indeed, the knowledge required to interpret the information is the key point of the survey. We will focus in this paper on several methods and techniques to automate partially the 3D acquisition and the morphological information acquisition. Compared to the previous works quoted in this section, we offer here different dedicated interfaces (e.g., user-friendly graphical user interface, implementation of photogrammetry algorithms, development of our own bundle adjustment, including specific constraints on distances on the quadrat to benefit from an optimal orientation quality...). The final survey should be done by the biologists and the main contribution of the approach described in this
paper is to be able to minimize the time spent underwater (data acquisition) and to be able to optimize
the analysis of the accurate and quite exhaustive 3D data (in several formats: 3D models, point clouds,
color map) in the laboratory.

3. Semi-Automatic Coral Measurements and Quadrat Sampling

In order to characterize population size structures, we developed a semi-automatic processing
chain based on photogrammetry. First, colonies are randomly sampled using a 20 × 20 cm plastic
quadrat. This quadrat will be used both to scale the future model and set up a coordinate system.
The whole process can be divided into five steps.

3.1. *Setting Up the Quadrat*

The quadrat is a homemade 3D printed plastic frame that can be easily made by anyone that
has the 3D mesh file. This mesh file is designed to comply with the standard used by biologists.
The frame's internal dimensions were fixed to 200 × 200 mm. The external dimensions were chosen to
ensure a good visibility of the targets: 230 × 270 × 15 mm. It incorporates two kinds of coded targets.
First, 18 classical circular targets used by Photoscan and Photomodeler software—second, 18 square
targets from the TagRuler system, which we have developed for our own photogrammetric platform
ARPENTEUR [7]. The presence of these two kinds of targets enables processing the photos both in the
commercial software and in our own photogrammetry modules. The 18 targets of each type are evenly
distributed on the frame. Their location is precisely known and they can be used as control points.
The shape of the quadrat itself, designed with CAD (Computer-Aided Design) techniques, takes into
consideration underwater constraints; no air pockets, thickness of 2 mm and cellular framework to
balance stiffness, weight and handling (see Figure 1).

<figure><img src="image_4.png" /><figcaption>**Figure 1.** Quadrat design and realization using a two-color 3D printer.</figcaption></figure>

3.2. *Shooting Methodology*

In practice, the printed quadrat is placed on the seabed (see Figure 2), or attached to buoys
(embedded in the quadrat) on a cave ceiling or sometimes held in place by another diver (see Figure 3)
who keeps it stationary with respect to the seabed. Then, between two and ten photographs are shot
(but more generally, three are enough), so that the common covered area and the base are maximized
(see Figure 4). Due to the difficult constraints of manned underwater acquisition (limited time, limited
visibility, instability, etc.), this approach is easy to carry out and minimizes the number of photos and
the processing time (see the next subsection).
<figure><img src="image_5.png" /><figcaption>**Figure 2.** Photograph using the quadrat in situ (Photo: Olivier Bianchimani).</figcaption></figure>

<figure><img src="image_6.png" /><figcaption>**Figure 3.** Diver holding the frame during the acquisition step (This approach is not recommended). The orientation phase is automatically performed using only the 36 coded targets printed on the quadrat. A post processing analysis of the residuals computed on measurements done manually will detect this problem (Photo: Olivier Bianchimani).</figcaption></figure>

<figure><img src="image_7.png" /><figcaption>**Figure 4.** Three photos of a quadrat containing colonies. They are shot with a convergent axis, maximizing the common covered area and the base (Photo: Olivier Bianchimani).</figcaption></figure>

### 3.3. *Automatic Orientation*

The photo orientation procedure is used to determine the pose (position and orientation) of the camera used for each capture (called external or exterior orientation). This is essential to be able to estimate and measure the 3D structure. In contrast to structure-from-motion (SfM) scenarios, no camera calibration is required in advance. However, we mention that using calibration parameters avoids
introducing more degrees of freedom in the optimization phase and leads to more robustness in our case. The camera parameters (referred to as interior or internal orientation) are estimated jointly with the exterior orientation in a closed-form manner, more accurately than self-calibration approaches.

Each single camera shot is enough to recover the interior and exterior orientation, and also the distortion coefficients, both radial and tangential. This is done based on Tsai's camera calibration scheme [39] adapted to planar targets as summarized below.

The photo orientation procedure is applied to each captured image separately at first, and then the recovered parameters are used as initial values in a further optimization phase. We consider that each photo contains the image of the quadrat with TagRuler or circular markers (with the region of focus at the center); here, five points are enough to recover all the parameters, but more points provide better accuracy. A robust TagRuler and circular marker detection are developed based on the AprilTags open-source project. The detected points together with their known relative position serve as input to the recovery procedure. Formally, using a pin-hole camera model, a 3D point $P_c = [X_c Y_c Z_c]^T$ in a camera coordinate system is projected to a 2D point $p_I = [x_I y_I]^T$ in image coordinates as

$$x_I = f \frac{X_c}{Z_c} \quad \text{and} \quad y_I = f \frac{Y_c}{Z_c}, \qquad (1)$$

where $f$ is the focus distance (in general, the ratio of the pixel dimensions parameter must be considered. It is not taken into account here as most modern cameras have square pixels). Using the ratio of both equations yields a result that is independent of the focus distance and distortion

$$\frac{x_I}{y_I} = \frac{X_c}{Y_c}. \qquad (2)$$

Transforming $P_c$ from world coordinates is achieved with the mean of a rotation matrix $R$ and translation vector $t = [t_x t_y t_z]^T$ (external parameters)

$$P_c = RP_W + t. \qquad (3)$$

By considering the origin of world coordinates at one corner of the quadrat, we can consider $Z_W = 0$ for all control points. After the expansion of the last equation and the substitution in Equation (2), we obtain:

$$\frac{x_I}{y_I} = \frac{r_{11}X_W + r_{21}Y_W + t_x}{r_{21}X_W + r_{22}Y_W + t_y}, \qquad (4)$$

where $r_{ij}$ is the element at row $i$ and column $j$ of $R$. This yields after simplification to a linear equation with six unknowns; $r_{11}, r_{12}, r_{21}, r_{22}, t_x, t_y$. For each 2D-3D correspondence point, we obtain one such equation. Using a least-squares approach to solve this equation, the obtained solution is up to scale. By applying the orthogonality constraint of a rotation matrix, we could simultaneously recover the five remaining elements using orthonormal equations. Hence, the focal distance and $t_z$ can be recovered using Equations (1) and (3).

The next step is to use the computed parameters as initial values in a nonlinear optimization phase. Since the origin of the coordinates is considered to be the same for all photos, the recovered external orientations are passed to the optimization step untouched. The focal distance, which is the only recovered internal parameter, is averaged for all photos. The principle point is initialized at the center of the image. Finally, we optimize the values for all external, internal, and distortion parameters using the Levenberg-Marquardt method and the reprojection error:

$$\sum_{N}^{i=1} (x_{li} - x_{pi})^2 + \sum_{N}^{i=1} (y_{li} - y_{pi})^2, \qquad (5)$$
where $(x_P, y_P)$ are the coordinates of the re-projected point using estimated parameters. Considering that the working distance is approximately 200 mm, the residuals obtained on the control points are generally less than 1 mm.

## 3.4. Red Coral Colonies' Measurements

We developed a fully functional graphical user interface (see Figure 5) that enables us to perform measurements on the colonies inside the quadrat. The interface allows the presentation of two photos side by side (chosen by the operator from the set of photographs taken on-site and regarding these colonies). The user selects homologous points on the photos, so that the resulting 3D point is computed by triangulation. The interface performs real-time controls to ensure the integrity of the data, for example, correlation to ensure the quality of the chosen pair of points. Finally, the user has the possibility to determine various information including the diameter of the base, maximum height and approximate volume as well as visual observations, such as the percentage of necrosis for each colony and the density of colonies and recruits.

A set of tools is available within this software and allows for constraining the homologous points determination, displaying the epipolar line, performing automatic cross correlation, and contrasting evaluation, local zoom, dynamic residual visualization, etc.

<figure><img src="image_8.png" /><figcaption>**Figure 5.** Application dedicated to red coral colonies' measurement. In this example, the viewpoint is skewed, but the measurements still obtain good accuracy, showing that the application is robust to acquisition approximations.</figcaption></figure>

## 3.5. Computation and Results

The results are calculated in a matter of a few seconds. The results are exported into standard widely used file formats (XML, Excel, PLY and VRML to visualize the point cloud resulting from the measurements) to fit the tools already used by biologists for their analysis.

Very recently, photogrammetric results and coral measurements have also been exported according to an ontology developed in the LIS laboratory in the framework of the ARPENTEUR project [40].

Below, we provide numerical results on an in situ study case to assess the quality of the measurements performed. We shot three photographs with a Nikon D300 (focal 20 mm, resolution: 4288 × 2848). The pictures are shown in Figure 6.
<figure><img src="image_9.png" /><figcaption>**Figure 6.** Three photos of a quadrat (test purpose) in an underwater environment, Scandola, Corsica. (Photo: Olivier Bianchimani).</figcaption></figure>

Two kinds of analysis are carried out: first, the reprojection residual (pixel unit) on the 36 coded targets present on the quadrat; next, the difference between the computed distances in mm between all the targets and the theoretical ones given by the coordinates of the control points (called an *accuracy*). For each case, the different error indicators computed are the average, the root mean square error (RMSE) and the median error. The results are shown in Table 1.

<table><caption>**Table 1.** Residual analysis.</caption><thead><tr><th></th><th>Reprojection Residual (Pixels) for 36 Ground Points</th><th>Residual on Distances (mm) for 544 Computed Distances</th></tr></thead><tbody><tr><td>Average</td><td>1.10</td><td>0.34</td></tr><tr><td>Median</td><td>1.07</td><td>0.41</td></tr><tr><td>RMSE</td><td>1.21</td><td>1.64</td></tr></tbody></table>

The reprojection residual, which we call a *precision*, are in the range of one pixel, which is a fairly good result, considering the additional constraints added underwater (mainly, distortion correction).

Concerning the distances, considering that the interior of the quadrat is a 20 cm square, the average distance error computed is less than 0.5 mm.

## 4. Estimation of the Temporal Changes and the Colony Spatial Distribution

The goal of this study is to determine the changes over time in the red coral populations. This implies both studying colony changes over time and also correlating genetic studies with the spatial position of colonies to better understand the mechanisms involved in the populations functioning (reproduction and recruitment).

Initially, we used a different process that relied on surveys of permanent plots, set up using PVC screws fastened to holes in the (rocky) substratum with epoxy putty. The pattern of this grid (number of transects and distances between the screws) varied, depending on the kind of substratum and microtopography. Then, areas were squared off using elastic cords to define normalized sampling areas. Then, transects were photographed, using different cameras (Nikon D70, Nikon D300) within housing and two electronic strobes. On each sampling area, the 200 × 200 mm quadrat was positioned (above and below the cord length) and two or three photos were shot using slightly different angles (about 30°). Notice that the covered area generally did not exceed 2 to 6 m². Nevertheless, it enabled us to monitor hundreds of colonies.

This acquisition process was heavy to carry, slow and limited, considering that it did not enable us to efficiently monitor spatial changes on the colonies.

Thus, we decided to use a different acquisition method that consists of positioning the quadrat at a given place. The first step is to shoot a large collection of photos of the area (generally, hundreds to thousands). These photos are taken by advanced divers, who efficiently cover the area and keep the same distance between them and the colonies, so that the light variations are minimized. Next, the set of photos can be oriented in a global system of coordinates given by the quadrat (here, global does not refer to a world geographic reference system, but rather means that we have a unique coordinate system for all the covered area, scaled thanks to the quadrat), which leads to a sparse point cloud. At the same
time, the scale is also applied. The next step computes the dense cloud. Alternatively, this dense cloud can be meshed, to produce a high-resolution 3D model. Finally, the most important step is building the orthophoto for each sampled area. From a suitable plane of projection, the previous 3D model is projected using an orthographic projection. This projection conserves distances, so the resulting orthophoto has the same scaled coordinate system attached as the one used by the mesh and the point cloud. During the monitoring process, biologists require estimating distances between colonies and quantifying temporal changes. Considering the goal of this study, orthophotos perfectly fit their needs. They offer a global view of the area (see Figure 7 for an example), with a very high resolution, while pixels are related to the true coordinates. Moreover, it is easier to work with orthophotos than with 3D models, which are not as handy to use.

<figure><img src="image_10.png" /><figcaption>**Figure 7.** Zoom on a global orthophoto. Original pixel dimensions are 21,000 × 15,000 (1 mm/pixel).</figcaption></figure>

## 4.1. Orthophoto Measurements

As a consequence, we developed an efficient dedicated software enabling us to load an orthophoto and pick points on it. It offers the possibility to perform measurements on selected points. The user can directly select interesting points on the orthophoto. The upper part of the application shows mouse coordinates in world units, while the bottom part summarizes the set of measurements already done. This point list can be managed (add, remove, edit, etc.). Moreover, in order to handle large-dimension orthophotos, they are divided into small blocks, which are merged in a transparent way for the final user. Finally, the distance map between all points can be computed. We provide an example of the user application in Figure 8.

We tried different methods to provide distance estimations. Due to the ruggedness of the seabed, Euclidean distances led to poor results, as they did not take into consideration the sea floor's topography.

Several methods have been developed to compute exact or approximations of the shortest path on graphs. Dijkstra developed a well-known algorithm to solve this problem in [41]. On a triangular mesh, this algorithm computes the optimal path that follows the ridges of the triangles, but it is slow on large meshes, as it needs to test every vertex present.

Dijkstra's algorithm enables a good approximation of a distance along the curve of the ground, provided the mesh is dense enough. However, there are two main problems. First, the computation time depends on the size of the mesh. Some improvements can be found in the literature to speed up the computation, like the A\* algorithm [42]. The counterpart is that the resulting path may not be the optimal one. In our case, we use a slight improvement of Dijkstra's algorithm: the list *S* is filled at runtime with the vertices discovered and not studied yet, instead of being totally filled at the beginning. The second problem concerns the computed distance, which can be over-estimated in some cases.
As the path follows the edges of the triangles, the meshing pattern can have a noticeable influence on the result. To provide the best fit of the desired distance, we decided to move to a geodesic-based algorithm (Algorithm 1). In this case, the estimated curve goes through the faces of the triangles, so that the meshing pattern and even the mesh density have a lower influence on the accuracy of the result. Several methods approximating the geodesic curve exist in the literature [43–45]. The one we have developed is based on the Darboux frame [46]. The method is initialized with the path estimated by the Dijkstra's algorithm. Then, the vertex positions are slightly moved to fit the geodesic curve. This is a good choice, as the path returned by Dijkstra is close to the geodesic one.

**Algorithm 1: Geodesic algorithm**

**Data:** Input data are the point list given by the computation of the Dijkstra path L, composed of a set of N points coinciding vertices $\mathbf{m}^{(i)}$ from the mesh, such as: $L: \{\mathbf{m}^{(i)}, 0 \le i \le N\}$.

**Result:** A list containing the geodesic path sampled.

**while** current iteration $\le$ max iteration **do**

Sampling of the Dijkstra path $L_s: \{\mathbf{d}_s^{(i)}, 0 \le i \le S\}$ to obtain S values interpolated using Spline curves;

**for** each sampled point $\mathbf{d}_s^{(i)}$ **do**

Computation of the interpolated normal $\mathbf{n}_s^{(i)}$;

Computation of the tangent vector to the path, $\mathbf{t}_s^{(i)}$;

Computation of the geodesic normal vector $\mathbf{g}_s^{(i)}$, which is both tangent to the path and orthogonal to the vector $\mathbf{t}_s^{(i)}$: $\mathbf{g}_s^{(i)} = \mathbf{n}_s^{(i)} \wedge \mathbf{t}_s^{(i)}$, where $\wedge$ is the cross product.

**end**

The Darboux formulas put in relation the previous three base vectors with three scalar quantities, called the normal curvature $\gamma_n$, the geodesic curvature $\gamma_g$ and the relative torsion $\tau_g$. This relation can be modeled with the following matrix form:

$$ \frac{d}{ds} \begin{pmatrix} \mathbf{t}_s \\ \mathbf{g}_s \\ \mathbf{n}_s \end{pmatrix} = \begin{pmatrix} 0 & \gamma_g & \gamma_n \\ -\gamma_g & 0 & -\tau_g \\ -\gamma_n & \tau_g & 0 \end{pmatrix} \begin{pmatrix} \mathbf{t}_s \\ \mathbf{g}_s \\ \mathbf{n}_s \end{pmatrix}. \quad (6) $$

We are interested in determining the geodesic curvature. Indeed, the geodesic path is given by [47]:

$$ \frac{dL_s}{dt} = -\gamma_g \mathbf{g}_s. \quad (7) $$

To determine $\gamma_g$, Equation (6) is rewritten:

$$ \frac{d}{ds} \begin{pmatrix} \mathbf{t}_s \\ \mathbf{g}_s \\ \mathbf{n}_s \end{pmatrix} = \begin{pmatrix} \mathbf{g}_s & \mathbf{n}_s & 0 \\ -\mathbf{t}_s & 0 & -\mathbf{n}_s \\ 0 & -\mathbf{t}_s & \mathbf{g}_s \end{pmatrix} \begin{pmatrix} \gamma_g \\ \gamma_n \\ \tau_g \end{pmatrix} = \mathbf{Cv}. \quad (8) $$

Inversion of matrix C leads to the determination of $\gamma_g$. The geodesic path is now computed for each sample point thanks to the following finite difference equation:

$$ \mathbf{d}_s^{(i+1)} = \mathbf{d}_s^{(i)} - c \operatorname{sign}(\gamma_g) \mathbf{g}_s, \quad (9) $$

where c is a low valued constant (e.g., c = 0.001) and sign(a) is a function that returns 1 if $a > 0$, -1 if $a < 0$ and 0 otherwise. The derivatives are numerically estimated.

**end**
To determine the 3D coordinates of the starting and ending points selected from the orthophoto, we have to do a few calculations. Knowing the plane of projection used to project the 3D model (generally the XY plane to simplify), we trace a ray from the (x, y) chosen point, with a direction equivalent to the normal of the plane (Z-axis if the plane is XY). The intersection between the ray and the triangles of the mesh corresponds to the requested 3D point. This search can be time-consuming, as it depends on the mesh size, so we decided to accelerate ray-triangles intersections using a primitive-based partitioning algorithm. While different methods exist in the literature, they are mainly divided into two categories: space partitioning algorithms, including octrees, kd-trees,... (see for example [48,49]) and primitive-based partitioning algorithms ([50,51]). This second category especially includes BVH (Bounding Volume Hierarchy). Recently, this approach has become increasingly popular in the ray-tracing community. It is, for example, used in the highly optimized library *Embree*, developed by Intel (Santa Clara, CA, USA) [52,53]. This is also the choice we made, as it is not really difficult to implement.

Finally, the key advantage of the presented tool is that it offers a handy graphical interface, enabling tracking of coral colonies and performing accurate 3D distance computations, while hiding the complex part of this computation step from the end user.

<figure><img src="image_11.png" /><figcaption>**Figure 8.** Presentation of the application that was developed to perform distance measurement from orthophotos. Notice that the quadrat on the left is an old version. The white transects are perfectly visible here.</figcaption></figure>

## 4.2. Visual Temporal Changes

Another useful tool was developed to visualize the evolution of the populations over time. It takes advantage of the fact that transects remain in place for the photogrammetric studies at different times. Thus, the orthophotos can be oriented in the same coordinate system. An example is provided in Figure 9. Moreover, the weighting value to blend the two images can be changed on the fly.
<figure><img src="image_12.png" /><figcaption>**Figure 9.** Blending of two orthophotos, with an equivalent weight. Notice that the transects are slightly shifted, as the strings move with the currents.</figcaption></figure>

**5. Inventory of Red Coral Colonies**

The third goal of the paper is to develop tools to characterize structurally complex red coral colonies in view to build a library of emblematic colonies. The tools developed modeling the colony and extracting from it relevant information, like morphological ones. Large colonies are particularly noticeable, as they exhibit complex shapes with multiple branches. Their structure and the acquisition context make this problem very challenging both for a human or an automatic-based analysis. We will show in this section the adopted methodology and solutions leading to appropriate and workable results. We rely on different kinds of tools to provide data from different sources. The global process can be divided into two steps: data acquisition and extraction of relevant information. The acquisition was carried out in situ, for the majority, but also in the lab on dead coral colonies that are held in public or private collections.

5.1. *Photogrammetric Acquisition*

The classical photogrammetric approach is achieved by shooting hundreds of photos, with a high
overlap. As said previously, the density of large colonies multiplies hidden areas. A photogrammetric
acquisition relies on the redundancy of data, ensured by the redundancy of the areas' shots. The more
complex the colony, the more photos we should take and the less missing data we get in the
reconstruction process. Moreover, red coral colonies grow in dark and narrow areas, like cave ceilings,
which present a complex topography. This makes the photogrammetric acquisition process challenging
and time-consuming. In a practical way, it is generally performed by a scuba diver. We tested both
HD (Olympus TG2 and Canon G1) and reflex cameras (Nikon D700 and Nikon D300). The whole
process enables computation of the 3D point cloud, and, next, the mesh based on the dense cloud. See
examples of acquired large colonies in Figures 10 and 11.

<figure><img src="image_13.png" /><figcaption>**Figure 10.** A large red coral colony (more than 50 cm in width) exposed for sale in a shop. Two days of photogrammetric processing on the 800 photos were required for its analysis.</figcaption></figure>

<figure><img src="image_14.png" /><figcaption>**Figure 11.** In situ acquisition of a dense colony; 50 shots, precision on targets close to 1 mm.</figcaption></figure>
## 5.2. Non-Photorealistic Rendering (NPR)

The next step consists of extracting relevant data from the mesh modeling the colonies. We first developed an automatic tool based on NPR (non-photorealistic rendering) [54]. NPR is a method that highlights contours, ridges or valleys and enables detection of noticeable points. The main advantage is that it uses information that is not available on images: the surface normals. Basically, the principle reads as following: consider a 3D model observed from a given point of view. For each vertex, we define the normalized vector **v** to be the vector from **t** to the viewpoint. We call **n** the normalized normal at vertex **t**. Then, we can compute the following scalar quantity resulting from the dot product:

$$I(\mathbf{t}) = \mathbf{v}(\mathbf{t}) \cdot \mathbf{n}(\mathbf{t}). \qquad (10)$$

Ideally, the set of vertices such as $I(\mathbf{t}) = 0$ represents the contour. To be less restrictive due to the acquisition noise and to extract more parts of the relief, it is a good idea to add a threshold, say *th*, so that we can keep suggestive contours. These suggestive contours are then defined by:

$$I(\mathbf{t}) < th. \qquad (11)$$

More formally, the suggestive contours are defined as the set of minima of $\mathbf{v} \cdot \mathbf{n}$ in the direction of $\mathbf{w}$, where $\mathbf{w}$ is the projection of the view vector onto the tangent plane. An NPR representation of a colony is given in Figure 12.

Hence, the NPR processing we use corresponds to the 3D contours from a given point of view. One way to obtain an accurate representation of these contours consists of oversampling the resulting image of contours, with a ray-tracing based method. For each contour pixel, a ray that passes through the pixel is cast from the camera. The intersection with the 3D mesh is tracked. By subdividing the pixels, an oversampled representation of the contours of the 3D model is obtained.

From four different viewpoints, the global shape of the model is extracted. This shape can be meshed again, so that we obtain a simplified representation of the colony, where the noise resulting from the acquisition has been smoothed. This process is important for the next step.

<figure><img src="image_15.png" /><figcaption>Figure 12. Example of an NPR representation of a large colony.</figcaption></figure>

## 5.3. Skeletonization

### 5.3.1. Getting the Skeleton

Manual extraction of morphological information on large colonies, such as the number of branches or diameter of each node, is a difficult task, even if it is not directly carried out in situ, as the populations show a really complex shape. Furthermore, working on the photos resulting from the acquisition process would not lead to good results, even in the case of an automatic process, because of the
occultations. This also motivated our choice to use an NPR-based approach, as it is the best way to determine the true contours, taking into account that some branches hide the rest of the geometry.

We propose in this section an automatic process for the extraction of morphological data. Now that we have an accurate and well-defined model of the coral colony, the idea is to apply skeletonization algorithms to extract the shape and obtain an analytic representation. With the skeleton, we will be able to localize the different parts of the colony and automatically count the number of branches or other associated dimensional data. Many methods for skeletonizing shapes exist in the literature. Few methods lead to a 3D skeleton, as the majority of the works are based on images [55–58]. Historically, the concept of medial skeletons was introduced for 2D shapes. It was primarily a way to reduce the amount of data necessary to describe the topology or symmetry of the shape. The concept was later extended to 3D shapes [59]. This yields new difficulties, as the desired skeletons are more rich and the shapes more complex, making this problem challenging.

In a general way, a medial skeleton meets several definitions. For a non-exhaustive list, we can write:

* The medial skeleton is the set of ball centers with the maximum radius inscribed in the shape.

* Setting a fire that propagates from the boundary to the center of a mesh, the medial skeleton is the set of points where the fire fronts meet.

* Reformulating the last definition in terms of distances, the medial skeleton is the set of locations with more than one closest point on the boundary of the shape.

For 3D shapes, the best families of skeletons that both extend the concept of wireframe of an image for volume and could be directly used on corals are the curve-based skeletons.

We decided to implement the method described in [60]. Liu et al. developed a robust algorithm to skeletonize noisy mesh that preserves connectivity. Their work is based on iteratively thinning on cell complexes. A cell complex is a closed set of k-cells. For example, a cube is a 3-cell, a rectangle is a 2-cell, its edges are 1-cell, and its vertices are 0-cell. An easy way to build the cell complex from a triangular mesh is to use classical methods to first turn it into a voxel representation. Then, 0-cells are all the vertices, 1-cells are the all the edges and so on. Each kind of cell enables characterizing a particular shape. 1-cells are useful at describing tubular shapes, which perfectly fits our needs. Regarding the definition of a medial skeleton in a thinning-based algorithm, it corresponds to the simulation of the way a fire propagates. The key in their method is to compute a global indicative value, the medial persistence, which characterizes how long a discrete cell persists in an isolated form. We briefly describe how the skeleton is created below:

The algorithm takes place in two steps. First, an iterative thinning process is performed and the medial persistence is computed for each cell $\delta$. Two scores quantify this indicator:

$$M_{abs}(\delta) = R(\delta) - I(\delta), \qquad (12)$$

$$M_{rel}(\delta) = 1 - \frac{I(\delta)}{R(\delta)}, \qquad (13)$$

where $I(\delta)$ is the number of the first iteration where the given cell ($\delta$) becomes isolated and $R(\delta)$ is the number of the iteration for which the cell is removed. In a second step, we keep the k-cells that score higher (now called medial cells) than the two defined thresholds $\epsilon_{abs}^k$ and $\epsilon_{rel}^k$. Practically, they were chosen to: $\epsilon_{abs}^k = 0.05L$, where L is the width of the bounding box and $\epsilon_{rel}^k = 0.5$. Finally, thinning is performed while keeping the medial cells. Examples of resulting skeletons are shown in Figures 13 and 14.
<figure><img src="image_16.png" /><figcaption>**Figure 13.** (left) superposition of a colony and its resulting skeleton, using the viewer from [61]; (right): advanced results, including skeleton smoothing and diameter estimation for the branches at node positions (sphere size).</figcaption></figure>

<figure><img src="image_17.png" /><figcaption>**Figure 14.** (left) superposition of another colony and its resulting skeleton, using the viewer from [61]; (right) advanced results, including skeleton smoothing and diameter estimation for the branches at node positions (sphere size).</figcaption></figure>

### 5.3.2. Extracting Relevant Information from the Skeleton

The next step is to build the tree that analytically describes the geometry of the polyps, from the set of lines of the skeleton. Each node of the tree corresponds to a junction between two branches. The distance between the nodes is also registered. Practically, from a starting point on the base, we iteratively read through the connected lines to this point. Each time three or more lines are connected to the same point, it is registered as a node and added to the tree. Once we have a structured representation of the skeleton, we can use it to extract new dimensional data. The first analysis that was performed was the estimation of the diameter of the branches. As branches are not perfect circles or ellipses and are also subject to noise, we decided that the best choice to evaluate the circumference would be the sphere that best fit the mesh centered on each node location. Practically, this is achieved
in an iterative way, by starting from a small sphere and making it grow until it reaches the mesh. As this latter is not perfectly circular, a criterion is added to stop the algorithm when a certain amount of faces has been intersected.

## 5.4. First Results

Using a Nikon D700, we took 118 photos of a dead coral colony placed on a desktop. Coded targets were positioned on the sides to scale the final model. After reconstruction, we obtained an accurate mesh containing about 450,000 faces. We first proceeded with the skeletonization as shown in Figure 15. Next, the graph describing the colony is built, as well as the spheres representing the local diameter at the node location. The result is visible in Figure 15. From these data, many extractions and interpretations are possible. For example, we decided to compute the total length of the branches and the approximate volume, using truncated cones between the arcs linking the spheres. Consequently, the resulting obtained length is $L = 0.53$ m and the approximated volume is $V = 1.28 \times 10^{-5}$ m³.

<figure><img src="image_18.png" /><figcaption>**Figure 15.** (top left) the resulting skeleton; (top right) original model; (bottom) Resulting simplified skeleton with local volume estimation.</figcaption></figure>

## 6. Plenoptic Approach

In traditional photogrammetry, it is necessary to completely scan the target area with many redundant photos taken from different points of view. This procedure is time-consuming and frustrating for large sites. Whereas many 3D reconstruction technologies are being developed for terrestrial sites relying mostly on active sensors such as laser scanners, solutions for underwater sites are much more limited. The idea of a plenoptic camera (also called light field camera) is to simulate a two-dimensional array of aligned tiny cameras. In practice, this can be achieved by placing an array of a micro-lenses between the image sensor and main lens as shown in Figure 16. This design allows the camera to capture the position and direction of each light beam in the image field, enabling the
scene to be refocused or to obtain different views from a single lens image after image capture. The first demonstration dates from 2004 by a research team [62] at Stanford University.

<figure><img src="image_19.png" /><figcaption>**Figure 16.** Diagram showing plenoptic camera internal design.</figcaption></figure>

In a plenoptic camera, each set of two close micro-lenses can be considered as a stereo pair if parts of the same scene appear in both acquired micro images. By applying a stereo matching method on the feature points of these micro images, it is possible to estimate a corresponding depth map. It is essential during image acquisition to respect a certain distance to the scene to remain within the working depth range. This working depth range can be enlarged but may affect the accuracy of the depth levels. In practice, the main lens focus should be set beyond the image sensor, in this case, each micro lens will produce a slightly different micro image for the object in regards to its neighbors. This allows the operator to have a depth estimation of the object. Recent advances in digital imaging resulted in the development and manufacturing of high-quality commercial plenoptic cameras. Main manufacturers are Raytrix (Kiel, Germany) [63] and Lytro (Mountain View, CA, USA) [64]. While the latter concentrates on image refocusing and offline enhancement, the former focuses on 3D reconstruction and modeling, which is our interest here. The plenoptic camera used in this project is a modified version of a Nikon D800 by Raytrix. A layer that is composed of around 18,000 micro-lenses is placed in front of the 36.3 mega pixels CMOS (Complementary Metal Oxide Semi-conductor) original image sensor. The micro-lenses are of three types, which differ in their focal distance. This helps them to enlarge the working depth range, by combining three zones that correspond to each lens type; namely; near, middle and far range lenses. The projection size of micro-lens (micro images) can be controlled by changing the aperture of the camera. In our case, the maximum diameter of a micro image is approximately 38–45 pixels depending on the lens type. Figure 17 illustrates two examples of captured plenoptic images showing micro images taken at two different aperture settings.
<figure><img src="image_20.png" /><figcaption>**Figure 17.** Two examples of plenoptic zoomed images showing micro images—taken at two different aperture settings (f/13 **left**, f/9 **right**) in the presence of a white filter.</figcaption></figure>

Our goal is to make this approach work in the marine environment so that only one shot is enough to obtain a properly scaled 3D model, which is a great advantage underwater, where diving time and battery capacity are limited. The procedure to work with a plenoptic camera is described in the following. Figure 18 shows examples plenoptic underwater images and reconstructed 3D models.

*   **Micro-lens Array (MLA) Calibration step:** This step computes the camera's intrinsic parameters, which include the position and the alignment of the MLA. The light field camera returns an uncalibrated raw image containing the pixels' intensity values read by the sensor during shooting. Individual micro-lenses can be distinguished in the raw image. To localize the exact position of each micro-lens in the raw image, a calibration image is taken using a white diffusive filter generating continuous illumination to highlight the edges and vignetting of each micro-lens as shown in Figure 17 (left). Using simple image processing techniques, it is possible to localize the exact position of each micro-lens, which corresponds to its optical center up to sub-pixel accuracy. In the same way, by taking another calibration image after changing the aperture of the camera so that the exterior edges of micro-lenses are touching (in Figure 17 (right)), we could use circle fitting to detect the outer edges of the micro-lenses with the help of the computed centers' positions. Hence, for any new image without the filter, it is possible to easily extract the micro images. There is no need to perform the MLA calibration underwater. The location and the size of micro lenses do not change as soon as the aperture and the focus of the camera is fixed. Finally, a metric calibration must also be performed to convert from pixels to metric units. Here, the use of a calibration grid enables determining very accurately the internal geometry of the micro-lens array. The calibration results provide many parameters: intrinsic parameters, orientation of MLA, and tilt of the MLA in regards to the image sensor as well as distortion parameters.

*   **MLA Imaging:** Each of these micro lenses capture a small fraction of the total image. These micro images can be synthesized to form the processed image in a focal plane. To do this, the depth (the distance from the object to the micro lens) must be taken for each pixel. For each point (x, y) in the object space, a set of micro lenses is selected that captures this point in their image plane. Then, the projection of this point (x, y) on the complete image plane by each micro lens is determined. The color value for the pixel in the full image plane is taken as the (weighted) average of the micro image values.

*   **Total focus processing:** As different types of micro lenses are present with a different focus for each type throughout the MLA, each group of similar micro lenses is considered as a sub-array. An image with "full focus" or an image in which a large depth of field is maintained can be synthesized using several sub-arrays of different focal lengths to compute a sharp image over a large distance. For more details on the total focusing process, we refer to Perwass and Wietzke [65].

*   **Depth map computation:** A brute force stereo points matching is performed among neighboring micro-lenses. The neighborhood is an input parameter that defines the maximum distance on the MLA between two neighboring micro-lenses. This process can be easily run in parallel. This is
achieved using RxLive software that relies totally on the GPU. The total processing time does not exceed 400 milliseconds per image. Next, using the computed calibration parameters, each matched point is triangulated in 3D in order to obtain a sparse point cloud.

* Dense point cloud and surface reconstruction: In this step, all pixel data in the image are used to reconstruct the surface in 3D using iterative filling and bilateral filtering for smoother depth intervals.

<figure><img src="image_21.png" /><figcaption>**Figure 18.** Underwater images taken using a plenoptic camera (first row); it shows also the repetition of an edge in the scene, which is essential to perform the depth estimation. The processed total focus image (second row) and reconstructed 3D models (third row and fourth row).</figcaption></figure>

## 7. Conclusions

This paper presents the results of more than 10 years collaboration between the computer science community and marine benthic biologists specialized in red coral. Through four kinds of the surveys, we developed different tools that can be applied to red coral study, our model species, to other species sharing the similar arborescents growth form, thus contributing to preserving it and protecting the marine coastal ecosystems.

These software programs will be available online (see http://perfect.medrecover.org) very soon.

Photogrammetry is the core of framework developed. Now widely used by many different scientific communities, the challenge is to extract relevant information from the reconstructed models.
First, we presented a graphical software enabling manual 3D measurements on couples of oriented
photos, showing red coral colonies by a quadrat, which is used as scale tool and surface standard.
The software provides reliable measurements (residuals generally lesser than 1 mm), for different
estimations of the characteristics of colonies.

Next, we described a set of tools especially designed to study temporal changes. Monitoring
the evolution of coral colonies over time and allowing for analyzing potential correlations between
spatial distribution and relevant population parameters such as parental relationships and recruitment
through genetic analysis in our case. The orthophotos of the studied sites are used to visualize the
red coral colonies and then to compute true geodesic distances. Other tools, such as those used for
orthophoto blending at different times, help advanced end users interpret the data.

Finally, the last survey presented was focused on archiving large red coral colonies and extracting
relevant data from them. New approaches, including NPR representations and skeletonization of the
3D models, were developed in order to automatically compute morphological information, like volume
and branch lengths, as well as obtaining simplified representations of these dense colonies. This density
can prevent or make difficult observations by knowledgeable people.

Much work remains in order to increase the accuracy of estimated measurements of large coral
colonies. Even if this study is focused on red coral, the different tools in this frame to better understand
and characterize the functioning of benthic communities could be easily extended to other species
with similar arborescent growth forms.

Beyond the geometrical and morphological survey and analysis we have done, we have started
managing the huge quantity of data collected since the beginning of this collaboration. This data is
mainly heterogeneous. We decided to use an ontology to manage all these data; indeed, LIS already
have a lot of experience in managing photogrammetric data with ontologies in an archaeological
context [66], and both the photogrammetric and the dimensional data of the coral colonies are already
stored using Ontology Web Language (OWL) formalism. We are currently working on an extension of
this approach in a 3D Information System dedicated to red coral survey based on photogrammetry
survey and knowledge representation for spatial reasoning.

**Acknowledgments:** This work is partially done in the framework of the PERFECT project, funded by the Foundation TOTAL, project 2014/257. The plenoptic camera was bought in the frame of the DGA RAPID LORI project (LOcalisation et Reconnaissance d'objets Immergés).

**Author Contributions:** Pierre Drap, Jean-Philip Royer, Mohamad Motassem Nawaf and Djamal Merad brought computer solutions to help manage Mediterranean red coral. They designed and implemented the photogrammetric and plenoptic approaches. Mauro Saccone played an important role in modeling the shape of the 3D quadrat and printing it, with the calibration grid. Olivier Bianchimani proceeded to in situ photo acquisitions. As biologists, Joaquim Garrabou, Jean-Baptiste Ledoux and Angel Lopez-Sanz determined the needs of the project and are processing the computer data.

**Conflicts of Interest:** The authors declare no conflict of interest.

References

1. Linnaeus, C. *Systema Naturae per Regna Tria Naturae: Secundum Classes, Ordines, Genera, Species, Cum Characteribus, Differentiis, Synonymis, Locis*; Impensis direct. Laurentii Salvii: London, UK, 1758; Volume 1, p. 881.

2. Ballesteros, E. *Coralligenous Assemblages: A Synthesis of Present Knowledge*. *Oceanogr. Marine Biol. Annu. Rev.* **2006**, *44*, 123–195.

3. Costantini, F.; Rossi, S.; Pintus, E.; Cerrano, C.; Gili, J.M.; Abbiati, M. Low connectivity and declining genetic variability along a depth gradient in *Corallium rubrum* populations. *Coral Reefs* **2011**, *30*, 991–1003.

4. Garrabou, J.; Harmelin, J.G. A 20-year study on life-history traits of a harvested long-lived temperate coral in the NW Mediterranean: Insights into conservation and management needs. *J. Anim. Ecol.* **2002**, *71*, 966–978.

5. Marschal, C.; Garrabou, J.; Harmelin, J.G.; Pichon, M. A new method for measuring growth and age in the precious red coral *Corallium rubrum* (L.). *Coral Reefs* **2004**, *23*, 423–432.
6. Montero-Serra, I.; Garrabou, J.; Doak, D.; Figuerola, L.; Hereu, B.; Jean-Baptiste, L.; Linares, C. Accounting for life-history strategies and timescales in marine restoration. *Conserv. Lett.* **2017**, doi:10.1111/conl.12341.

7. Drap, P.; Merad, D.; Mahiddine, A.; Seinturier, J.; Gerenton, P.; Peloso, D.; Boï, J.M.; Bianchimani, O.; Garrabou, J. In situ underwater measurements of red coral: Non-intrusive approach based on coded targets and photogrammetry. *Int. J. Heritage Digit. Era* **2014**, doi:10.1260/2047-4970.3.1.123.

8. Rogers, C.S.; Gilnack, M.; Fitz, H.C. Monitoring of Coral Reefs with Linear Transects: A Study of Storm Damage. *J. Exp. Mar. Biol. Ecol.* **1983**, *66*, 286-300.

9. Dahl, A.L. Monitoring Coral Reef for Urban Impact. *Bull. Mar. Sci.* **1981**, *33*, 544–551.

10. Mayor, A.G. Structure and Ecology of Samoan Reef. *Carnegie Inst. Wash. Publ.* **1924**, *340*, 1–25.

11. Dahl, A.L.; Lamberts, A.E. Environmental Impact of a Samoan Coral Reef: A Resurvey of Mayor’s 1917 Transect. *Pac. Sci.* **1977**, *31*, 309–319.

12. Sanchez, F.; Serrano, A.; Ballesteros, M.G. Photogrammetric quantitative study of habitat and benthic communities of deep Cantabrian Sea hard grounds. *Cont. Shelf Res.* **2009**, *29*, 1174–1188.

13. Singh, H.; Weyer, F.; Howland, J.; Duester, A.; Yoerger, D.; Bradley, A. Quantitative stereo imaging from the autonomous benthic explorer (ABE). In Proceedings of the OCEANS ‘99 MTS/IEEE—Riding the Crest into the 21st Century, Washington, DC, USA, 13–16 September 1999; Volume 1, pp. 52–57.

14. Hollister, C.D. Photogrammetry and hebble (High-Energy Benthic Boundary Layer Experiment). In *Underwater Photography: Scientific and Engineering Applications*; Van Nostrand Reinhold Co.: New York, NY, USA, 1984; p. 413.

15. Fryer, J.G. Stereoscopic Coral Maps from Underwater Photogrammetry. *Cartogr. J.* **1983**, *20*, 23–25, doi:10.1179/caj.1983.20.1.23.

16. Abdo, D.A.; Seager, J.W.; Harvey, E.; McDonald, J.I.; Kendrick, G.A.; Shortis, M. Efficiently measuring complex sessile epibenthic organisms using a novel photogrammetric technique. *J. Exp. Mar. Biol. Ecol.* **2006**, *339*, 120–133.

17. Madjidi, H.; Nagahdaripour, S. 3-D photo-mosaicking of benthic environments. In Proceedings of the OCEANS 2003, San Diego, CA, USA, 22–26 September 2003; Volume 4, pp. 2317–2318.

18. Shortis, M.; Seager, J.W.; Williams, A.; Barker, B.A.; Sherlock, M. A towed body stereo-video system for deep water benthic habitat surveys. In Proceedings of the Eighth Conference on Optical 3-D Measurement Techniques, Zurich, Switzerland, 2007; Volume II, pp. 150–157.

19. Leon, J.; Roelfsema, C.M.; Saunders, M.I.; Phinn, S.R. Measuring coral reef terrain roughness using ‘Structure-from-Motion’ close-range photogrammetry. *Geomorphology* **2015**, *242*, 21–28.

20. Zapata-Ramírez, P.; Scaradozzi, D.; Sorbi, L.; Palma, M.; Pantaleo, U.; Ponti, M.; Cerrano, C. Innovative study methods for the Mediterranean coralligenous habitats. *Adv. Oceanogr. Limnol.* **2013**, *4*, 102–119, doi:10.1080/19475721.2013.849758.

21. Shortis, M.; Harvey, E.; Abdo, D. A Review of Underwater Stereo-Image Measurement for Marine Biology and Ecology Applications. In *Oceanography and Marine Biology: An Annual Review*; Crc Press-Taylor and Francis Group: Boca Raton, FL, USA, 2009; Volume 47, pp. 257–292.

22. Canese, S.; Giusti, M.; Rossi, L.; Angiolillo, M.; Salvati, E.; Cardinali, A.; Bo, M.; Greco, S. Integration of different technologies into a Geographic Information System (GIS) to study coralligenous biocenosis in Calabrian coastal waters (South Italy). In Proceedings of the 1st Symposium on the Coralligenous and Other Calcareous Bio-Concretions of the Mediterranean Sea, Tabarka, Tunisia, 15–16 January 2009; pp. 63–67.

23. Giusti, M.; Canese, S.; Angiolillo, M.; Bo, M.; Salvati, E.; Greco, S. Three-dimentional distribution of Gerardia Savaglia in relation to depth, orientation and slope of the substratta in the south Tyrrhenian Sea. In Proceedings of the 1st Symposium on the Coralligenous and Other Calcareous Bio-Concretions of the Mediterranean Sea, Tabarka, Tunisia, 15–16 January 2009; pp. 63–67.

24. Ledoux, J.-B.; Garrabou, J.; Bianchimani, O.; Drap, P.; Féral, J.-P.; Aurelle, D. Fine-scale genetic structure and inferences on population biology in the threatened Mediterranean red coral, Corallium rubrum. *Mol. Ecol.* **2010**, *19*, 4204–4216.

25. Bythell, J.; Pan, P.; Lee, J. Three-dimensional morphometric measurements of reef corals using underwater photogrammetry techniques. *Coral Reefs* **2001**, *20*, 193–199.

26. Agisoft. Agisoft PhotoScan, 2016. Available online: http://www.agisoft.com (accessed on 22 January 2017).
27. Burns, J.; Delparte, D.; Gates, R.; Takabayashi, M. Integrating structure-from-motion photogrammetry with geospatial software as a novel technique for quantifying 3D ecological characteristics of coral reefs. *Coral Reefs* **2015**, doi:10.7717/peerj.1077.

28. Burns, J.H.R.; Delparte, D.; Gates, R.D.; Takabayashi, M. Utilizing Underwater Three-Dimensional Modeling to Enhance Ecological and Biological Studies of Coral Reefs. In Proceedings of the International Archives of the Photogrammetry, Remote Sensing and Spatial Information Sciences, Piano di Sorrento, Italy, 16–17 April 2015; pp. 61–66.

29. Burns, J.; Delparte, D.; Kapono, L.; Belt, M.; Gates, R.; Takabayashi, M. Assessing the impact of acute disturbances on the structure and composition of a coral community using innovative 3D reconstruction techniques. *Methods Oceanogr.* **2016**, *15–16*, 49–59.

30. Ferrari, R.; McKinnon, D.; He, H.; Smith, R.N.; Corke, P.; González-Rivero, M.; Mumby, P.J.; Upcroft, B. Quantifying Multiscale Habitat Structural Complexity: A Cost-Effective Framework for Underwater 3D Modelling. *Remote Sens.* **2016**, *8*, 113.

31. Figueira, W.; Ferrari, R.; Weatherby, E.; Porter, A.; Hawes, S.; Byrne, M. Accuracy and Precision of Habitat Structural Complexity Metrics Derived from Underwater Photogrammetry. *Remote Sens.* **2015**, *7*, 16883–16900.

32. Bryson, M.; Ferrari, R.; Figueira, W.; Pizarro, O.; Madin, J.; Williams, S.; Byrne, M. Characterization of measurement errors using structure-from-motion and photogrammetry to measure marine habitat structural complexity. *Ecol. Evol.* **2017**, *7*, 5669–5681.

33. Bennecke, S.; Kwasnitschka, T.; Metaxas, A.; Dullo, W.C. In situ growth rates of deep-water octocorals determined from 3D photogrammetric reconstructions. *Coral Reefs* **2016**, *35*, 1227–1239.

34. Bryson, M.; Johnson-Roberson, M.; Pizarro, O.; Williams, S. Automated registration for multi-year robotic surveys of marine benthic habitats. In Proceedings of the 2013 IEEE/RSJ International Conference on Intelligent Robots and Systems, Tokyo, Japan, 3–7 November 2013; pp. 3344–3349.

35. Johnson-Roberson, M.; Bryson, M.; Friedman, A.; Pizarro, O.; Troni, G.; Ozog, P.; Henderson, J.C. High-Resolution Underwater Robotic Vision-Based Mapping and Three-Dimensional Reconstruction for Archaeology. *J. Field Robot.* **2017**, *34*, 625–643.

36. Johnson-Roberson, M.; Pizarro, O.; Williams, S.B.; Mahon, I. Generation and visualization of large-scale three-dimensional reconstructions from underwater robotic surveys. *J. Field Robot.* **2010**, *27*, 21–51.

37. Nawaf, M.; Boi, J.M.; Merad, D.; Drap, P. Low Cost Embedded Stereo System for Underwater Surveys. In Proceedings of the 5th International Workshop on LowCost 3D, ISPRS—Technical Commission II, Hamburg, Germany, 28–29 November 2017.

38. Casella, E.; Collin, A.; Harris, D.; Ferse, S.; Bejarano, S.; Parravicini, V.; Hench, J.L.; Rovere, A. Mapping coral reefs using consumer-grade drones and structure from motion photogrammetry techniques. *Coral Reefs* **2017**, *36*, 269–275.

39. Tsai, R. A versatile camera calibration technique for high-accuracy 3D machine vision metrology using off-the-shelf TV cameras and lenses. *IEEE J. Robot. Autom.* **1987**, *3*, 323–344.

40. Drap, P. The Arpenteur Project. Available online: http://www.arpenteur.org (accessed on 20 September 2017).

41. Dijkstra, E.W. A note on two problems in connexion with graphs. *Nume. Math.* **1959**, *1*, 269–271.

42. Hart, P.E.; Nilsson, N.J.; Raphael, B. A Formal Basis for the Heuristic Determination of Minimum Cost Paths. *IEEE Trans. Syst. Sci. Cybern.* **1968**, *4*, 100–107.

43. Kimmel, R.; Sethian, J.A. Computing Geodesic Paths on Manifolds. *Proc. Natl. Acad. Sci. USA* **1998**, *95*, 8431–8435.

44. Mitchell, J.S.B.; Mount, D.M.; Papadimitriou, C.H. The Discrete Geodesic Problem. *SIAM J. Comput.* **1987**, *16*, 647–668.

45. Surazhsky, V.; Surazhsky, T.; Kirsanov, D.; Gortler, S.J.; Hoppe, H. Fast Exact and Approximate Geodesics on Meshes. *ACM Trans. Graph.* **2005**, *24*, 553–560.

46. Seinturier, J.; Riedinger, C.; Mahiddine, A.; Peloso, D.; Boï, J.M.; Merad, D.; Drap, P. Towards a 3d Based Platform for Cultural Heritage Site Survey and Virtual Exploration. In Proceedings of the ISPRS—International Archives of the Photogrammetry, Remote Sensing and Spatial Information Sciences, Strasbourg, France, 2–6 September 2013; pp. 573–578.
47. Lai, R.; Shi, Y.; Sicotte, N.; Toga, A.W. Automated corpus callosum extraction via Laplace-Beltrami nodal parcellation and intrinsic geodesic curvature flows on surfaces. In Proceedings of the IEEE International Conference on Computer Vision, Barcelona, Spain, 6–13 November 2011; pp. 2034–2040.

48. Laboratory, R.P.I.I.P.; Meagher, D. *Octree Encoding: A New Technique for the Representation, Manipulation and Display of Arbitrary 3-D Objects by Computer*; Technical Report; Rensselaer Polytechnic Institute: Troy, NY, USA, 1980.

49. Bentley, J.L. Multidimensional Binary Search Trees Used for Associative Searching. Commun. ACM **1975**, *18*, 509–517.

50. Wald, I. On fast Construction of SAH-based Bounding Volume Hierarchies. In Proceedings of the 2007 IEEE Symposium on Interactive Ray Tracing, Ulm, Germany, 10–12 September 2007; pp. 33–40.

51. Havran, V. Heuristic Ray Shooting Algorithms. Ph.D. Thesis, Department of Computer Science and Engineering, Faculty of Electrical Engineering, Czech Technical University in Prague, Prague, Czech Republic, 2000.

52. Intel. Embree, High Performance Ray Tracing Kernels. Available online: https://embree.github.io/ (accessed on 20 May 2017).

53. Wald, I.; Woop, S.; Benthin, C.; Johnson, G.S.; Ernst, M. Embree: A Kernel Framework for Efficient CPU Ray Tracing. *ACM Trans. Graph.* **2014**, *33*, 143:1–143:8.

54. DeCarlo, D.; Finkelstein, A.; Rusinkiewicz, S.; Santella, A. Suggestive Contours for Conveying Shape. *ACM Trans. Graph.* **2003**, *22*, 848–855.

55. Blum, H. A Transformation for Extracting New Descriptors of Shape. In *Models for the Perception of Speech and Visual Form*; MIT Press: Cambridge, MA, USA, 1967; pp. 362–380.

56. Serra, J. Image analysis and mathematical morphology. In *Image Analysis and Mathematical Morphology; Academic Press: Cambridge, MA, USA, 1982; Volume 1*.

57. Siddiqi, K.; Shokoufandeh, A.; Dickinson, S.J.; Zucker, S.W. Shock Graphs and Shape Matching. Int. J. Comput. Vis. **1999**, *35*, 13–32.

58. Bai, X.; Latecki, L.J.; Liu, W. Skeleton Pruning by Contour Partitioning with Discrete Curve Evolution. *IEEE Trans. Pattern Anal. Mach. Intell.* **2007**, *29*, 449–462.

59. Tagliasacchi, A.; Delame, T.; Spagnuolo, M.; Amenta, N.; Telea, A. 3D Skeletons: A State-of-the-Art Report. *Comput. Graph. Forum* **2016**, *35*, 573–597.

60. Liu, L.; Chambers, E.W.; Letscher, D.; Ju, T. A simple and robust thinning algorithm on cell complexes. *Comput. Graph. Forum* **2010**, *29*, 2253–2260.

61. Liu. 2010. Computer Software. Available online: http://www.cs.wustl.edu/~taoju/lliu/paper/pgcc/pgcc.html (accessed on 20 May 2017).

62. Ng, R.; Levoy, M.; Brédif, M.; Duval, G.; Horowitz, M.; Hanrahan, P. Light field photography with a hand-held plenoptic camera. *Comput. Sci. Tech. Rep. CSTR* **2005**, *2*, 1–11.

63. Raytrix. Raytrix 3D Light Field Camera Technology, 2016. Available online: https://www.raytrix.de/ (accessed on 22 January 2017).

64. Lytro. Lytro—Home, 2016. Available online: https://www.lytro.com/ (accessed on 22 January 2017).

65. Perwass, C.; Wietzke, L. Single lens 3D-camera with extended depth-of-field. In Proceedings of the SPIE, Human Vision and Electronic Imaging XVII, Volume 8291, Burlingame, CA, USA, 22 January 2012.

66. Drap, P.; Papini, O.; Sourisseau, J.c.; Gambin, T. Ontology-based photogrammetry survey in underwater archaeology. In Proceedings of the 14th ESWC, Portorož, Slovenia, 28 May–1 June 2017; Blomqvist, E., Hose, K., Paulheim, H., Ławrynowicz, A., Ciravegna, F., Hartig, O., Eds.; Springer: Cham, Switzerland, 2017.

<figure><img src="image_22.png" /></figure>

© 2018 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).