# gdal_calc.py

Command line raster calculator with numpy syntax.

## Synopsis

    gdal_calc.py --calc=expression --outfile=out_filename [-A filename]
                 [--A_band=n] [-B...-Z filename] [other_options]

DESCRIPTION

Command line raster calculator with numpy syntax. Use any basic arithmetic supported by numpy arrays such as `+`, `-`, `*`, and `\` along with logical operators such as `>`. Note that all files must have the same dimensions, but no projection checking is performed.

`--help```

  Show this help message and exit

`-h```

  The same as `--help`.

`--calc``=expression`

  Calculation in gdalnumeric syntax using `+`, `-`, `/`, `*`, or any numpy array functions (i.e. `log10()`).

`-A``<filename>`

  Input gdal raster file, you can use any letter (A-Z).

`--A_band``=<n>`

  Number of raster band for file A (default 1).

`--outfile``=<filename>`

  Output file to generate or fill.

`--NoDataValue``=<value>`

  Output nodata value (default datatype specific value).

`--type``=<datatype>`

  Output datatype, must be one of [`Int32`, `Int16`, `Float64`, `UInt16`, `Byte`, `UInt32`, `Float32`].

`--format``=<gdal_format>`

  GDAL format for output file.

`--creation-option``=<option>`

  Passes a creation option to the output format driver. Multiple options may be listed. See format specific documentation for legal creation options for each format.

`--co``=<option>`

  The same as creation-option.

`--allBands``=[A-Z]`

  Process all bands of given raster (A-Z).

`--overwrite```

  Overwrite output file if it already exists.

`--debug```

  Print debugging information.

`--quiet```

  Suppress progress messages.

## Example

Add two files together:

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

Average of two layers:

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

Set values of zero and below to null:

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