Slide Menu
Search Help Show/Hide Menu
Short Citation:
National Geodetic Survey, 2024: 2020 NOAA NGS Topobathy Lidar DEM: Alamagan, CNMI, https://www.fisheries.noaa.gov/inport/item/64808.

Item Identification

Title: 2020 NOAA NGS Topobathy Lidar DEM: Alamagan, CNMI
Status: Completed
Publication Date: 2021
Abstract:

Woolpert, Inc. was contracted to acquire and process topographic-bathymetric lidar for the islands of Anatahan, Alamagan, Guguan, and Sarigan in response to Hurricane Yutu for Quantum Spatial, Inc. (QSI). Woolpert collected lidar using their Leica HawkEye 4X (HE4X) topo-bathy lidar sensor that consists of a Chiroptera 4X (CH4X) sensor, with an additional Leica 40kHz deep bathymetric channel to provide high density topo lidar. The HE4X is a latest generation topographic and bathymetric lidar sensor. The system provides denser data than previous traditional bathymetric lidar systems. It is unique in its ability to acquire bathymetric lidar, topographic lidar and 4-band digital camera imagery simultaneously. The HE4X provided 300 kHz topographic data, an effective 140 kHz shallow bathymetric data and 40 kHz deep bathymetric data. 4-band 80 MP digital camera imagery was also collected simultaneously with the sensor’s RCD-30 camera.

The bathymetric and topographic lasers are independent and do not share an optical chain or receivers, so they are optimized for their specific function. As with any bathymetric lidar, maximum depth penetration is a function of water clarity and seabed reflectivity. The HE4X is designed to penetrate to 3 times the secchi depth. This is also represented as Dmax = 4/K, where K is the diffuse attenuation coefficient, and assuming K is between 0.1 and 0.3, a normal sea state and 15% seabed reflectance.

Both the topographic and bathymetric sub-systems use a palmer scanner to produce an elliptical scan pattern of laser points with a degree of incidence ranging from +/-14 degrees (front and back) to +/-20 degrees (sides), providing a 40 degrees field of view. This has the benefit of providing multiple look angles on a single pass and helps to eliminate shadowing effects. This can be of particular use in urban areas, where all sides of a building are illuminated, or for bathymetric features such as the sides of narrow water channels or features on the seafloor such as smaller objects and wrecks. It also assists with penetration in the surf zone where the back scan passes the same ground location a couple of seconds after the front scan, allowing the areas of whitewater to shift.

All topo lidar data for this project were collected simultaneous to meet United States Geological Survey, Quality Level 1 (USGS QL1) with a minimum of 8 pts per square meter at an accuracy of 10cm RMSEz. A minimum of 2 points per square meter were acquired for bathymetric lidar data. For practical purposes the survey area was divided into survey blocks in each island, allowing acquisition to be conducted in the most efficient and consistent manner possible.

The data includes topobathy data in an LAS 1.4 format file along with associated bare earth digital elevation models (DEM). This file is the project specified 1 meter bare earth DEM dataset. The dataset was derived from topobathymetric data in a LAS format 1.4, point data record format 6, with the following classifications in accordance with project specifications and the American Society for Photogrammetry and Remote Sensing (ASPRS) classification standards:

1 - unclassified

2 - ground

7 - noise

40 - bathymetric bottom or submerged topography

41 - water surface

42 - derived water surface

43 - manmade submerged feature

45 - water column

1 Overlap - edge clip

1 Withheld - bathy land

User data values differentiates between NIR and green lasers. A value of of 1 indicates the point is from the NIR laser, and a values of 2-5 indicate the green laser.

This dataset is the 1m orthometric NAVD88 (using Geoid12b) DEM.

Purpose:

This lidar data was required by National Oceanic and Atmospheric Administration (NOAA) and the National Geodetic Survey (NGS), Remote Sensing Division Coastal Mapping Program (CMP) to enable accurate and consistent measurement of the national shoreline. The CMP works to provide a regularly updated and consistent national shoreline to define America's marine territorial limits and manage coastal resources.

