The estimation of depth in optically shallow waters using satellite imagery can be efficient and cost-effective. Active sensors measure the distance traveled by an emitted laser pulse propagating through the water with high precision and accuracy if the bottom peak intensity of the waveform is greater than the noise level. However, passive optical imaging of optically shallow water involves measuring the radiance after the sunlight undergoes downward attenuation on the way to the sea floor, and the reflected light is then attenuated while moving back upward to the water surface. The difficulty of satellite-derived bathymetry (SDB) arises from the fact that the measured radiance is a result of a complex association of physical elements, mainly the optical properties of the water, bottom reflectance, and depth. In this research, we attempt to apply physics-based algorithms to solve this complex problem as accurately as possible to overcome the limitation of having only a few known values from a multispectral sensor. Major analysis components are atmospheric correction, the estimation of water optical properties from optically deep water, and the optimization of bottom reflectance as well as the water depth. Specular reflection of the sky radiance from the water surface is modeled in addition to the typical atmospheric correction. The physical modeling of optically dominant components such as dissolved organic matter, phytoplankton, and suspended particulates allows the inversion of water attenuation coefficients from optically deep pixels. The atmospheric correction and water attenuation results are used in the ocean optical reflectance equation to solve for the bottom reflectance and water depth. At each stage of the solution, physics-based models and a physically valid, constrained Levenberg–Marquardt numerical optimization technique are used. The physics-based algorithm is applied to Landsat Operational Land Imager (OLI) imagery over the shallow coastal zone of Guam, Key West, and Puerto Rico. The SDB depths are compared to airborne lidar depths, and the root mean squared error (RMSE) is mostly less than 2 m over water as deep as 30 m. As the initial choice of bottom reflectance is critical, along with the bottom reflectance library, we describe a pure bottom unmixing method based on eigenvector analysis to estimate unknown site-specific bottom reflectance.

This report addresses system characterization of the Maxar WorldView-3 satellite and is part of a series of system characterization reports produced and delivered by the U.S. Geological Survey Earth Resources Observation and Science Cal/Val Center of Excellence in 2020. These reports present and detail the methodology and procedures for characterization; present technical and operational information about the specific sensing system being evaluated; and provide a summary of test measurements, data retention practices, data analysis results, and conclusions.WorldView-3 is a high-resolution multispectral satellite launched in 2014 by Maxar Technologies on an Atlas V launch vehicle from Vandenberg Air Force Base in California for Earth resources monitoring. WorldView-3 provides substantial technical improvements to previous WorldView satellites, including spectral bands, ground sample distance, and swath. The WorldView-3 satellite was designed and built by Lockheed Martin for Maxar Technologies using the BCP–5000 bus with the WorldView-3 Imager and the Clouds, Aerosols, Vapors, Ice, and Snow sensor. The high-resolution WorldView-3 Imager is the main instrument, and the Clouds, Aerosols, Vapors, Ice, and Snow sensor provides additional data on obscurants and other atmospheric effects used in data production. More information on Maxar WorldView satellites and sensors is available within the “2020 Joint Agency Commercial Imagery Evaluation—Remote Sensing Satellite Compendium” and from the manufacturer at https://www.maxar.com/.The Earth Resources Observation and Science Cal/Val Center of Excellence system characterization team completed data analyses to characterize the geometric (interior and exterior), radiometric, and spatial performances. Results of these analyses indicate that WorldView-3 has a range of interior geometric performance of −0.09 (−0.07 pixel) to 0.24 meter (0.19 pixel) in band-to-band registration; an exterior geometric performance in the range of a −21.10- (−2.11 pixels) to 28.23-meter (2.82 pixels) offset in comparison to Sentinel-2; a radiometric performance in the range of 0.007 to 0.874 (offset and slope); and a spatial performance in the range of 1.2 to 1.7 pixels at full width at half maximum with a modulation transfer function at a Nyquist frequency in the range of 0.093 to 0.185.

The Leica Geosystems CountryMapper hybrid system has the potential to collect data that satisfy the U.S. Geological Survey (USGS) National Geospatial Program (NGP) and 3D Elevation Program (3DEP) and the U.S. Department of Agriculture (USDA) National Agriculture Imagery Program (NAIP) requirements in a single collection. This research will help 3DEP determine if this sensor has the potential to meet current and future 3DEP topographic lidar collection requirements. We performed an accuracy analysis and assessment on the lidar point cloud produced from CountryMapper. The boresighting calibration and co-registration by georeferencing correction based on ground control points are assumed to be performed by the data provider. The scope of the accuracy assessment is to apply the following variety of ways to measure the accuracy of the delivered point cloud to obtain the error statistics. Intraswath uncertainty from a flat surface was computed to evaluate the point cloud precision. Intraswath difference between opposite scan directions and the interswath overlap difference were evaluated to find boresighting or any systematic errors. Absolute vertical accuracy over vegetated and non-vegetated areas were also assessed. Both horizontal and vertical absolute errors were assessed using the 3D absolute error analysis methodology of comparing conjugate points derived from geometric features. A three-plane feature makes a single unique intersection point. Intersection points were computed from ground-based lidar and airborne lidar point clouds for comparison. The difference between two intersection points form one error vector. The geometric feature-based error analysis was applied to intraswath, interswath, and absolute error analysis. The CountryMapper pilot data appear to satisfy the accuracy requirements suggested by the USGS lidar specification, based upon the error analysis results. The focus of this research was to demonstrate various conventional accuracy measures and novel 3D accuracy techniques using two different error computation methods on the CountryMapper airborne lidar point cloud.

