High Resolution Digital Elevation Model (DEM) of Looe Key, Florida, 2022

Frequently anticipated questions:

• What does this data set describe?
  1. How might this data set be cited?
  2. What geographic area does the data set cover?
  3. What does it look like?
  4. Does the data set describe conditions during a particular time period?
  5. What is the general form of this data set?
  6. How does the data set represent geographic features?
  7. How does the data set describe geographic features?
• Who produced the data set?
  1. Who are the originators of the data set?
  2. Who also contributed to the data set?
  3. To whom should users address questions about the data?
• Why was the data set created?
• How was the data set created?
  1. From what previous works were the data drawn?
  2. How were the data generated, processed, and modified?
  3. What similar or related data should the user be aware of?
• How reliable are the data; what problems remain in the data set?
  1. How well have the observations been checked?
  2. How accurate are the geographic locations?
  3. How accurate are the heights or depths?
  4. Where are the gaps in the data? What is missing?
  5. How consistent are the relationships among the data, including topology?
• How can someone get a copy of the data set?
  1. Are there legal restrictions on access or use of the data?
  2. Who distributes the data?
  3. What's the catalog number I need to order this data set?
  4. What legal disclaimers am I supposed to read?
  5. How can I download or order the data?
• Who wrote the metadata?

What does this data set describe?

1. How might this data set be cited?
   Kranenburg, Christine J., Hatcher, Gerald A., and Warrick, Jonathan A., 20240628, High Resolution Digital Elevation Model (DEM) of Looe Key, Florida, 2022: data release DOI:10.5066/P1QRS3SK, U.S. Geological Survey - St. Petersburg Coastal and Marine Science Center, St. Petersburg, FL.

   Online Links:

   • https://doi.org/10.5066/P1QRS3SK

   This is part of the following larger work.
   Kranenburg, Christine J., Hatcher, Gerald A., Zawada, David G., Warrick, Jonathan A., Yates, Kimberly K., and Johnson, Selena A., 20240628, Underwater Photogrammetry Products of Looe Key, Florida From Images Acquired Using the SQUID-5 System in July 2022: data release DOI:10.5066/P1QRS3SK, U.S. Geological Survey - St. Petersburg Coastal and Marine Science Center, St. Petersburg, FL.

   Online Links:

   • https://doi.org/10.5066/P1QRS3SK

2. What geographic area does the data set cover?

   West_Bounding_Coordinate: -81.41016166
   East_Bounding_Coordinate: -81.40065243
   North_Bounding_Coordinate: 24.54779751
   South_Bounding_Coordinate: 24.54442068
   

3. What does it look like?

   LKR22_dem_orthoavg_hirez_webview.png (PNG)

   Sample view of the high resolution (1-cm) DEM and corresponding orthoimage (tile B-1) in Portable Network Graphics (PNG) file format. The image is available on the webpage for this data release.

4. Does the data set describe conditions during a particular time period?

   Beginning_Date: 12-Jul-2022

   Ending_Date: 16-Jul-2022

   Currentness_Reference:

   ground condition

5. What is the general form of this data set?

   Geospatial_Data_Presentation_Form: raster digital data
   

6. How does the data set represent geographic features?
   1. How are geographic features stored in the data set?
      This is a Raster data set. It contains the following raster data types:
      • Dimensions, type Pixel
   2. What coordinate system is used to represent geographic features?
      

      Grid_Coordinate_System_Name: Universal Transverse Mercator
      Universal_Transverse_Mercator:

      UTM_Zone_Number: 17N
      Transverse_Mercator:

      Scale_Factor_at_Central_Meridian: 0.9996
      Longitude_of_Central_Meridian: -81
      Latitude_of_Projection_Origin: 0.0
      False_Easting: 500000.0
      False_Northing: 0.0
      

      Planar coordinates are encoded using coordinate pair
      Abscissae (x-coordinates) are specified to the nearest 0.01
      Ordinates (y-coordinates) are specified to the nearest 0.01
      Planar coordinates are specified in Meters
      The horizontal datum used is North American Datum of 1983 (2011).
      The ellipsoid used is GRS 1980.
      The semi-major axis of the ellipsoid used is 6378137.000000.
      The flattening of the ellipsoid used is 1/298.257222101.
      

      Vertical_Coordinate_System_Definition:

      Altitude_System_Definition:

      Altitude_Datum_Name: North American Vertical Datum 1988 (NAVD88)
      Altitude_Resolution: 0.001
      Altitude_Distance_Units: meters
      Altitude_Encoding_Method: Attribute values
      

