Слияние спутниковых снимков и сохранение координат - PullRequest
0 голосов
/ 04 сентября 2018

Спасибо, что заглянули сюда. В настоящее время я работаю над проектом, и я еще не так силен с Python. Поэтому я надеялся на конструктивную обратную связь по этому вопросу.

У меня есть набор данных, содержащий выборки керна, все они хранятся с идентификатором выборки, широтой, долготой, содержанием и другими данными, не относящимися к этому вопросу.

Теперь я импортировал этот набор данных и нарезал его так, как я хочу. Для изображений я использую модуль растерио, чтобы открыть 2 спутниковых изображения, которые покрывают регион. Я использую модуль utm для преобразования туда и обратно между значениями latlong-> UTM-> Pixel (что также, кажется, вызывает некоторые странные координаты в некоторых точках). Достаточно досадно, что два изображения Sentinel-2 вырезаны прямо по центру карты. Так как я делаю ограничивающие рамки поверх тех мест, где берутся образцы, это проблема, так как мне нужно извлечь 10x10 пикселей в этой области. Это приводит к тому, что многие образцы не получают должного среза.

Поэтому я подумал, почему бы не объединить два изображения в один большой прямоугольный бит. Но мне все еще нужно сохранить метаданные с координатами UTM.

Как бы вы предложили мне продолжить. Можно ли сделать это простым способом? Есть ли какой-то другой аспект в этом аспекте, который я упустил?

Спасибо, что уделили время.

enter image description here enter image description here

1 Ответ

0 голосов
/ 05 сентября 2018

Я не уверен, что полностью понимаю вопрос, но если вы просто пытаетесь объединить 2 изображения, вы смотрели на инструмент командной строки gdal_merge.py ?

Очень простой пример:

gdal_merge.py -o merged_image.tif image1.tif image2.tif

...