The accuracy assessment of airborne lidar point cloud typically estimates vertical accuracy by computing RMSEz (root mean square error of the z coordinate) from ground check points (GCPs). Due to the low point density of the airborne lidar point cloud, there is often not enough accurate semantic context to find an accurate conjugate point. To advance the accuracy assessment in full three-dimensional (3D) context, geometric features, such as the three-plane intersection point or two-line intersection point, are often used. Although the point density is still low, geometric features are mathematically modeled from many points. Thus, geometric features provide a robust determination of the intersection point, and the point is considered as a GCP. When no regular built objects are available, we describe the process of utilizing features of irregular shape called amorphous natural objects, such as a tree or a rock. When scanned to a high-density point cloud, an amorphous natural object can be used as ground truth reference data to estimate 3D georeferencing errors of the airborne lidar point cloud. The algorithm to estimate 3D accuracy is the optimization that minimizes the sum of the distance between the airborne lidar points to the ground scanned data. The search volume partitioning was the most important procedure to improve the computational efficiency. We also performed an extensive study to address the external uncertainty associated with the amorphous object method. We describe an accuracy assessment using amorphous objects (108 trees) spread over the project area. The accuracy results for ∆x, ∆y, and ∆z obtained using the amorphous object method were 3.1 cm, 3.6 cm, and 1.7 cm RMSE, along with a mean error of 0.1 cm, 0.1 cm, and 4.5 cm, respectively, satisfying the accuracy requirement of U.S. Geological Survey lidar base specification. This approach shows strong promise as an alternative to geometric feature methods when artificial targets are scarce. The relative convenience and advantages of using amorphous targets, along with its good performance shown here, make this amorphous object method a practical way to perform 3D accuracy assessment.

This report addresses system characterization of Planet’s SuperDove and is part of a series of system characterization reports produced and delivered by the U.S. Geological Survey Earth Resources Observation and Science Cal/Val Center of Excellence. These reports present and detail the methodology and procedures for characterization; present technical and operational information about the specific sensing system being evaluated; and provide a summary of test measurements, data retention practices, data analysis results, and conclusions.Since 2013, Planet has launched more than 360 Dove 3U CubeSats, where U stands for 10-centimeter (cm) x 10-cm x 10-cm stowed dimensions, each weighing about 5.8 kilograms. Since 2015, all Dove satellites have had four-band imagers with about a 3-meter (m) pixel ground sample distance. Since 2016, all Doves have been launched into Sun-synchronous orbits varying from 474 to 524 kilometers, with inclinations between 97 and 98 degrees. The Dove series satellites do not have orbit maintenance...

The traditional practice to assess accuracy in lidar data involves calculating RMSEz (root mean square error of the vertical component). Accuracy assessment of lidar point clouds in full 3D (three dimension) is not routinely performed. The main challenge in assessing accuracy in full 3D is how to identify a conjugate point of a ground-surveyed checkpoint in the lidar point cloud with the smallest possible uncertainty value. Relatively coarse point-spacing in airborne lidar data makes it challenging to determine a conjugate point accurately. As a result, a substantial unwanted error is added to the inherent positional uncertainty of the lidar data. Unless we keep this additional error small enough, the 3D accuracy assessment result will not properly represent the inherent uncertainty. We call this added error “external uncertainty,” which is associated with conjugate point identification. This research developed a general external uncertainty model using three-plane intersections and accounts for several factors (sensor precision, feature dimension, and point density). This method can be used for lidar point cloud data from a wide range of sensor qualities, point densities, and sizes of the features of interest. The external uncertainty model was derived as a semi-analytical function that takes the number of points on a plane as an input. It is a normalized general function that can be scaled by smooth surface precision (SSP) of a lidar system. This general uncertainty model provides a quantitative guideline on the required conditions for the conjugate point based on the geometric features. Applications of the external uncertainty model were demonstrated using various lidar point cloud data from the U.S. Geological Survey (USGS) 3D Elevation Program (3DEP) library to determine the valid conditions for a conjugate point from three-plane modeling.