Orthomosaic imagery from the UAS survey of the coral reef off Waiakane, Molokai, Hawaii, 24 June 2018

Aerial imagery was collected using a Department of Interior-owned 3DR Solo quadcopter fitted with a Ricoh GR II digital camera featuring a global shutter. The UAS was flown on pre-programmed autonomous flight lines which were oriented roughly shore-normal and were spaced to provide approximately 75 percent overlap between images from adjacent lines, at an approximate altitude of 100 meters above ground level (AGL). The camera was triggered at 1 Hz using a built-in intervalometer. Before each flight, the camera’s digital ISO, aperture, and shutter speed were adjusted for ambient light conditions. A total of five flights were conducted for the survey between 16:40 and 17:45 UTC (06:40 and 07:45 HST). Flight F01 was a reconnaissance flight, and no mapping imagery was collected during the flight. Flights F02 and F03 were conducted at an approximate altitude of 100 meters above ground level (AGL), resulting in complete coverage of the mapping area with a nominal ground-sample-distance (GSD) of approximately 2.5 centimeters per pixel. Flights F04 and F05 were conducted using the same flight lines and altitudes of F02 and F03, but the camera was fitted with a circular polarizing filter to reduced reflections and provide improved imaging of the seafloor through the water surface. Only the images collected without the polarizing filter (F02 and F03) were used for the creation of the orthomosaic presented in this data release.

Survey control was established using twenty temporary ground control points (GCPs) distributed throughout the survey area. The GCPs were placed using a combination of kayaking, wading, and snorkeling. The GCPs consisted of: nine submerged targets consisting of small (80 centimeter X 80 centimeter) square tarps with black-and-white cross patterns anchored to the shallow (less than 1.5 meters deep) seafloor using 0.9 kilogram fishing weights; nine sub-aerial targets consisting of orange plastic five-gallon bucket lids (32 centimeter diameter) painted with a black “X” pattern and affixed in a horizontal orientation to vertical rebar stakes placed in areas of reef rubble to provide the targets with sufficient elevation to remain above the water surface during the survey; and two sub-aerial ground targets consisting of small (80 centimeter X 80 centimeter) square tarps with black-and-white cross patterns placed in the sand at the shoreline. Two of the submerged targets were disturbed by waves or currents during the survey and were not used for SfM processing. All GCP positions were measured using post-processed kinematic (PPK) GPS, using corrections from a GPS base station on a temporary benchmark (MK02) located approximately 1 kilometer away from the study area. Reference coordinates for MK02 were established using the mean position derived from four static GPS occupations with durations greater than 4 hours submitted to the National Geodetic Survey Online Positioning User Service (NGS OPUS).

The image files were renamed using a custom python script. The file names were formed using the following pattern Fx-YYYYMMDDThhmmssZ_Ryz.*, where: - Fx = Flight number - YYYYMMDDThhmmssZ = date and time in the ISO 8601 standard, where 'T' separates the date from the time, and 'Z' denotes UTC ('Zulu') time. - Ry = RA or RB to distinguish camera 'RicohA' from 'RicohB' - z = original image name assigned by camera during acquisition - * = file extension (JPG or DNG)
The approximate image acquisition coordinates were added to the image metadata (EXIF) ('geotagged') using the image timestamp and the telemetry logs from the UAS onboard single-frequency 1-Hz autonomous GPS. The geotagging process was done using the GeoSetter software package. To improve timestamp accuracy, the image acquisition times were adjusted to true ('corrected') UTC time by comparing the image timestamps with several images taken of a smartphone app ('Emerald Time') showing accurate time from Network Time Protocol (NTP) servers. For this survey, + 00:00:02 (2 seconds) were added to the image timestamp to synchronize with UTC time. The positions stored in the EXIF are in geographic coordinates referenced to the WGS84(G1150) coordinate reference system (EPSG:7660), with elevation in meters relative to the WGS84 ellipsoid.
Additional information was added to the EXIF using the command-line 'exiftool' software with the following command: exiftool ^ -P ^ -Copyright="Public Domain. Please credit U.S. Geological Survey." ^ -CopyrightNotice="Public Domain. Please credit U.S. Geological Survey." ^ -ImageDescription="Low-altitude aerial image of the shallow coral reef off Waiakane, Molokai, HI, USA, from USGS Unmanned Aircraft System (UAS) survey 2018-617-FA." ^ -Caption-Abstract="Coral reef off Waiakane, Molokai, HI, USA, from USGS survey 2018-617-FA." ^ -Caption="Aerial image of the shallow coral reef off Waiakane, Molokai, HI, USA, from USGS Unmanned Aircraft System (UAS) survey 2018-617-FA." ^ -sep ", " ^ -keywords="coral reef, Molokai, Maui County, Hawaii, Waiakane, 2018-617-FA, Unmanned Aircraft System, UAS, drone, aerial imagery, U.S. Geological Survey, USGS, Pacific Coastal and Marine Science Center" ^ -comment="Low-altitude aerial image from USGS Unmanned Aircraft System (UAS) survey 2018-617-FA." ^ -Credit="U.S. Geological Survey" ^ -Contact="pcmsc_data@usgs.gov" ^ -Artist="U.S. Geological Survey, Pacific Coastal and Marine Science Center"

Structure-from-motion (SfM) processing techniques were used to create the point clouds in the Agisoft Photoscan/Metashape software package using the following workflow with the JPG images from F02 and F03 (images collected without a polarizing filter): 1. Preliminary image alignment and sparse point cloud error reduction was performed using shoreline imagery to develop an a priori camera lens model. The lens model was then fixed for the subsequent SfM processing steps. 2. Additional image alignment for all the images was performed with the following parameters - Accuracy: 'high'; Pair selection: 'reference', 'generic'; Key point limit: 0 (unlimited); Tie point limit: 0 (unlimited). 3. Sparse point cloud error reduction was performed using an iterative gradual selection and camera optimization process in which all sparse points exceeding a Reconstruction Uncertainty of 10 were removed from the sparse point cloud. Additional sparse points obviously above or below the true surface were manually deleted after the last error reduction iteration, and a final camera optimization was performed. 4. Ground control points (GCPs) were manually marked for all GCPs. All submerged GCPs were disabled in subsequent processing steps to serve as validation check points. 5. A final camera optimization was performed, and a dense point cloud was created using the 'high' accuracy setting, with 'aggressive' depth filtering. 6. A preliminary Digital Surface Model (DSM) with a native resolution of 5.3 centimeters per pixel was created using the dense point cloud. This DSM was used only for orthorectification of imagery. 7. An RGB orthomosaic with a native resolution of 2.5 centimeters per pixel was created using the DSM as the orthorectification surface. 8. The orthomosaic was clipped using a manually digitized clipping mask to exclude areas on the perimeter of the survey which suffered from obvious off-nadir edge effects. 9. The orthomosaic was exported to a GeoTIFF format with a 2.5-centimeter pixel resolution and converted to a cloud-optimized GeoTIFF formats (using internal JPEG compression with a quality of 90) for compatibility with cloud storage services using the GDAL software package.