7. How does the data set describe geographic features?

   Entity_and_Attribute_Overview:

   The DEMs are presented as tiled 32-bit floating point GeoTIFFs where pixels represent elevation in meters relative to NAVD88 (GEOID 12B) and the no-data value is -32,767. The horizontal projection is NAD83(2011) UTM Zone 17N. These DEM GeoTIFFs are available in 1-cm resolution for eight blocks (or tiles) of the survey area. Each tile measures 200 meters on a side and are labeled SQUID5_LKR_2022_DEM-col-row.tif where col represents the column name and can have a value of A-E, and row is the row number and can have a value of 0-1.

   Entity_and_Attribute_Detail_Citation: U.S. Geological Survey
   

Who produced the data set?

1. Who are the originators of the data set? (may include formal authors, digital compilers, and editors)
   • Christine J. Kranenburg
   • Gerald A. Hatcher
   • Jonathan A. Warrick
2. Who also contributed to the data set?
   Data collection was funded by the U.S. Geological Survey Pacific Coastal Marine Science Center and the U.S. Geological Survey Saint Petersburg Coastal and Marine Science Center. The authors would like to thank Dr. Jason Spadaro, Assistant Professor, Marine Science and Technology, College of the Florida Keys for installing calibration targets on the reef, Lisa Symons, Regional Response Coordinator, and the staff of the Eastern Region, Office of National Marine Sanctuaries, NOAA Florida Keys National Marine Sanctuary, for coordination efforts.
3. To whom should users address questions about the data?
   U.S. Geological Survey, St. Petersburg Coastal and Marine Science Center

   Attn: SPCMSC Data Management Group

   600 4th St South

   St. Petersburg, FL
   

   727-502-8000 (voice)

   gs-g-spcmsc_data_inquiries@usgs.gov

Why was the data set created?

How was the data set created?

1. From what previous works were the data drawn?

   raw images (source 1 of 2)

   Kranenburg, Christine J., Hatcher, Gerald A., Zawada, David G., Warrick, Jonathan A., and Yates, Kimberly K., 20231222, Overlapping seabed images collected at Looe Key coral reef, Florida, 2022: U.S. Geological Survey - St. Petersburg Coastal and Marine Science Center, St. Petersburg, Florida.

   Online Links:

   • https://doi.org/10.5066/P9M3NYWI
   • https://cmgds.marine.usgs.gov/data-services/rscc/SQUID/FL_2022_LooeKey/images/

   Type_of_Source_Media: TIFF
   Source_Contribution: Raw images to which SfM techniques were applied.
   

   GNSS antenna positions (source 2 of 2)

   Kranenburg, Christine J., Hatcher, Gerald A., Zawada, David G., Warrick, Jonathan A., and Yates, Kimberly K., 20231222, GNSS locations of seabed images collected at Looe Key coral reef, Florida, 2022: U.S. Geological Survey - St. Petersburg Coastal and Marine Science Center, St. Petersburg, Florida.

   Online Links:

   • https://doi.org/10.5066/P9M3NYWI
   • https://cmgds.marine.usgs.gov/data-services/rscc/SQUID/FL_2022_LooeKey/image_positions/FL_2022_LooeKey_Image_Locations.txt

   Type_of_Source_Media: comma-delimited text file
   Source_Contribution:

   Location data for the raw images to which SfM techniques were applied.

