coincident.io.download.fetch_usgs_lpc_tiles#

coincident.io.download.fetch_usgs_lpc_tiles(aoi, project, tnmdataset='Lidar Point Cloud (LPC)', output_dir=None)#

Query USGS TNMAccess for LPC products intersecting an AOI, assemble their bounding boxes into a GeoDataFrame (EPSG:4326), and optionally write to disk.

Steps:

  1. Reproject AOI to EPSG:4326 and explode to a single geometry.

  2. Convert that geometry into the TNM APIs polygon string format.

  3. Call query_tnm_api to retrieve all LPC products (LAS/LAZ) in the AOI.

  4. Filter the returned items by matching project via _filter_items_by_project.

  5. Build a Shapely box from each items boundingBox coordinates.

  6. Assemble into a GeoDataFrame in EPSG:4326.

  7. If output_dir is provided, save as usgs_lpc_products.geojson.

Parameters:

  • aoi (gpd.GeoDataFrame) – Area of interest.

  • project (str) – A substring or identifier used to filter product titles.

  • tnmdataset (str) – Dataset name for TNM API (default: “Lidar Point Cloud (LPC)”, changing this will break the function).

  • output_dir (Optional[str]) – Directory to write out the GeoJSON. If None, no file is written.

Returns:

Contains columns from the API plus a geometry column

with each products bounding box in WGS84.

Return type:

gpd.GeoDataFrame