Example: Equalization

Functions > Image Processing > Image Manipulation > Example: Equalization

equalize

Use the equalize function to perform a histogram equalization. Equalization can be used as an advanced form of level mapping to enhance features of a picture, or to use a particular display device more efficiently.

For information on using this example, refer to About Image Processing Examples.

1. Read in an image and plot the following histogram.

<region id="ID0EL1VX" top="403.20000000000005" width="600" left="19.200000000000003" height="312" xmlns="http://schemas.mathsoft.com/worksheet50">
  <plot background-type="white">
    <xyPlot>
      <title class="- topic/title " wwtype:type="Paragraph" xmlns:wwtype="urn:WebWorks-Type-Schema" />
      <legend />
      <traces>
        <trace resultRef="76">
          <traceStyle color="#FFFF0000" symbol="none" line-weight="1" line-style="Solid">lines</traceStyle>
        </trace>
      </traces>
      <graph-size width="500" height="254.4" />
      <axes>
        <xAxis rank="1" legend-position="PlotBoundaryBottom" start="0" end="275">
          <axisLine position="ticknumberlock" positionticmark="0" legendWidth="107.796666666667" />
          <axisGrid>
            <gridFrequency>12</gridFrequency>
            <gridLabels display="true" />
            <gridLines />
            <tickMarks display="true" />
          </axisGrid>
          <axisLabel />
          <markers />
          <numberFormat>
            <decimal precision="3" show-trailing-zeros="false" radix="dec" imaginary-value="i" />
          </numberFormat>
          <plotEquations>
            <plotEquation>
              <math resultRef="77">
                <id xml:space="preserve" xmlns="http://schemas.mathsoft.com/math50">k_range</id>
              </math>
              <math resultRef="78">
                <placeholder xmlns="http://schemas.mathsoft.com/math50" />
              </math>
            </plotEquation>
          </plotEquations>
          <xyDomain scale-type="linear" auto-scale="true">
            <startValue>
              <placeholder xmlns="http://schemas.mathsoft.com/math50" />
            </startValue>
            <endValue>
              <placeholder xmlns="http://schemas.mathsoft.com/math50" />
            </endValue>
          </xyDomain>
        </xAxis>
        <yAxis rank="1" legend-position="PlotBoundaryLeft" start="0" end="8800">
          <axisLine position="ticknumberlock" positionticmark="0" legendWidth="76.77" />
          <axisGrid>
            <gridFrequency>12</gridFrequency>
            <gridLabels display="true" />
            <gridLines />
            <tickMarks display="true" />
          </axisGrid>
          <axisLabel />
          <markers />
          <numberFormat>
            <decimal precision="3" show-trailing-zeros="false" radix="dec" imaginary-value="i" />
          </numberFormat>
          <plotEquations>
            <plotEquation>
              <math resultRef="79">
                <id xml:space="preserve" xmlns="http://schemas.mathsoft.com/math50">HK</id>
              </math>
              <math resultRef="80">
                <placeholder xmlns="http://schemas.mathsoft.com/math50" />
              </math>
            </plotEquation>
          </plotEquations>
          <xyDomain scale-type="linear" auto-scale="true">
            <startValue>
              <placeholder xmlns="http://schemas.mathsoft.com/math50" />
            </startValue>
            <endValue>
              <placeholder xmlns="http://schemas.mathsoft.com/math50" />
            </endValue>
          </xyDomain>
        </yAxis>
      </axes>
    </xyPlot>
  </plot>
</region>

This image is very dark. You can spread the intensities over a larger portion of the available spectrum.

2. Define the cumulative histogram by the difference equation for CK.

