gis - gdalwarp 太慢(与 gdal_merge 相比)

Question

我有 70 多个 TIFF 格式的光栅图像，我正在尝试合并它们。

原件可以在这里找到:
http://www.faa.gov/air_traffic/flight_info/aeronav/digital_products/vfr/

在预处理(pct2rgb、gdalwarp 个人图表、gdal_translate 以切割衣领)之后，我尝试通过 gdalwarp 运行它们以使用如下命令对它们进行镶嵌:

gdalwarp --config GDAL_CACHEMAX 3000 -overwrite -wm 3000 -r bilinear -srcnodata 0 -dstnodata 0 -wo "NUM_THREADS=3" /data/aeronav/sec/c/Albuquerque_c.tif .....70 other file names ...master.tif

处理 12 小时后:

创建 321521P x 125647L 的输出文件。
处理输入文件/data/aeronav/sec/c/Albuquerque_c.tif。
0...10...20...30...40...

这意味着 gdalwarp 永远不会结束。

相比之下。像这样的 gdal_merge 命令:

gdal_merge.py -n 0 -a_nodata 0 -o /data/aeronav/sec/master.tif /data/aeronav/sec/c/Albuquerque_c.tif ......70 plus files.....

几个小时内完成。

由于“平均”采样，gdal_merge 的问题是质量较差的输出。如果可能，我想至少使用“双线性”和“三次”采样，为此需要 gdalwarp。

为什么两者的性能会有如此大的差异？为什么 gdalwarp 不想完成？是否有任何其他命令行选项可以在 gadalwarp 中加快速度，或者有没有办法将采样选项添加到 gdal_merge ？

Answer 1

好像gdalwarp不是合并这些 GeoTiffs 的理想命令(因为我对再次变形不感兴趣)。相反我用

gdalbuildvrt /data/aeronav/sec/master.virt .... 70+ files in order

构建虚拟马赛克。然后我用了 gdal_translate将 virt 文件转换为 GeoTiff:

gdal_translate -of GTiff  /data/aeronav/sec/master.virt /data/aeronav/sec/master.tif

就是这样 - 这花了不到一个小时(甚至比 gdal_merge 还要快，并且保留了原始文件的质量)。