Supplemental Information:

This project consisted of two separate mobilizations to acquire all lidar data. The HE4X system was installed in a Reims-Cessna F406 (ZK-XLF) owned and operated by Kiwi Air (New Zealand) in June 2019. Calibration flights were collected over the Niue International Airport (IUE) on January 9, 2020. After transiting, 19 survey flights were performed from January 20 to February 15, 2020, with eight weather days during this period. The second mobilization took place from June 19 to July 13, 2020, which included six survey flights, ten weather days, and three airplane transit days. A second calibration flight was collected over Garapan, Saipan on July 14, 2020, and results from this calibration were used in processing the lidar data from the second mobilization. Alamagan was acquired 20200621- 20200702.

NOAA NGS Project: GU1903-TB-C

CONTRACTOR: Woolpert, Inc

Keywords

Theme Keywords

Thesaurus Keyword
Global Change Master Directory (GCMD) Science Keywords
EARTH SCIENCE > LAND SURFACE > TOPOGRAPHY > TERRAIN ELEVATION
Global Change Master Directory (GCMD) Science Keywords
EARTH SCIENCE > OCEANS > BATHYMETRY/SEAFLOOR TOPOGRAPHY > BATHYMETRY > COASTAL BATHYMETRY
Global Change Master Directory (GCMD) Science Keywords
EARTH SCIENCE > OCEANS > COASTAL PROCESSES > COASTAL ELEVATION
ISO 19115 Topic Category
elevation
UNCONTROLLED
None beach
None DEM
None erosion

Spatial Keywords

Thesaurus Keyword
Global Change Master Directory (GCMD) Location Keywords
VERTICAL LOCATION > LAND SURFACE
Global Change Master Directory (GCMD) Location Keywords
VERTICAL LOCATION > SEA FLOOR

Instrument Keywords

Thesaurus Keyword
Global Change Master Directory (GCMD) Instrument Keywords
LIDAR > Light Detection and Ranging

Platform Keywords

Thesaurus Keyword
Global Change Master Directory (GCMD) Platform Keywords
Airplane > Airplane

Physical Location

Organization: Office for Coastal Management
City: Charleston
State/Province: SC

Data Set Information

Data Set Scope Code: Data Set
Data Set Type: Elevation
Maintenance Frequency: None Planned
Data Presentation Form: Model (digital)
Distribution Liability:

Any conclusions drawn from the analysis of this information are not the responsibility of NOAA, the Office for Coastal Management or its partners

Data Set Credit: National Oceanic and Atmospheric Administration (NOAA), National Geodetic Survey (NGS), Remote Sensing Division (RSD), Coastal Mapping Program (CMP)

Support Roles

Data Steward

CC ID: 1070601
Date Effective From: 2021
Date Effective To:
Contact (Organization): NOAA Office for Coastal Management (NOAA/OCM)
Address: 2234 South Hobson Ave
Charleston, SC 29405-2413
Email Address: coastal.info@noaa.gov
Phone: (843) 740-1202
URL: https://coast.noaa.gov

Distributor

CC ID: 1070599
Date Effective From: 2021
Date Effective To:
Contact (Organization): NOAA Office for Coastal Management (NOAA/OCM)
Address: 2234 South Hobson Ave
Charleston, SC 29405-2413
Email Address: coastal.info@noaa.gov
Phone: (843) 740-1202
URL: https://coast.noaa.gov

Metadata Contact

CC ID: 1070600
Date Effective From: 2021
Date Effective To:
Contact (Organization): NOAA Office for Coastal Management (NOAA/OCM)
Address: 2234 South Hobson Ave
Charleston, SC 29405-2413
Email Address: coastal.info@noaa.gov
Phone: (843) 740-1202
URL: https://coast.noaa.gov

Originator

CC ID: 1070602
Date Effective From: 2021
Date Effective To:
Contact (Organization): National Geodetic Survey (NGS)
Address: 1315 East-West Hwy
Silver Spring, MD 20910
URL: https://geodesy.noaa.gov/

