Metadata-Version: 2.4
Name: seamless-3dep
Version: 0.5.1
Summary: Easy retrieval of seamless DEM from 3DEP across US.
Project-URL: Changelog, https://docs.hyriver.io/seamless-3dep/latest/CHANGELOG/
Project-URL: CI, https://github.com/hyriver/seamless-3dep/actions
Project-URL: Documentation, https://docs.hyriver.io/seamless-3dep
Project-URL: Homepage, https://docs.hyriver.io/seamless-3dep
Project-URL: Issues, https://github.com/hyriver/seamless-3dep/issues
Author-email: Taher Chegini <cheginit@gmail.com>
License: MIT
License-File: AUTHORS.md
License-File: LICENSE
Classifier: Development Status :: 4 - Beta
Classifier: Intended Audience :: Science/Research
Classifier: License :: OSI Approved :: MIT License
Classifier: Operating System :: OS Independent
Classifier: Programming Language :: Python
Classifier: Programming Language :: Python :: 3 :: Only
Classifier: Programming Language :: Python :: 3.12
Classifier: Programming Language :: Python :: 3.13
Classifier: Programming Language :: Python :: 3.14
Classifier: Topic :: Scientific/Engineering
Classifier: Topic :: Scientific/Engineering :: GIS
Classifier: Typing :: Typed
Requires-Python: >=3.12
Requires-Dist: numpy>=2
Requires-Dist: rasterio>=1.5
Requires-Dist: tiny-retriever>=0.3
Provides-Extra: dev
Requires-Dist: geopandas>=1; extra == 'dev'
Requires-Dist: ipykernel; extra == 'dev'
Requires-Dist: ipywidgets; extra == 'dev'
Requires-Dist: jupytext; extra == 'dev'
Requires-Dist: matplotlib; extra == 'dev'
Requires-Dist: nbconvert; extra == 'dev'
Requires-Dist: pytest; extra == 'dev'
Requires-Dist: pywbt>=0.2.8; extra == 'dev'
Requires-Dist: rioxarray; extra == 'dev'
Provides-Extra: docs
Requires-Dist: jupytext; extra == 'docs'
Requires-Dist: mike>=2; extra == 'docs'
Requires-Dist: mkdocs-jupyter; extra == 'docs'
Requires-Dist: mkdocs-materialx; extra == 'docs'
Requires-Dist: mkdocs>=1.6; extra == 'docs'
Requires-Dist: mkdocstrings[python]>=0.27; extra == 'docs'
Requires-Dist: ruff; extra == 'docs'
Provides-Extra: lint
Requires-Dist: codespell; extra == 'lint'
Requires-Dist: pre-commit; extra == 'lint'
Provides-Extra: test
Requires-Dist: coverage[toml]; extra == 'test'
Requires-Dist: pytest-cov; extra == 'test'
Requires-Dist: pytest-sugar; extra == 'test'
Requires-Dist: rioxarray; extra == 'test'
Requires-Dist: shapely; extra == 'test'
Provides-Extra: typecheck
Requires-Dist: pandas-stubs; extra == 'typecheck'
Requires-Dist: pyproj; extra == 'typecheck'
Requires-Dist: pyright; extra == 'typecheck'
Requires-Dist: rioxarray; extra == 'typecheck'
Requires-Dist: scipy-stubs; extra == 'typecheck'
Requires-Dist: shapely; extra == 'typecheck'
Requires-Dist: types-geopandas; extra == 'typecheck'
Requires-Dist: types-pyyaml; extra == 'typecheck'
Requires-Dist: types-shapely; extra == 'typecheck'
Requires-Dist: xarray; extra == 'typecheck'
Description-Content-Type: text/markdown

# Seamless3DEP: Streamlined Access to USGS 3DEP Topographic Data

