Слияние спутниковых снимков и сохранение координат
Спасибо, что заглянули сюда. В настоящее время я работаю над проектом, и я еще не так силен с Python. Поэтому я надеялся на конструктивную обратную связь по этому вопросу.
У меня есть набор данных, содержащий образцы керна, все они хранятся с идентификатором образца, широтой, долготой, содержимым и другими данными, не относящимися к этому вопросу.
Теперь я импортировал этот набор данных и нарезал его так, как я хочу. Для изображений я использую модуль растерио, чтобы открыть 2 спутниковых изображения, которые покрывают регион. Я использую модуль utm для преобразования туда и обратно между значениями latlong->UTM->Pixel (что также, кажется, вызывает некоторые странные координаты в некоторых точках). Достаточно досадно, что два изображения Sentinel-2 вырезаны прямо по центру карты. Поскольку я делаю ограничивающие рамки поверх тех мест, где берутся образцы, это проблема, так как мне нужно извлечь 10x10 пикселей в этой области. Это приводит к тому, что многие образцы не получают надлежащего вырезания.
Поэтому я подумал, почему бы не объединить два изображения в один большой прямоугольный бит. Но мне все еще нужно сохранить метаданные с координатами UTM.
Как бы вы предложили мне продолжить. Можно ли сделать это простым способом? Есть ли какой-то другой аспект этого вопроса, который я упустил?
Спасибо за ваше время.
1 ответ
Я не уверен, что полностью понимаю вопрос, но если вы просто пытаетесь объединить 2 изображения, вы смотрели на инструмент командной строки gdal_merge.py?
Очень простой пример:
gdal_merge.py -o merged_image.tif image1.tif image2.tif