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