Jun 12, 2024

Public workspaceThree-dimensional modeling of coral skeletal fragments for morphometrics

DOI
dx.doi.org/10.17504/protocols.io.yxmvmeep6g3p/v1
  • 1Marine Laboratory, University of Guam;
  • 2Climate Change Cluster, University of Technology Sydney
Colin J Anthony
Icon indicating open access to content
QR code linking to this content
DOI: dx.doi.org/10.17504/protocols.io.yxmvmeep6g3p/v1
Protocol Citation: Colin J Anthony, Colin Lock, Renee Crisostomo, Bastian Bentlage 2024. Three-dimensional modeling of coral skeletal fragments for morphometrics. protocols.io https://dx.doi.org/10.17504/protocols.io.yxmvmeep6g3p/v1
License: This is an open access protocol distributed under the terms of the Creative Commons Attribution License,  which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited
Protocol status: Working
We use this protocol and it's working
Created: June 01, 2024
Last Modified: June 12, 2024
Protocol Integer ID: 101060
Keywords: MeshLab, phenomics, 3D scanner, surface area, coral reefs, scleractinia, coral, morphology
Funders Acknowledgements:
NSF Guam EPSCoR
Grant ID: OIA-1946352
Abstract
This protocol details a comprehensive method for scanning and modeling coral fragments to accurately calculate surface area, volume, and other key morphological features essential for phenomics and ecophysiology studies. The procedure encompasses high-resolution 3D scanning, digital reconstruction, and data generation, thus providing precise quantification of coral morphologies. This protocol offers a standardized approach for morphological analyses, supporting diverse research applications in coral biology and ecology.
Materials
Equipment
Airbrush and compressor (example)
Drying oven (examples)
Balance
Jewelry scanner (D3D-S)
Computer (Intel CORE i7 preferred)
Portable Incremental Counter

