Сравнение значений NDVI внутри полигонов ⇐ Python

[Anonymous]

Я создал три изображения NDVI на основе изображений Landsat5. Теперь я хотел бы посмотреть, различаются ли значения NDVI (средние) внутри трех разных областей на каждом изображении. Если вы посмотрите на изображение по ссылке, вы увидите два многоугольника (красный и синий). Я хотел бы использовать шейп-файл с полигонами и извлечь средние значения NDVI (из растра) внутри чтения, синего цвета, а также средние значения за пределами двух полигонов. Изображение с многоугольниками

Я понятия не имею, какие инструменты использовать, чтобы получить значения NDVI изнутри и непосредственно внутри полигонов. Я думаю, мне понадобится третий многоугольник для «внешних» значений, однако из него придется исключить красные и синие значения.

В конце концов я хочу сделать это программно, используя Python, однако мне нужно знать, какие шаги предпринять и какие инструменты использовать.