Размещение векторного слоя поверх растрового слоя в python
Я пытаюсь классифицировать кадрирование по спутниковым снимкам (изображения Sentinel-2 в формате jp2) из шейп-файлов, содержащих несколько полигонов с классификаторами. Я могу сделать это в Qgis, и это выглядит следующим образом: изображение вектора поверх растра.
Я хочу сделать это на Python. Я использую растерио для открытия файлов jp2 и pyShp для чтения шейп-файла. Теперь я не уверен, как добавить эти 2 изображения поверх друг друга, чтобы определить значения пикселей определенного многоугольника (фигуры), я что-то упускаю из виду совершенно очевидно? Я полагаю, что можно было бы как-то объединить данные из двух файлов, но я не смог найти решение для этого.
import shapefile
import rasterio
satImage = rasterio.open("my/path.jp2")
shp = shapefile.Reader("my/shape.shp")