Expendables
Filtered Seawater
RUST-OLEUM Chalked, Linen White spray paint (Internet # 307244837, Model # 339822)
ULTRAGEL CONTROL Loctite super glue (#43903)
standard 26/6 staples

Software
MeshLab v2020.04
Skeletal Preparation
Skeletal Preparation
12h
12h
Airbrush tissue from coral skeleton with filtered seawater until only a bare skeleton remains.

Example Acropora pulchra fragment
after removing tissue.

Note
Here we describe a protocol designed with smaller fragments (<5g), but depending on your scanner and study requirements, this protocol is easily scalable to larger samples.

Place fragments in drying oven overnight (50°C) to evaporate remaining water from fragment.
12h
Weigh dried fragments.

This will be used later to calculate skeletal density.
Coat coral fragment in white, chalked spray paint.

We use RUST-OLEUM Chalked, Linen White spray paint (Internet # 307244837, Model # 339822). This reduces light scatter and allows the fragment to be detected by our 3D-scanner.
Depending on species, mount coral fragment on a metal pin to prevent scanner clip from covering up parts of the skeleton that are of interest.

For example, we mount staghorn Acropora branches on standard 26/6 staples with ULTRAGEL CONTROL Loctite super glue (#43903). This makes the model cleaning process much easier, as it provides a clear, easily patchable cutoff point for our coral model.
Scanning fragments and generating point clouds
Scanning fragments and generating point clouds
Depending on your research question, there are two routes you can take to model your fragment: 1) simplified morphometric modeling or 2) corallite-focused morphometric modeling.

The simplified workflow (SW) uses less detailed models quickly calculate a few informative metrics, such as skeletal density or complexity. The simplified workflow is also often paired with eco-physiological workflows, such as Symbiodiniaceae cell counting, as it inherently calculates skeletal surface area, which can be used to normalize cell counts. The corallite-focused workflow (CW) allows you to quantify characteristics traditionally used in taxonomy, but also requires more computing power and time to complete. Also note, the corallite-focused workflow can generate all of the same characters as those described in the simplified workflow.

Critical steps indicate steps that differ between the two workflows. Within these steps, differences will be indicated by wither SW or CW.

AB
Skeletal surface area & volumeOnly used to normalize other measurements
Skeletal densityRatio of fragment weight to volume
Skeletal complexityRatio of fragment area to volume
Corallite densityNumber of non-emersed corallites per cm^2
Table 1. Simple characters measured in simplified morphometric modeling

AB
Radial crowdingAverage number of emersed radial corallites /3 transects
Radial wall thicknessWidth of radial wall
Radial corallite profile lengthMaximum distance from base to outer edge of radial corallite
Radial corallite diameterMaximum diameter of radial corallite from inner to outer wall
Radial calice diameterMaximum diameter of radial calice from inner to outer wall
Axial tip diameterMaximum diameter of axial tip
Table 2. Traditional taxonomic characters (Wallace, 1999) measured in corallite-focused modeling

Mount fragment on 3D scanner. We use a D3D-s jewelry scanner, but many scanners can be used, with local protocol modification.
Name project and set directory for point cloud outputs.
Set the scanning length and diameter to ~1.5x the length and diameter of the coral fragment.

This accounts for variation, and ensures portions of the fragment are not missed during the scanning process.
Set the scanning depth to 1.3-1.5x the diameter used in the previous step.
Set horizontal turns, verticle turns, and/or horizontal angles to ones that sufficiently capture your coral.

For Acropora fragments (20x30mm), we do 12 vertical turns, and scan from 0 and 70 degrees (SW) or 4 horizontal turns and 12 vertical turns (CW). This will require local modification and fine tuning depending on coral morphology and fragment size. The goal is to sufficiently scan the entire surface area in a minimal amount of time. Generally, more turns produce a better model
Critical
SW - Set resolution to medium or high, depending on your coral's corallite complexity.

CW - Set resolution to maximum
Set camera gain to 3. (May be adjusted locally).
Start scanning.
Modeling in MeshLab
Modeling in MeshLab
Import point clouds into MeshLab

Open MeshLab > "File" > "Import Mesh" > Highlight all "*.STL" files generated by scanner > "Open"


Example point clouds directly after import (SW)

Align point clouds

Click the "Align" tool (Yellow circle with black "A" on the top tool bar) > "Glue Here Visible Meshes" > edit default ICP parameters to sample 6000 points (CW) > "Process" > Close the "Align Tool" window

I like to make sure my Avg Err is below 0.04.


Example of aligned point clouds (SW)

Critical
Reconstruct fragment surface

"Filters" > "Remeshing, Simplification and Reconstruction" > "Surface Reconstruction: Screened Poisson" > Check box next to "Merge all visible layers" > Set reconstruction depth to 9 (SW) or 11 (CW) > Check box next to "Pre-Clean" > "Apply" > Close "Surface Reconstruction: Screened Poisson" window > Hold "Ctrl" and click on the eyeball next to "Poisson mesh" on the right side of the screen > right click "Poisson mesh" > "Duplicate Current layer" (MeshLab does not have an undo function, so duplicating your mesh serves as a good backup) > Close the eyeball next to the non-selected Poisson mesh
Critical
Remove mount from surface reconstruction

"Select Faces in a rectangular Region" > highlight and delete all portions of the reconstruction that are not part of the coral by dragging the tool over unwanted regions and pressing "delete"
Close holes in reconstruction

"Filters" > "Remeshing, Simplification and Reconstruction" > "Close Holes" > set "Max size to be closed" to 5000 > "Apply"

Sometimes, it doesn't work properly, so you may need to redo or revisit step 18 with mesh copy generated in step 17 until it can effectively close the hole.


Example of final surface reconstruction (SW)

Calculate skeletal volume

"Filters" > "Quality Measure and Computations" > "Compute Geometric Measures" > Record "Mesh Volume"
Calculate skeletal surface area

If normalizing cell counts or corallites to skeletal surface area, delete all faces that did not contain live tissue (see Step 18). Recalculate surface area (see Step 20).
Calculate corallite density (CS) and radial crowding (CW)

Select "Z-painting" paint brush > Simultaneously use a Portable Incremental Counter and the paint brush tool to count non-emersed corallites for corallite density (CS) and emersed corallites for radial crowding > Record number of corallites, normalize to surface area calculated in step 21.

Emersed corallites require a higher resolution model for consistent detection, hence their inclusion in the CW workflow and not the CS workflow.
Critical
Calculate radial wall thickness, radial corallite profile length, radial corallite diameter, radial calice diameter, and axial tip diameter (CW) using the built-in "Measuring Tool" within MeshLab.

For Acropora fragments, we collect a subsample of ~30 corallites per fragment per character. Radial wall thickness, radial corallite profile length, radial corallite diameter, and radial calice diameter measurements.

The first 2 cm of each fragment were not sampled to ensure that measurements of immature corallites did not skew measurements of matured corallites. Axial tip diameter only produced one measurement per fragment.



Example Acropora aspera model (A) with radial corallite profile length (B - yellow arrow), radial corallite diameter (B - red arrow), radial calice diameter (B - pink arrow), and radial wall thickness (C - green arrow) marked with ruler tool in MeshLab.


Critical
Export final Mesh and close MeshLab

File > "Export Mesh As..." > Name file > "Save" > Close MeshLab
Protocol references
Wallace CC (1999) Staghorn Corals of the World: A Revision of the Coral Genus Acropora. Collingwood: CSIRO Publishing. 421 p.

Cignoni P, Callieri M, Corsini M, Dellepiane M, Ganovelli F, Ranzuglia G (2008) MeshLab: an Open-Source Mesh Processing Tool. Sixth Eurographics Italian Chapter Conference. 129-136.

Kazhdan M, Hoppe H (2013) Screened Poisson surface reconstruction. ACM Transactions on Graphics (TOG). 32(3), 29.