[GIS] combine raster with another raster mask file – gdal

cclipgdalmaskingraster

I have 2 rasters in GeoTiff format. One has all of the land and water and the other is a mask file that shows either 0 for water or 1 for land. I basically want my output to just have the land values and all the areas that were water in the first image to just be 0,0,0. Is there a way to do this with any of the gdal command line tools?

Basically I want to do what gdal_rasterize does with a shp file mask but with a raster instead.

Best Answer

You may want to check gdal_translate (see http://www.gdal.org/gdal_translate.html) in which you can define a mask band. But first, you'll have to stack your two raster in a multiband image.

Related Solutions

QGIS Raster Masking – Clipping Raster Layer with Raster Mask Layer in QGIS

The simple way to do this in QGIS is to use the Raster Calculator (Raster->Raster Calculator). You have a couple of options. The easiest to explain/understand is to make a unitary raster from your mask (all data set to either 1 or NoData) and then multiply your clip layer by the unitary mask layer.

To ensure the extents match the mask layer, in the raster calculator window select the mask layer in the 'Raster bands' list on the left and then click the 'Current layer extent' button on the right.

You can create a unitary mask on the fly by using a conditional statement (see the link) something like this:

(maskLayer@1 >= 0) * clipLayer@1

This statement basically says: treat everything in my mask layer that is not NoData as being equal to 1 (NoData stays as NoData). Just be sure to remember to set the extent (see above)!

Raster Calculator and paletted data:

Any output from the raster calculator will be just values and not carry over any information contained within a colour pallet. You have a couple of options to 'get your palette back':

Go to the style tab of the original layer's properties and (making sure it is paletted), click on Save Style button at the bottom. You can then apply this same style to your new layer and if you have only clipped it as per the above instructions, it will appear the same as before. This is easy but your raster is not paletted in a persistent way outside of the QGIS project.
To make the presentation persistent, you can right-click the layer and 'Save As'. Check the 'Rendered Image' radio button at the top of the save-dialog box. This will make a new raster with the exact colour scheme as your palatted layer BUT it will no longer be a single band raster but a four band RGBA raster.
If you MUST have an actual stand-alone single-band raster with a persistent palette (and not RGBA), then I don't know of a way of doing this in QGIS but you could take the output of step 2 above and open it in GIMP (or Photoshop) and change the mode from RGB to Indexed Colour. But you will need to create a world file as saving the image in GIMP/Photoshop will destroy the georeferencing.

[GIS] Mask image based on no-data pixels in another image with different projection – GDAL or Python

You could create a mask from the image, polygonize this and then use the same gdalwarp command you normally use. If you wanted to do this in Python you could use the bandmath function in RSGISLib (rsgislib.org). For example:

# Import RSGISLib
import rsgislib
from rsgislib import imagecalc
from rsgislib.imagecalc import BandDefn
import subprocess

# 1) Create mask
inImage = 'input_image.tiff'
outMask = 'image_mask.tiff'

gdalFormat = 'GTiff'
dataType = rsgislib.TYPE_8UINT
expression = 'b1 > 0?1:0' # Only include values greater than zero.
bandDefns = []
bandDefns.append(BandDefn('b1', inImage, 1)) # Use band one (variable b1) of in image
imagecalc.bandMath(outMask, expression, gdalFormat, dataType, bandDefns)

# 2) Generate polygon from mask
subprocess.call("gdal_polygonize.py -mask image_mask.tiff image_mask.tiff -b 1 -f 'ESRI Shapefile' image_Footprint.shp",shell=True)

There are other ways you could create an image mask if you didn't want to use RSGISLib, you could use RIOS (https://bitbucket.org/chchrsc/rios/wiki/applier/simpleexample) to accomplish the same thing.

Best Answer

Related Solutions

QGIS Raster Masking – Clipping Raster Layer with Raster Mask Layer in QGIS

[GIS] Mask image based on no-data pixels in another image with different projection – GDAL or Python

Related Question