Extents

Currentness Reference: Ground Condition

Extent Group 1

Extent Group 1 / Geographic Area 1

CC ID: 1295401
W° Bound: 145.80744
E° Bound: 145.855674
N° Bound: 17.623287
S° Bound: 17.576499

Extent Group 1 / Time Frame 1

CC ID: 1295400
Time Frame Type: Range
Start: 2020-06-21
End: 2020-07-02

Spatial Information

Spatial Representation

Representations Used

Grid: Yes
Vector: No
Text / Table: No
TIN: No
Stereo Model: No
Video: No

Reference Systems

Reference System 1

CC ID: 1295398

Coordinate Reference System

CRS Type: Projected
EPSG Code: EPSG:8693
EPSG Name: NAD83(MA11) / UTM zone 55N
See Full Coordinate Reference System Information

Access Information

Security Class: Unclassified
Data Access Procedure:

Data is available online for bulk or custom downloads

Data Access Constraints:

None

Data Use Constraints:

Users should be aware that temporal changes may have occurred since this data set was collected and some parts of this data may no longer represent actual surface conditions. Users should not use this data for critical applications without a full awareness of its limitations.

Distribution Information

Distribution 1

CC ID: 1070603
Start Date: 2021-05
End Date: Present
Download URL: https://coast.noaa.gov/dataviewer/#/lidar/search/where:ID=9338
Distributor: NOAA Office for Coastal Management (NOAA/OCM) (2021 - Present)
File Name: Customized Download
Description:

Create custom data files by choosing data area, product type, map projection, file format, datum, etc. A new metadata will be produced to reflect your request using this record as a base. Change to an orthometric vertical datum is one of the many options.

File Type (Deprecated): Zip
Compression: Zip

Distribution 2

CC ID: 1070604
Start Date: 2020-05
End Date: Present
Download URL: https://noaa-nos-coastal-lidar-pds.s3.us-east-1.amazonaws.com/dem/NGS_Alamagan_CNMI_Topobathy_DEM_2020_9338/index.html
Distributor: NOAA Office for Coastal Management (NOAA/OCM) (2021 - Present)
File Name: Bulk Download
Description:

Bulk download of data files in LAZ format, geographic coordinates, orthometric heights. Note that the vertical datum (hence elevations) of the files here are different than described in this document. They will be in an orthometric datum.

File Type (Deprecated): GeoTIFF
Distribution Format: GeoTIFF

URLs

URL 1

CC ID: 1070605
URL: https://noaa-nos-coastal-lidar-pds.s3.amazonaws.com/laz/geoid12b/9331/supplemental/extent_ngs_alamagan_topobathy_2020_m9331.kmz
Name: Browse graphic
URL Type:
Browse Graphic
File Resource Format: KML
Description:

This graphic displays the footprint for this lidar data set.

URL 2

CC ID: 1070606
URL: https://coast.noaa.gov/dataviewer/
Name: NOAA's Office for Coastal Management (OCM) Data Access Viewer (DAV)
URL Type:
Online Resource
File Resource Format: HTML
Description:

The Data Access Viewer (DAV) allows a user to search for and download elevation, imagery, and land cover data for the coastal U.S. and its territories. The data, hosted by the NOAA Office for Coastal Management, can be customized and requested for free download through a checkout interface. An email provides a link to the customized data, while the original data set is available through a link within the viewer.

URL 3

CC ID: 1071715
URL: https://noaa-nos-coastal-lidar-pds.s3.amazonaws.com/laz/geoid12b/9331/supplemental/NOAA_Typhoon_Yutu_Topobathymetric_Lidar_Data_Report.pdf
Name: Project Report
URL Type:
Online Resource
File Resource Format: pdf
Description:

Project Report

Technical Environment

Description:

OS Independent

Data Quality

Horizontal Positional Accuracy:

