NOAA's NOS/NCCOS/CCMA Biogeography Branch, in collaboration with NOAA vessel Nancy Foster and territory, federal, and private sector partners, acquired multibeam bathymetry data in Puerto Rico from 2/25/08 to 3/8/08. Data was acquired with a hull-mounted Kongsberg Simrad EM 1002 multibeam echosounder (95 kHz) in 2008. It was processed by a NOAA contractor using CARIS HIPS software. Data has all correctors applied (attitude, sound velocity) and has been reduced to mean lower low water (MLLW) using final approved tides and zoning from NOAA COOPS. Data is in UTM zone 19 north, datum NAD83. The processed CARIS data was used to generate a CARIS BASE surface based on swath angle. An ASCII XYZ file was exported from the BASE surface and imported to ESRI ArcMap 9.2 using the Convert XYZ to Raster tool developed by CCMA for this purpose. The Benthic Terrain Modeller (BTM) tool, a collection of ArcGIS terrain visualization tools developed by the Oregon State University (OSU) Department of Geosciences and the National Oceanic and Atmospheric Administration (NOAA) Coastal Services Center, was used to calculate the rugosity of that bathymetric grid. More information on the specific algorithms used can be found in the BTM's documentation.
The project was conducted to meet IHO Order 1 and 2 accuracy standards, dependant on the project area and depth. All users should individually evaluate the suitability of this data according to their own needs and standards.
These data were prepared by an agency of the United States Government. Neither the United States Government nor any agency thereof, nor any of their employees, make any warranty, expressed or implied, or assumes any legal liability or responsibility for the accuracy, completeness, or usefulness of any information, apparatus, product, or process disclosed in this report, or represents that its use would not infringe privately owned rights. Reference therein to any specific commercial product, process, or service by trade name, trademark, manufacturer, or otherwise does not necessarily constitute or imply its endorsement, recommendation, or favoring by the United States Government or any agency thereof. Any views and opinions of authors expressed herein do not necessarily state or reflect those of the United States Government or any agency thereof. Although all data have been used by NOAA, no warranty, expressed or implied, is made by NOAA as to the accuracy of the data and/or related materials. The act of distribution shall not constitute any such warranty, and no responsibility is assumed by NOAA in the use of these data or related materials.
NOAA - NMAO personnel and R/V NANCY FOSTER
"...rugosity derivation relies, in part, on a neighborhood analysis using a 3 grid cell by 3 grid cell neighborhood. An algorithm is passed through the Raster Map Algebra Operation object within Spatial Analyst that calculates the planar distance between the center point of the center cell and of each of the eight surrounding cells in the neighborhood. Next, using the Pythagorean Theorem, the surface distance is calculated for each planar distance using the difference in elevation between the cells. The result of this function is sixteen separate grid data sets with each cell value equal to this surface distance. The next step in the process is to calculate the area formed by three adjacent sides. The result is eight triangular surface area grids. These grid datasets are combined to obtain a surface area data set for the input bathymetric data set. The final step in the process is to create a data set that represents the ratio of surface area to planar area. This final data set represents rugosity for the study area."
More information on the specific algorithms used, and contact points for questions, can be found in the BTM's documentation.
Rugosity values near 1 represent flat smooth terrain, while higher values reflect increasing rugosity or terrain roughness. For rugosity grid interpretation, is recommended that the grid be reclassified according to standard deviation divisions.