2. How were the data generated, processed, and modified?

   Date: 27-Dec-2022 (process 1 of 3)

   IMAGERY COLOR CORRECTION Because of the strong color modifications caused by light absorption and scattering in underwater imaging, a color correction process was conducted on the raw images. The color correction was a twofold process. First, images were corrected for the high absorption (and low color values) in the red band using the color balancing techniques of Ancuti and others (2017). For this, the red channel was modified using the color compensation equations of Ancuti and others (2017, see equation 4 on page 383) that use both image-wide and pixel-by-pixel comparisons of red brightness with respect to green brightness. After compensation, the images were white balanced using the "greyworld" assumption that is summarized in Ancuti and others (2017). Combined, these techniques ensured that each color band histogram was centered on similar values and had similar spread of values. The remaining techniques of Ancuti and others (2017), which include sharpening techniques and a multi-product fusion, were not employed. The resulting images utilized only about a quarter to a half of the complete 0-255 dynamic range of the three-color bands. Thus, the brightness values of each band were stretched linearly over the complete range while allowing the brightest and darkest 0.05 percent of the original image pixels (that is, 2506 of the 5.013 million pixels) to be excluded from the histogram stretch. This final element was included to ensure that light or dark spots in the images, which often occurred from water column particles or image noise, did not exert undo control on the final brightness values. Color-corrected images were output with the same file names and file types as the originals to make replacement within the SfM photogrammetry project easy. As a courtesy, the script used to implement this procedure is provided as a supplemental support file (OrthoImage_Color_Correction_Procedure.m), included in this related data release https://cmgds.marine.usgs.gov/data-releases/datarelease/10.5066-P93RIIG9/. Person who carried out this activity:
   

   Jonathan A. Warrick

   U.S. Geological Survey, Pacific Coastal and Marine Science Center

   Research Geologist

   2885 Mission St.

   Santa Cruz, CA
   

   USA

   831-460-7569 (voice)

   jwarrick@usgs.gov

   Data sources used in this process:

   • raw images

   Data sources produced in this process:

   • color-corrected images

   Date: 28-Feb-2023 (process 2 of 3)

   SfM PHOTOGRAMMETRY Digital imagery and position data recorded by the SQUID-5 system were processed using SfM photogrammetry techniques that generally follow the workflow outlined by Hatcher and others (2020 and 2023). These techniques are detailed here and include specific references to parameter settings and processing workflow. The primary software used for SfM processing was Agisoft Metashape Professional, version 1.6.6, build 11715, which will be referred to as "Metashape" in the discussion herein. Because of the large number of images in this dataset, processing was conducted on a 792-CPU-core linux-based High-Performance Computing (HPC) cluster at the USGS Advanced Research Computing (ARC) group (https://doi.org/10.5066/P9XE7ROJ). First, the raw images collected during the five mission days were added to a new project in Metashape. Raw images were used over the color-corrected images, owing to their larger dynamic range, which generally resulted in more SfM tie points. The images from each camera were assigned a unique camera calibration group in the Camera Calibration settings. Within the Camera Calibration settings, the camera parameters were also entered as 0.00345 x 0.00345 millimeter (mm) pixel sizes for all camera sensors, 8 mm focal length for the central camera (CAM13), and 6 mm focal lengths for the remaining cameras (CAM01, CAM39, CAM75, CAM82). These different focal lengths represented different lenses chosen for each camera. Additionally, the cameras required offsets to transform the GNSS positions to each camera's entrance pupil (that is, optical center). Initial measurements of these offsets were obtained using a separate SfM technique, outlined in Hatcher and others (2020), which found the offsets to be: Camera X(m) Y(m) Z(m) CAM01 -0.320 -0.205 0.823 CAM13 0.033 0.036 0.739 CAM39 0.170 -0.280 0.838 CAM75 0.047 0.396 0.698 CAM82 -0.110 -0.690 0.675 Where X and Y are the camera sensor parallel offsets, and Z is the sensor normal offset. The accuracy settings were chosen to be 0.01 m for CAM13 and 0.025 m for the other 4 cameras. Lastly, these offsets were allowed to be adjusted using the "Adjust GPS/INS offset" option, because slight camera shifts may occur with each rebuild and use of the SQUID-5 system. The SQUID-5 GNSS antenna positions were then imported into the project and matched with each image by time. The coordinates were converted in Metashape to the North American Datum of 1983 (NAD83 [2011]) Universal Transverse Mercator (UTM) Zone 17 North (17N) projected coordinate system, and altitudes were converted to the North American Vertical Datum of 1988 (NAVD88) orthometric heights (in meters). Prior to aligning the images, the Metashape reference settings were assigned. The coordinate system was "NAD83(2011) / UTM zone 17N". The camera accuracy was set to 0.10 m in the horizontal and 0.15 m in vertical dimensions, following an examination of the source GNSS data. Tie point accuracy was set at 1.0 pixels. The remaining reference settings were not relevant, because there were no camera orientation measurements, marker points, or scale bars in the SfM project. The data were then aligned in Metashape using the "Align Photos" workflow tool. Settings for the alignment included "High" accuracy, Generic preselection turned OFF and "Reference" preselection turned "ON" and using the "Source" information. This last setting allowed the camera position information to assist with the alignment process. Additionally, the key point limit was set to 50,000 and the tie point limit was assigned a value of zero, which allows for the generation of the maximum number of points for each image. Lastly, neither the "Guided image matching" nor the "Adaptive camera model fitting" options were used. This process resulted in over 205 million tie points. The total positional errors for the cameras were reported to be 0.042 m, 0.027 m, and 0.044 m in the east, north and altitude directions, respectively. Thus, the total positional error was 0.067 m. To improve upon the camera calibration parameters and computed camera positions, an optimization process was conducted that was consistent with the techniques of Hatcher and others (2020), which are based on the general principles provided in Over and others (2021). First, a duplicate of the aligned data was created in case the optimization process eliminated too much data using the "Duplicate Chunk" tool. Within the new chunk, the least valid tie points were removed using the "Gradual Selection" tools. As noted in Hatcher and others (2020), these tools are used less aggressively for the underwater imagery of SQUID-5 than commonly used for aerial imagery owing to the differences in image quality. First, all points with a "Reconstruction Uncertainty" greater than 20 were selected and deleted. Then, all points with a "Projection Accuracy" greater than 8 were selected and deleted. The camera parameters were then recalibrated with the "Optimize Cameras" tool. Throughout this process the only camera parameters that were adjusted were f, k1, k2, k3, cx, cy, p1, and p2. Once the camera parameters were adjusted, all points with "Reprojection Errors" greater than 0.4 were deleted, and the "Optimize Cameras" tool was used one final time. This optimization process resulted in slightly over 146.5 million tie points, a reduction of roughly 28 percent of the original tie points. The camera positional errors were reported to be 0.040 m, 0.025 m, and 0.044 m in the east, north and altitude directions, respectively, and the total positional error was 0.064 m. The final computed arm offsets were found to be: Camera X(m) Y(m) Z(m) CAM01 -0.317 -0.200 0.847 CAM13 0.019 -0.113 0.768 CAM39 0.168 -0.259 0.861 CAM75 0.049 0.417 0.717 CAM82 -0.108 -0.675 0.699 Following the alignment and optimization of the SQUID-5 data, mapped SfM products were generated in Metashape. For these steps, the original raw images were replaced with color-corrected images. This replacement was conducted by resetting each image path from the raw image to the color-corrected image. First, a three-dimensional dense point cloud was generated using the "Build Dense Cloud" workflow tool. This was run with the "High" quality setting and the "Moderate" depth filtering, and the tool was set to calculate both point colors and confidence. The resulting dense cloud was over 10.5 billion points over the 0.13 square kilometer survey area, or roughly 81,000 points per square meter (8.1 points per square centimeter). The dense points were classified by thresholding Metashape-computed confidence values, which are equivalent to the number of image depth maps that were integrated to make each point. Values of one were assigned "low noise", and values of two and greater were assigned "unclassified". The final Dense cloud was partitioned into blocks (also referred to as tiles) measuring 200 meters on a side and exported with point colors and classification as a LAZ file type. Person who carried out this activity:
   

   Christine J. Kranenburg

   U.S. Geological Survey, St. Petersburg Coastal and Marine Science Center

   Cartographer

   600 4th Street South

   St. Petersburg, FL
   

   USA

   727-502-8000 (voice)

   ckranenburg@usgs.gov

   Data sources used in this process:

   • raw images
   • color-corrected images
   • GNSS antenna positions

   Data sources produced in this process:

   • point cloud

   Date: 28-Feb-2023 (process 3 of 3)

   GENERATION OF DIGITAL ELEVATION MODEL (DEM) A digital elevation model (DEM), which is a x,y raster of elevation values, was generated from the point cloud using the Metashape "Build DEM" workflow tool using a geographic projection, dense cloud source data, interpolation disabled, and the default output resolution of 0.00511 meters. The resulting DEMs are available in a 1-cm resolution pixel size, as a 5-by-2 tile index grid (8 tiles total, due to 2 tiles containing no data). A "quicklook" version of this DEM was also created by exporting the data at a 1-m pixel resolution. Person who carried out this activity:
   

   Christine J. Kranenburg

   U.S. Geological Survey, St. Petersburg Coastal and Marine Science Center

   Cartographer

   600 4th Street South

   St. Petersburg, FL
   

   USA

   727-502-8000 (voice)

   ckranenburg@usgs.gov

   Data sources used in this process:

   • point cloud

   Data sources produced in this process:

   • DEM

3. What similar or related data should the user be aware of?
   Hatcher, Gerald A., Warrick, Jonathan A., Kranenburg, Christine J., and Ritchie, Andrew C., 20230726, Accurate Maps of Reef-scale Bathymetry with Synchronized Underwater Cameras and GNSS: Remote Sensing 15(15), 3727.

   Online Links:

   • https://doi.org/10.3390/rs15153727

   Hatcher, Gerald A., Warrick, Jonathan A., Ritchie, Andrew C., Dailey, Evan T., Zawada, David G., Kranenburg, Christine, and Yates, Kimberly K., 20200626, Accurate Bathymetric Maps From Underwater Digital Imagery Without Ground Control: Frontiers in Marine Science Volume 7, Article 525.

   Online Links:

   • https://doi.org/10.3389/fmars.2020.00525

   Ancuti, Codruta O., Ancuti, Cosmin, Vleeschouwer, Christophe De, and Bekaert, Philippe, 2017, Color Balance and Fusion for Underwater Image Enhancement: IEEE Transactions on Image Processing Volume 27, Number 1.

   Online Links:

   • https://doi.org/10.1109/TIP.2017.2759252

   Over, Jin-Si R., Ritchie, Andrew C., Kranenburg, Christine J., Brown, Jenna A., Buscombe, Daniel, Noble, Tom, Sherwood, Christopher R., Warrick, Jonathan A., and Wernette, Philippe A., 2021, Processing Coastal Imagery with Agisoft Metashape Professional Edition, Version 1.6--Structure from Motion Workflow Documentation: U.S. Geological Survey Open-File Report 2021-1039.

   Online Links:

   • https://doi.org/10.3133/ofr20211039

How reliable are the data; what problems remain in the data set?

1. How well have the observations been checked?
   The accuracy of the position data used for SfM data processing is based on the accuracy and precision of the Global Navigation Satellite System (GNSS) equipment and camera timing. The post-processed GNSS navigation data produced a 10-hertz (Hz) vehicle trajectory with an estimated 2-sigma accuracy of 10 cm horizontal and 15 cm vertical. The horizontal and vertical accuracies of the surface models generated by SfM were assessed with positional error assessments of the cameras and found to be less than 5 cm in both dimensions.
2. How accurate are the geographic locations?
   Previous SfM-based measurements of the field-based Sediment Elevation Table (SET) stations at USGS field sites in the Florida Keys were within 3 cm of the total uncertainty of the field-based GNSS measurements. Additionally, the average horizontal scaling of the models was found to be between 0.016 percent and 0.024 percent of water depth (Hatcher and others, 2020). No independent assessment of horizontal accuracy was possible from the Looe Key field site.
3. How accurate are the heights or depths?
   Previous SfM-based measurements of the field-based Sediment Elevation Table (SET) stations at USGS field sites in the Florida Keys were within 3 cm of the total uncertainty of the field-based GNSS measurements. The average vertical scaling of the models is between 0.016 percent and 0.024 percent of water depth (Hatcher and others, 2020). No independent assessment of vertical accuracy was possible from the Looe Key field site.
4. Where are the gaps in the data? What is missing?
   Gaps in the data coverage are coded with a NoData value of -32,767 in the DEM. Gaps primarily occur either because the line spacing briefly widened such that there was insufficient sidelap for reconstruction, or over sand patches which have moving sandwaves and / or lack texture, which is necessary for image correlation. The largest of these data voids is an approximately 450 square meter void in the south-central section of the dataset. Additionally, small gaps may be present at vertical or concave surfaces due to the lack of visibility of these surfaces when viewed from above. These high-relief areas are common around the perimeter of reef spurs. No data were collected in 2 of the 10 potential tiles covered by the 5-by-2 index map. Data in those tiles were either too deep to be captured with an optical system or too shallow to access with the vessel used during data collection.
5. How consistent are the relationships among the observations, including topology?
   Fifteen submerged marker buoys are present and visible in the data products. These buoys are anchored to the seafloor with ropes and resemble sea mounts in the DEM. They are generally 3-5 meters tall, less than 1 meter in diameter and are located at or near the following coordinates: Easting(m), Northing(m) 458539.4, 2714590.8 458564.1, 2714612.2 458565.5, 2714614.8 458765.1, 2714620.4 458787.2, 2714697.4 458849.8, 2714717.5 458871.1, 2714674.9 458905.5, 2714717.1 458958.7, 2714733.4 458998.4, 2714711.2 459033.9, 2714742.6 459057.8, 2714704.1 459130.6, 2714791.0 459184.2, 2714776.0 459185.5, 2714832.4

How can someone get a copy of the data set?

1. Who distributes the data set? (Distributor 1 of 1)
   
   U.S. Geological Survey, St. Petersburg Coastal and Marine Science Center

   600 4th St. South

   St. Petersburg, FL
   

   727-502-8000 (voice)

   gs-g-spcmsc_data_inquiries@usgs.gov
2. What's the catalog number I need to order this data set? SQUID5_LKR_2022_DEM-col-row.tif*
3. What legal disclaimers am I supposed to read?
   Unless otherwise stated, all data, metadata and related materials are considered to satisfy the quality standards relative to the purpose for which the data were collected. Although these data and associated metadata have been reviewed for accuracy and completeness and approved for release by the U.S. Geological Survey (USGS), no warranty expressed or implied is made regarding the display or utility of the data on any other system or for general or scientific purposes, nor shall the act of distribution constitute any such warranty. Any use of trade, product, or firm names is for descriptive purposes only and does not imply endorsement by the U.S. Government.
4. How can I download or order the data?
   • Availability in digital form:

     Data format:	TIFF images can be opened directly with any TIFF-compatible image viewer. in format TIFF Size: 5729
     Network links:	https://doi.org/10.5066/P1QRS3SK

   • Cost to order the data: None

Who wrote the metadata?

Dates:

Last modified: 14-Aug-2024

Metadata author:

U.S. Geological Survey, St. Petersburg Coastal and Marine Science Center

Attn: SPCMSC Data Management Group

600 4th St. South

St. Petersburg, FL

727-502-8000 (voice)

gs-g-spcmsc_data_inquiries@usgs.gov

Metadata standard:

Content Standard for Digital Geospatial Metadata (FGDC-STD-001-1998)