In [8]:
Copied!
import geopandas as gpd
import rioxarray as rxr
import matplotlib.pyplot as plt
from streamkit import get_huc_data
import geopandas as gpd
import rioxarray as rxr
import matplotlib.pyplot as plt
from streamkit import get_huc_data
In [20]:
Copied!
dem_30, flowlines_mr = get_huc_data("180101020502", nhd_layer="medium",
crs="EPSG:3310", dem_resolution=30)
dem_30, flowlines_mr = get_huc_data("180101020502", nhd_layer="medium",
crs="EPSG:3310", dem_resolution=30)
In [21]:
Copied!
dem_10, flowlines_hr = get_huc_data("180101020502", nhd_layer="high",
crs="EPSG:3310", dem_resolution=10)
dem_10, flowlines_hr = get_huc_data("180101020502", nhd_layer="high",
crs="EPSG:3310", dem_resolution=10)
In [38]:
Copied!
fig, axes = plt.subplots(1,2, figsize=(12,8))
dem_30.plot(ax=axes[0], add_colorbar=False)
flowlines_mr.plot(ax=axes[0], column='OBJECTID', cmap='gist_ncar')
dem_10.plot(ax=axes[1], add_colorbar=False)
flowlines_hr.plot(ax=axes[1], column='OBJECTID', cmap='gist_ncar')
fig, axes = plt.subplots(1,2, figsize=(12,8))
dem_30.plot(ax=axes[0], add_colorbar=False)
flowlines_mr.plot(ax=axes[0], column='OBJECTID', cmap='gist_ncar')
dem_10.plot(ax=axes[1], add_colorbar=False)
flowlines_hr.plot(ax=axes[1], column='OBJECTID', cmap='gist_ncar')
Out[38]:
<Axes: title={'center': 'band = 1, spatial_ref = 0'}, xlabel='x coordinate of projection\n[metre]', ylabel='y coordinate of projection\n[metre]'>