>Lidar horizontal accuracy is a function of Global Navigation Satellite System GNSS derived positional error, flying altitude, and INS derived attitude error. The obtained RMSEr value is multiplied by a conversion factor of 1.7308 to yield the horizontal component of the National Standards for Spatial Data Accuracy NSSDA reporting standard where a theoretical point will fall within the obtained radius 95 percent of the time ACCr. This project was collected using a Leica HawkEye 4X(HE4X) which is a Chiroptera 4X (CH4X) sensor with an additional Leica 40kHz deep bathymetric channel. The sensor combination was dually mounted resulting in two seperate GNSS/IMU combinations. For the Leica Chiroptera 4x data, based on a maximum flying altitude of 600 meters, an IMU error of 0.007 decimal degrees, and a GNSS positional error of 0.015 meters, the RMSEr value for the CNMI area is 0.132 meters, with a ACCr of 0.23 meters at the 95% confidence level. For the additional mounted Leica 40kHz deep bathy channel data, based on a maximum flying altitude of 600 meters, an IMU error of 0.007 decimal degrees, and a GNSS positional error of 0.015 meters, the RMSEr value for the CNMI area is 0.132 meters, with a ACCr of 0.23 meters at the 95% confidence level. The project specification requires horizontal positions to be accurate to 1.0mRMSE.

Vertical Positional Accuracy:

Ground survey points were not able to be collected in the CNMI due to the remoteness of the islands. Thus, the absolute accuracy was not able to be assessed for the Commonwealth Northern Marina Islands (CNMI) of Anatahan, Alamagan, Guguan, or Sarigan. Data were produced to meet the accuracy specifications in planning the project. Additionally, the acquisition of the island of Guam was planned in parallel to the CNMI area acquisition and can be a good reflection of the quality of data for the CNMI area. Below are the absolute vertical accuracy statistics that were calculated for the island of Guam.

Absolute accuracy was assessed using both Non-Vegetated Vertical Accuracy (NVA) and Vegetated Vertical Accuracy (VVA) survey methods. Survey checkpoints were evenly distributed as feasible throughout the project area. NVA compares known ground check point data that were withheld from the calibration and post-processing of the lidar point cloud to the derived gridded bare earth DEM. NVA is a measure of the accuracy of lidar point data in open areas with level slope (less than 20 degrees) where the lidar system has a high probability of measuring the ground surface and is evaluated at the 95% confidence interval (1.96*RMSE). In the Guam area, 29 survey checkpoints were used to assess the non-vegetated vertical accuracy. Project specifications require NVA meet 0.196 m accuracy at the 95% confidence interval. Ground check points located in land cover categories other than bare earth, urban, or submerged topography were used to compute the Vegetated Vertical Accuracy (VVA). In Guam area, 13 survey checkpoints were used to assess the vegetated vertical accuracy. Project specifications require VVA meet 0.36 m based on the 95th percentile tested against the derived bare earth DEM. Please refer to the NOAA Typhoon Yutu Topobathymetric lidar final data report for final accuracies, to be provided upon project completion.

The Guam area dataset Non-Vegetated Vertical Accuracy tested 0.092 m vertical accuracy at 95% confidence level against the derived bare earth DEM in open terrain using 29 ground check points, based on RMSEz (0.047 m) x 1.9600.

The Guam dataset Vegetated Vertical Accuracy tested 0.284 m vertical accuracy at the 95th percentile against the derived bare earth DEM using 13 landclass points.

Completeness Report:

Data covers the project boundary.

Conceptual Consistency:

Not applicable

Data Management

Have Resources for Management of these Data Been Identified?: Yes
Approximate Percentage of Budget for these Data Devoted to Data Management: Unknown
Do these Data Comply with the Data Access Directive?: Yes
Actual or Planned Long-Term Data Archive Location: NCEI-CO
How Will the Data Be Protected from Accidental or Malicious Modification or Deletion Prior to Receipt by the Archive?:

Data is backed up to tape and to cloud storage.

Lineage

Sources

Acquisition and Processing

