开源GIS数据处理工具GDAL：gdal_calc.py用法

gdal_calc.py --calc=expression --outfile=out_filename [-A filename]

             [--A_band=n] [-B...-Z filename] [other_options]

Options:

  -h, --help            显示此帮助消息并退出

  --calc=expression     使用+ - / *或任何的gdalnumeric语法进行计算

                         numpy数组函数（即log10（））

  -A filename           输入gdal光栅文件，可以使用任何字母（A-Z）

  --A_band=n            文件A的光栅带数（默认为1）

  --outfile=filename    输出文件生成或填充

  --NoDataValue=value   输出节点值（默认数据类型特定值）

  --type=datatype       输出数据类型，必须是['Int32'，'Int16'，'Float64'，

                        'UInt16'，'Byte'，'UInt32'，'Float32'）之一

  --format=gdal_format  GDAL格式的输出文件（默认'GTiff'）

  --creation-option=option, --co=option

                        将创建选项传递给输出格式驱动程序。

                         可以列出多个选项。 请参阅具体格式

                         每个法律创建选项的文档

                         格式。

  --allBands=[A-Z]      处理给定栅格（A-Z）的所有频带

  --overwrite           覆盖输出文件（如果已经存在）

  --debug               打印调试信息

  --quiet               抑制进度消息

gdal_calc.py -A input1.tif -B input2.tif --outfile=result.tif --calc="A+B"

gdal_calc.py -A input.tif -B input2.tif --outfile=result.tif --calc="(A+B)/2"

gdal_calc.py -A input.tif --outfile=result.tif --calc="A*(A>0)" --NoDataValue=0