[![PyPi](https://img.shields.io/pypi/v/seamless-3dep.svg)](https://pypi.python.org/pypi/seamless-3dep)
[![Conda Version](https://img.shields.io/conda/vn/conda-forge/seamless-3dep.svg)](https://anaconda.org/conda-forge/seamless-3dep)
[![CodeCov](https://codecov.io/gh/hyriver/seamless-3dep/branch/main/graph/badge.svg)](https://codecov.io/gh/hyriver/seamless-3dep)
[![Python Versions](https://img.shields.io/pypi/pyversions/seamless-3dep.svg)](https://pypi.python.org/pypi/seamless-3dep)
[![Downloads](https://static.pepy.tech/badge/seamless-3dep)](https://pepy.tech/project/seamless-3dep)

[![CodeFactor](https://www.codefactor.io/repository/github/hyriver/seamless-3dep/badge)](https://www.codefactor.io/repository/github/hyriver/seamless-3dep)
[![Ruff](https://img.shields.io/endpoint?url=https://raw.githubusercontent.com/astral-sh/ruff/main/assets/badge/v2.json)](https://github.com/astral-sh/ruff)
[![pre-commit](https://img.shields.io/badge/pre--commit-enabled-brightgreen?logo=pre-commit&logoColor=white)](https://github.com/pre-commit/pre-commit)
[![Binder](https://mybinder.org/badge_logo.svg)](https://mybinder.org/v2/gh/hyriver/seamless-3dep/HEAD?labpath=docs%2Fexamples)

Seamless3DEP is a lightweight Python package that simplifies access to topographic data
from the USGS
[3D Elevation Program (3DEP)](https://www.usgs.gov/core-science-systems/ngp/3dep).
Whether you need elevation data or its derivatives, Seamless3DEP provides an efficient
interface to both static and dynamic 3DEP products. For global coverage or
bathymetry-inclusive workflows, Seamless3DEP also exposes the
[NOAA global DEM mosaic](https://gis.ngdc.noaa.gov/arcgis/rest/services/DEM_mosaics/DEM_global_mosaic/ImageServer)
and a generic helper for any ArcGIS `ImageServer/exportImage` endpoint.

Seamless3DEP utilizes connection pooling across threads (safely) to optimize service
calls and minimize redundant connections. This reduces both the service load and the
time required to retrieve data, making it an ideal tool for handling large-scale
topographic data requests.

📚 Full documentation is available [here](https://docs.hyriver.io/seamless-3dep).

## Available Products

### Static DEMs

- 1/3 arc-second (10 meters)
- 1 arc-second (30 meters)
- 2 arc-second (60 meters)

### Dynamic Products

- Digital Elevation Model (DEM)
- Hillshade Derivatives:
    - Gray Hillshade
    - Multidirectional Hillshade
    - GreyHillshade with Elevation Fill
    - Hillshade with Elevation Tint
- Terrain Analysis:
    - Aspect (Degrees and Map)
    - Slope (Degrees and Map)
    - Height (Ellipsoidal)
- Contours:
    - Contour 25 (Dynamically generates 25 contours for the area of interest)
    - Contour Smoothed 25 (Smoothed version of the 25 contours)

## Core Functions

Seamless3DEP offers eight main functions designed for efficient data retrieval and
processing:

- `get_dem`: Retrieves static DEMs within a specified bounding box. The function
    automatically splits large areas into manageable tiles, downloads data as GeoTIFF
    files in EPSG:4326, and supports resolutions of 10m, 30m, or 60m. A `buff_npixels`
    option produces overlapping tiles to avoid NoData strips along seams when mosaicking
    with `tiffs_to_da`.
- `get_map`: Fetches any 3DEP product (including DEMs and terrain derivatives) for areas
    within the US. Works with all available product types, allows custom resolution
    settings, and downloads in EPSG:3857. Also accepts `buff_npixels` for overlapping
    tiles.
- `get_global_dem`: Fetches the NOAA global DEM mosaic (SRTM + GEBCO + ICESat) for any
    bounding box on Earth. Returns GeoTIFFs in EPSG:3857. Unlike 3DEP, it includes
    sub-zero bathymetric values (down to roughly -9982 m), making it the right choice
    outside the US or for coastal and marine workflows. Native resolution is
    approximately 30 m (1 arc-second).
- `get_image_server`: General-purpose downloader for any ArcGIS
    `ImageServer/exportImage` endpoint. Accepts any output CRS, an optional Esri
    rendering rule, and handles URL validation, bbox tiling, and resume-on-failure
    automatically. `get_map` and `get_global_dem` are thin wrappers over this function.
- `elevation_bygrid`: Samples elevation values from the 10 m seamless DEM at a grid of
    longitude/latitude coordinates. Reads directly from the USGS Cloud-Optimized
    GeoTIFFs (EPSG:4269) and supports configurable resampling methods including nearest,
    bilinear, cubic, cubic spline, and Lanczos.
- `decompose_bbox`: Handles large area requests by breaking down extensive bounding
    boxes into optimal sizes based on resolution and maximum pixel count, ensuring
    efficient data retrieval.
- `build_vrt`: Creates virtual raster datasets by combining multiple GeoTIFF files.
    Requires `libgdal-core` installation and supports efficient data handling for large
    areas. Note that `libgdal-core` is an optional dependency and is not installed when
    `seamless-3dep` is installed from PyPI. However, it is installed as a dependency
    when `seamless-3dep` is installed from `conda-forge`.
- `tiffs_to_da`: Converts a list of GeoTIFF files to an `xarray.DataArray` object. This
    function is useful for combining multiple GeoTIFF files that `get_map`, `get_dem`,
    and `get_global_dem` produce into a single `xarray.DataArray` for further analysis.
    Note that `shapely` and `rioxarray` need to be installed to use this function.

## Important Notes

- Bounding box coordinates should be in decimal degrees (WGS84) format: (west, south,
    east, north)
- Default projection for requesting maps is EPSG:3857
- `get_map` is restricted to US coverage; use `get_global_dem` or `get_image_server` for
    data outside the US

## Installation

Choose your preferred installation method:

### Using `pip`

```console
pip install seamless-3dep
```

### Using `micromamba` (recommended)

```console
micromamba install -c conda-forge seamless-3dep
```

Alternatively, you can use `conda` or `mamba`.

## Quick Start Guide

We can retrieve topographic data using Seamless3DEP in just a few lines of code. Then,
we can visualize or even reproject the data using `rioxarray`.

### Retrieving a DEM

```python
from pathlib import Path
import seamless_3dep as s3dep
import rioxarray as rxr

# Define area of interest (west, south, east, north)
bbox = (-105.7006276, 39.8472777, -104.869054, 40.298293)
data_dir = Path("data")

# Download DEM
tiff_files = s3dep.get_dem(bbox, data_dir)

# Convert to xarray.DataArray
dem = s3dep.tiffs_to_da(tiff_files, bbox, crs=4326)
```

![DEM Example](https://raw.githubusercontent.com/hyriver/seamless-3dep/main/docs/examples/images/dem.png)

### Retrieving a Slope Map

```python
slope_files = s3dep.get_map("Slope Degrees", bbox, data_dir)
dem = s3dep.tiffs_to_da(slope_files, bbox)
```

![Slope Example](https://raw.githubusercontent.com/hyriver/seamless-3dep/main/docs/examples/images/slope_dynamic.png)

## Contributing

We welcome contributions! Please see the
[contributing](https://docs.hyriver.io/seamless-3dep/latest/CONTRIBUTING/) section for
guidelines and instructions.