CC ID: 1295389
Contact Role Type: Originator
Contact Type: Organization
Contact Name: Woolpert

Process Steps

Process Step 1

CC ID: 1295390
Description:

The CH4X sensor was mounted in a Leica PAV100 gyro-stabilized mount integrated with a NovAtel SPAN GNSS and LCI-100C IMU. Real time navigation and GNSS/IMU data logging was provided by Leica FlightPro software. Lidar data were logged on the Airborne Hydrography, AB (AHAB) operator console.

Physical mounting offsets between the GNSS antenna, IMU and gyro-stabilized mount were determined through a combination of manual measurements and iterative processing in NovAtel Inertial Explorer software. Manual measurements were taken from the GNSS antenna to the reference point on the IMU in the CH4X sensor head. These measurements are added to the known offset between the IMU reference point and the rotation center of the gyro-stabilized mount to calculate the preliminary offset between the GNSS antenna and sensor reference point. This preliminary value was then used to seed the post-processing software which, through an iterative computation, uses the dynamic accelerations and rotations during flight to refine the offsets. Once the solution converges, the final offsets are entered into the flight management software and used in subsequent post-processing of the GNSS/IMU data for final trajectories.

Process Date/Time: 2020-05-04 00:00:00
Process Contact: NGS Communications and Outreach Branch
Phone (Voice): (301) 713-3242
Email Address: ngs.infocenter@noaa.gov

Process Step 2

CC ID: 1295391
Description:

Field calibration of the HE4X system is carried out to eliminate systematic errors by calculating corrections for boresight errors, scanner angle errors, remaining IMU angle errors and any necessary internal timing errors.

a. 2 x Line A over mixed terrain with flat or gentle slopes and features such as peaked roof buildings (1 x each direction)

b. 1 x Line B offset +50% from Line A in one direction

c. 1 x Line C offset -50% from Line A in the same direction as Line B

d. 2 x Line D orthogonal to previous lines (1 x each direction)

A set of calibration lines were acquired at 500m and 400m altitude for the Niue calibration; and at 600m, 500m, and 400m for the Saipan calibration. Persistent high clouds in the area prevented collection of calibration lines at higher altitudes. All sets of lines are used to calibrate and verify the topographic lidar, while the 500m and 400m lines are used for the bathymetric lidar.

Calibration values are calculated using the automatic calibration routine within the Leica Lidar Survey Studio (LSS) software. This utility first identifies patches or areas of gentle slope within the overlap region of all the lines to use for calibration. Patch selection prevents areas of vegetation, side of cars or buildings, from being used in the calibration process. Next, the utility compares the front side and back side of the elliptical scan within the same line, as well as comparing all lines to each other, to identify suitable calibration parameters such that data within the patches match. The procedure is iterative and continues until the best possible solution is computed.

Calibration for each channel (topo and shallow) is done independently. Topo channel calibration was computed using 500m altitude lines. The 600m and 400m lines were then used for verification. Calibration of the shallow channel were computed using 500m altitude. Any lower altitude data were used for verification.

At each step of the calibration process, quality assurance is conducted to ensure values being calculated are valid. This is done using the Leica LSS Quality Control Utility. Two types of checks are done; firstly, the front scan is compared to the back scan for every line. Secondly, a single line is chosen as a baseline and is compared to every other line. We would expect the average errors from both of these checks to be small; less than. In addition, the data is visually reviewed. In particular, features are studied to ensure lines from different directions show structures in the same position, in other words, verifying horizontal accuracy is maintained. These tests all provide assurance of relative accuracy.

For this project, calibration lines were acquired over Niue International Airport on January 9, 2020, and Garapan, Saipan on July 14, 2020. Results are good and indicate that both calibrations were successful. Values computed were used for all data collected during a given mobilization period.

Woolpert acquired a set of ground truth points with Real Time Kinematic (RTK) GNSS within the calibration areas. Ground truth is not used within the automatic calibration routine; however, a comparison to the lidar data was used to verify absolute accuracy. Results show data is well within required accuracy specifications.

