Crop Sentinel 2 Product/ Изображение с использованием Python
Я пытаюсь сделать мониторинг земли, используя Sentinel 2 и Python. Для этого я собираюсь применить некоторые показатели, такие как NDVI,NDWI, к продуктам Sentinel 2. На данный момент у меня есть скрипт, использующий утилиту sentinelsat python для поиска и получения всей информации о продуктах между двумя датами и двумя процентами облачного покрытия. Я также сделал еще один скрипт для загрузки этих продуктов и их хранения в MongoDB. Тогда я могу использовать показатели, которые я сказал на них.
Проблема заключается в том, что, как вы знаете, размер продуктов Sentinel 2 составляет 100 x100 км, в то время как размер области, которую я хочу, намного меньше. Я получаю продукты с использованием файла geojson, поэтому идея состоит в том, чтобы разработать еще один сценарий, который может обрезать только пиксели в каждой полосе, чтобы получить новые полосы только с этой небольшой областью для эффективной работы.
Я провел несколько дней, исследуя структуру продуктов и читая статьи о них, но не могу найти ничего, с чего начать. какие-либо предложения? заранее спасибо