USGS ScienceSearch

USGS · 70279336

Satellite embeddings for crop type classification: A comparative examination

Abstract

Embedding datasets encode complex relationships among multiple sources of Earth observation data into a compact format. Here, we evaluated the utility of a 10-m global Satellite Embedding product (SE) for classifying crop types in central California for the year 2020. We compared the classification accuracy of a random forest model based exclusively on the SE layer to an existing random forest model with multiple imagery inputs. Our results showed the SE-based classification had higher agreement with the reference dataset (California Department of Water Resources crop map) than the classification based on Landsat and National Agricultural Imagery Program inputs (94.7% versus. 91.9% overall accuracy, respectively). The performance of individual crop types was consistent across models, ranging from high agreement for rice (98.4% versus 98% accuracy) to lower agreement for pasture, grain, and fallow/young perennial classes (< 65% accuracy in both models). The SE-based workflow used three times less cloud-based computational resources and represented substantial savings of predictor development time. Geospatial embedding products can aid classification efforts by reducing predictor development time and processing demands while maintaining classification accuracy.

Explore related subjects

90° N90° S · 180° W ← longitude → 180° E
Source-reported bounding extent: 34.4221364° to 40.7017767° latitude; -124.1894753° to -117.889804° longitude. This indicates report coverage, not an exact sampling location. View area on OpenStreetMap.

Keep this discovery

Explore connections, maps & timelines

BibTeXRIS

Britt Windsor Smith, Jessica J. Walker, Christopher E. Soulard. 2026-08-25. Satellite embeddings for crop type classification: A comparative examination. https://doi.org/10.1016/j.jag.2026.105531

Cite the original work for its findings. Save a collection to share your selection of sources.

KEEP EXPLORING

Related USGS reports

The impacts of cover crop biomass on satellite-based detectability of cover crops in Maryland

Cover crop adoption in the U.S. has increased over the past decades, increasing the need to quantify their performance and environmental benefits. While remote sensing (RS)-based approaches for detecting cover crop presence have been developed, there has been limited research on how cover crops with varied biomass and management practices influence detectability. Using unique field-level cover crop presence and biomass datasets in Maryland, U.S., we investigated how RS-based detectability changes for cover crops with varied aboveground biomass, planting, and termination dates. Specifically, we proposed a time-integrated satellite-based greenness feature from Harmonized Landsat-8 and Sentinel-2 (HLS) time series from 2017 to 2021 to estimate biomass of cover crops and evaluated their detectability using a phenology-based cover crop detection framework. The impacts of cover crop planting and termination dates on cover crop biomass were also analyzed. Our results demonstrated that Normalized Difference Vegetation Index (NDVI) estimated biomass with higher accuracy compared to other vegetation indices, and the time-integrated model estimated biomass with higher accuracy than the single-date “snapshot” linear model (R2 from 0.53 to 0.66 and RMSE from 922 kg/ha to 747 kg/ha). While the snapshot models were species sensitive, the time-integrated models showed strong robustness across different species. Detectability increased with cover crop biomass, as detected cover crops averaged 963.3 ±719.5 kg/ha compared to 297.2 ± 209.0 kg/ha for non-detected cover crops. Detection accuracy reached 96.1% for fields exceeding 500 kg/ha, compared with an overall accuracy of 62.7%. Earlier planting and later termination increased biomass and detectability, with biomass rising by 4.14 kg/ha/day (p < 0.01). This study demonstrates how management practices affect cover crop biomass and detectability via satellite time series and provides insights that can inform management of cover crops and monitoring of their effects on agroecosystems.

Maryland

Aligning legacy NLCD land cover maps based on Landsat Collection 1 to Collection 2

The transition from Landsat Collection 1 to Collection 2 introduced significant improvements in radiometric and geometric accuracy. However, the improvements cause location misalignment between the existing Landsat-derived land cover products and the new collection. The legacy National Land Cover Database (NLCD) has been used as a cornerstone land cover source for a variety of research. Therefore, a method aligning the legacy NLCD product to Collection 2 is required to ensure its continuity and consistency of service. We developed a strategy to not only align legacy NLCD to match new Collection 2 geometric locations but also improve land cover labeling in the region that was affected by the geometric shifts. The method identifies boundary pixels of homogeneous land cover patches as potential problem areas that are likely impacted by geometric shifts and generates candidate labels from 3 × 3 window with the target pixel at the center and segmentation-derived majority label. Standard phenology patterns of each candidate land cover type are established based on the random samples except boundary pixels within a 1000-pixels × 1000-pixels processing window region. The phenological distance to each standard land cover type pattern is calculated through a penalty dynamic time warping (DTW) method for each target pixel in the boundary region. Finally, the method determines the most suitable label based on the phenological distance from the candidate labels. Both visual and accuracy assessment results demonstrate that the alignment preserves the overall land cover patterns in the original legacy NLCD product while reducing the spatial discrepancies between the Landsat Collection 2 and land cover. In addition, it enhances the accuracy of land cover labeling of boundary pixels. The overall accuracy (OA) was increased by 7% in the land cover boundary regions after alignment. The quality and confusion matrix comparison between the alignment results and the original legacy NLCD confirm the reliability of the method. Our alignment method has the potential to serve as a framework for aligning other Landsat-derived land cover products to future collections.

conterminous United States

An automated compositing method for producing annual clear images from Landsat Collection 2 for annual NLCD production

Quality image input is fundamental to the quality of derived land cover products. Substantial time and effort are usually required to prepare images. Here, we present a novel and streamlined compositing algorithm that ingests Landsat Collection 2 Analysis Ready Data (ARD) and outputs cloud-free and gap-free composite imagery, which can be directly used for classification. This method leverages and improves the previous National Land Cover Database (NLCD) Virtual Median Value Point (VMVP) compositing method, the first part of the image preparation for NLCD 2019 operational production. The NLCD 2019 image preparation approach includes a second part, a residual cloud and cloud shadow detection and gap-filling method, to produce final cloud-free and gap-free composite imagery. The second part requires one clear reference image for each target year. Additional reference images are needed for producing reasonable observations for perennial ice/snow areas because Pixel QA (Quality Assessment) from ARD has difficulties differentiating ice/snow areas from clouds. Unlike the NLCD 2019 image preparation approach, our new compositing method, which is referred to as Automated VMVP (AVMVP), uses Landsat ARD as the only input and does not require reference images and extra steps. In this method, we developed new spectral filter criteria coupled with counts of clear observations using Pixel QA to identify potential cloud and cloud shadow observations on initially selected observations from the NLCD VMVP compositing algorithm. We also automate “gap-filling” using clear observations retrieved from a maximum of ±2 years around the target year when needed. Finally, a percentile-filtered compositing method was developed for the perennial ice/snow areas. All these steps are streamlined, pixel-based, and directly run on Landsat Collection 2 ARD. We have run successful tests on the conterminous United States (CONUS). Composite images derived from our innovative method were used to produce the CONUS Annual NLCD Collection 1 product suite that covers the period from 1985 to 2023.

conterminous United States