zoukankan      html  css  js  c++  java
  • 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%in;%PATH%

    SET PYTHONHOME=%OSGEO4W_ROOT%appsPython27

    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

  • 相关阅读:
    dxCalloutPopup 简单使用教程
    Delphi INI文件保存与读取
    AlertWindowManager 弹出提示窗口使用帮助(下)
    AlertWindowManager 弹出提示窗口使用帮助(上)
    可输入弹出窗口-[POPUP_GET_VALUES_USER_HELP]
    [BAPI]采购申请PR审批-BAPI_REQUISITION_RELEASE_GEN
    如何取域值 (当一些业务需要的值只有数字或者字母时 ,汉字描述在域里面)
    采购订单、采购申请审批策略相关表
    [BAPI]如何修改工单状态-BAPI_ALM_ORDER_MAINTAIN
    [函数]读取采购订单、采购申请更改历史-ME_CHANGEDOC_READ2
  • 原文地址:https://www.cnblogs.com/Sonet-life/p/5865168.html
Copyright © 2011-2022 走看看