<region id="ID0E24VX" top="940.80000000000018" width="599" left="19.200000000000003" height="312" xmlns="http://schemas.mathsoft.com/worksheet50">
  <plot background-type="white">
    <xyPlot>
      <title class="- topic/title " wwtype:type="Paragraph" xmlns:wwtype="urn:WebWorks-Type-Schema" />
      <legend />
      <traces>
        <trace resultRef="93">
          <traceStyle color="#FFFF0000" symbol="none" line-weight="1" line-style="Solid">lines</traceStyle>
        </trace>
      </traces>
      <graph-size width="499" height="254.4" />
      <axes>
        <xAxis rank="1" legend-position="PlotBoundaryBottom" start="0" end="275">
          <axisLine position="ticknumberlock" positionticmark="0" legendWidth="107.796666666667" />
          <axisGrid>
            <gridFrequency>12</gridFrequency>
            <gridLabels display="true" />
            <gridLines />
            <tickMarks display="true" />
          </axisGrid>
          <axisLabel />
          <markers />
          <numberFormat>
            <decimal precision="3" show-trailing-zeros="false" radix="dec" imaginary-value="i" />
          </numberFormat>
          <plotEquations>
            <plotEquation>
              <math resultRef="94">
                <id xml:space="preserve" xmlns="http://schemas.mathsoft.com/math50">k_range</id>
              </math>
              <math resultRef="95">
                <placeholder xmlns="http://schemas.mathsoft.com/math50" />
              </math>
            </plotEquation>
          </plotEquations>
          <xyDomain scale-type="linear" auto-scale="true">
            <startValue>
              <placeholder xmlns="http://schemas.mathsoft.com/math50" />
            </startValue>
            <endValue>
              <placeholder xmlns="http://schemas.mathsoft.com/math50" />
            </endValue>
          </xyDomain>
        </xAxis>
        <yAxis rank="1" legend-position="PlotBoundaryLeft" start="0" end="71500">
          <axisLine position="ticknumberlock" positionticmark="0" legendWidth="75.54" />
          <axisGrid>
            <gridFrequency>12</gridFrequency>
            <gridLabels display="true" />
            <gridLines />
            <tickMarks display="true" />
          </axisGrid>
          <axisLabel />
          <markers />
          <numberFormat>
            <decimal precision="3" show-trailing-zeros="false" radix="dec" imaginary-value="i" />
          </numberFormat>
          <plotEquations>
            <plotEquation>
              <math resultRef="96">
                <id xml:space="preserve" xmlns="http://schemas.mathsoft.com/math50">CK</id>
              </math>
              <math resultRef="97">
                <placeholder xmlns="http://schemas.mathsoft.com/math50" />
              </math>
            </plotEquation>
          </plotEquations>
          <xyDomain scale-type="linear" auto-scale="true">
            <startValue>
              <placeholder xmlns="http://schemas.mathsoft.com/math50" />
            </startValue>
            <endValue>
              <placeholder xmlns="http://schemas.mathsoft.com/math50" />
            </endValue>
          </xyDomain>
        </yAxis>
      </axes>
    </xyPlot>
  </plot>
</region>

3. Linearize the cumulative histogram.

4. Plot the new image histogram.

<region id="ID0EFCWX" top="1747.2000000000003" width="596" left="19.200000000000003" height="312" xmlns="http://schemas.mathsoft.com/worksheet50">
  <plot background-type="white">
    <xyPlot>
      <title class="- topic/title " wwtype:type="Paragraph" xmlns:wwtype="urn:WebWorks-Type-Schema" />
      <legend />
      <traces>
        <trace resultRef="110">
          <traceStyle color="#FFFF0000" symbol="none" line-weight="1" line-style="Solid">lines</traceStyle>
        </trace>
      </traces>
      <graph-size width="496" height="254.4" />
      <axes>
        <xAxis rank="1" legend-position="PlotBoundaryBottom" start="0" end="275">
          <axisLine position="ticknumberlock" positionticmark="0" legendWidth="107.796666666667" />
          <axisGrid>
            <gridFrequency>12</gridFrequency>
            <gridLabels display="true" />
            <gridLines />
            <tickMarks display="true" />
          </axisGrid>
          <axisLabel />
          <markers />
          <numberFormat>
            <decimal precision="3" show-trailing-zeros="false" radix="dec" imaginary-value="i" />
          </numberFormat>
          <plotEquations>
            <plotEquation>
              <math resultRef="111">
                <id xml:space="preserve" xmlns="http://schemas.mathsoft.com/math50">k_range</id>
              </math>
              <math resultRef="112">
                <placeholder xmlns="http://schemas.mathsoft.com/math50" />
              </math>
            </plotEquation>
          </plotEquations>
          <xyDomain scale-type="linear" auto-scale="true">
            <startValue>
              <placeholder xmlns="http://schemas.mathsoft.com/math50" />
            </startValue>
            <endValue>
              <placeholder xmlns="http://schemas.mathsoft.com/math50" />
            </endValue>
          </xyDomain>
        </xAxis>
        <yAxis rank="1" legend-position="PlotBoundaryLeft" start="0" end="8800">
          <axisLine position="ticknumberlock" positionticmark="0" legendWidth="71.93" />
          <axisGrid>
            <gridFrequency>12</gridFrequency>
            <gridLabels display="true" />
            <gridLines />
            <tickMarks display="true" />
          </axisGrid>
          <axisLabel />
          <markers />
          <numberFormat>
            <decimal precision="3" show-trailing-zeros="false" radix="dec" imaginary-value="i" />
          </numberFormat>
          <plotEquations>
            <plotEquation>
              <math resultRef="113">
                <id xml:space="preserve" xmlns="http://schemas.mathsoft.com/math50">JK</id>
              </math>
              <math resultRef="114">
                <placeholder xmlns="http://schemas.mathsoft.com/math50" />
              </math>
            </plotEquation>
          </plotEquations>
          <xyDomain scale-type="linear" auto-scale="true">
            <startValue>
              <placeholder xmlns="http://schemas.mathsoft.com/math50" />
            </startValue>
            <endValue>
              <placeholder xmlns="http://schemas.mathsoft.com/math50" />
            </endValue>
          </xyDomain>
        </yAxis>
      </axes>
    </xyPlot>
  </plot>