Process Date/Time: 2021-05-04 00:00:00

Process Step 3

CC ID: 1295392
Description:

For this project, the flight parameters were used to provide 100% coverage using a 15% sidelap for the topo and 20% sidelap for the bathymetry.

During acquisition, flight lines are shown on a pilot’s display, and the aircraft is controlled by the pilot at all times. The HE4X system includes a NovAtel SPAN GNSS system with an LCI-100C IMU for aircraft position and orientation. One IMU is in the main Chiroptera sensor head, which includes the topo channel, shallow channel and RCD30 camera. Information from this IMU are also used in real-time by the PAV100 gyro-stabilized mount to compensate for deviations in pitch and roll. A second IMU is contained within the deep channel sensor head, installed over a second hatch in the aircraft. This head does not include a gyro-stabilized mount. Aircraft bank angles were restricted to 20º to avoid any potential GNSS dropouts. No flights were planned if the PDOP was expected to go above 3.0.

Data were monitored for quality during acquisition using the Operators Console running on the AHAB collection computer. The operator monitored system status of the scanners and receivers, waveforms, camera images, data coverage, flight lines and the health of the navigation system.

All data were recorded to a removable solid-state hard disk. At the end of the flight the hard disk was removed and taken to the field office where data was copied on to backup disks for transmittal back to the main processing office. Data was reviewed daily in the field for quality and coverage.

second hatch in the aircraft. This head does not include a gyro-stabilized mount. Aircraft bank angles were restricted to 20º to avoid any potential GNSS dropouts. No flights were planned if the PDOP was expected to go above 3.0.

Data were monitored for quality during acquisition using the Operators Console running on the AHAB collection computer. The operator monitored system status of the scanners and receivers, waveforms, camera images, data coverage, flight lines and the health of the navigation system.

All data were recorded to a removable solid-state hard disk. At the end of the flight the hard disk was removed and taken to the field office where data was copied on to backup disks for transmittal back to the main processing office. Data was reviewed daily in the field for quality and coverage.

Process Date/Time: 2020-05-04 00:00:00

Process Step 4

CC ID: 1295393
Description:

Position and orientation data were acquired in the aircraft using a NovAtel SPAN with LCI-100C IMU. All data were post-processed using NovAtel Inertial Explorer software to provide a tightly coupled position and orientation solution.

Due to the distance of Alamagan, Anatahan, Guguan, and Sarigan from the single base station on Guam and their remoteness a precise point positioning (PPP) solution was used for them on ITRF2014.

Initial data coverage analysis and quality checks to ensure there were no potential system issues were carried out in the field prior to demobilization of the sensor. Final processing was conducted in Woolpert’s offices.

In general, data were initially processed in Leica’s Lidar Survey Studio (LSS) using final processed trajectory information. LAS files from LSS were then imported to a Terrascan project where spatial algorithms were used to remove noise and classify bare earth/ground. Manual review was conducted in both Terrascan and LP360 prior to product creation.

Final trajectory data were post processed in NovAtel Inertial Explorer. Inertial Explorer accounts for the fixed offset between the reference point and IMU and uses a multi-pass algorithm to compute a tightly coupled solution. A GNSS base station or precise point positioning was used for processing. Average Forward and Reverse Separation RMS for the project was 0.019m in Easting and Northing, and 0.042m in Height.

Process Date/Time: 2020-05-04 00:00:00

Process Step 5

CC ID: 1295394
Description:

Lidar processing was conducted using the Leica Lidar Survey Studio (LSS) software. Calibration information, along with processed trajectory information were combined with the raw laser data to create an accurately georeferenced lidar point cloud for the entire survey in LAS v1.4 format. All points from the topographic and bathymetric laser include 16-bit intensity values.

