Histogram matching

This tool alters the cumulative distribution function (CDF) of a raster image to match, as closely as possible, the CDF of a reference histogram. Histogram matching works by first calculating the histogram of the input image. This input histogram and reference histograms are each then converted into CDFs. Each grid cell value in the input image is then mapped to the corresponding value in the reference CDF that has an equivalent (or as close as possible) cumulative probability value. Histogram matching provides the most flexible means of performing image contrast adjustment.

The reference histogram must be specified to the tool in the form of a text file (.txt). This file must contain two columns (delimited by a tab, space, comma, colon, or semicolon) where the first column contains the x value (i.e. the values that will be assigned to the grid cells in the output image) and the second column contains the frequency or probability. Note that 1) the file must not contain a header, 2) each x value/frequency pair must be on a separate row, and 3) the frequency/probability must not be cumulative (i.e. the file must contain the histogram and not the CDF). The CDF will be computed for the reference histogram automatically by the tool. It is possible to create this type of histogram using the wide range of distribution tools available in most spreadsheet programs (e.g. LibreOffice's Calc program). You simply must remember to save the file as a text-only (ASCII) file.

Histogram matching is related to the Histogram Matching (Two Images tool, which can be used when a reference CDF can be derived from a reference image. Histogram Equalization and Gaussian Contrast Stretch are similarly related tools frequently used for image contrast adjustment, where the reference CDFs are uniform and Gaussian (normal) respectively.


See Also:


The following is an example of a Python script that uses this tool:

wd = pluginHost.getWorkingDirectory()
inputRaster = wd + "input.dep"
inputHisto = wd = "histo.txt"
outputFile = wd + "output.dep"
args = [inputFile, inputHisto, outputFile]
pluginHost.runPlugin("HistogramMatching", args, False)

This is a Groovy script also using this tool:

def wd = pluginHost.getWorkingDirectory()
def inputFile = wd + "input.dep"
def inputHisto = wd = "histo.txt"
def outputFile = wd + "output.dep"
String[] args = [inputFile, inputHisto, outputFile]
pluginHost.runPlugin("HistogramMatching", args, false)