DEM_filepath – path to DEM raster 
sample_points_filepath – path to sample points shapefile
In [ ]:
    
DEM_filepath = ""
    
In [ ]:
    
sample_points_filepath = ""
    
In [3]:
    
import matplotlib.pylab as plt
    
In [4]:
    
%matplotlib inline
    
In [1]:
    
import rasterio
import fiona
    
In [5]:
    
plt.plot([1,2,3,4])
plt.ylabel('some numbers')
plt.show()
    
    
In [17]:
    
with rasterio.drivers():
    with rasterio.open(DEM_filepath) as source_dem:
        array_dem = source_dem.read(1)
        source_dem.close()
    
In [20]:
    
plt.imshow(array_dem)
plt.ylabel("pixels")
    
    Out[20]:
    
In [6]:
    
with fiona.open(sample_points_filepath, 'r') as source_points:
    points = [f['geometry']['coordinates'] for f in source_points]
    
    #plt.figure()
    for f in source_points:
        x, y = f['geometry']['coordinates']
        plt.plot(x, y, 'ro')
    plt.show()
    
    source_points.close()