gdal_merge.py

1 gdal_merge.py: 合并（Merge）/镶嵌（Mosaic）工具。要求图像必须是相同坐标系统、具有相同的波段数；可以不同分辨率，可以有重叠区域（后加入图像覆盖先加入的图像）。

注意：只能使用nearest重采样方法。

gdal_merge.py [-o out_filename] [-of out_format] [-co NAME=VALUE]*

              [-ps pixelsize_x pixelsize_y] [-tap] [-separate] [-v] [-pct]

              [-ul_lr ulx uly lrx lry] [-n nodata_value] [-init "value [value...]"]

              [-ot datatype] [-createonly] input_files

参考：http://www.gdal.org/gdal_merge.html

http://osdir.com/ml/gdal-development-gis-osgeo/2006-03/msg00087.html

使用时需要python、gdal(python库)支持，设置方法：

set OSGEO4W_ROOT=D:\Program Files (x86)\QGIS Chugiak

set PATH=%OSGEO4W_ROOT%\bin;%PATH%

SET PYTHONHOME=%OSGEO4W_ROOT%\apps\Python27

2 应用

2.1 图像合并（Merge）/镶嵌（Mosaic）

目的：将多个图像合并为一幅图像。

原理：读入源图像（多个），复制（如果分辨率不一致，需要进行重采样）源图像到目标图像的相应位置。

方法：gdal_merge.py

-n 源图像中无值时赋此值。

-a_nodata 同-n(>1.9)。

-co 参数。格式：key=value。COMPRESS=JPEG：使用JPEG压缩方式。JPEG_QUALITY=75。JPEG压缩率。

-of 输出格式。默认为GTiff(GeoTiff）。

-o 输出图像。默认（out.tif）。

示例：

gdal_merge.bat -n 0 -a_nodata 0 -coCOMPRESS=JPEG -co JPEG_QUALITY=75 -of GTiff -o D:/qgis-ppt/mosaic/merge4.tifD:/qgis-ppt/mosaic/J1.tif D:/qgis-ppt/mosaic/J2.tif D:/qgis-ppt/mosaic/J3.tifD:/qgis-ppt/mosaic/J4.tif

2.2 保持输入图像波段不变，

目的：将多个图像合并为一幅图像，并保持在目标图像中保留原始图像的所有波段。

原理：读入源图像（多个），复制（如果分辨率不一致，需要进行重采样）源图像到目标图像的相应位置。

方法：gdal_merge.py

-separate 将为每个输入图像的通道保留在目标图像的单独通道中。

示例：

gdal_merge.bat -separate -of GTiff -oD:/qgis-ppt/mosaic/merge4-layerstack.tif D:/qgis-ppt/mosaic/J1.tifD:/qgis-ppt/mosaic/J2.tif D:/qgis-ppt/mosaic/J3.tif D:/qgis-ppt/mosaic/J4.tif

可以看到结果中，源图像的3个通道都被依次加入到目标图像中。

2.3 保存PCT信息

目的：将多个图像合并为一幅图像，并保存PCT颜色表。

原理：读入源图像（多个），复制（如果分辨率不一致，需要进行重采样）源图像到目标图像的相应位置，将PCT信息保存在Tiff等输出文件的文件头。

方法：gdal_merge.py

-pct 读入第一幅图像的PCT信息，并写入输出文件（假设所有输入文件使用相同的PCT）。

示例：

gdal_merge.bat -pct -of GTiff -oD:/qgis-ppt/mosaic/merge4-pct.tif D:/qgis-ppt/mosaic/J1.tifD:/qgis-ppt/mosaic/J2.tif D:/qgis-ppt/mosaic/J3.tif D:/qgis-ppt/mosaic/J4.tif

打开图像的16进制格式，可以看到，PCT文件比无PCT的文件，多了PCT文件头信息，其它部分都一样。

2.4 输出指定部分数据（将结果进行裁剪）

目的：将多个图像合并为一幅图像，输出指定范围数据。

原理：读入源图像（多个），复制（如果分辨率不一致，需要进行重采样）源图像到目标图像的相应位置，将结果数据进行范围裁剪。

方法：gdal_merge.py

-ul_lr upperleft_lowerright,左上角到右下角范围。默认使用全部范围。

示例：

gdal_merge.bat -ul_lr 114.01 39.99 114.03 39.97 -ofGTiff -o D:/qgis-ppt/mosaic/merge4-ul_lr.tif D:/qgis-ppt/mosaic/J1.tifD:/qgis-ppt/mosaic/merge4.tif

2.5 输出指定分辨率数据

目的：将多个图像合并为一幅图像，并输出指定分辨率数据。

原理：读入源图像（多个），进行重采样，将结果得到到目标图像的相应位置。

方法：gdal_merge.py

-ps xpixelsize ypixelsize,xy两个方向的输出分辨率。默认使用第一幅图像的分辨率。注意Y的方向。

示例：

gdal_merge.bat -of GTiff -ps 2.14577e-05-2.14577e-05 -o D:/qgis-ppt/mosaic/merge4-res.tif D:/qgis-ppt/mosaic/J1.tifD:/qgis-ppt/mosaic/J2.tif D:/qgis-ppt/mosaic/J3.tif D:/qgis-ppt/mosaic/J4.tif