Я пытаюсь объединить несколько геотифов в один большой геотиф с обзорами, однако последний объединенный геотиф показывает ряд горизонтальных артефактов по краям исходных объединенных геотифов ( см. Здесь пример ).
Я создаю объединенный файл, используя следующий код:
'''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, так как ошибки происходят наребра оригинальных геотифов, и в этом случае это может быть связано или похоже на проблему, найденную в этом посте .
Этот код использует VRT для объединения геотифов в целях скорости, но, возможно, было бы лучше просто объединить их с gdalwarp?
Редактировать: Я сократил количество флагов и пропустил обзоры в приведенном выше коде, как это предлагается в комментарии ниже Бенджамином.Ошибка, кажется, выдается в коде выше.Я думаю, что проблема может заключаться в процессе маскировки.Я думаю, что в какой-то момент в процессе укладки полос входы искажаются.Как правило, нецелесообразно gdal_translate VRT?