</region>

5. Plot the cumulative histogram of the equalized image.

<region id="ID0EREWX" top="2208.0000000000005" width="597" left="19.200000000000003" height="312" xmlns="http://schemas.mathsoft.com/worksheet50">
  <plot background-type="white">
    <xyPlot>
      <title class="- topic/title " wwtype:type="Paragraph" xmlns:wwtype="urn:WebWorks-Type-Schema" />
      <legend />
      <traces>
        <trace resultRef="123">
          <traceStyle color="#FFFF0000" symbol="none" line-weight="1" line-style="Solid">lines</traceStyle>
        </trace>
      </traces>
      <graph-size width="497" height="254.4" />
      <axes>
        <xAxis rank="1" legend-position="PlotBoundaryBottom" start="0" end="275">
          <axisLine position="ticknumberlock" positionticmark="0" legendWidth="107.796666666667" />
          <axisGrid>
            <gridFrequency>12</gridFrequency>
            <gridLabels display="true" />
            <gridLines />
            <tickMarks display="true" />
          </axisGrid>
          <axisLabel />
          <markers />
          <numberFormat>
            <decimal precision="3" show-trailing-zeros="false" radix="dec" imaginary-value="i" />
          </numberFormat>
          <plotEquations>
            <plotEquation>
              <math resultRef="124">
                <id xml:space="preserve" xmlns="http://schemas.mathsoft.com/math50">k_range</id>
              </math>
              <math resultRef="125">
                <placeholder xmlns="http://schemas.mathsoft.com/math50" />
              </math>
            </plotEquation>
          </plotEquations>
          <xyDomain scale-type="linear" auto-scale="true">
            <startValue>
              <placeholder xmlns="http://schemas.mathsoft.com/math50" />
            </startValue>
            <endValue>
              <placeholder xmlns="http://schemas.mathsoft.com/math50" />
            </endValue>
          </xyDomain>
        </xAxis>
        <yAxis rank="1" legend-position="PlotBoundaryLeft" start="0" end="71500">
          <axisLine position="ticknumberlock" positionticmark="0" legendWidth="75.8766666666667" />
          <axisGrid>
            <gridFrequency>12</gridFrequency>
            <gridLabels display="true" />
            <gridLines />
            <tickMarks display="true" />
          </axisGrid>
          <axisLabel />
          <markers />
          <numberFormat>
            <decimal precision="3" show-trailing-zeros="false" radix="dec" imaginary-value="i" />
          </numberFormat>
          <plotEquations>
            <plotEquation>
              <math resultRef="126">
                <id xml:space="preserve" xmlns="http://schemas.mathsoft.com/math50">DK</id>
              </math>
              <math resultRef="127">
                <placeholder xmlns="http://schemas.mathsoft.com/math50" />
              </math>
            </plotEquation>
          </plotEquations>
          <xyDomain scale-type="linear" auto-scale="true">
            <startValue>
              <placeholder xmlns="http://schemas.mathsoft.com/math50" />
            </startValue>
            <endValue>
              <placeholder xmlns="http://schemas.mathsoft.com/math50" />
            </endValue>
          </xyDomain>
        </yAxis>
      </axes>
    </xyPlot>
  </plot>
</region>

6. Compare the original and equalized images.

(arcs.bmp)

(arcs_eq.bmp)

The equalized image contains details that are not visible in the original image.

The above technique can be used in the following situations:

• The image contains mainly highlights and insufficient shadows.

• The full dynamic range of the image is not used. For example, the image contains large gaps where no intensities are used.

Some highlights or shadows can become over exaggerated in equalized photographic images due to the sensitivity of the human eye and PC monitor. You can perform level mapping and attempt to correct the problem.

Scaling and clipping are other useful techniques for adjusting the dynamic range of an image. You can use these techniques only for images with intensity spectrum that is narrower than 0 to 255. Equalization can be used even when intensity levels 0 and 255 are present.

Was this helpful?