#gdal #geotiff
#gdal #geotiff
Вопрос:
Я пытаюсь объединить несколько геотиффов в один большой geotiff с обзорами, однако окончательный объединенный geotiff показывает несколько горизонтальных артефактов по краям исходных объединенных геотиффов (смотрите здесь пример).
Я создаю объединенный файл, используя следующий код:
'''Produce Combined VRT'''
string ='gdalbuildvrt -srcnodata "0 0 0 0" -hidenodata -r bilinear %s -overwrite %s' %(tmp_vrt, GDal_merge_string)
os.system(string)
'''Convert VRT to Geotiff'''
string ='gdal_translate -b 1 -b 2 -b 3 -mask 4 --config GDAL_TIFF_INTERNAL_MASK YES -of GTiff %s %s' %(tmp_vrt,tmp_fname)
os.system(string)
У меня есть подозрение, что это может быть связано с использованием gdal_translate в vrt, поскольку ошибки возникают на краях исходных геотиффов, и в этом случае это может быть связано или похоже на проблему, описанную в этом сообщении.
Этот код использует VRTS для объединения геотиффов в целях ускорения, но, возможно, было бы лучше просто объединить их с gdalwarp?
Редактировать: я уменьшил количество флагов и опустил обзоры в приведенном выше коде, как предложено в комментарии ниже Бенджамином. Ошибка, похоже, возникает в приведенном выше коде. Я думаю, проблема может заключаться в процессе маскирования. Я предполагаю, что в какой-то момент в процессе укладки полос входные данные искажаются. нецелесообразно ли вообще использовать gdal_translate VRTS?
Комментарии:
1. Это помогло бы удалить некоторые ненужные теги, чтобы упростить устранение неполадок (например, просмотреть результирующее изображение без разбиения на плитки, сжатия, обзорных изображений и т.д.).
2. Спасибо за ваш комментарий @Benjamin , я отредактировал приведенный выше код, чтобы уменьшить его только до флагов, необходимых для создания изображения.