During this LSS processing stage, an automatic land/water discrimination is made for the bathymetric waveforms. This allows the bathymetric (green) pulses over water to be automatically refracted for the pulse hitting the water surface and travelling through the water column, producing the correct depth. Another advantage of the automatic land/water discrimination is that it permits calculation of an accurate water surface over smaller areas, allowing simple bathymetric processing of smaller, narrower streams and drainage channels. Sloping water surfaces are also handled correctly.

Prior to processing, the hydrographer can adjust waveform sensitivity settings dependent on the environment encountered and enter a value for the refraction index to be used for bathymetry. The index of refraction is an indication of the water type. A value 1.3423 was used for the index of refraction, indicating salty water.

In the field, default waveform sensitivity settings were used for processing. In order to determine the optimal waveform sensitivity settings for final processing, sample areas were selected and processed with multiple different settings, to iteratively converge on the best possible settings. This is done by reviewing the processed point cloud and waveforms within sample areas. Settings affect which waveform peaks are classified as valid seabed, and which peaks are classified as noise. Optimal settings strike a balance between the amount of valid data that is classified as seabed bottom, and the amount of noise that is incorrectly classified due to peaks in the waveforms. Ideally all valid data is selected, while only a small amount of noise remains to be edited out. Once optimal threshold settings were chosen, these were used for the entire project.

It is important to note that all digitized waveform peaks are available to be reviewed by the hydrographer; both valid seabed bottom and peaks classed as noise. This allows the hydrographer to review data during Terrascan and LP360 editing for valid data such as objects that may have been misclassified as noise.

LSS processing produced LAS files in 1.4 format.

Additional QC steps were performed prior to import to Terrascan. Firstly, the derived water surface was reviewed to ensure a water surface was correctly calculated for all bathymetry channels. No significant issues were apparent. Spot checks were also made on the data to ensure the front and back of the scans remained in alignment and no calibration or system issues were apparent prior to further data editing in Terrascan.

LSS stores data in multiple LAS files for a single flight line. Each file corresponds to a single .dat file from the raw airborne data. Woolpert merged these multiple files into a single file per flight line and moved data into a standard class definition in preparation for data editing using proprietary scripts within SAFE’s FME software.

Process Date/Time: 2020-05-04 00:00:00

Process Step 6

CC ID: 1295395
Description:

After data were processed in LSS and the data integrity reviewed. Anatahan, Alamagan, Guguan, and Sarigan were transformed from the ITRF2014 ellipsoid to the NAD83(2011) Epoch 2010 ellipsoid using VDatum. With the entire project now on the correct ellipsoid, data were organized into tiles within a Terrascan project. The tile layout is the same as that used for the imagery and is provided with the project deliverables.

Data classification and spatial algorithms were applied in Terrasolid’s Terrascan software. Customized spatial algorithms, such as isolated points and low point filters, were run to remove gross fliers in the topographic and bathymetric data. A grounding algorithm was also run on the topographic data to distinguish between points representing the bare earth, and other valid topo lidar points representing features such as vegetation and buildings. Algorithms were run on the entire dataset.

Data were reviewed manually to reclassify any valid bathy points incorrectly identified by the automated routines in LSS as invalid, and vice versa. In addition, any topo points over the water were reclassified to Water to correct the ground representation. Manual editing was conducted both in Terrascan and LP360. Steps for manual editing included:

• Re-class any topo unclassified laser data and bathy seabed data from the water surface to a water surface class

• Review bathymetry in cross section.

o Re-class suitable data to Seabed (Class 40).

o Re-class any noise in the bathy ground class to bathy noise (Class 45).

o Re-class man made submerged objects to Submerged Object (Class 43).

• Review topo ground points in areas of gaps or spikes.

o Add points to ground (Class 2) if points are available to fill gaps in the ground model.

o Re-class any noise in the ground class to Topo Unclassified (Class 1) if valid vegetation or other feature, or Noise if the point is not valid Low Noise (Class 7).

Once editing was completed in TerraScan the data was vertically transformed to the GUVD04 datum using GEOID12B.

