使用gdalbuildvrt和gdaltranslate在geotiff边缘之间的水平伪像

Question

我正在尝试将多个Geotiff合并为一个具有概述的大型Geotiff，但是最终的合并Geotiff在原始合并Geotiff（see here for an example）的边缘周围显示了许多水平伪像。

我使用以下代码创建合并文件：

'''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)

我有一个预感，这可能与在vrt上使用gdal_translate有关，因为错误发生在原始geotiff的边缘，在这种情况下，它可能与the issue found in this post相关或相似。

此代码使用VRT组合geotiff以达到速度目的，但是将它们与gdalwarp合并可能会更好？

编辑：正如本杰明在下面的注释中所建议的，我减少了标志的数量，并省略了上面代码中的概述。该错误似乎是在上面的代码中产生的。我认为问题可能在于掩盖过程。我想在乐队堆叠过程中的某些时候，输入会失真。 gdal_translate VRT通常不建议吗？