Я использую некоторые спутниковые изображения Sentinel-2 в python.Теперь у меня нет проблем с использованием более новых (в прошлом 2016 году).Но мне нужно использовать некоторые из 2016 года. Они не обрабатываются таким же образом Европейским космическим агентством!
Обычно, когда вы загружаете плитку, вы обычно получаете файл .jp2 для каждой из полос спутников.Но в более новой версии они предварительно обрабатывают RGB-версию для вас вместе с обычными группами.Эта версия работает отлично для меня в Python.Однако, чтобы создать RGB-версию старых изображений, мне нужно объединить три полосы (4,3,2 - R, G, B) в один файл.Сначала Gdal_merge справляется с этим довольно хорошо.Когда я открываю изображение, оно выглядит великолепно!Но, прочитав его на python, я сразу заметил, что что-то не так.Изображение выглядит как чисто белое изображение с голубыми полосами на нем.Теперь я обратился к gdalinfo как с новым рабочим примером, так и со старой версией, которую я сам сшил, и это вывод.
Как вы можете заметитьНа первый взгляд размеры кажутся нормальными.Но полосы не имеют ни правильного типа, ни правильного цвета.Поэтому я явно делаю что-то не так, когда объединяю файлы.
Вот эта команда, которую я использую, чтобы объединить 3 полосы в 1 файл .jp2.
gdal_merge.py -o outfile.jp2 -separate B04.jp2 BO3.jp2 BO2.jp2
Теперь, как уже упоминалось.Это создает файл, и файл выглядит красиво, когда я открываю его в QGIS.Но это бесполезно для меня в python.
Вот скриншот импорта из python.
img - это предварительно обработанное изображение ESA.
img1 - мой убогий импорт gdal_merge.
А вот изображение сбоя: p 
Теперь мне кажется, что мне не хватаетнекоторое базовое понимание такого рода манипуляций с изображениями.Итак, благодаря объединенной мудрости Stackoverflow - что я могу сделать, чтобы правильно соединить мои группы в сексуальный RGB, который может быть правильно прочитан модулем растерио.
Заранее спасибо :)