Digital Elevation Models (DEMs) were then created using TerraScan at 1m resolution using the Topo Ground, Seabed, and Submerged Object Classes (Class 2, 40, 43) at a 1m resolution on a 5000m x 5000m tile layout provided in the project deliverables.

Process Date/Time: 2020-05-04 00:00:00

Process Step 7

CC ID: 1295396
Description:

Although the bathymetry data includes intensity values, these are raw values. For intensity (reflectance) to correctly represent the reflectance of the seabed, the intensities must be normalized for any losses in signal as the light travels through the water column, so that the intensity value better reflects the intensity of the seabed itself.

One of the fundamental issues that exists with reflectance imagery is the variance in return due to water clarity differences occurring spatially along line, and temporally from day to day. This is challenging for any bathymetric lidar sensor.

If water clarity is relatively consistent along a line, then it is possible to achieve an overall homogenous reflectance image for an area. To a certain extent, variation in reflectivity intensity can be minimized by limiting the size of flight blocks and trying to ensure similar environmental parameters exist within a single flight block. In other words, where changes in water clarity or environment may be expected, flight blocks should be split to allow different normalization parameters to be used per block for the reflectance processing. Where this is not possible, and water clarity varies significantly along a line, variation in reflective intensity will be seen in the output imagery. While this imagery can still be analyzed and used for manual seabed classification, it prohibits the use of unsupervised, or semi-automated classification.

For this survey, cloud shadows (ambient light) had an effect on the resulting reflectance images.

Woolpert used proprietary in-house scripts developed in MATLAB to compute project specific correction parameters and normalize the raw intensity data for depth. This provides intensities that more closely represent the reflectance of the actual seabed. Corrected values were used to create 1m reflectance images per flightline using Applied Imagery’s QT Modeler software. Individual flightline reflectance images were then used in Trimble’s OrthoVista software to create a final reflectance image for the entire area.

Process Date/Time: 2020-05-04 00:00:00

Process Step 8

CC ID: 1295397
Description:

Quality control is carried out through every phase of the project. Several checks were used to ensure data integrity and quality was maintained. Specific statistics were generated from lidar check points and image air targets.

The airborne operator monitored system status of the scanners and receivers, waveforms, camera images, data coverage, flight lines and health of the navigation system during data acquisition. Flight logs are maintained during data acquisition. Logs not only track lines acquired, but also any relevant information on weather or water clarity, instances when sensor issues occur, and so on. These logs are a valuable resource during processing.

During acquisition, aircraft bank angles were restricted to 20º to avoid any potential GNSS dropouts. No flights were planned if the PDOP was expected to go above 3.0. Separation plots and additional statistics were reviewed for each flight trajectory processed.

Throughout data editing adjacent survey lines of data are compared to ensure there are no data busts, or system artifacts. During processing Terrasolid’s TMatch software is run to examine the Delta Z differences between overlapping lines. If differences are greater than 0.02m, then a simple Z correction is applied per flight line to remove any vertical differences between flight lines. TMatch can then be run again once all corrections are applied to ensure adjacent lines agree within specification. This provides a measure of inter-swath accuracy.

Due to the rugosity and canopy cover of Anatahan, Alamagan, Guguan, and Sarigan a reliable comparison to adjacent lines was not obtained.

Check points were not obtained on Anatahan, Alamagan, Guguan, or Sarigan due to the remoteness of the islands.

Process Date/Time: 2020-05-04 00:00:00

Catalog Details

Catalog Item ID: 64808
GUID: gov.noaa.nmfs.inport:64808
Metadata Record Created By: Maryellen Sault
Metadata Record Created: 2021-06-23 11:19+0000
Metadata Record Last Modified By: Kirk Waters
Metadata Record Last Modified: 2024-01-10 19:17+0000
Metadata Record Published: 2024-01-10
Owner Org: NGS
Metadata Publication Status: Published Externally
Do Not Publish?: N
Metadata Last Review Date: 2022-03-16
Metadata Review Frequency: 1 Year
Metadata Next Review Date: 2023-03-16