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

posted @ 2016-09-12 16:10  Sonet  阅读(3550)  评论(0编辑  收藏  举报