Plot a dataset overlain on imagery

PlotOnImagery.py

#
# First, Convert the ground track to imagery coordinate system
#

from pyproj import Transformer
transformer = Transformer.from_crs("EPSG:4326", "EPSG:32617")

#
# Make a plot for the imagery
#

import rasterio
from rasterio.plot import show

fig,ax = plt.subplots(figsize=(10,10))
sat_data = rasterio.open('SomeTifFile.tif')
show(sat_data,cmap="gray")

#
# Convert our dataset from lat-lon to projected coordinates, then plot
#
[x,y] = transformer.transform( lat , lon )
ax.scatter(x,y)

plt.show()