Я хотел бы использовать Zonal Statistics для извлечения значений RGB из tif в файл shp. Это можно легко сделать с помощью плагина в QGIS, но у меня есть почти 40 отдельных файлов shp и связанный tiff.
В идеале я хотел бы немного простого кода, который бы добавил файл Shp с новым столбцом для средних значений R, G и B.
Любые идеи, где I go из кода ниже, который получает R для каждого шейп-файла. Я не могу понять, как вывести его в качестве дополнительного столбца.
statsr=zonal_stats('vector.shp','raster.tif', band= 1, nodata=0, stats=['mean'],geojon_out=True)