"""
Copernicus DEM around Mount Rainier
===================================

The 30 m Copernicus digital surface model, loaded from the Planetary Computer
and shaded by elevation. ``source="earth-search"`` reads the same 2021 release
from the unsigned AWS Open Data bucket instead.
"""

import matplotlib.pyplot as plt

import easysnowdata as esd

aoi = (-121.94, 46.72, -121.54, 46.99)  # Mount Rainier, WA

dem = esd.terrain.dem.load(aoi, resolution=30)

fig, ax = plt.subplots(figsize=(7, 5))
dem.plot.imshow(ax=ax, cmap="terrain", cbar_kwargs={"label": "elevation [m]"})
ax.set_title("Copernicus DEM GLO-30")
ax.set_aspect("equal")
fig.tight_layout()