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