import topotoolbox
import matplotlib.pyplot as plt
from rasterio import CRS
dem = topotoolbox.load_opentopography(south=50, north=50.1, west=14.35,
    east=14.6, dem_type="SRTMGL3")
dem = dem.reproject(CRS.from_epsg(32633), resolution=90)
_= dem.plot(cmap="terrain")
plt.show()