import topotoolbox
import matplotlib.pyplot as plt
from rasterio import CRS
dem = topotoolbox.load_dem('kunashiri')
dem = dem.reproject(CRS.from_epsg(32655), resolution = 30)
_= dem.plot(cmap="terrain")
plt.show()