Planet Imagery Geometric Assessment
Explore the source record for details and available documents.
Engineering topics
Publications and source records attributed to Guoqing (Gary) Lin.
Explore the source record for details and available documents.
Atmospheric motion vectors (AMVs), derived by tracking patterns, represent the winds in a layer characteristic of the pattern. AMV height (or pressure), important for applications in atmospheric research and operational meteorology, is usually assigned using observed IR brightness temperatures with a modeled atmosphere and can be inaccurate. Stereoscopic tracking provides a direct geometric height measurement of the pattern that an AMV represents. We extend our previous work with multi-angle imaging spectro–radiometer (MISR) and GOES to moderate resolution imaging spectroradiometer (MODIS) and the GOES-R series advanced baseline imager (ABI). MISR is a unique satellite instrument for stereoscopy with nine angular views along track, but its images have a narrow (380 km) swath and no thermal IR channels. MODIS provides a much wider (2330 km) swath and eight thermal IR channels that pair well with all but two ABI channels, offering a rich set of potential applications. Given the similarities between MODIS and VIIRS, our methods should also yield similar performance with VIIRS. Our methods, as enabled by advanced sensors like MODIS and ABI, require high-accuracy geographic registration in both systems but no synchronization of observations. AMVs are retrieved jointly with their heights from the disparities between triplets of ABI scenes and the paired MODIS granule. We validate our retrievals against MISR-GOES retrievals, operational GOES wind products, and by tracking clear-sky terrain. We demonstrate that the 3D-wind algorithm can produce high-quality AMV and height measurements for applications from the planetary boundary layer (PBL) to the upper troposphere, including cold-air outbreaks, wildfire smoke plumes, and hurricanes.
Explore the source record for details and available documents.
Sensor spatial resolution (SSR) is a determining factor in developing remote sensing instruments for earth observations. The SSR is defined as the full width at half maximum (FWHM) of sensor point/line spread function (PSF/LSF). This presentation first discusses large variations in image absolute SSR and SSR ratios to ground sampling distance (GSD) or grid pixel size from multiple instruments including AVHRR (Advanced Very-High Resolution Radiometer), MODIS (Moderate Resolution Imaging Spectroradiometer), VIIRS (Visible Infrared Imaging Radiometer Suite), WorldView and Planet. When the physics based SSR is more than twice the nominal grid pixel size, a more appropriate grid pixel size is proposed.
Commercial companies such as Maxar, PlanetScope, and BlackSky have launched many satellites with revisit times ranging from two hours to one day and have built large archives of high resolution (1-3 m) Earth observing data. The high temporal resolution and global coverage of these satellites makes commercial satellite images ideal for scientists studying rapidly changing processes such as flooding and fires. A key step in assessing changes in these images is co-registration of various images to each other. High geolocation accuracy (at most 0.5 pixels of offset) allows for easy co-registration of images across different times and sensors. Here, we assess the latest PlanetScope imagery and evaluate it for geolocation accuracy. For our assessment, we compare the target image to a reference image with known geolocation accuracy (WorldView imagery) and determine the offset between these images by shifting them to maximize their Pearson Cross-Correlation (PCC) value. Offsets required to maximize the PCC give the geolocation accuracy of the target image relative to the reference image. Previously, our global assessment of Planet data revealed large variability in geolocation accuracy from one continent to another. The measured root mean squared errors (RMSEs) range from 5.4 m in North America to 14.9 m in Africa. We have updated this previous assessment to include both their newest SuperDove series as well as a more robust assessment of the temporal stability of PlanetScope's geolocation accuracy.
When end-users are using multiple satellite images for a study the first step is often the assessment and rectification of geolocation offsets. We have developed an automated algorithm to calculate the geometric offsets between two satellite images. It splits both images into subset images (chips) and resamples them to a common resolution. Offsets are imposed on matching chips and the Pearson Cross Correlation (PCC) value calculated for each offset, where the best PCC gives the true offsets. With the true offsets corrected for, a quality of image registration is calculated, measurement uncertainty (MU). The MU threshold value is non-uniform across sensors because it depends on spectral characteristics of each sensor. Previously, the MU threshold was determined by manual inspection for each new sensor assessed. In this study, we automated MU threshold determination. 1/MU vs calculated offset approximates a Gaussian function in shape, see abstract image. A Gaussian function is fit to this dataset, and the root of the third derivative is calculated and applied as a quality threshold, removing the data of lower quality from the final assessment of image offsets. We apply this algorithm to commercial satellite data for proof of concept.
Explore the source record for details and available documents.
Science users of commercial satellite data build their studies on the properties of the satellite data they work with, pixel size being one of the determining factors. However, pixel size and nominal image footprint size can be different, impacting the scale of features discernable in the satellite images. Here, we assess geometric properties such as nominal image footprint size for Planet Labs' SuperDove series, MAXAR's WorldView series, and potentially BlackSky. Additionally, we assess temporal change in nominal footprint size for the SuperDove series. Nominal Image footprint size is assessed over the CalVal sites in Baotou, China; Shadnagar, India; and Big Spring, TX. Edge spread functions are constructed along the black/white transitions, from which the line spread functions are constructed and modulation transfer functions of the scene are estimated. Scene full width half maximum (FWHM), which represents sensor footprint size, is estimated from the line spread function. The average nominal footprint size is 3.3 pixels for SuperDove series, 1.5 pixels for WorldView-2, and 1.3 pixels for WorldView-3. The SuperDove series nominal footprint size improves with time in orbit, from an average of 3.4 pixels soon after launch to an average of 3.2 pixels one or more years after launch.
Explore the source record for details and available documents.
NASA Terra satellite has provided earth observations from its 5 onboard instruments for 25 years since its launch in December 1999. All these instruments are still mostly healthy, and observations are expected to extend 2 more years of life into late 2026. The onboard fuel is preserved by ceasing inclination adjust maneuvers (IAMs, stopped 2020) and drag make-up maneuvers (DMUs, stopped 2022). Without IAMs, the orbital plane is drifting, changing the sunlight conditions by late 2026 from nominal 10:30 AM to 8:30 AM local time at descending node (LTDN). The ceasing of DMUs combined with a constellation exit maneuver in October 2022 result in orbit altitude drifting down. Altitude will drift from a nominal mean of 715 km to 702 km by late 2026. A companion earth observing satellite, Aqua, is going through a similar orbit drift process. This presentation shows the details of the drifts, that have potentials to impact data products, especially those products from the Moderate Resolution Imaging Spectroradiometer (MODIS) instruments on Terra and Aqua satellites.
Some satellite data is delivered in images with gridded pixels. This gridded pixel size is often assumed to be the spatial resolution of the satellite sensor; however, this is not always the case. An image can be grided to any arbitrary pixel size, but the sensor resolution will remain constant. For example, an image with a pixel grid size much smaller than the sensor resolution will appear blurry along what should be sharp transitions. This discrepancy between an image’s pixel size and true sensor spatial resolution can be the source of much confusion and even misinformation among data users, which may lead them to waste time and resources on using images that do not suit their spatial resolution needs. This presentation will highlight our evaluation of the true spatial resolution of various government and commercial images in the pixel size range of 0.3 m to 60 m. Images evaluated include ESA’s Sentinel-2 (60 m, 20 m, 10 m pixels), USGS’s Landsat 8/9 (30 m & 15 m pixels), Planet’s SuperDoves (3 m pixels), BlackSky’s Globals (~1 m pixels), and the optical bands of Maxar’s WorldView-2 (2.4 m – 0.41 m pixels) and WorldView-3 (1.38 m – 0.31 m pixels). Our evaluation of true sensor spatial resolution, or ‘footprint size’ is based on the sensor’s line spread function (LSF). We calculate the width at half the height of the LSF to find the full width at half maximum (FWHM). The FWHM is how we report sensor spatial resolution. Different objects are examined for constructing the LSF depending on the sensor spatial resolution. Coarser resolution sensors in this evaluation such as Sentinel-2 and Landsat 8/9 are examined at bridges over a dark water background. The bright bridge acts as a line impulse, giving a sensor’s line spread function (LSF) in one direction. Additionally, we simulate the impacts of bridge width on the apparent LSF to obtain a true LSF without the effects of bridge width for these sensors. Finer resolution sensors will image the irregularities in bridges such as trusses, sidewalks, and in some cases painted lines, interfering with the LSF construction. Instead, these sensors are evaluated at large (60 m – 140 m) black and white checkerboards known as Cal/Val sites. At these locations, the image’s transition from black to white is extracted as an edge spread function (ESF). We calculate the derivative of this ESF to obtain the sensor’s LSF. From there, we find the FWHM as we do for the coarser resolution images. With the FWHM and pixel size, we determine how over- or under-sampled the images are. When the ratio of a sensor’s spatial resolution and the gridded image’s pixel size is less than 1, the image is considered under-sampled. In this case, each pixel’s information is unique but only a portion of that pixel’s ground area has been measured. On the other side, if the ratio is greater than 1, the image is considered over-sampled. That is, each pixel’s information is sourced from within the ground extent of the pixel and some extent outside additionally. We will show the true spatial resolution and the extent of over-/under-sampling in the imagery from ESA’s Sentinel-2 (60 m, 20 m, 10 m pixels), USGS’s Landsat 8/9 (30 m & 15 m pixels), Planet’s SuperDoves (3 m pixels), BlackSky’s Globals (~1 m pixels), and the optical bands of Maxar’s WorldView-2 (2.4 m – 0.41 m pixels) and WorldView-3 (1.38 m – 0.31 m pixels).